Математическая модель изотермической внутренней эрозии

УДК 519.87

Математическая модель изотермической внутренней эрозии*

А.А. Папин1, В.А. Вайгант2, А.Н. Сибин1

1 Алтайский государственный университет

2 Боннский университет

Mathematical Model of Isothermal Internal Erosion

A.A. Papin1, W.A. Weigant2, A.N. Sibin1

1 Altai State University

2 University of Bonn

Рассматривается математическая модель изотермической внутренней эрозии без учета деформации пористой среды. При достижении определенной величины скорости фильтрации происходит вынос частиц грунта из области течения. В качестве математической модели используются уравнения сохранения массы для воды, подвижных твердых частиц и неподвижного пористого скелета, а также закон Дарси для воды и подвижных твердых частиц и соотношение для интенсивности суффозионного потока. Дается постановка задачи и проводится преобразование системы уравнений. В результате преобразований для насыщенности водной фазы возникает вырождающееся на решении параболическое уравнение, для давления — эллиптическое уравнение и уравнение первого порядка для пористости грунта. Имеется аналогия с классической моделью Маскета — Леверетта. Описаны гипотезы, которые определяют интенсивности фазового перехода. Кроме того, приведен краткий обзор моделей внутренней суффозии. Рассматривается случай автомодельного движения без учета силы тяжести и скорости твердого скелета. Получено уравнение для концентрации подвижных твердых частиц грунта.

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

БМ 10.14258/izvasu(2015)1.1-16

In this paper, a mathematical model of isothermal internal erosion without deformation of a porous medium is studied. Removal of soil particles from a flow occurs at a certain magnitude of filtration velocity. Equations of mass conservation of water, moving solids, and stationary porous skeleton along with Darcy's law for water and moving solid particles and equation for the intensity of suffusion aquifer are used as a mathematical model of the problem. Section 1 of this paper describes the problem formulation and development of the system of equations. Results of the development are a parabolic equation for water phase saturation, elliptical equation for pressure, and equation of the first order for soil porosity. The developed model demonstrates similarity with the classical Masket — Leverett model. Assumptions describing the phase transition intensity are considered in section 2. Among other things, a brief review of internal erosion models is provided. A case of self-similar motion without consideration of gravity forces and a velocity of the solid skeleton is also investigated. An equation for concentration of soil flexible solid particles is deduced.

Key words: multiphase flow, porous medium,

suffusion, phase transition, saturation.

1. Постановка задачи. Рассматриваются процессы фильтрации подземных вод и внутренней суффозии. Грунт моделируется как трехфазная сплошная пористая среда. Поры полностью заполнены смесью воды = 1) и подвижных твердых частиц (^ = 2). Доля пор в грунте (^ = 3) опре-

* Работа выполнена при финансовой поддержке государственного задания Министерства №2014/2 и гранта РФФИ №13-08-01097.

деляется пористостью ф = (VI + У2)/У, где V = V+У2+У3 — общий объем грунта; У, У2, Уз — соответственно объемы воды, подвижных твердых частиц и скелета грунта.

Уравнения сохранения массы для каждой из фаз с учетом фазового перехода имеют вид [1]:

^ + У-

др2 т

дрз

т

+ V • (р2«2) = т; + V • (рэ«э) = -т,

и скелета грунта; V = - опера-

нде т — интенсивность фазового перехода (суф-фозионный поток); г?1, г?2, из — соответственно истинные скорости воды, подвижных твердых частиц грунта и скелета грунта; р1 = фв 1 р0, р2 = = фв2р°, рэ = (1 — Ф)рЭ — приведенные плотности воды, подвижных твердых частиц грунта и скелета; вх = VI/(VI + V2), в2 = V2/(VI + V2) — концентрации воды (насыщенность) и подвижных твердых частиц в порах; р0, р°,р3 — истинные плотности воды, подвижных твердых частиц грунта

д д д дх2 ' дх3

тор градиента, х = (ж1,ж2,жэ). В рассматриваемом случае р3 = р°, так как подвижные частицы захватываются суффозионным потоком из грунта.

Уравнения сохранения импульса для воды и подвижных твердых частиц грунта берем в виде [2-4]

8гф{щ-йъ) = -Ko(ф) — (^7Pi + p0ig), г = 1,2. (4) Мг

Здесь Ко(ф) — симметрический тензор фильтрации пористой среды; к^ — относительные фазовые проницаемости (к0г = к0г(вг) > 0, =о = 0, 0 < вг < 1); Мг — коэффициенты динамической вязкости; д — ускорение силы тяжести; Р1,Р2 — соответственно давления первой и второй фаз.

Пусть в = вх, тогда 1 — в = в2. Разность фазовых давлений удовлетворяет соотношению вида [2; 4]

Р2 — Р1 = Рс(в,х) > 0,

где рс — заданная функция, обладающая свойствами [5; 6]

Рс(х, в) = ро(х)^'(в), ро(х) > 0, з(в) > 0,

дз

3(0) = 0, з (1) = 1,

дв

< 0.

В дальнейшем предполагается, что скелет грунта неподвижен (иэ = 0), истинные плотности р0 постоянны. В этом случае система приводится к эллиптико-параболической системе [3] и уравнению кинетики. Действительно, используя сделанные предположения вместо (1)-(4), получим

д{зф) дЬ

+ V• (вфих) =0;

д(1 — в)ф „ .. . , . т

Л-^ + 1 -8)фй2) = —-

р02

д(1 — ф) т

р0 рэ

ког

я1фй1 = -К0{ф)к

Мг

Сложив уравнения (6), (7) и (8), выводим

V-

Положим,

V = вфмх + (1 — в)фг?2-

Используя (9) и (5) для V получим следующее представление:

—V = Ко(ф)(ко1 (УР1 + р?д)+ +к02 (у(Р1 + Рс) + р2д)) =

т>- / 7 / , ко2 дрс , ,

= Ко(к(\уР1 + V «)+

+ко2 V Рс + д(р°ко2 + р0ко1)) = 1

^02(0 дРс

= К0(к V - У '

к(С) де

¿е)+

в

+к02 V Рс + д(р°к02 + р0к01)) =

= ВД)к

где р — так называемое «приведенное» давление [3]:

1

[ к02(О дрс

в

а также введены обозначения:

к(в) = &01 + к02, 1

/ ^ и ^ ,/ С к02

/ = Ко(ко2 Х/х Рс + к J ^ V-щ-

в

+(к01р0 + к02 р°)д).

Здесь символ Vх применяется только по переменной х, входящей явно, например:

▽ хРс(в, х) =

/дРс

у дх1 ' дх2 ' дхэ С учетом (12) для «1 = вфЙ1 имеем

—V! = К0а V в + К>к01 V Р + /0,

где

&01&02 дРс

к дв

/"0 = К0^01(

ко2(0 ^ дРс(£, х) о

к

де

Используя (11), находим

Коки у Р = —ко1 (V + /)/к. Тогда вместо (13) получим

где

—VI = Коа у в — Ь« + р,

Подставляя (13) в (6), привлекая (8) и (10), приходим к системе уравнений относительно в, р и ф:

д(эф) дЬ

= у

у (Ко(ф)к(в) у Р + Л = 0;

о <9

Используя (14), приходим к следующей эквивалентной системе относительно в, р, ф и V:

У • V =0;

рд

= ~т>

—V = Ко(ф)к(в) у р + р.

Замечание 1. Решение начально краевой задачи для системы (8), (15), (16), в которой т является функцией ф, в, VI и «2, можно построить следующим образом: подставляя в (15), (16), вместо ф известную функцию фо, приходим к эллиптико-параболической системе для р и в (данная система исследовалась в [3; 5]). Найденные значения р и в позволяют определить «1, «2 и тем самым из (8) найти новое значение ф.

Особенностью рассматриваемой задачи является возможное вырождение на решение уравнения (15), поскольку а(0) = а(1) = 0, а коэффициент фильтрации, как правило, задается следующим образом: Ко = Коф3/(1 — ф)2 [6]. Кроме того, пористость и насыщенность должны удовлетворять условиям 0 < в < 1, 0 < ф< 1.

Замечание 2. Частный случай модели (6)-(9) при р1 = р2 рассматривается в работах [7-14].

2. Интенсивность фазового перехода. На основе обработки экспериментальных данных были предложены различные формулы для интенсивности фазового перехода. В работе [7] суф-фозионный поток задается следующим образом:

где Л — определяемая экспериментально функция (отвечает за устойчивость грунта суффозионному воздействию); «2 = ф(1 — в)(г/2 — из) — поток подвижных частиц грунта.

Обобщение формулы (17) предложено в серии работ [6-8]:

: т = тег — т^ер,

где тег — поток твердых частиц (процесс суффозии) ; т¿ер — поток осевших твердых частиц (процесс кольматации).

Для определения т¿ер в работе [7] предлагается использовать соотношение

т<1ер = Р2А(1 - ф)-111-

вСГ

Здесь всг — критическое значение концентрации подвижных твердых частиц грунта, при достижении которого (в2 = всг) процессы суффозии и кольматации уравновешивают друг друга. Подставив (17) и (19) в (18), получим соотношение [7]

в сг

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

V* =4.2 фуД^-

1

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

В работах [6; 11] для т используется зависимость

В работах Б. БопвШ (см. например: [16, с. 187]) движение воды, подвижных частиц и отрыв частиц от скелета моделируются на основе подходов, развитых в задачах с неизвестной границей. Вода и подвижные частицы грунта рассматриваются как однородная смесь (несжимаемая вязкая жидкость со стандартной реологией), которая движется со скоростью и и имеет плотность р = фр? + (1 — ф)р2. Неизвестная граница Г между областями, занятыми смесью и твердым скелетом, определяется из уравнения переноса вида

дф _ ^ ,

0,

где функция ф имеет следующие свойства: ф = 0 на Г, ф > 0 — в твердом скелете, а в области

2

2

фильтрации ф < 0; с = сгп, п — вектор нормали к границе Г. Скорость движения последней равна

сг =

kd(r - тс),т > тс;

0,т < тс.

где — коэффициент пропорциональности; т — модуль касательного напряжения

\/(Тп)2 — (пТп)2;

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

Тензор напряжений и тензор скоростей деформации имеют вид:

T = -PI + D(U);

Д«) = J(v« + (v«)T)-

В данном подходе интенсивность фазового перехода определяется формулой

m = р(сг — и  n).

В настоящей работе предлагается использовать следующее соотношение для определения суффозионного потока:

m = ¿(s)Ä(^)max{|vi| — , 0}.

Здесь

Ф)

0, s > 1;

1 — s, 0 < s < 1;

1, s< 0.

1, ф > 1; Д(ф) = { ф(1 — ф), 0 < в < 1; 0, ф < 0.

3. Автомодельный случай. Рассмотрим одномерное движение при условиях д = 0, «3 = 0. Система (6)-(8) преобразуется к виду:

д(яф) dvi

dt dx

д (1 — dv2

dt dx

0;

p2

д(1 — ф) m p2

dt

Решение системы (21)-(23) ищется в области (—то, сй) в предположении, что все искомые функции зависят лишь от переменной е = х — ct (с — неизвестная постоянная). Вместо (21)-(23) получим

д (вф) 1 д«1

д(1 — s^ dv2

0;

зе

ö£

p2

д (1 — ф) m

зе

P20'

Здесь искомыми являются функции в(е), ф(е), «1, «2 и постоянная с. Из (24) следует

«1 — свф = А1.

Складывая уравнения (25) и (26), получим

«2 — с

где А1 и А2 — постоянные. Определим А1 и А2 из следующих краевых условий: «1|^=0 = , «21^=0 = , в|5=0 = в+, ф|5=0 = ф+, =

ф-.

Из (27) и(28) следует

v- — cs^- = Ai; v- — c(1 — s-ф-) = A2;

v+ — cs+ ф+ = Ai; v+ — c(1 — s+ф^) = A2. Рассмотрим случай, когда

s+ ф+ = s ф , v + = v + —, = «2

Положим,

_ v+s-ф

2

^ ~V1

Й+Ф+ '

s+ ф+

±--vt.

тогда

A 1 =0, A2 =0, с

(v+ + v + )ф+

При сделанных ранее предположениях из уравнения (9) выводим

tv- , dp 1

г>1 = -Ä0fc0i —;

de

jv- , Ф2

de

Вычитая из равенства (30) и (29), предварительно домножив (29) на к02, а (30) на к0 1, получим

дРс ¿в

/г01с

Разделив (31) на к и К0, получим уравнение для в:

¿в к01с свф

=

' ¿е кК0 К0 '

На множестве А = {е|0 <в< 1; 0 <ф< 1; |«11 > «к} имеем «1 = |с|вф > 0 и т = 0. Уравнение (8) с учетом (20) на множестве А имеет вид:

с^ = А(1-5)ф(1-ф)(|фф-^).

Последняя система уравнений определяет искомые насыщенность в и пористость ф.

+

v

+

v

2

+

+

v

v

i

i

+

Заключение. В работе дан краткий литературный обзор моделей внутренней эрозии почвы. Построена новая изотермическая модель суффо-

зионного выноса грунта и рассмотрено автомодельное решение (типа бегущей волны).

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

1. Korobkin A., Khabakhpasheva T., Papin A. Waves Propagating along a Channel with Ice Cover // European Journal of Mechanics B/Fluids. - 2014. - V. 47.

2. Хабиров В.В., Хабиров C.B. Разработка газогидратов современными технологиями / / Труды Института механики УНЦ РАН. — Уфа, 2010.

3. Антонцев С.Н., Кажихов A.B., Монахов В.Н. Краевые задачи механики неоднородных жидкостей. — Новосибирск, 1983.

4. Gard S.K., Pritchett J.W. Dynamics of Gas — Fluidized Beds. Journal of Applied Phisics // Journal of Applied Phisics. —1975. — Vol. 46, №10.

5. Папин A.A. Краевые задачи двухфазной фильтрации. — Барнаул, 2009.

6. Wang J., Walters D.A., Settari A., Wan R.G. Simulation of Cold Heavy Oil Production Using an Integrated Modular Approach with Emphasis on Foamy Oil Flow and Sand Production Effects // 1st Heavy Oil Conference — 2006.

7. Vardoulakis I., Stavropoulou M., Papanas-tasiou R. Hydro-Mechanical Aspects of the Sand Production Problem // Transport in Porous Media. — 1996. — №22.

8. Vardoulakis I., Stavropoulou M., Papanas-tasiou R. Sand Erosion in Axial Flow Conditions // Transport in Porous Media. — 2001. — №45.

9. Vardoulakis I. Sand-production and sand internal erosion: Continuum modeling // Alert

School: Geomechanical and Structural Issues in Energy Production. — 2006.

10. Папин А.А., Гагарин Л.А., Шепелев В.В., Сибин А.Н., Хворых Д.П. Математическая модель фильтрации грунтовых вод, контактирующих с многолетнемерзлыми породами / / Известия Алтайского гос. ун-та. — 2013. — №1/2 (77).

11. Кузиков С.С., Папин А.А., Сибин А.Н. Численное моделирование процесса суффозион-ного выноса грунта // Сб. тр. 17-й регион. конф. по математике «МАК-2014». — Барнаул, 2014.

12. Сибин А.Н. Математическая модель деформации мерзлого грунта вблизи термокарстовых озер // Анализ, геометрия и топология : c6. тр. Всерос. молодеж. школы-семинара. — Барнаул, 2013.

13. Папин А.А., Сибин А.Н., Хворых Д.П. Об одной задаче фильтрации в условиях вечной мерзлоты // Сб. тр. 16-й регион. конф. по математике «МАК-2013». — Барнаул, 2013.

14. Кузиков С.С., Папин А.А., Сибин А.Н. Численное исследование профильной задачи внутренней эрозии в межмерзлотном водоносном слое // Известия Алтайского гос. ун-та. — 2014. — №1/2 (85).

15. Рекомендации по методике лабораторных испытаний грунтов на водопроницаемость и суф-фозионную устойчивость. — Ленинград, 1983.

16. Bonelli S. Erosion of Geomaterials. — UK, 2012.