CC BY
89
УДК 539.375
АСИМПТОТИЧЕСКИЙ АНАЛИЗ ПРОЦЕССА РАЗВИТИЯ ТРЕЩИНЫ ГИДРАВЛИЧЕСКОГО РАЗРЫВА ПЛАСТА
В. И. Астафьев
Самарский государственный аэрокосмический университет им. ак. С. П. Королёва (национальный исследовательский университет), факультет летательных аппаратов (№ 1),
443086, Самара, Московское шоссе, 34.
E-mail: [email protected]
Представлено автомодельное решение задачи о 'развитии трещины гидроразрыва пласта. Предполагается, что трещина ориентирована вертикально и исследование процесса ее 'развития проводится в рамках связанной задачи теории гидродинамики (течение вязких жидкостей) и линейной теории упругости (движение хрупких трещин). Найдены асимптотики автомодельного решения при больших и малых расстояниях от вершины трещины, предложены простые ап-прокимации решения во всем диапазоне изменения переменных.
Ключевые слова: трещина гидроразрыва пласта, автомодельное решение,
асимптотический анализ.
Введение. В настоящее время большинство нефтегазовых месторождений России находится на стадии завершающей разработки. Новые месторождения, вовлекаемые в разработку, содержат трудноизвлекаемые запасы углеводородов, коллекторы характеризуются низкой проницаемостью и слабым дренированием. В течение длительного срока эксплуатации скважин параметры призабойной зоны значительно ухудшаются. Это связано с изменением проницаемости, выпадением парафинов и асфальтенов и значительной обводненностью скважин. Эффективность работы таких скважин за время эксплуатации значительно уменьшается. Одним из основных методов интенсификации разработки сложных и проблемных нефтегазовых месторождений является гидравлический разрыв пласта (ГРП). В настоящее время около трети запасов углеводородов можно извлечь только с использованием этой технологии.
Гидравлический разрыв пласта (ГРП) представляет собой процесс закачки жидкости в пласт при давлении, превосходящем предел прочности породы, при котором происходит ее разрушение и образование трещины. При повышении давления в породах пласта образуются новые и открываются уже имеющиеся трещины. Этот метод применяется для освоения скважин, для повышения продуктивности нефтяных и газовых месторождений, для повышения приемистости нагнетательных скважин, при изоляции пластовых вод и т.д.
1. Постановка задачи. Впервые модель вертикальной трещины ГРП в условиях плоской деформации (рис. 1) была предложена в работе Ю. П. Желтова и С. А. Христиановича [1].
Процесс развития такой трещины описывается двумя группами уравнений. Первая группа описывает движение вязкой несжимаемой жидкости внут-
Владимир Иванович Астафьев (д.ф.-м.н., проф.), профессор, каф. теоретической механики.
ри трещины и имеет вид [2]
dw д q
dt дх ’ др _ _ Я_
дх w3’
(1)
(2)
где и>(х, I) — раскрытие трещины, д(х, £) — скорость потока жидкости в трещине, р(х, I) —давление жидкости в трещине, ¡л — вязкость жидкости.
Вторая группа описывает деформирование упругого пласта при воздействии на него давления р(х, I) со стороны трещины и сводится к уравнению [3]
р{х, t) — ро = —
Е
№
2тг(1 - г/2) Jо
dw(s, t) sds ds s2 — x2 '
(3)
где po — внешнее пластовое давление, E и v — упругие модули, l(t) — текущая длина трещины.
Таким образом, задача о развитии трещины ГРП описывается тремя уравнениями (1)-(3) для трёх неизвестных функций q(x, t), w(x, t) и p(x, t). Для уравнений (1) и (2) необходимо ещё поставить начальные и граничные условия. В качестве начального условия будем считать, что при t = 0 нам известно начальное распределение давления ро{х) в трещине (или начальное смещение берегов трещины Wq(x), связанное с ро(х) соотношением (3)), т. е.
w(x, i)|t=0 = Wo(x), 0 ^ X ^ ¿0-
(4)
Одним из граничных условий для уравнений (1) и (2) является условие закачки жидкости в трещину:
*)|.г=0 =
Qo(t)
t > о,
(5)
где Qo(t) — заданная скорость закачки.
Вторым является условие непротекания жидкости через вершину трещины х = l(t):
t>0■ ^
Учитывая, что текущая длина трещины ¿(¿) также является неизвестной функцией, определяемой в процессе решения задачи, необходимо поставить ещё одно дополнительное условие для нахождения этой зависимости. В линейной механике разрушения таким условием является критерий разрушения [4]
где Kj(t) —коэффициент интенсивности напряжений в упругом пласте с трещиной длины l(t), Кic—трещиностойкость пласта (характеристика способности материала пласта сопротивляться развитию в нем трещины нормального отрыва).
Величина Ki(t) определяется по характеру раскрытия трещины w (х, t) в её вершине как [4]
Таким образом, четыре неизвестные функции — закон распространения трещины ¿(¿), величина раскрытия трещины и}(х, ¿), характер распределения давления жидкости в трещине р(х, ¿) и скорость потока жидкости в трещине д(х, ¿) при £ > 0 находятся из трёх интегро-дифференциальных уравнений (1)—(3) в области 0 < х < ¿(¿), начального условия (4), граничных условий (5), (6) и критерия распространения трещины (8).
2. Автомодельная постановка задачи. Полное решение начально-краевой задачи с неизвестной границей (1)—(8) может быть выполнено только численно [5]. При попытках аналитического решения сложных начально-краевых задач часто прибегают к использованию автомодельных переменных и поиску автомодельных решений [6].
Подробный анализ автомодельных решений в задаче о развитии трещины ГРП был представлен в работе [7]. Ряд автомодельных решений был получен в работах [8-10]. В настоящей работе автомодельное решение строится для случая полубесконечной трещины. Модель полубесконечной трещины обычно используется при исследовании особенностей поведения решения в окрестности вершины трещины, т. е. при X —>■ ¿(¿).
В случае полубесконечной трещины начало координат удобно поместить в вершину трещины (точка х = ¿(¿), у = 0), а ось £ направить вдоль трещины в противоположном к оси х направлении. Тогда подвижная система координат (£, г]) будет связана с неподвижной системой координат соотношениями £ = -х + Щ, г] = у.
Предположим, что решение задачи (1)—(8) для полубесконечной трещины является автомодельным, т. е. д(х, ¿) = <?({), и}(х, ¿) = ги(£) и р(х, ¿) = р(£). В этом случае соотношения (1)—(3) примут следующий вид:
KT(t) = К1С
(7)
di dC
dq ydw
(9)
(10)
(И)
где V — скорость движения трещины.
Критерий распространения трещины (8) запишется как дополнительное условие на характер поведения функции ги(£) в окрестности вершины трещи-
Таким образом, автомодельная задача о развитии полубесконечной трещины ГРП сводится к решению уравнений (10) и (11) для функций р(£) и «;(£) с условием ограниченности функции р(£) на бесконечности и с дополнительным условием (12) на поведение функции ги(£) в окрестности вершины трещины.
Как видим, решение автомодельной задачи р(£) и ги(£) зависит от четырех размерных параметров: рУ, Е' = Е/{ 1 — г/2), Кіс и ро- Поэтому предварительно приведем уравнения (10)—(12) к безразмерному виду. Выберем в качестве масштаба длины величину Ь = 12цУЕ'2/р$, масштаба давлений — величину Р = ро и масштаба смещений — величину \¥ = роЬ/Е1. Тогда безразмерные функции П = (р — ро)/Р и О, = IV/}¥ будут зависеть от безразмерной переменной £ = £/Ь и безразмерной трещиностойкости к = 4 \Z2Kjc /ро л/кЬ следующим образом [8]:
где за безразмерной переменной £ сохранено прежнее обозначение £.
Таким образом, в безразмерном виде автомодельное решение П(£) и П(£) зависит лишь от одного автомодельного параметра задачи — величины безразмерной трещиностойкости
3. Решение автомодельной задачи.
3.1. Инверсия соотношения (15). При решении уравнений (14) и (15) при дополнительном условии (16) и условии ограниченности функции П(£) на бесконечности удобно обратить соотношение (15), т. е. представить его в виде интегральной зависимости функции П(£) от П(£).
В механике хрупкого разрушения инверсия соотношения (15) записывается в виде [4, 11]
ны при £ —>■ 0:
(12)
Из уравнения (9) следует, что
(13)
(Ш 1
(14)
(15)
(16)
(17)
Сходимость интеграла (18) накладывает на поведение функции П(г?) на бесконечности определённое ограничение. Пусть П(г?) ведёт себя на бесконечности как г]а. Тогда подынтегральное выражение в (18) ??1^2П(г?)/({ — г?)
будет вести себя как Для сходимости иа бескоиечиости
необходимо, чтобы а < —1/2. Следовательно, инверсию (18) можно использовать лишь в том случае, когда П(г?) ведёт себя на бесконечности как г]а при а < —1/2. Если в процессе решения (14)—(16) обнаружится, что характер убывания П(г?) на бесконечности более слабый, чем г?-1/2, инверсия (18) будет приводить к расходящемуся интегралу. На это было обращено внимание еще в работе [4], где отмечалось, что инверсия (18) может использоваться, если функция П(г?) задана лишь на некотором конечном интервале (О, (I) и равна нулю при (1 < г) < оо.
Инверсное соотношение с более слабыми требованиями, наложенными на поведение функции П(г?) на бесконечности, было приведено в работе [5] и имело следующий вид:
с 4 Г( + ог,1у/тш,1
где 0 ^ И ^ оо.
При I) —у оо из (19) вытекает формула (18), а при И = 0 соотношение (19) принимает вид
2\/£ к] о \Г)) С-Г]
Вновь предполагая, что П(г?) ведёт себя на бесконечности как г]а, из усло-[ П(г?) <Ьг)
вия сходимости / ----- на бесконечности получаем, что а < 1/2. Как
У у/Щ С-11
видим, соотношение (21) накладывает на характер поведения функции П(г?) менее слабые условия, чем (18), и даже допускает возможность возрастания П(г?) на бесконечности, но не быстрее, чем г?1/2. Из условия (16) следует, что константа С в соотношении (20) равна безразмерной трещиностойкости к, т.е. С = к.
Таким образом, искомую задачу (14)—(16) с учётом соотношения (20) можно записать как
— =______ (21)
# п2(£Г ( )
ад *(22)
^ 2^ё тг Уо Ы (-Г1 1 К >
Соотношение (22) можно проинтегрировать и записать в виде
ВД = к^+ - Г ^(£, т])Щт])(1т], (23)
к .'о
лД + л/П
где ^(£, г?) = 1п
— 2л/—.
V
Граничным условием для функции П(£) будет условие её затухания на бесконечности:
lim П(£) = 0. (24)
£—)-оо
3.2. Асимптотическое поведение решения при £ —> 0 и £ —> оо. Определим характер асимптотического поведения решения при £ —>■ 0 и £ —>■ оо. Поведение функции П(£) при £ —>■ 0 заложено в дополнительном условии (16), которое учтено в (22) после замены константы С на безразмерный параметр трещиностойкости к. Таким образом, при £ —>■ 0 можно записать, что
ВД = <Ш + о(£1/2), (25)
где По(£) = кл/£, а о(£1/2) является малой более высокого порядка, чем £^2 при £ —> 0.
Для поиска асимптотики П(£) при £ —>■ оо перепишем (22) в виде
k d i'°° Я?
П'(0 = А - - / П(£г) f (26)
2\/{ 7Г 7о V^(l - Z)
Пусть Q(£) —)• f2oo(£) = с£" при £ —> оо. Тогда из (21) следует, что П(£) —> Поо(£) = £1_2"/с2(1 — 2а). Условие (24) накладывает на а ограничение а > 1/2, т. е. на бесконечности функция П(£) растёт быстрее, чем
По(£) = клД,- Для нахождения с и а разделим левую и правую части (26) на
Qoc/ = ca;£a_1 и найдём предел при £ —>■ оо. Получим, что
Um —Ц Г = !■ (27)
£^оо са£а-1 \2\/{ Wo
Заменяя в (27) П(£г) при £ —> оо функцией Поо(£,г) = £1_2az1_2"/c2(l — 2а) и учитывая, что
к
lim ------ —-г- = 0,
2са£а_1/2
условие (27) перепишем в виде
ГП oo(£z) ^ 7Г
^соУо ПооЧО ^(1-^) 4‘ 1
Для того чтобы
lim Mg) = lim (W1"-
П'ооШ с3ск(1 - 2ск)£“-1
был конечен, необходимо положить а = 2/3. В этом случае
и Поо(^) =__________9___
iU£) 2с W
Следовательно, для константы с из (28) получается значение
18 Г°° dz
с3 = - / -щ--г = 18л/3. (29)
7Г Jо Z5/b{ 1 - г)
Таким образом, при £ —> оо можно записать, что
ВД = гМО + о(е2/3), (зо)
где Qoo = се2/3, с = (18л/3), а о(£2/3) является бесконечно малой более высокого порядка, чем £2/3 при £ —>■ оо.
Учитывая, что Поо(£) = о^1/2) при £ —>• 0, а Г2о(£) = <э(£2/3) при £ —>■ оо, из соотношений (25) и (30) можно записать простейшее приближенное представление решения системы (22):
Q = fcV£ + c£2/3. (31)
Это соотношение полностью описывает асимптотический характер поведения решения при £ —>■ 0 и £ —>■ оо. Функция П(£), соответствующая представлению (31), находится из уравнения (21) и может быть записана в виде
dri 6 /.at 1 \
= Усо Л2г?(1 + да?1/6)2 = p(nai + l + crf + J’ ( ^
где a = c/k] £ = t6.
Асимптотическое поведение функции П(£), представленной в виде (32), получается следующим:
Inf я
П0(£) = -#, Поо(0 = -
к2 > °° С2Є1/3 '
3.3. Промежуточная асимптотика решения при 0 < £ < оо. Сравнение приближенного асимптотического представления решения (31), (32) с результатами численного решения, представленного в работе [8], показывает, что выражение (31) не описывает поведение П(£) в промежуточной области 0 < £ < оо, где численное решение не соответствует ни асимптотике По(£), ни асимптотике Qoo(£). С целью нахождения аналитического характера поведения решения П(£) в промежуточной области будем искать решение П(£) и П(£) при £ —>■ 0 в виде двухчленного представления
т = клД+ш), (зз)
где Qi(£) = о(Є1/2), Т.е. lim ЩО. = 0;
V£
П(0 = ^|+Пі(0, (34)
где Пі(£) = о(1п£), т.е. lim = 0.
§—s-o in £
Подставляя двухчленные представления (33) и (34) в основные уравнения (21), (22), для функций Qi(£) и Щ(£) получаем уравнения:
dlli _ _2Qi(0
d( kVC ’ { }
= -4 Г,Д1п" + *2п‘(")^ (36)
^ кк2 Уо у V С - V
Предположим, что второй член асимптотического разложения (33) (промежуточная асимптотика) имеет вид ^1(£) = а£". Тогда соотношение (36) будет выглядеть как
а« 1 4 Г00 /?1п, + №(,)
=-^1 уч ^ (37)
Сделав замену г? = £,г(0<,г< оо), из соотношения (37) для величины а
получаем
7г£;2 1 /'00 1п£ + 1пг +/г2П1(£г)
—— аа = - Иш £1-а / —------=---------------------------ч ^ 38
Учитывая, что
[°° <1г
из соотношения (38) получаем, что а = 1, а = 43/(ттк2), где
Г°° \nzdz 2
5 = / —=------- = 7Г .
Уо \/ф - !)
Таким образом, в задаче о развитии трещины ГРП промежуточная асимптотика для функции П(£) будет иметь вид ^1(£) = а£, что отличает её от традиционной асимптотики ^1(£) = 0(£3^2) линейной механики разрушения [4, 11]. Величина П1(£) в этом случае находится из (35) и имеет вид 1Ш = -4 а^/к3.
Построим приближенную аналитическую аппроксимацию решения (22), имеющего асимптотики (33) и (34) при £ —>■ 0 и Поо(£), Поо(£) при £ —>■ оо. В этом случае приём построения (31) и (32) не проходит, так как функция Ло(£) + ^1(0 уже не является о(£2/3) при £ —> оо.
Будем искать в этом случае функцию П(£) в виде П(£) = ку/£(1 + Сс>(£)). Из условий асимптотического поведения функции Г2(£) при £ —>■ 0 и £ —>■ оо на функцию ш(^) накладываются следующие условия:
Нт ш(^) = 0, Ит -^-г4г = у-
£->ос £!/6 к
Учитывая промежуточную асимптотику (33) для функции П(£), условие
Ит = 0 можно заменить на условие Ит = —. В этом случае функ-
^0 А1/2 к сь с
ЦИЮ Сс>(£), имеющую асимптотики (л)о = и СОоо = те/в, можно аппрок-
гь К
симировать следующим образом:
= ^о(0^оо(0 = ас С2/3
{и ио(0+^оо(0 ка^ + с^/6' 1 ;
Таким образом, аппроксимация Г2(£) будет иметь вид
т = к-Л+ ас* (40)
с +
где а = Аж/к2, с = (18л/3)1/3.
Приближенное выражение для П найдём из (21) и представим его в виде
П=['^=[*___________________________________________(41)
П2(г))
асг]
а?у1/3_|_с
Графики зависимостей Г2(£) и П(£) представлены на рис. 2 и 3. Сравнение этих зависимостей с результатами численного решения [5] показывает их полное соответствие, т.е. представление (40), (41) достаточно точно аппроксимирует поведение решения системы (21), (22) во всем диапазоне 0 < £ < оо.
Рис. 2. График зависимости Й(£)
Рис. 3. График зависимости П(£)
4. Анализ решения. Из полученного решения (40), (41) следует, что давление р(х,1) = ро(1 + П(£)) начиная с некоторого расстояния х = Хо(1) становится отрицательным, и р(х,1) —> — оо при х —> 1(1). Однако отрицательные давления в жидкости физически недопустимы, т. е. полученное решение (40), (41) в некоторой окрестности вершины распространяющейся трещины
ИЗ
О < £ < d становится некорректным. Для устранения этой некорректности необходимо уточнить математическую модель (1)—(3).
Предположим, что жидкость заполняет не всю трещину 0 < £ < оо, а только лишь часть её d < £ < оо, и перед вершиной трещины 0 < £ < d находится сухая зона, не заполненная жидкостью разрыва (lag в англоязычной научной литературе [8 и др.]). Давление р(х, t) в этой сухой зоне равно нулю, граница сухой зоны определяется из условия обращения в нуль давления р(х, t) на границе «жидкость — сухая зона». Таким образом, уравнение (14) для безразмерной функции П(£) в автомодельной переменной £ = l(t) — х следует переписать в виде
Уравнение (15) и условия (16), (17) разрушения линейного упругого тела при этом не меняются, и их можно сохранить в прежнем виде.
Безразмерный параметр к является единственным параметром данной автомодельной задачи, а неизвестная величина d находится из условия непрерывности функции П(£) при £ = d, т. е.
Если, подобно тому как это было сделано в п. 3.2, провести асимптотический анализ поведения решения при £ —^ 0 и ( ->■ оо, то обнаружится, что наличие сухой зоны не меняет асимптотик П (£) при £ —>■ 0 и £ —>■ оо, т.е. вновь
В этом случае аппроксимация П(£) во всём диапазоне 0 < £ < оо функцией
где а = с/к] £ = ¿6.
Таким образом, приближенная оценка размера сухой зоны d будет находиться из следующего уравнения:
d < £ < оо.
(42)
(43)
т = клД+о(е/2), е^о, ^(0 = с£2/3 + °(С2/3), С ->■ с»-
(44)
(45)
т = клД+се/з
позволяет записать для П(£) следующее представление:
(46)
(47)
аТ 1 к2
111 ------------ ----------------- —-----------
аТ + 1 аТ + 1 6
Уравнение (48) является трансцендентным уравнением, описывающим зависимость с1(к), и может быть решено только численно. Простые асимптотические решения этого уравнения легко находятся при к —> 0 и к —> оо. Так, при к —> 0, разлагая левую часть (48) в ряд по малой величине (аТ)-1, получаем, что (аТ)~1 = к/\/3, т. е. Т = л/З/с, или с? = (л/3/( 18\/с?)1/3)6 = 1/36. Следовательно, при к —> 0 размер сухой зоны становится слабо зависящим от величины к и может быть оценен как с! ~ 0,03.
Аналогично, при к —> оо величина аТ становится малой и уравнение (48) будет иметь решение Т = (к/с)е~к^6 или с? = (к/с)6е~к'~ ¡=з 1,03 • 10~3к6е~к'\ Как видим, с возрастанием к величина с!(к) убывает и при больших значениях к становится пренебрежимо малой, т. е. наличием сухой зоны у вершины трещины можно пренебречь.
Полученная зависимость с1(к), определяемая из уравнения (48), изображена на рис. 4 и соответствует результатам численного анализа, выполненного в работе [8].
Рис. 4. Зависимость размера сухой зоны (1 от безразмерной трещиностойкости к
Заключение. Выполненный асимптотический анализ решения задачи о росте трещины ГРП показал, что рассмотрение полубесконечной трещины позволяет установить особенности поведения решения у вершины трещины, но не позволяет оценить влияние скорости закачки <5о(^) и давления ро(Ь) на характер роста трещин, так как в этом случае необходимо рассматривать трещину конечной длины. Из полученного решения для функции и>(х, I) видно, что особенности постановки данной задачи (по существу, это связанная задача «жидкость — твёрдое тело») уточняет форму раскрытия трещины по сравнению с традиционными подходами линейной механики разрушения.
Связанность постановки задачи проявилась ещё и в том, что в любой трещине ГРП обнаруживается наличие «сухой зоны», примыкающей к вершине трещины. При больших значениях параметра к её влияние несущественно, но по мере уменьшения к её значимость возрастает, и наличием этой «сухой зоны» нельзя пренебрегать в постановке задач о росте трещин ГРП.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Желтое Ю.П. О гидравлическом разрыве нефтеносного пласта// Изв. АН СССР. ОТН., 1955. - №5. - С. 3-41.
2. Лойцянский А. Г. Механика жидкости и газа. — М.: Дрофа, 2003. — 840 с.
3. Партон В. 3., Перлин П. И. Интегральные уравнения теории упругости. — М.: Наука,
1977. - 312 с.
4. Черепанов Г. 77. Механика хрупкого разрушения. — М.: Наука, 1974. — 640 с.
5. Garagash D., Detoumay Е. The Tip Region of a Fluid-Driven Fracture in an Elastic Medium // J. Appl. Mech., 2000. — Vol. 67, No. 1. — P. 183-192.
6. Баренблатт Г. 77. Подобие, автомодельность и промежуточная асимптотика. — JL: Гид-рометеоиздат, 1982. — 256 с.
7. Пергамент А.Х., Улькин Д. А. Автомодельные асимптотики в задаче о распространении трещины гидроразрыва в плоско-деформированной среде: Препринт института прикладной математики им. М. В. Келдыша РАН. — М., 2007. — 31 с.
8. Garagash D., Détournât Е. Similarity solution of a semi-infinite fluid-driven fracture in a linear elastic solid // Comptes Rendus de I Académie des Sciences - Series IIB - Mecha,nics-Physics-Chemistry-Astronomy, 1998. — Vol. 326, No. 5. — P. 285-292.
9. Астафьев В. П., Федорченко Г. Д. Автомодельное решение задачи о развитии трещины гидроразрыва пласта// Вестн. СамГУ. Естественнонаучн. сер., 2007. — №4(54). — С. 34-41.
10. Зазовский А. А., Одишария М. Г., Песляк 70. А. Автомодельные решения задачи о распространении трещины гидроразрыва в непроницаемой горной породе // Изв. АН СССР. МТТ, 1986. - №5. - С. 92-100.
11. Райс Д. Ж. Математические методы в механике разрушения / В сб.: Разрушение: Т. 2: Математические основы теории разрушения; ред. Г. Либовид. — М.: Мир, 1975. — С. 204-335.
Поступила в редакцию 19/IV/2010; в окончательном варианте — 19/VIII/2010.
MSC: 74R10, 74F10, 74B05, 74G70
ASYMPTOTIC ANALYSIS OF HYDRAULIC FRACTURE CRACK PROCESS GROWTH
V. I. Astafiev
S. P. Korolyov Samara State Aerospace University
(National Research University), Faculty of Aircraft Construction,
34, Moskovskoe shosse, Samara, 443086, Russia.
E-mail: vlastassu.samara.ru
In this work the self-similar solution of a problem of propagation of hydraulic fracture crack is presented. It is proposed that a crack is focused vertically and research of of its development process is analyzed within the coupled problem of the hydrodynamics theory (flow of viscous fluids) and the linear theory of elasticity (fracture crack growth). Asymptotics of the self-similar solution at the long and small distances from, tip of a crack were obtained,. Simple approximations of the solution in all length range are offered.
Key words: hydraulic fracture crack, self-similar solution, asymptotical analysis.
Original article submitted 19/IV/2010; revision submitted 19/VIII/2010.
Vladimir I. Astafiev (Dr. Sci. (Phys. & Math.)), Professor, Dept, of Theoretical Mechanics. 116
