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

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

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

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

Н.Н. Хлапук

1. Построение системы дифференциальных уравнений для совместного описания фильтрации в зернистой среде и суффозионных процессов, вызванных действием градиентов напора

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

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

Рассмотрим грунт, который состоит из относительно недвижимого скелета и движущихся суффозионных частиц, или которые могут перейти в состояние движения при определенных гидродинамических условиях в конкретной точке (элементарном объеме). Между частицами есть пространство, формирующееся в поры сложной структуры [24… 27]. Пусть n - пористость; m_cк, m_c, соответственно, доли объемов частиц, формирующих скелет, и суффозионных частиц в единице объема грунта. Нами принимается, что

m_cк+ m_c(x, y, t)+ n(x, y, t, m_c) = 1, (1)

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

(2)

где U_r, U , U_z - составляющие действительной средней скорости движения воды в порах; V_cr, V_c , V_cz - составляющие средней скорости движения суффозионных частиц.

Для осесимметричной фильтрации (U_r = U, V_cr = V_c) имеем

. (3)

Поскольку взаимодействие фильтрационного потока с суффозионными частицами происходит на уровне действительных скоростей в довольно сложной форме и без фазовых преобразований, то переход на фиктивные скорости движения частиц, когда они считаются ассоциированными с жидкостью, становится невозможным. При наличии в пористой среде рассматриваемого двухкомпонентного конгломерата, преобразование массы одной компоненты в массу другой не происходит. Содержимое каждой из этих компонент внутри фиксированного элементарного объема может изменяться лишь за счет входа и выхода самой массы этой компоненты через ограничивающую поверхность [18]. Поэтому уравнение изменения общей массы (3) распадается на два независимых уравнения изменения массы каждой компоненты, в результате получим систему

 (4)

Для описания движения жидкой компоненты, как известно, используется закон Дарси. Однако в данном случае, когда вследствие изменения пористости грунта n изменяется его коэффициент фильтрации k, следует использовать закон Дарси в форме (аналогичной, предлагаемой нами в [20])

nU=V=- k(r, t)h(r, t) /r. (5)

Тогда из первого уравнения системы (4) имеем объединенное уравнение движения жидкости и сохранение массы жидкой компоненты

n(r, t) / t - (r k(r, t) h(r, t) / r) / rr = 0. (6)

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

Сначала, для обоснования общей формы такого уравнения движения частиц, рассмотрим течение Куэтта [17, 23], модифицированного нами до вида, при котором подвижная плоскость заменяется подвижной пористой стенкой некоторой толщины. Очевидно, что движение стенки по недвижимой плоскости станет возможным лишь при достижении действующей на нее со стороны потока касательной силы определенного (критического) значения. Значение этой критической силы, отнесенной к единице площади поверхности (касательное напряжение), можно определить в соответствии с законом Ньютона

_c.кр = -dU/dy_y=h,

где h - глубина потока (расстояние от верхней недвижимой плоскости к подвижной стенке); y - вертикальная координата, нормальная к поверхности стенки. Из решения уравнения , где f_kp - критическое значение силы трения, которое возникает между стенкой и неподвижной нижней плоскостью, находим критический градиент напора I_kp = 2f_kp / gh и, соответственно, критическую среднюю скорость потока = f_kph/6. Как видим, для модифицированного течения Куэтта при увеличении h (для грунта это аналогично увеличению размеров пор) критический градиент уменьшается, а скорость возрастает. Пусть пористая стенка под действием потока равномерно движется со скоростью V_c, что в определенной мере можно рассматривать как модель движения массы суффозионных частиц. Тогда в результате решения уравнения , то есть du/dy_y=h =0,5 ghI_kp (значение производной при y = h найдем путем интегрирования уравнения dh/dx = ( /g)d²u/dy² при предельных условиях U(o) = 0, U(h) = V_c, приняв градиент напора в направлении движения потока x постоянным), получим

. (7)

Если рассмотреть такое же течение между коаксиальными цилиндрами радиусами r₁ и r₂ , когда один из них, внутренний (r = r₁), движется со скоростью V_c, то по аналогии с модифицированным течением Куэтта получим

, (8)

или при замене градиентов на средние скорости

. (9)

В целом для данных двух случаев можно записать

V_c = A(I - I_kp), V_c = B(), (10)

где A, B - параметры, значения которых зависит только от геометрии потока и физических характеристик жидкости. Как видим, движение твердого тела начинается лишь при достижении потоком определенной скорости (градиента), а скорость движения V_c пропорциональна разности .

Рассмотрим теперь фильтрацию в зернистой среде, которая для упрощения выкладок однозначно характеризуется средним диаметром частиц. Выделим из массива фильтрации довольно малый объем грунта W, состоящего из несуффозионных частиц со средним диаметром d. На частицы грунта действует объемная фильтрационная сила в направлении линии тока s [8, 22]

. (11)

Фильтрационная сила, действующая на все частицы N в объеме W, _w = W, а сила, соотнесенная к одной частице, ₁ = _w /N. Определив количество частиц по формуле N = 6mW/ d³, находим силу ₁

₁ = ( g d³/6m)I . (12)

Для того, чтобы выразить градиент напора І через скорость фильтрации V, воспользуемся формулой сопротивления для ламинарной фильтрации, выведенной Д.М. Минцем [13] на основании значительного экспериментального материала по фильтрации жидкостей и газов сквозь слои шарообразных частиц и справедливой для пористых сред разнообразной геометрической структуры

= 5,1 ²/Re, (13)

где - коэффициент формы частиц грунта; , Re - коэффициент сопротивления и число Рейнольдса, которые определяются по формулам

= gIn³d/6(1 - n)_sV² , Re = _sVd/6(1 - n). (14)

Из формул (12)…(14) получим

, (15)

где k_n = 131 - коэффициент пропорциональности.

Если в грунте между частицами скелета диаметром d находятся суффозионные частицы диаметром d_с, то, поскольку d_с << d, их можно рассматривать как новый грунт с такими параметрами: n_н = n/(n + m_c), m_cн = m_c/(n + m_c). Тогда с заменой n и m на n_н і m_н имеем формулу для определения гидродинамической силы, которая действует на суффозионную частицу,

= f(n, m_c , d_c). (16)

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

_с1кр = f (n, m_c, d_c), (17)

_с1р = f (n, m_c, d_c)(). (18)

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

Уравнение (17) можно записать с учетом n_o= V_кр=k_o I_кр в виде

_с1кр = f (n, m_c, d_c) , (19)

где k_o, n_o - соответственно, коэффициент фильтрации и пористость среды до начала движения частиц.

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

Уравнение (18) описывает межфазовое взаимодействие жидкости с твердой движущейся суффозионной частицей, и по форме аналогично тем уравнениям, которые используются, например, в задачах динамики водонасыщенной среды, в теории движения наносов, в моделях движения инерционных частиц в ламинарном потоке и т.п. [3, 4, 5, 7, 11, 12, 19].

Скорость движения суффозионной частицы V_c теперь можно найти из уравнения равномерного движения

_с1кр =_с1р. (20)

Обоснованием возможности применения уравнения равномерного движения (20) является то, что ускорение суффозионной частицы можно принять равным нулю. Как показала П.Я. Полубаринова-Кочина [14], в уравнениях движения отдельных частиц жидкости в недеформированной среде для реальных значений коэффициента фильтрации (k =1...100 м/сут) членами, которые содержат дифференцирование по времени, можно пренебречь. В нашем случае скорости суффозионных частиц будут меньшими скорости фильтрационного потока, поэтому их изменением во времени также можно пренебречь. Тогда из уравнений (17)…(20) получим зависимость

V_c = , (21)

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

Подставив (21) в (4), с использованием (5), получим

(22)

В уравнениях (22) неизвестными есть величины h, k, n, m_c. Целесообразно выразить пористость n и количество суффозионных частиц m_c посредством коэффициента фильтрации k с помощью соответствующих уравнений состояния.

Величину m_c можно найти из уравнения (1) как функцию пористости

m_c = а_n - n, (23)

где а_n = 1 - m_cк = const - объем, который занимают поры и суффозионные частицы в единице объема грунта.

Коэффициент фильтрации выразим через пористость, воспользовавшись, например, формулой В. Мацкрле [6]

k = (k_o /n_o) n = a_кn, (24)

где k_o, n_o - соответственно, начальный коэффициент фильтрации и пористость грунта.

Исследовательские данные Д.М. Минца и В.З. Мельцера [6] зависимости изменения коэффициента фильтрации в процессе кольматажа можно также описать линейным уравнением

k = a_кз n. (25)

Если в качестве модели потока жидкости в пористой среде рассматривать модель капиллярной трубки [2, 9], в которой изменять просвет (то есть площадь пор), то для этого случая нами получено

k = a_кт n. (26)

В формулах (25), (26) a_кз и a_кт - коэффициенты пропорциональности, которые характеризуют интенсивность изменения k в зависимости от n.

Нами показано (результаты здесь не приводятся), что более общие зависимости k(n), полученные из формулы Козени - Кармана [6] (например, кубическая зависимость), также могут быть аппроксимированы линейной зависимостью и использованы для данного вида деформаций из-за того, что диапазон изменения пористости находится в узких границах.

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

 (27)

с двумя неизвестными функциями h = h(r, t), k = k(r, t), где [a] = а, если а > 0, [a] = 0, если а 0, то есть второе уравнение вырождается до k /t = 0 при k (r, t) a_к а_n или при h/t I_кр. При этом: k(r, t) = k_o, если r_c(t) r R_o; k(r, t) = k_в, если r_o r r_в(t), где r_в(t), r_c(t) - соответственно, решения уравнений k(r, t) = k_в и I(r, t) = I_кр; k_в = a_к а_n - предельное значение коэффициента фильтрации, который отвечает полному вымыву суффозионных частиц; r_o - радиус буровой скважины (дрены); R_o - радиус области влияния.

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

2. Постановка и решение задач фильтрации при стабилизации суффозионных процессов

Рассмотрим, так сказать, «конечный» случай, что целесообразно с практической точки зрения, когда время t довольно большое, процесс суффозии (кольматажа) близок к завершению, общий процесс стабилизируется, а разность напоров H = H_o - h_o такая, при которой вокруг дрены (скважины) возникает локальная область, в которой действующие градиенты превышают критические значения. Тогда в уравнениях (27) можно поставить k/t = 0. В результате решения уравнения (rk(r)h/r)/r = 0, где k(r) = k_o, если r_c < r R_o, k(r) = f(r), если r_o r <r_c (f(r) - неизвестная функция, f(r_c) = k_o, при предельных условиях h(r_o) = h_o, h(R_o) = H_o и условиях сопряжения h(r_c + 0) = h(r_c - 0), h(r_c + 0) = h(r_c - 0), имеем

(28)

(29)

. (30)

В результате решения уравнения имеем

. (31)

Неизвестную функцию f(r) ищем путем решения следующей задачи:

(32)

которая эквивалентна интегральному уравнению

. (33)

Подставив (30) в (31), получаем трансцендентное уравнение для нахождения неизвестной величины

. (34)

Тогда будем имеем

, (35)

то есть

(36)

(37)

Как видим из (38), при стабилизации процесса действующий градиент во всей зоне r_o r <r_c постоянен и равен I_кр. Это значит, что при принятом в уравнении (27) постоянном значении критического градиента I_кр деформации прекращаются только тогда, когда действующие градиенты спадают до уровня I_кр. В действительности же, в процессе вымыва суффозионных частиц пористость и, соответственно, коэффициент фильтрации возрастает и уже их дальнейший вымыв должен происходить при другом значении градиента, то есть критический градиент в зоне суффозии изменяется в зависимости от изменения коэффициента фильтрации - I_кр= I_кр(k(r)). В данном случае, когда диапазон изменения пористости (коэффициента фильтрации) грунта небольшой, можно принять следующее аналитическое выражение для критического градиента (что согласовывается с лабораторными исследованиями)

I_кр(k) = I_кро, (38)

где , - коэффициенты, характеризующие влияние изменения коэффициента фильтрации на критический градиент, + = 1; I_кро - критический градиент при исходном значении коэффициента фильтрации k_о.

В результате, аналогично предыдущему решению задачи при k/t = 0, для случая изменения I_кр по формуле (39) получим:

f(r) = k_о, (39)

то есть

k(r) = (40)

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

I(r) = (41)

которые при замене и на а и b принимают вид

I(r) = ,

где а 1/, b = -/ - коэффициенты, характеризующие свойство грунта изменять коэффициент фильтрации под действием градиентов напора.

Приведенные решения задач фильтрации при стабилизации суффозионных процессов показывают, что для большинства практических задач фильтрации, когда нужно установить лишь конечный эффект суффозионных деформаций, достаточно построить локально нелинейные модели процесса фильтрации, учитывающие изменение коэффициента фильтрации грунта под действием больших, чем критические градиентов напора. При этом нами предлагается оставить локальную пространственно-временную линейность закона Дарси, несмотря на то, что новая модель будет базироваться на нелинейном уравнении div(k(I_д(r, t), I_кр)I) = 0, где I_д(r, t) - действующий градиент, при котором происходят деформационные процессы.

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

Необходимые коэффициенты a, b (или , ), которые являются характеристиками грунта, могут быть получены из простых опытов на установке Дарси [15] по разработанной методике [21] или вычислены с использованием [16].

фильтрация градиент стабилизация суффозия

Библиографический список

1. Аравин В.И., Нумеров С.Н. Теория движения жидкостей и газов в недеформируемой пористой среде. - М.: ГИТТЛ, 1953. 616 с.

2. Бэр Я., Заславски Д., Ирмей С. Физико-математические основы фильтрации воды: /Пер. с англ. - М.: Мир, 1971. 452 с.

3. Васильев О.Ф. Динамика потока гидросмеси на деформируемом ложе //Изд. АН СССР, ОТН. 1957. № 10. С. 34-39.

4. Ватажин А.Б., Клименко А.Ю. Континуальные модели движения инерционных частиц в ламинарном и турбулентном потоках, основанные на уравнениях Фоккера-Планка. //МЖГ. 1994. № 2. С. 27-35.

5. Векслер А.Б. Основные уравнения одноразмерного руслового потока в размываемом русле. //Изв. ВНИИГ. 1969. Т. 90. С. 169-179.

6. Гаврилко В.М., Алексеев В.С. Фильтры буровых скважин. - М.: Недра, 1985. 334 с.

7. Динамика сплошных сред в расчетах гидротехнических сооружений / Б.И. Дидух, В.Л. Лобысев, В.М. Лятхер и др. /Под ред. В.М. Лятреха и Ю.С. Яковлева. - М.: Энергия, 1976.

8. Иванов П.Л. Грунты и основания гидротехнических сооружений. - М.: Высшая школа, 1985. 352 с.

9. Кондратьев В.Н. Фильтрация и механическая суффозия в несвязных грунтах. - Симферополь: Крымиздат, 1958. 74 с.

10. Корбанова В.Н. Физические свойства горных пород. - М.: Гостоптехиздат, 1962. 490 с.

11. Лятхер В.М., Дидух Б.И. Движение водонасыщенной грунтовой среды при сейсмических воздействиях // Труды координационных совещаний по гидротехнике. Сейсмостойкость больших плотин. Дополнительные материалы. 1973. С. 111-123.

12. Лятхер В.М., Дидух Б.И. Одномерные краевые задачи водонасыщенной среды. //Труды Гидропроекта. 1971. Сб. 20. С.167-196.

13. Минц Д.М., Шуберт С.А. Гидравлика зернистых материалов. - М.: Изд-во Минкоммунхоза РСФСР, 1955. 112 с.

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

15. Рекомендации по методике лабораторных испытаний грунтов на водопроницаемость и суффозионную устойчивость. П 12-83 / ВНИИГ. - Л.: Энергия, 1983. 64 с.

16. Рекомендации по проектированию обратных фильтров гидротехнических сооружений. П 92-80 / ВНИИГ. - Л.: Энергия, 1981. 105 с.

17. Романков П.Г., Курочкина М.И. Гидромеханические процессы химической технологии. - Л.: Химия, 1982. 288 с.

18. Слезкин Н.А. О дифференциальных уравнениях фильтрации. //Докл. АН. СССР. 1951. Т. 79. № 5. С. 755-758.

19. Франкль Ф.И. К теории движения взвешенных наносов. //Докл. АН СССР. 1953. Т. 92. № 2. С. 247-250.

20. Хлапук М.М., Бомба А.Я. Особливості закону Дарсі при моделюванні процесів фільтрації в середовищах, що деформуються // Зб. статей за матеріалами III наук.-техн. конференції. Ч.2. Гідротехнічне будівництво. - Рівне, 1997. С.70-73.

21. Хлапук М.М., Яцик А.В., Іващенко А.П. Фізичне та математичне моделювання деформацій грунту в навколодренній зоні //Водне господарство України. 1999. № 1-2. С. 36-38.

22. Чугаев Р.Р. О фильтрационных силах. //Изв. ВНИИГ. 1960. Т. 63. С. 115-141.

23. Шлихтинг Г. Теория пограничного слоя: Пер. с нем. - М.: Наука, 1974. 712 с.

24. Карнаухов А.П. Модели пористых систем // Моделирование пористых систем /Сиб. отд. АН СССР. 1976. С. 42-59.

25. Карнаухов А.П. Некоторые общие принципы моделирования пористых систем // Моделирование пористых материалов. /Сиб. отд. АН СССР. 1976. С. 31-41.

26. Карнаухов А.П. Текстура и классификация пористых материалов. //Моделирование пористых материалов. /Сиб. отд. АН СССР. 1976. С. 5-30.

27. Хейфец Л.И., Неймарк А.В. Многофазные процессы в пористых средах. - М.: Химия, 1982. 320 с.

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

...