The quasi-one-dimensional hyperbolic model of hydraulic fracturing

Full Text

Введение. В настоящей работе моделируется режим развития трещины гидроразрыва пласта (ГРП), в котором коэффициент интенсивности напряжений намного превосходит коэффициент трещиностойкости породы. В результате происходит разрушение породы на масштабах порядка размера пор в виде трещины-предвестника, фронт которой намного опережает фронт трещины ГРП, которая образуется вследствие нагнетания жидкости разрыва. В этом случае раскрытие вертикальной трещины ГРП происходит вследствие превышения давлением жидкости минимального горизонтального напряжения в породе и последующего заполнения уже существующей трещины© 2016 Самарский государственный технический университет. Образец для цитирования И л ь я с о в А. М., Б у л г а к о в а Г. Т. Квазиодномерная модель гиперболического типа гидроразрыва пласта // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2016. T. 20, № 4. С. 739-754. doi: 10.14498/vsgtu1522. Сведения об авторах 739 И л ь я с о в А. М., Б у л г а к о в а Г. Т. предвестника, образовавшейся в результате разрушения породы. Основы теории развития вертикальных трещин в пористой среде заложены в работе [1]. В последние два десятилетия как в отечественной, так и в зарубежной литературе появилось более двух десятков тысяч теоретических и прикладных работ по моделированию образования и развития трещин ГРП. В данном обзоре рассматриваются только квазиодномерные модели развития трещины ГРП, применение которых широко распространено в нефтедобывающей промышленности. В этих моделях жидкости гидроразрыва предполагаются ньютоновскими и несжимаемыми, а окружающий пористый массив считается линейно-упругим. Вследствие этого в одномерных моделях возникает линейная связь между давлением в трещине и ее раскрытием. По характеру этой связи различают локальные и нелокальные модели. В локальных моделях между раскрытием и давлением имеется конечная связь, а в нелокальных моделях эта связь представляет собой сингулярное интегральное уравнение [2]. Теоретические исследования этих моделей в основном посвящены нахождению автомодельных решений. В работах [3-6] изучаются локальные модели, а в работах [9-13] - нелокальные модели. Модель распространения вертикальных трещин в пористой среде в случае ньютоновской жидкости гидроразрыва, когда длина трещины много больше ее постоянной высоты и ширины, разработана в работах [3, 4] - локальная модель Перкинса-Керна-Нордгрена (модель PKN). В модели PKN течение в пласте не рассматривается и вводится эмпирический коэффициент фильтрационной утечки жидкости в пласт через стенки трещины. В рамках этой модели в работе [4] получены аналитические автомодельные решения в случаях малой и большой утечек, а также численные решения в общем случае. В работе [5] рассмотрено обобщение классической модели PKN с зоной пропитки с прямолинейным фильтрационным течением, зависящим только от времени в несжимаемом пласте. Получены автомодельные решения в предельных случаях малых и больших фильтрационных утечек при заданном постоянном расходе или постоянном давлении на входе в трещину. В этих предельных случаях степенные зависимости искомых функций от времени близки к аналогичным зависимостям, полученным в [4]. Также в [5] показано, что в общем случае автомодельное решение возможно только при линейной зависимости давления от времени на входе в трещину. В работе [7] предложена двухзональная модель PKN, которая обобщает модель, представленную в [5]. В работе [7] кроме зоны пропитки также учитывается изменение порового давления в дальней зоне пласта, где предполагается упругий режим фильтрации пластовой жидкости [6]. Методом последовательных приближений с применением многочленов Эрмита получено автомодельное решение при постоянном давлении на входе в трещину. В работе [8] рассмотрена та же модель PKN с зоной пропитки, что и в [5]. Получены автомодельные решения степенного и экспоненциального вида в предельных случаях малых и больших фильтрационных утечек при расходе или давлении на входе в трещину такого же вида. В [8] показано, что в общем случае конечной фильтрационной утечки автомодельное степенное решение возможно только при линейной зависимости давления от времени на входе в трещину или квадратичной зависимости расхода от времени. В работе [9] рассмотрена модель PKN с нелокальной связью раскрытия 740 Квазиодномерная модель гиперболического типа гидроразрыва пласта и давления в трещине. Получено автомодельное решение для раскрытия трещины и давления в виде разложения по полиномам Чебышева второго рода при постоянном расходе на входе и отсутствии фильтрационной утечки. Показано, что при умеренных расходах закачки форма трещины близка к эллиптической. При больших расходах градиент давления на конце трещины увеличивается, а длина трещины - уменьшается. В рамках нелокальной модели с геометрией трещины, аналогичной геометрии классической модели PKN, методами асимптотических разложений получены полуаналитические автомодельные решения для различных предельных режимов течения в трещине. В работе [10] получено полуаналитическое решение для режима течения в трещине с непроницаемыми стенками в пренебрежении коэффициентом трещиностойкости породы. В работе [11] получено решение для режима течения в трещине с непроницаемыми стенками с высоким коэффициентом трещиностойкости породы. В отличие от работ [10, 11], в работе [12] представлен численный алгоритм для режима течения с фильтрационной утечкой, но с малой трещиностойкостью породы. Было проведено сравнение полученного численного решения на малых и больших временах с соответствующими асимптотическими решениями. Решение для наиболее общего случая, когда учитывается влияние всех трех механизмов на развитие трещины ГРП - коэффициента трещиностойкости, вязкости жидкости и фильтрационной утечки, получено в работе [13], в которой использован многопараметрический метод сингулярных возмущений. Для полученного решения рассмотрены предельные режимы малой и большой вязкости жидкости гидроразрыва. Во всех указанных работах используется приближение теории смазки [14]. Как следствие, локальные модели представляют собой системы уравнений параболического типа. В данной работе рассматривается обобщение локальной модели PKN, учитывающее конвективные и нестационарные слагаемые в уравнении движения жидкости. 1. Физическое обоснование модели. В представленном выше обзоре рассмотрены только работы, в которых жидкости разрыва являются ньютоновскими. Это связано с тем, что во многих технологических процессах, связанных с ГРП, в качестве жидкости разрыва используется вода. К таким процессам относятся: 1) первая часть тестового (калибровочного) ГРП, в которой при скоростях закачки жидкости гидроразрыва пласта из скважины вытесняется пластовая вода линейным (несшитым) гелем; в этом случае в течение нескольких минут вытесняемой водой создается трещина ГРП; 2) ГРП в низкопроницаемых нефтенасыщенных пластах; 3) явление авто-ГРП - самопроизвольное распространение магистральной трещины от нагнетательной скважины. Во всех перечисленных процессах вследствие малой вязкости воды в трещине ГРП могут существовать достаточно протяженные области, в которых силы инерции имеют один порядок с вязкими силами или превосходят их. Выполним оценку числа Рейнольдса Re. Предположим, что расход жидкости при ГРП равен 2 ÷ 4 м3 /мин. Тогда при высоте трещины h порядка 10 ÷ 30 м и ширине трещины w у скважины (порядка диаметра перфорационных отверстий), равной 0.01 м, скорость течения в начале трещины v будет составлять 741 И л ь я с о в А. М., Б у л г а к о в а Г. Т. 0.11 ÷ 0.66 м/с. При этом Re = 1100 ÷ 6600. Если конечная длина трещины L (полудлина, если трещина симметрична) равна 10 ÷ 50 м, то соотношение инерционных и вязких сил для трещины в целом оценивается приведенным числом Рейнольдса ReL = w Re /L, равным 0.2 ÷ 6.6. Эта оценка показывает, что в случае, когда жидкостью разрыва является вода, нельзя пренебрегать силами инерции. Характерное время T тестового ГРП (наименее продолжительного из процессов 1-3) составляет 180 с. Если характерное расстояние L0 от стенок трещины до границ пласта порядка 200 м, а скорость продольной волны в упругом пласте a0 равна 4000 м/с, то характерное время пробега звуковой волны оценивается величиной t0 = L0 /a0 , равной 0.05 с. Отсюда следует оценка T t0 . Следовательно, процесс упругого деформирования пласта при закачке жидкости в трещину можно считать квазистационарным. 2. Квазиодномерная модель течения жидкости в вертикальной трещине. Рассматривается процесс распространения вертикальной трещины постоянной высоты h в насыщенном пористом пласте с коэффициентом проницаемости k при нагнетании несжимаемой ньютоновской жидкости в пористый массив через вертикальную скважину. Предполагается, что трещина расположена симметрично относительно скважины. Как и в модели PKN, предполагается, что высота трещины много больше ее ширины δ и много меньше ее длины L, т.е. выполняется неравенство δ h L. Кроме этого, предполагается, что пласт является линейно-упругим, а по длине трещины условия напряженного состояния изменяются незначительно, что выполняется, если изменение давления жидкости внутри трещины много меньше характерного напряжения в пласте. Это дает возможность считать, что в плоскостях, перпендикулярных направлению трещины, приближенно выполняется условие плоской деформации, а сечения не влияют друг на друга (гипотеза плоских сечений [3, 4]). Принимается, что ось 0x направлена вдоль трещины, ось 0y направлена перпендикулярно оси 0x в направлении изменения ширины трещины δ, а ось 0z направлена по высоте трещины. Далее считается, что трещина расположена симметрично относительно плоскости x0y, следовательно, верхний край трещины имеет координату z = h/2, а нижний край - координату z = -h/2. В предыдущем пункте показано, что процесс упругого деформирования можно считать квазистационарным. Следовательно, можно использовать известное стационарное решение для ширины δ симметричной трещины, находящейся в условиях плоской деформации в полубесконечном пространстве, на контур которой действует постоянное давление [15]. При введенных выше допущениях для каждого поперечного сечения трещины ГРП можно приближенно написать следующее: δ(x, z, t) = 2(1 - ν) p(x, t) - σ G h2 - z2, 4 (1) где ν, G - соответственно коэффициент Пуассона и модуль сдвига породы; p - давление в трещине; σ - постоянное минимальное горизонтальное напряжение в пласте, в направлении которого трещина раскрывается при гидравлическом разрыве. 742 Квазиодномерная модель гиперболического типа гидроразрыва пласта Интегрирование (1) по высоте пласта позволяет найти площадь поперечного сечения как функцию координаты x и времени t: S(x, t) = π(1 - ν)h2 p(x, t) - σ . 4G (2) Если S(x, t) - переменное поперечное сечение трещины, то с учетом несжимаемости жидкости можно ввести осредненные по сечению скорость, квадрат скорости и давление согласно выражениям v = 1 S vds, (1 + β) v 2 = S 1 S v 2 ds, p = S 1 S pds. (3) S В подынтегральных выражениях в (3) записаны истинные скорости и давления, а β - поправка Кориолиса на неравномерное распределение истинной скорости по сечению [16]. Для стационарных ламинарных течений с параболическим профилем скорости β = 1/3, для квазистационарных турбулентных течений β = 0. При неустановившихся течениях поправка Кориолиса является функцией времени. С учетом (3) уравнения изменения массы и количества движения несжимаемой жидкости с постоянной плотностью ρ в одномерном потоке в канале переменного сечения S(x, t) с малой кривизной (необязательно вертикального) и проницаемыми стенками имеют вид ∂ (ρS) ∂ (ρ v S) + = -ρuL Π, ∂t ∂x 2 ∂ p ∂ (ρ v S) ∂ (1 + β) ρ v S + = -S - Πτ - ρgS sin α, ∂t ∂x ∂x (4) (5) где uL - нормальная стенкам канала скорость фильтрационной утечки; Π - смоченный периметр; α - угол наклона канала к горизонту; τ - напряжение трения на стенке канала; g - ускорение силы тяжести. Как и в работе [5], можно заменить произвольное сечение S(x, t) прямоугольным сечением с площадью, задаваемой соотношением (2): S(x, t) = w(x, t)h. (6) Из (2) и (6) следует локальная связь между давлением в трещине и ее раскрытием w - осредненной шириной трещины δ по ее высоте h: p - σ = bw, b= 4G > 0, πh(1 - ν) (7) где b - жесткость трещины. С учетом (7) система уравнений (4), (5) запишется в дивергентном виде: ∂ (w) ∂ (wv) + = -2uL , ∂t ∂x ∂ (wv) ∂ (1 + β) wv 2 + bw2 /(2ρ) τ + = -2 , ∂t ∂x ρ (8) (9) 743 И л ь я с о в А. М., Б у л г а к о в а Г. Т. где учтено, что трещина вертикальна, а смоченный периметр, в силу допущений модели, следующий: Π(x, t) = 2w(x, t) + 2h ≈ 2h. Скобки осреднения в уравнениях опущены. Систему (8), (9) также можно представить в неконсервативном виде: ∂U ∂U 1 0 , +B = d, U = (w, v) , A = v w ∂t ∂x v w = B, d = -2uL , -2τ /ρ 2 (1 + β) v + bw/ρ 2 (1 + β) wv (10) A . (11) Отсюда следует, что данная система уравнений является квазилинейной. Для замыкания системы уравнений (8), (9) или (10), (11) необходимо задать скорость фильтрации жидкости через ее стенки и напряжение трения на стенках трещины. Далее предполагается следующее: 1) пласт бесконечен и в начальный момент времени пластовое (поровое) давление постоянно и равно pk ; 2) течение в пласте квазистационарно и подчиняется уравнениям упругого режима фильтрации; 3) в каждом сечении трещины фильтрационная утечка жидкости через стенки трещины не зависит от утечек в соседних сечениях и представляет собой прямолинейно-параллельный поток в пласте. При допущениях 1-3 применение автомодельного решения для прямолинейного течения в бесконечном линейно-упругом пласте при заданном постоянном давлении на границе полупространства (течение от нагнетательной галереи) [6] приводит к двухмерному распределению давления p(x, y, t) в пласте: y , (12) p(x, y, t) = p(x, t) + pk - p(x, t) erf √ 2 χt где χ = k/(µβ∗ ) - коэффициент пъезопроводности породы; β∗ = βl m + βs - коэффициент упругой сжимаемости пласта; βl m, βs - соответственно коэффициенты упругой сжимаемости пластовой жидкости и скелета породы. С учетом (12) скорость фильтрационной утечки жидкости через стенки трещины принимается в виде uL = - k ∂p(x, y, t) µ ∂y y=0 = k p(x, t) - pk , µ πχ (t - tL ) (13) где k - проницаемость пласта; µ - вязкость закачиваемой жидкости. В формуле (13) учтено запаздывание утечки вследствие распространения фронта трещины. В формуле (13) обозначено tL = f -1 (L), а L = f (t) - закон движения фронта трещины. Нужно отметить, что допущения 1-3 являются модельными и вместо (13) можно использовать другие эмпирические или аналитические выражения для 744 Квазиодномерная модель гиперболического типа гидроразрыва пласта скорости фильтрационной утечки. При известной поправке Кориолиса β полученная система уравнений (8), (9) замыкается заданием напряжения трения τ и коэффициента трения Cw при ламинарном режиме (Re 2000) течения [17]: ρv 2 λ k0 · 64 24 τ = Cw , Cw = = = , (14) 2 4 4 Re Re где k0 = 3/2 - поправочный коэффициент для каналов прямоугольного сечения, когда одна сторона прямоугольника много больше другой его стороны; λ - коэффициент гидравлического сопротивления. В формуле (14) число Рейнольдса Re определяется через гидравлический диаметр DH сечения трещины: 4S 4wh ρvDH , DH = = ≈ 2w. Re = µ Π 2(w + h) При турбулентном режиме (Re > 2000) течения [17]: τ = Cw ρv 2 , 2 Cw = λ1 λk0 k0 · 0.3164 0.08701 = = , 0.25 = 4 4 4 Re Re0.25 (15) где пренебрегаем шероховатостью стенок трещины. Здесь λ1 - коэффициент гидравлического сопротивления при турбулентном режиме течения в канале прямоугольного сечения; λ - коэффициент гидравлического сопротивления при турбулентном режиме течения в круглой трубе (формула Блазиуса); k0 = 1.1 - поправочный коэффициент для турбулентных потоков в каналах прямоугольного сечения, когда одна сторона прямоугольника много больше другой его стороны. 3. Характеристики, соотношения на характеристиках и инварианты Римана. Матрица A в (10) не вырождена при условии w = 0. Умножение векторного уравнения (10) слева на обратную матрицу A-1 дает новое представление данного уравнения: ∂U ∂U +C = f, ∂t ∂x v w , (βv 2 + ρ-1 bw)/w (2β + 1)v C = A-1 B = f = A-1 d = -2uL , 2(uL v - τ /ρ)/w , U = (w, v) . (16) (17) Согласно общей теории квазилинейных уравнений [18-21], вычисляются собственные значения матрицы C, что при условии w = 0 дает два вещественных собственных значения на решении U : λ± = (β + 1) v ± (β 2 + β) v 2 + bw/ρ. Дифференциальные уравнения характеристических кривых имеют следующий вид: C+ : C- : dx = (β + 1) v + dt dx = (β + 1) v - dt (β 2 + β) v 2 + bw/ρ, (18) (β 2 + β) v 2 + bw/ρ. (19) 745 И л ь я с о в А. М., Б у л г а к о в а Г. Т. Находятся левые собственные векторы: 1 -βv + w 1 l2 = -βv - w l1 = (β 2 + β) v 2 + bw/ρ , 1 , (β 2 + β) v 2 + bw/ρ , 1 . Матрица ΩL , составленная из левых собственных векторов, не вырождена при условии w = 0, поскольку det ΩL = det l11 l12 l21 l22 = 2 w (β 2 + β) v 2 + bw/ρ = 0, (20) где номер строки обозначает номер собственного вектора, а номер столбца - его координату. Таким образом, система уравнений (16), (17) является строго гиперболической системой [18]. Соотношения на характеристиках имеют следующий вид: 1 -βv + w (β 2 + β) v 2 + bw/ρ D+ w + D+ v = =- 1 -βv - w 2 -βv + w (β 2 + β) v 2 + bw/ρ uL + 2 (uL v - τ /ρ) , (21) w (β 2 + β) v 2 + bw/ρ D- w + D- v = = 2 βv + w (β 2 + β) v 2 + bw/ρ uL + 2 (uL v - τ /ρ) . (22) w Здесь введены операторы внутреннего дифференцирования вдоль характеристик [20]: D+ = ∂t + λ+ ∂x , D- = ∂t + λ- ∂x , а в правых частях (21) и (22) скорость фильтрации в пласт и пристеночное трение выражаются через искомые функции по формулам (13)-(15). Как известно [19], квазилинейную гиперболическую систему из двух уравнений можно записать в инвариантах Римана. Однако при любом β > 0 инварианты Римана системы (21), (22) не выражаются в элементарных функциях. Это удается сделать, если пренебречь поправкой Кориолиса, т.е. положить β = 0. В этом случае соотношения вдоль характеристик запишутся в виде 2 2 bw/ρuL + (uL v - τ /ρ) , w w 2 2 D- I2 = bw/ρuL + (uL v - τ /ρ) , w w D+ I1 = - (23) (24) где инварианты Римана суть I1,2 = v ± 2 746 bw/ρ. (25) Квазиодномерная модель гиперболического типа гидроразрыва пласта Правые части соотношений (23) и (24) выражаются через инварианты Римана с учетом (7), (13)-(15) и следующих из (25) представлений искомых функций через эти инварианты: w= ρ (I1 - I2 )2 , 16b 1 v = (I1 + I2 ). 2 При β = 0 дифференциальные уравнения характеристик (18), (19) имеют вид dx dx C+ : C- : (26) = v + bw/ρ, = v - bw/ρ, dt dt а величину a = bw/ρ = (p - σ)/ρ (27) можно трактовать как скорость малых возмущений раскрытия трещины или давления (играет роль скорости звука). Отметим, что (27) совпадает со скоростью малых возмущений давления в трубке с упругими стенками, вычисленной для линейной связи (7) через растяжимость стенок трубки [22]. Для наглядности на рисунке показаны зависимости местной скорости a малых возмущений давления (27) несжимаемой жидкости с плотностью ρ = = 1000 кг/м3 от раскрытия w трещины для различных значений высоты трещины h для породы с модулем сдвига G = 25 ГПа и коэффициентом Пуассона ν = 0.13 (песчаник). Зависимость местной скорости малых возмущений давления a от раскрытия трещины w для различных значений высоты трещины h [The dependence of local speed of the pressure disturbance on the crack-tip opening displacement for different values of fracture height] Очевидно, что уравнения модели развития трещины гидроразрыва и все предыдущие соотношения можно легко переписать в терминах скорости и давления жидкости. 4. Граничные условия для трещины ГРП. Эволюционность границ трещины. Постановка начально-краевой задачи. Далее рассматриваются вопросы о граничных и начальных условиях уравнений представленной модели развития трещины гидроразрыва пласта. Вопросы корректности начально-краевых 747 И л ь я с о в А. М., Б у л г а к о в а Г. Т. задач для систем квазилинейных гиперболических уравнений и уравнений газовой динамики рассмотрены в [19]. Вопросы о количестве граничных условий на подвижных и неподвижных границах для уравнений газовой динамики рассмотрены в [19, 21, 23], а для квазилинейных уравнений - в [18, 19, 23]. Трещина ГРП образуется в области перфорационных отверстий вертикальной скважины, т.е. входная граница является неподвижной. Обычно на практике при инициации трещины ГРП задается постоянный расход Q0 закачиваемой жидкости, поэтому граничным условием на входе трещины ГРП является следующее нелинейное граничное условие: (28) Q0 - hwv = 0. Рассмотрим вопрос об эволюционности входной границы трещины ГРП. Пусть в области x > 0 в начальный момент времени заданы начальные условия. Из области течения x > 0 к неподвижной границе x = 0 приходит одна характеристика с характеристической скоростью λ- = v - a < 0. От входной границы в область течения уходит характеристика второго семейства λ+ = = v + a > 0. Входная граница будет эволюционной, если на ней задано одно граничное условие (28), так как в этом случае на границе определяются раскрытие w и скорость течения v. Предположим, что фронт трещины представляет собой сильный разрыв ее раскрытия (давления). В случае пренебрежения неоднородностью профиля скорости (β = 0) применение интегральных законов сохранения массы и количества движения [21] к установившемуся квазиодномерному течению в канале с упругими стенками с учетом сил давления, возникающих на скачке в системе координат, связанной с разрывом, приводит к следующим условиям на скачке: 2 S2 v2 = S1 v1 - ΠuL dx, 1 2 S2 p + ρv22 = S1 p + ρv12 + 2 p(S)dS - 1 2 Πτ dx - ρ 1 ΠvuL dx. 1 Здесь считается, что скачок раскрытия имеет некоторую ширину. Индексом 1 обозначены параметры перед разрывом, а индексом 2 - параметры за разрывом. Для простоты предполагается, что плотности закачиваемой и пластовой жидкостей совпадают. При пренебрежении фильтрационной утечкой и вязкостью жидкости эти соотношения совпадают с условиями на ударной волне в идеальной несжимаемой жидкости [22] в трубке с упругими стенками, если связь p = p(S) имеет место на ударной волне. Поскольку на практике выполняется неравенство uL v, последними интегралами в правых частях законов сохранения можно пренебречь. Если, кроме того, второй интеграл в правой части закона сохранения количества движения много меньше первого, то в системе координат наблюдателя в случае линейной связи (7) условия на скачке раскрытия примут вид w1 (v1 - D) = w2 (v2 - D) = m, b b ρw1 (v1 - D)2 + w12 = ρw2 (v2 - D)2 + w22 . 2 2 748 (29) (30) Квазиодномерная модель гиперболического типа гидроразрыва пласта В формулах (29), (30) D - скорость скачка; m - поток массы нагнетаемой жидкости через разрыв. Далее рассматривается вопрос об эволюционности фронта трещины ГРП, который, по предположению, является сильным разрывом раскрытия трещины. В этом случае m = 0 и жидкость течет через разрыв. Рассмотрим течение жидкости разрыва в области x 0 с подвижной внутренней границей x = x(t) - сильным скачком раскрытия трещины ГРП, если заданы начальные условия и краевое условие (28). На скачке заданы пять неизвестных величин w1 , v1 , w2 , v2 , D. Для того чтобы разрыв раскрытия трещины был эволюционным, число заданных граничных условий на нем должно составлять N = 5 - s, где s - общее число приходящих на разрыв характеристик из областей течения слева и справа от разрыва. В [22] показано, что в ударных волнах в упругой трубке, аналогично ударным волнам в газовой динамике, скорость течения жидкости относительно скачка перед скачком сверхкритическая (|v1 - D| > a1 ), а за скачком - докритическая (|v2 - D| < a2 ). Это означает, что из области слева от разрыва на скачок приходит одна характеристика со скоростью λ+ = v2 - D + a > 0, а из области справа от разрыва на скачок приходят две характеристики со скоростями λ+ = v1 - D + a < 0 и λ- = v1 - D - a < 0. Итак, s = 3, N = 2, и задания двух законов сохранения (29) и (30) на скачке достаточно для его эволюционности. В случае m = 0 жидкость не течет через разрыв. Из (29) и (30) следует, что в этом случае течение на фронте трещины ГРП является непрерывным, поскольку выполняются равенства w1 = w2 , v1 = v2 . Скорость скачка в этом случае является скоростью распространения фронта трещины. Поскольку в данной модели развитие трещины разрыва породы не моделируется, предлагается следующая идеализированная постановка задачи о развитии трещины ГРП: при t > 0 в области x 0 найти обобщенное решение [19] гиперболической системы уравнений (16), (17) с граничным условием (28) на входной границе, условиями на бесконечности lim w(x, t) = w0 = 0, x→∞ lim v(x, t) = 0 x→∞ и начальными условиями w (x, 0) = w0 = 0, v (x, 0) = 0. Требование неравенства нулю раскрытия трещины в начальный момент времени следует из условия (20) гиперболичности системы уравнений (16), (17). 5. Аналогия с газовой динамикой. Для численного решения полученной системы уравнений методами типа Годунова [18] необходимо решить задачу о распаде произвольного начального разрыва в пренебрежении диссипативными слагаемыми - трением в трещине ГРП и фильтрационной утечкой. В этом случае при β = 0 система уравнений (16), (17) с учетом (27) примет вид ∂w ∂v ∂w +v +w = 0, ∂t ∂x ∂x ∂v a2 ∂w ∂v + +v = 0. ∂t w ∂x ∂x (31) (32) 749 И л ь я с о в А. М., Б у л г а к о в а Г. Т. Уравнения (31) и (32) аналогичны уравнениям газовой динамики в гидравлическом приближении [23], если в последних уравнениях газ считать несжимаемой жидкостью. Можно также заметить, что существует аналогия с уравнениями одномерной газовой динамики с изоэнтропическими плоскими волнами [19-21], если учесть, что в уравнениях (31) и (32) раскрытие трещины играет одновременно как роль плотности, так и роль давления. Применительно к системе уравнений (31), (32) можно ввести определения в терминах раскрытия трещины w и скорости жидкости v, инварианты Римана I1 , I2 [22], а также получить решения типа простых волн (ПВ): волн Римана (ВР) и центрированных волн и доказать утверждения, аналогичные утверждениям в одномерной газовой динамике с изоэнтропическими плоскими волнами [20]. 1. Вырожденная ПВ системы (31), (32) есть постоянное течение. Ее линиями уровня являются траектории частиц жидкости, покрывающие плоскость событий. Линиями уровня невырожденной ПВ являются прямолинейные характеристики (26), покрывающие плоскость событий. 2. Течение, удовлетворяющее системе уравнений (31), (32), является невырожденной ПВ тогда и только тогда, когда один из инвариантов Римана тождественно постоянен. Если I1 ≡ const, то линии уровня ПВ - прямолинейные характеристики C- . Если I2 ≡ const, то линиями уровня ПВ являются прямолинейные характеристики C+ . 3. (Достаточный признак ПВ) Если в непрерывном течении имеется характеристика C+ (C- ), вдоль которой постоянны w и v, причем w = 0, то с каждой ее стороны течение является либо постоянным, либо ВР. 4. Невырожденная ВР является волной сжатия (разрежения) тогда и только тогда, когда ее прямолинейные характеристики с ростом времени сходятся (расходятся) на плоскости событий. В волне сжатия частицы разгоняются, а в волне разрежения - тормозятся. Заключение. В работе предложена математическая модель развития трещины ГРП гиперболического типа, являющаяся обобщением модели PKN. Исследована эволюционность границ трещины ГРП. Предложена постановка начально-краевой задачи для представленных уравнений модели. Показано, что при пренебрежении диссипативными слагаемыми в представленной модели можно построить теорию простых волн, аналогичную теории одномерной газовой динамики с изоэнтропическими плоскими волнами. Декларация о финансовых и других взаимоотношениях.

About the authors

Aidar M Ilyasov

Ufa State Aviation Technical University

Email: amilyasov67@gmail.com
(Cand. Phys. & Math. Sci.; amilyasov67@gmail.com), Senior Researcher, Dept. of Mathematics 12, K. Marx st., Ufa, 450000, Russian Federation

Guzel T Bulgakova

Ufa State Aviation Technical University

Email: bulgakova.guzel@mail.ru
(Dr. Phys. & Math. Sci.; bulgakova.guzel@mail.ru; Corresponding Author), Professor, Dept. of Mathematics 12, K. Marx st., Ufa, 450000, Russian Federation

References

Supplementary files