Движение проппанта в раскрывающейся трещине гидроразрыва пласта Текст научной статьи по специальности «Физика»

УДК 501.1

ДВИЖЕНИЕ ПРОППАНТА В РАСКРЫВАЮЩЕЙСЯ ТРЕЩИНЕ

ГИДРОРАЗРЫВА ПЛАСТА

А. В. Татосов, А. С. Шляпкин

Татосов Алексей Викторович, доктор физико-математических наук, профессор кафедры фундаментальной математики и механики, Тюменский государственный университет, Россия, 625003, Тюмень, Семакова, 10, atatosov@utmn.ru

Шляпкин Алексей Сергеевич, аспирант, Тюменский государственный университет, Россия, 625003, Тюмень, Семакова, 10

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

Ключевые слова: трещина гидроразрыва, пористая среда, вязкая жидкость. РО!: 10.18500/1816-9791 -2018-18-2-217-226

ВВЕДЕНИЕ

Для повышения дебита нефтяных скважин иногда применяется технология гидроразрыва пласта. В основополагающих работах [1-3] обоснованы формы и построены математические теории развития трещины гидроразрыва. Качественный обзор проведен в [4]. С целью формирования остаточной трещины в скважину подается вязкая жидкость (жидкость гидроразрыва) с примесью твердых частиц (проппант). В зарубежной литературе тема о транспорте проппанта по трещине широко освещается [5,6]. Наиболее полно изучен рост осадка, обусловленный движением вниз взвешенных частиц. В данной работе предложена модель подачи проппанта в раскрывающуюся трещину гидроразрыва пласта с учетом потерь жидкости на просачивание в пористую среду и падения взвешенных частиц под действием силы тяжести. Акцент сделан на изучении процесса осаждения частиц за счет протекания жидкости гидроразрыва в грунт. Несмотря на малую скорость фильтрации, жидкость просачивается в пористую среду через боковую поверхность полости, площадь которой значительно превосходит площадь сечения Христиановича. Закачиваемая вязкая смесь, надавливая на берега начальной трещины, расширяет ее. Ввиду просачивания жидкости гидроразрыва в пласт у берегов трещины скапливается дисперсная фаза и быстро стекает вниз. Цель данного исследования — определить влияние примеси частиц в смеси на динамику раскрытия трещины и ее конечную форму.

1. МАТЕМАТИЧЕСКАЯ ПОСТАНОВКА ЗАДАЧИ

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

Рис. 1. Схематическое представление вертикальной трещины (а); сечение в горизонтальной

плоскости (Христиановича) (б) Fig. 1. Schematic representation of a vertical crack (a); cross-section in the horizontal plane

(Khristianovicha) (b)

Форма трещины. На основании гипотезы Перкинса [1] с учетом формы плоских трещин в упругой среде введем связь избыточного давления р со средней шириной 5 трещины [7,8]

Р = Ь5, Ь = ^ , (1)

п(1 - V)Н

где — коэффициент Пуассона, — модуль сдвига материала. Геометрическая форма трещины такова, что

5 < Н < Ь.

Здесь Н, Ь — соответственно высота и длина трещины. На носике трещины, в приближении Перкинса, избыточное давление равно нулю

х = Ь(Ь) : р = 0.

В процессе подачи проппанта происходит осаждение частиц. Площадь поперечного сечения полости 5 = 5Н. Выделим часть сечения свободного от осадка, пусть 5 = еН площадь свободной области.

Уравнения неразрывности. Скорость движения жидкости и взвешенных частиц в продольном направлении будем считать одинаковой. В отличие от [8,9] поступательный поток вещества происходит только по свободной от осадка области. Уравнение неразрывности, отражающее сохранение общего объема смеси, имеет вид

I + « = (2)

дЬ дх

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

а + а2 = 1; а* — объемное содержание дисперсной фазы в осевшем слое. Условие сохранения объема частиц есть

д(а* (5 - в) + в) + д(а2ви) = ^ дЬ дх

Введем обозначение а = а2/а*, тогда

д(5 - в + ав) + д(ави) = 0 (3)

дЬ дх

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

/ дРг П

Уп = МхЛ — = — -куп,

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

к Р /„ч

Уп = - (4)

п У

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

дУ

— = Уг, Уп = в Уг, (5)

где в — объемное содержание жидкости гидроразрыва в пласте, близкое к величине пористости.

Ввиду потери жидкости через обе боковые поверхности трещины интенсивность утечки есть

I = 2кУп. (6)

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

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

а2 (т9 + Ws) = а* Ws.

Здесь ws — скорость движения контактной поверхности П относительно близлежащих к ней частиц осадка. Скорость роста объема осадка за счет падения взвешенных частиц, отнесенная к единице длины трещины, есть 5ws, или (а/(1 — а))5wg.

Интенсивность потерь жидкости гидроразрыва за счет протекания в грунт определена как / = 2Нуп. Единице объема смеси соответствует а объем осадка и (1 — а) объем жидкости, просочившейся в грунт. Таким образом, скорость роста объема осадка (отнесенная к единице длины трещины) за счет просачивания жидкости гидроразрыва в пористую среду будет определяться выражением (а/(1 — а))/.

Учитывая влияние обоих процессов на динамику роста осадка, получим:

д(S — s) dt

а

1 — а

(2hvn + 5wg).

(7)

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

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

др 12п

дХ = ^

-u.

(8)

где па — вязкость смеси.

Основная система уравнений. Все приведенные выше уравнения движения имеют место лишь при условии существования свободной от осадка области течения. Будем считать, что поток смеси вдоль трещины возможен только при s > 0. Если же в некотором сечении s = 0, а S > 0, то трещина перекрыта. В таком сечении движение смеси отсутствует, избыточное давление падает до нуля, просачивание приостанавливается. Раскрытая трещина при этом удерживается слоем частиц, находящимся в напряженном состоянии. Положим в дальнейшем h, b = const. Исключая часть переменных, систему уравнений (1)- (8) приведем к виду

дб д (eu) dt дх д (ае) + д (aeu)

—2v, 2а

де д (eu) дЪ дх

2

а

а

1 — а b

v

1а

w,

дЪ

дх

1а

v

1а

w,

u=

б2 ^

12па дх

(9)

в

где

vn, е > 0, 0, е = 0,

Величина vn определена выражением

w

i(5/h)wg, е > 0, 0, е = 0.

б

П Y

Аппроксимируем wg согласно [5]

(1 — а* а)п

Wg = -Wg0 ,

na

wg0

2 (Pp— pi )gr 9 n

Па = n(1 — а)"

где рр, р1 — плотности частиц проппанта и жидкости гидроразрыва; д — ускорение свободного падения; г — радиус частиц; а* = 0.66, т = 1.82.

Система (9) дополняется граничным условием на носике трещины

х = 1(г) : 8 = 0

(10)

и граничным условием на входе в трещину.

В сравнении с [10] предложенная модель дополнительно учитывает влияние силы тяжести и технологическую неоднородность подаваемой смеси.

2. ЗНАЧИТЕЛЬНОЕ ПРОТЕКАНИЕ В ГРУНТ

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

- < 1, V

(11)

считая, что для основного объема взвеси а и (1 — а) одного порядка. Введем безразмерные переменные:

х

х 1*

—, £

- £ = -

г* ' 8* '

£П

£

£

-, £

* '

£

£* ^ V

Учитывая вид уравнений, примем в качестве характерных величин

1 ПРо

У_

у *

(12)

1

1* = -

1

р0

4 V 6/Зк Ъ2' /вкРо,

V 6 ф1

г* =

4в кЪ2'

у * = Ро 2в Ъ '

Р0

8 =£ = у

V* = 2вкЪ. П

(13)

Здесь р0 — параметр, в качестве которого удобно принять некоторое избыточное давление на входе в трещину, такой величине давления будет соответствовать £(0,£) = 1.

Выбранный масштаб соответствует тому, что в течение времени г* в пористую среду просочилось порядка половины количества жидкости, закаченной в трещину 2v*г* = 2вУ* = 8*. Это приведет к заполнению осадком а/(1 — а) доли объема всей трещины. Следовательно, характерное время Tf заполнения всей трещины частицами за счет фильтрации можно оценить, как Tf ~ ((1 — а)/а)г* ~ г*, характерное же время падения взвешенных частиц тд ~ Н/шд. При этом ш ~ (8*/Н)тд, V ~ V*, так что условие (11) есть Tf /тд ^ 1.

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

др 12цс

дх £2

-и.

Безразмерная форма уравнений движения. Система уравнений движения в переменных (12), (13) с учетом условия (11) принимает следующую безразмерную форму (черта опущена):

д8 2 дц дг дх

де дц_ дг + дх

ц = —(1 — а)т е8 где ц = еи, V = 8/У при е > 0, V = 0 при е = 0.

1 д(ае)

V, —а дг

о ду

дх ж =1

+ 2

д (ач)

дх

а

1а

(14)

В дальнейшем, используются только безразмерные величины.

Начальные и краевые условия. Процесс гидроразрыва будем проводить при постоянном давлении на входе в трещину. Предварительно до момента т1 подается жидкость гидроразрыва без примеси частиц. Затем в течение времени т2 — проппант. Полное время закачки смеси т = т1 + т2, после чего подача смеси прекращается. Граничные условия имеют вид

5 (0,t) = 1, 0 <t < т; q (0,t) = 0, t>T; a(0,t) = i0, 0 <* ^ Tl' (15)

[.A, t>Ti.

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

5 (0,t) = 1, a(0,t) = {0; 0<t ^' (16)

I A, t > t1 .

В этом случае подача смеси поддерживается максимально продолжительное время. Время закачки т определится моментом остановки потока во входном сечении: q(0, t) = 0, при t > т.

Начальные условия соответствуют состоянию покоя

5(x, 0) = £(x, 0) = q(x, 0) = Y(x, 0) = 0, x > 0.

Анализ расчетов. Система уравнений в безразмерном виде (14) исследовалась численно. Использовалась явная разностная схема по аналогии с [11]. Схема тестировалась для случая подачи в пласт жидкости гидроразрыва без примеси частиц. При выбранном граничном условии постоянного давления на входе в трещину решение задачи автомодельно и приведено в [8]. В обратном предельном случае отсутствия оторочки воспроизводится решение из [10].

На рис. 2 приведен пример развития трещины в виде распределения величин 5 и е по ее длине в моменты времени t = 1, 2,6,10 — соответственно кривые 1-4.

а / a б / b

Рис. 2. Распределение безразмерных величин 6 (а) и е (б) по длине трещины в моменты

времени t =1, 2, 6,10 — соответственно кривые 1-4; A = 0.2, т1 = 1 Fig. 2. The distribution of the dimensionless quantities 6 (a) and е (b) along the length of the crack at the instants t = 1,2,6,10 are, respectively, curves 1-4; A = 0.2, т1 = 1

Граничные условия соответствуют случаю (16) при A = 0.2, т1 = 1. До момента t = т1 функции 5 и е совпадают (кривые 1). В дальнейшем при подаче проппанта полость трещины заполняется частицами (кривые 2). Уменьшение свободного пространства особенно заметно в начальном сечении. Здесь е меньше единицы. Контактный разрыв, отделяющий слой частиц от первоначальной оторочки, приближается к носику трещины. В дальнейшем наличие примеси частиц в жидкости гидроразрыва играет существенную роль в судьбе трещины. Начинается более интенсивное выпадение осадка в раскрывающейся ее части. В некоторый момент происходит закупоривание носика трещины и ее дальнейший рост прекращается. Ввиду продолжающегося просачивания жидкости гидроразрыва в грунт свободное пространство сокращается (кривые 3-4) и за конечный промежуток времени полностью исчезает. По мере сужения свободного пространства происходит незначительное уменьшение средней ширины трещины.

При постановке граничного условия (16) длина трещины не может превысить некоторой максимальной величины. Для различных значений параметров A и т1 наблюдается остановка роста трещины, которая происходит раньше прекращения закачки смеси. Обозначим предельную длину трещины и момент ее достижения как Lm и tm соответственно.

Исследуем влияние продолжительности предварительной подачи в скважину жидкости гидроразрыва без примеси частиц (оторочки) на процесс формирования трещины. Поставив граничное условие (16) с постоянным значением A = 0.2, варьируем параметр т1. На рис. 3 приведены зависимости предельной длины трещины, времени подачи смеси и момента остановки роста от длительности предварительной закачки. Подача оторочки увеличивает предельную длину трещины и снижает среднюю скорость ее роста до предельной длины. Отношение времени заполнения всей полости трещины частицами т к моменту достижения предельной длины tm изменяется незначительно.

Ограничим время закачки, поставив граничное условие (15). Положим A = 0.2, т = const = 10 или т2 = const = 10; параметр т1 — варьируем. В первом случае — постоянно полное время подачи смеси, во втором — время подачи проппанта. В обоих случаях меняем время закачивания оторочки. На рис. 4 представлена окончательная форма трещины; значениям т1 = 0,2, 5 соответствуют кривые 1-3.

При т1 = 0 к моменту сброса давления т трещина почти целиком заполнена частицами. Для двух других значений т1 изменения более значительны. После принудительного прекращения закачки поток смеси в начальном сечении останавливается. Из полости трещины жидкость гидроразрыва протекает в грунт. Присутствующие в ней частицы пополняют осадок. С течением времени трещина сужается по длине

Рис. 3. Зависимость предельной длины трещины, времени подачи смеси и момента остановки роста от длительности закачивания оторочки: 1 — Lm, 2 — 0.05т, пунктир — 0.5т/tm, A = 0.2

Fig. 3. Dependence of the limiting length of the crack, the time of the mixture feeding and the stopping time of growth from the duration of pumping of the rim: 1 — Lm, 2 — 0.05т, dotted line —

0.5т/tm, A = 0.2

5

1.0 —I

\ 1 1

0.5- 0.5

\ -T-1-T- 1 ' У-1-Г-1— 3 \ -г-1-T- -Г"1-T-1-T 1-1

-'-1-'-Г—'-1--Г^ -'-1-'-Г—'-1-^-1 ^

0.0 0.4 0.8 1.2 1.6 х 0.0 0.4 0.8 1.2 1.6 x

а / a б / b

Рис. 4. Окончательная форма трещины с остановкой подачи проппанта; значениям т\ = 0, 2, 5 соответствуют кривые 1-3, A = 0.2, т = const = 10 (а), т2 = const = 10 (б) Fig. 4. The final shape of the crack with stopping proppant flow; The values т\ = 0,2, 5 correspond to the curves 1-3, A = 0.2, т = const = 10 (a), т2 = const = 10 (b)

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

Исследуем влияние объемного содержания частиц в смеси на процесс формирования трещины. Подадим однородную смесь проп-панта без оторочки при различных значениях а. Поставив граничное условие (16), положим т\ = 0. На рис. 5 представлены зависимости предельной длины трещины, времени подачи смеси и момента остановки роста от доли дисперсной фазы. Уменьшение объемного содержания частиц в смеси приводит к замедлению выпадения осадка и более длительному росту трещины. В результате достигается большая длина Lm. Время закачки существенно возрастает вследствие непрерывного уменьшения потока подаваемой смеси.

ЗАКЛЮЧЕНИЕ

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

Рис. 5. Зависимость предельной длины трещины, времени подачи смеси и момента остановки роста от объемного содержания частиц: 1 — Lm, 2 — 0.05т, пунктир — 0.5т/tm, т\ = 0

Fig. 5. Dependence of the limiting length of the crack, the time of the mixture feeding and the stopping time of growth on the volume content of the particles: 1 — Lm, 2 — 0.05т, dotted line — 0.5т/tm, T1 = 0

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

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

1. Perkins T. K., Kern L. R. Widths of hydraulic fractures // Journal of Petroleum Technology. 1961. Vol. 13, iss. 09. Paper SPE 89. P. 937-949. DOI: 10.2118/89-PA

2. Nordgren R. P. Propagation of a vertical hydraulic fracture // Society of Petroleum Engineers. 1972. Vol. 12, iss. 04. Paper 7834. P. 306-314. DOI: 10.2118/3009-PA

3. Желтов Ю. П., Христианович С. А. О гидравлическом разрыве нефтеносного пласта // Изв. АН СССР. Отд-ние техн. наук. 1955. № 5. С. 3-41.

4. Есипов Д. В., Куранаков Д. С., Лапин В. Н., Черный С. Г. Математические модели гидроразрыва пласта // Вычислительные технологии. 2014. Т. 19, № 2. С. 33-61.

5. Mobbs A. T., Hammond P. S. Computer Simulations of Proppant Transport in a Hydraulic Fracture // SPE Production and Facilities. 2001. Vol. 16, № 2. P. 112-121. DOI: 10.2118/69212-PA

6. Dontsov E. V., Peirce A. P. Slurry flow, gravitational settling and a proppant transport model for hydraulic fractures // Journal of Fluid Mechanics. 2014. Vol. 760. P. 567-590. DOI: 10.1017/jfm.2014.606

7. Новацкий В. Теория упругости. М. : Мир, 1975. 256 с.

8. Ивашнев О. Е., Смирнов Н. Н. Формирование трещины гидроразрыва в пористой среде // Вестн. Моск. ун-та. Сер. 1, Математика. Механика. 2003. № 6. С. 28-36.

9. Смирнов Н. Н., Тагирова В. П. Анализ степенных автомодельных решений задачи о формировании трещины гидроразрыва // Вестн. Моск. ун-та. Сер. 1: Математика. Механика. 2007. № 1. С. 48-54.

10. Татосов А. В. Модель закачки проппанта в трещину гидроразрыва // Вычислительные технологии. 2005. Т. 10, № 6. С. 91-101.

11. Самарский А. А., Галактионов В. А., Курдюмов С. П., Михайлов А. П. Режимы с обострением в задачах для квазилинейных параболических уравнений. М. : Наука, 1987.

Образец для цитирования:

Татосов А. В., Шляпкин А. С. Движение проппанта в раскрывающейся трещине гидроразрыва пласта // Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2018. Т. 18, вып. 2. С. 217-226. DOI: 10.18500/1816-9791-2018-18-2-217-226

The Motion of Propping Agent in an Opening Crack in Hydraulic Fracturing Plast

A. V. Tatosov, A. S. Shlyapkin

Alexey V. Tatosov, https://orcid.org/0000-0002-4229-3812, Tyumen State University, 10, Semakova, Tyumen, 625003, Russia, atatosov@utmn.ru

Alexey S. Shlyapkin, Tyumen State University, 10, Semakova, Tyumen, 625003, Russia

In the present study the process of hydraulic fracture formation when pumping a viscous fluid with an admixture of particles into a well is considered. A model of a crack propagation taking account of the loss of liquid on seepage into a porous medium and the sedimentation of suspended particles under the action of gravity is developed. Detailed analysis of the sedimentation growth caused by leakage of hydraulic fracturing fluid into a porous medium is carried out. It is shown that the presence of particles has a significant effect on the pro-

cess of crack opening. The crack growth in the presence of particles is limited, its final state depends on the composition of the mixture, injection method, inlet pressure, the volume content of the particles, the volume of the rim (pure fracturing fluid without admixture). All these factors are taken into consideration in the proposed model based on special dimensionless forms of the equations of motion. The results make it possible to estimate the crack residual and choose the technological parameters to provide the desired state of the ruptured formation.

Key words: hydraulic fracture, porous medium, viscous liquid. References

1. Perkins T. K., Kern L. R. Widths of Hydraulic Fractures. Journal of Petroleum Technology, 1961, vol. 13, iss. 09, paper SPE 89, pp. 937-949. DOI: 10.2118/89-PA

2. Nordgren R. P. Propagation of a Vertical Hydraulic Fracture. Society of Petroleum Engineers, 1972, vol. 12, iss. 04, paper 7834, pp. 306-314. DOI: 10.2118/3009-PA

3. Zheltov Yu. P., Khristianovich S. A. O gidravlicheskom razryve neftenosnogo plasta [On hydraulic fracturing of oil reservoirs]. Izv. AN SSSR, Otd-nie tekhn. nauk [Proc. USSR Acad. Sci. Sect. Tech. Sci.], 1955, no. 5, pp. 3-41 (in Russian).

4. Esipov D. V., Kuranakov D. S., Lapin V. N., Chernyi S. G. Mathematical models of hydraulic fracturing. Computational technologies, 2014, vol. 19, no. 2, pp. 33-61 (in Russian).

5. Mobbs A. T., Hammond P. S. Computer Simulations of Proppant Transport in a Hydraulic Fracture. SPE Production & Facilities, 2001, vol. 16, no. 2, pp. 112-121. DOI: 10.2118/69212-PA

6. Dontsov E. V., Peirce A. P. Slurry flow, gravitational settling and a proppant transport model for hydraulic fractures. Journal of Fluid Mechanics, 2014, vol. 760, pp. 567-590. DOI: 10.1017/jfm.2014.606

7. Novatsky B. Teorija uprugosti [Theory of Elasticity]. Moscow, Mir, 1975. 256 p. (in Russian).

8. Ivashnev O. E., Smirnov N. N. Formirovanie treshhiny gidrorazryva v poristoj srede [Formation of a hydraulic fracture in a porous medium]. Vestnik Moskovskogo Universiteta. Ser. 1, Matematika. Mekhanika, 2003, no. 6, pp. 28-36 (in Russian).

9. Smirnov N. N., Tagirova V. P. Analiz stepennyh avtomodel'nyh reshenij zadachi o formirovanii treshhiny gidrorazryva [Analysis of power-law self-similar solutions of the problem of fracture formation]. Vestnik Moskovskogo Universiteta. Ser. 1, Matematika. Mekhanika, 2007, no. 1, pp. 48-54 (in Russian).

10. Tatosov A. V. Model of crack filling in hydraulic fracturing. Computational technologies, 2005, vol. 10, no. 6, pp. 91-101 (in Russian).

11. Samarskii A. A., Galaktionov V. A., Kurdyumov S. P., Mikhailov A. P. Rezhimy s obostre-niem v zadachah dlja kvazilinejnyh parabolicheskih uravnenij [Regimes with peaking in problems for quasilinear parabolic equations]. Moscow, Nauka, 1987. 480 p. (in Russian).

Cite this article as:

Tatosov A. V., Shlyapkin A. S. The Motion of Propping Agent in an Opening Crack in Hydraulic Fracturing Plast. Izv. Saratov Univ. (N.S.), Ser. Math. Mech. Inform., 2018, vol. 18, iss. 2, pp. 217-226 (in Russian). DOI: 10.18500/1816-9791-2018-18-2-217-226