Сопряженные задачи распространения трещин гидроразрыва тема автореферата и диссертации по механике, 01.02.05 ВАК РФ

Защита состоится "¿У" А/а Я _ 1998г. в часов на заседании Диссертационного совета Д.002.87.01 при ИПМ РАН по адресу: 137526, Москва, пр-т Вернадского, 101.

С диссертацией можно проблем механики РАН.

ознакомиться в библиотеке Института

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

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

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

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

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

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

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

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

Основные задачи и цели исследования :

- Постановка и исследование задач импульсного гидравлического разрыва пласта, описывающих распространение трещин, расклиниваемых высокоскоростными потоками сжимаемого газа

- Развитие методов решения сопряжеппых задач гидравлического разрыва пласта и исследование специфических особенностей решений.

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

- Учет влияния анизотропии упругих характеристик пороуп-ругого пласта, на гидравлический разрыв;

- Построение точных опорных решений сопряженных задач распространения трещин гидроразрыва в рамках двумерных и

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

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

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

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

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

Рассмотрена возможность понижения размерности задачи о распространении трещины гидроразрыва в упругом пласте. Методами теории возмущений для псевдо-трехмерных моделей массированного гидроразрыва (сильно вытянутых трещин) по-

лучены поправки к нулевому приближению (модели трещины Перкинса-Керна), которые позволяют учесть нелокальные аффекты, связанные с конечностью длины трещины и неоднородностью распределений давления и закрепителя по ее длине.

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

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

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

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

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

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

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

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

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

Публикации: Основные результаты диссертации изложены в 31 печатных работах.

Апробация работы: Результаты исследований докладывались па всесоюзных семинарах «Современные проблемы математике и методы фильтрации» в Москве 1984г., Ш-ем семинаре «Современные проблемы теории фильтрации» 1989г., всесоюзной И-ой школе- семинаре «Физические основы прогнозирования разрушения горных пород» в Фрунзе 1985г., Х-ой школе -семинаре «Теоретические и прикладные проблемы вычислительной математике и математической физике» Риге 1985г., У11-ом всесоюзном съезде по теоретической и прикладной механике в Москве 1991г., семинарах Института Механики при МГУ в 1988г., Института проблем механики РАН в 1996г., Государственной академии нефти и газа им. И.М.Губкина в 1997г., Институт Проблем нефти и газа РАН в 1995,1996,1997г., МИФИ в 1996, 1997г..

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

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

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

( й — скорость газа; а, Ъ — функции раскрытия трещины п числа Рейиольдса; р— плотность газа.)

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

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

Анализ найденных инвариантных решений позволил сформулировать следующие утверждения.

При использовании двучленного закона сопротивления имеет место конечность скорости распространения фронта газа, как и при описании течения газа (жидкости) в канале линейным законом сопротивления. Для ламинарных течений (линейный закон сопротивления) инвариантные решения были получены и исследованы Г.И.Баренблаттом 1952, 1953г.. В этих работах, в част-

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

В другом предельном случае (больших скоростей течения) двучленный закон сопротивления переходит в квадратичный закон (турбулентный режим течения). Инвариантные решения для такого закона течения были впервые найдены А.Г.Бондаренко, Н.А.Кудряпювым, В.М.Колобашкпным 1980. Ими было показано, что скорость распространения возмущений при таком законе сопротивления бесконечна.

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

В § 2 в рамках поршневого приближения найдены инвариантные решения задачи о вытеснении одного газа другим при квадратичном законе сопротивления. Для двух жидкостей аналогичная задача была рассмотрена Л.М.Гаджиевым, И.Р.Мурадовым, С.У.Агаевым 1983. Здесь, в частности, было показано, что при больших числах Рейнольдса для поршневого вытеснения одного газа другим максимальная скорость течения достигается на границе раздела между газами. Чем больше плотность вытесняемого газа, тем с большей скоростью необходимо закачивать газ, чтобы вывести границу раздела между газами на то же расстояние.

При импульсном гидравлическом разрыве, когда скорости движеиия газа (жидкости) в трещинах велики, становится существенным влияние на поток конвективных членов в уравнении движения. Такие высокоскоростные течения рассматриваются в §3. Построенные в этом параграфе инвариантные решения позволили сформулировать следующие утверждения.

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

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

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

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

Заметим, что система уравнений, описывающая течение газа в ледеформируемых трещинах, полностью совпадает с системой уравнений, описывающей течение слабопроводящего газа в поперечных магнитных полях. Явление сильного возрастания температуры проводящего газа при его движении в магнитных полях наблюдалось экспериментально и было получено в численных экспериментах и получило название Т-слоя (Зайцев С.Г., Фаворский И.К. 1970; Васильев Р.В., Ерофеев A.B., Зуев A.A. 1987). В этих же работах было установлено и возникновение ударных волн торможения. Можно полагать, что математически обнаруженный эффект максимума имеет туже природу.

Для расчета произвольных (не описываемых инвариантными решениями) высокоскоростных потоков газа в канале с силами сопротивления в §5 предложен численный метод, основанный па методе коррекции потоков. Построенный алгоритм является алгоритмом сквозного счета и позволяет получить все не выделенные заранее поверхности разрыва. Он позволяет добиться вы-

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

Во второй главе рассматриваются задачи о распространении трещин гидроразрыва в непроницаемых пластах. При этом, разрыв пласта несжимаемой жидкостью рассматривается в рамках модели трещины Перкинса-Керна, а разрыв пласта газом -в рамках модели трещины Желтова-Христиановича. Отмстим, что эти модели используются для описания гидравлического разрыва пласта вертикальными трещинами прямоугольной формы. В обоих моделях, для того чтобы остаться в рамках двухмерных постановок задач, несмотря на трехмерную геометрию трещины, используется приближенный принцип независимости между собой соответстующих сечений трещины. Так если в модели Перкинса-Керна в каждом вертикальном сечении трещины ее раскрытие дается формулой Снеддона (решение плоской двумерной задачи упругости), то в модели Желтова-Христиановича формула Снеддона используется для определения раскрытия тре

щины в горизонтальных сечениях.

В §1 рассматриваются задачи о гидравлической! разрыве пласта в рамках модели Перкинса-Керна, когда в результате разрыва образуется вертикальная трещина постоянной высоты, сильно вытянутая в горизонтальном направлении. Движение жидкости разрыва в трещине описывается степенным законом сопротивления, причем нелинейность обусловлена либо реологией жидкости разрыва, либо турбулентным характером течения в трещине. Дано обобщение модели вертикальной трещины гидроразрыва на случай степенного закона сопротивления при движении жидкости разрыва в трещине. При этом можно моделировать разрыв ньютоновской жидкостью, неныотоновской жидкостью со степенным реологическим соотношением и разрыв при турбулентном движении ньютоновской жидкости в трещине в соответствии с квадратичным законом сопротивления. Показано, что рассматриваемая задача сводится к решению уравнения

Это уравнение относится к типу нелинейных уравнений диффу-

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

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

1(1 + акт)У\ Ь[Т) 1 ехр(ат), '

n + l-(3(m + l) „ п + 1

к =-—Ш. ^ pQ = l/(m + n + l)

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

Режимы гидравлического разрыва пласта

N Режимы (¿0 ч «0 Р/к А Х/к

1 -1</?<0 0 + — 0 п т+п -\ 2 _ «'±"±2 п -1

IV р = 0 + 0 - 0 -1 /п -(п+1)

VI Р = Ря + + 0 т + п + 2 0 0

VIII Р= Ре + + + оо Щ±11+2 " т+Г оо

р < — 1 - излив 0 = —1 - закрытая скважина 0

0 - постоянное давление закачки Р = Рд - постоянный расход

0 = 0е - экспоненциальное распространение трещины (предельное автомодельное решение) 0 > Д. - режим с обострением

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

2. п = 1, /? = & = -2 ЩХ)

Показано, что существует предельное максимальное значение скорости отбора жидкости из трещины, при котором устье трещины закрывается /?* - при Р < Р* не существует решения.

Найденные в § 1-2 автомодельные решения дают возможность решать обратную задачу определения длины трещины и параметров пласта по управляемым параметрам, контролируемым в устье скважины: давление нагнетания ро(0> Расх°Д 15 объем закачанной жидкости, который равен объему трещины.

В §3 предложен способ восстановления параметров пласта: высоты вертикальной трещины 2Я, бокового горного давления о- и длины трещины как функции времени по измерениям давления нагнетания, расхода и объема закачанной жидкости разрыва. Эти параметры не поддаются прямым измерениям, но знание их необходимо для проведения массированного гидравлического разрыва пласта и дальнейшей эксплуатации скважины. Широкий спектр найденных автомодельных решений и предложенная методика дают возможность, используя различные режимы закачки жидкости разрыва в трещину, получить надежные данные об этих параметрах пласта.

В §4 приводится постановка задач о развитии магистральных трещин в непроницаемых породах под действием расклинивающих потоков газа и анализируются физические факторы влияющие на распространение трещин. Рассмотрение проводится в рамках квазистационарной модели трещины Желтова-Христиановича. Строятся автомодельные решения, получены асимптотические решения автомодельных задач и исследован выход решений неавтомодельных задач на автомодельные асимптотики.

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

В §5 проводится исследование влияния теплообмена на течение газа в деформируемой трещине, которая может распространяться в среде.

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

В §6 рассмотрено расклинивание растущей и деформирующейся отдельной магистральной трещины высокоскоростным потоком газа. Исследование таких течений проведено как на основании построенных инвариантных решений, так и численных решений неавтомодельных задач.

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

Для новых технологий гидравлического разрыва пласта - импульсного гидроразрыва характерны не только высокие скорости закачки газа в трещины, но и высокие скорости распространения самих трещин, часто сравнимые со скоростями упругих волн в среде. При этом раскрытие трещин должно определяться из решения системы уравнений теории динамической упругости.

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

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

Основная трудность, возникающая при математическом моделировании импульснщго гидроразрыва, заключается в создапии эффективных методов решения связанной системы уравнений газовой динамики и динамической теории упругости:

1. Течение флюида в трещине

/э-^г? + Vp = —F, F = ай+ Ьр\и\й

шр— le + — = —div(wpu) — 2 wqr — p-^w

2. Динамические уравнения деформации упругого тела (а = 1,2)

Wn = Ua + Va, Ql«a = 0, 02va -

д д д д ОХ 2 ¿/i'l С/Ж1 £7Ж2

□ = + il _ 2 ^ с2а dt2

Отметим, что при взрывных нагрузках, в течение длительного промежутка времени распространение трещины происходит с постоянной скоростью, которая затем медленно убывает. В данной главе будет рассмотрена только первая стадия (v сошг). В этом случае одним из эффективных методов их решения является метод функционально-инвариантных решений. В §1 построены инвариантные решения задач о импульсном гидроразрыве пласта при постоянном давлении закачки жидкости в трещину.

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

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

механическими свойствами среды. Впервые задача о звездообразной трещине в условиях антиплоской деформации была рассмотрена Шером E.II., Мартышоком П,Л. 1976.

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

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

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

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

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

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

В §1 рассматривается распространение вертикальной трещины гидроразрыва постоянной высоты Н с раскрытием ги под действием закачиваемой в нее вязкой жидкости в рамках модели Перкинса-Керна. Вязкость жидкости разрыва считается отличной от вязкости пластовой жидкости /iq (у/ - граница раздела между вытесняющей и вытесняемой жидкостями)

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

к „ д ( д2 дг ^ М<ОД> I у l> z//; Ы>о

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

Найдено точное автомодельное решение задачи, соответствующее режиму нагнетания при постоянном давлении:

У; (а:, 0 )=2//(аз), p(0,t)=Po 18

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

В §2 рассматривается распространение трещины гидроразрыва в рамках модели Желтова-Христиановича.

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

Проанализированы асимптотики найденных автомодельных решений вблизи скважины и вершины трещины. Исследована возможность с помощью полученных решений по параметрам на скважине: давлении закачки, расхода и объема закачанной в трещину жидкости разрыва (р, д, V) определять длину трещины I, объем трещины Ус и объем жидкости V/ отфильтровавпюйся в пласт при создании соответствующих режимов закачки жидкости в пласт.

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

В общем случае решение такого рода задач проводится прямым численным моделированием и сопряжено с большими вы-

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

характерная скорость распространения трещины; «о -характерная скорость движения газа в трещине; А — константа, определяемая из решения внешней задачи; Ь

длина трещины в момент времени £.

Эта асимптотика корректна только па стадии формирования трещины гидроразрыва, когда скорость распространения трещины отлична от нуля. Отметим, что давление жидкости разрыва вершине трещины меньше чем начальное давление поровой жидкости на бесконечности. Причем давление жидкости падает по линейному закону вдоль трещины в некоторой окрестности ее вершины [Ь — X) < е. Вне рассматриваемой окрестности давление в трещине определяется гидродинамическим переносом жидкости вдоль трещины. Величина падения давления в вершине трещины определяется параметром аА, который связан со скоростью распространения трещины. При больших скоростях распространения трещины

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

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

нагнетания в них суспензии, содержащей закрепитель во взвешенном состоянии , когда раскрытие трещины и> определяется из интегрального уравнения

Цель исследования - построение псевдо-трсхмерной модели вертикальной трещины гидроразрыва переменной высоты ir большой протяженности в горизонтальном направлении при неоднородном распределении по глубине массива бокового горного давления с улетом эффектов расклинивания трещины при ее закреплении.

Общая характеристика основных стадий гидравлического разрыва неоднородно проницаемого пласта, включая закрепление трещины, рассматривается в §1.

Профиль раскрытия сильно вытянутых трещин в неоднородном пласте получен в §2 аналитическими методами двумерной теории упругости. Асимптотическим методом (Гольдштейп Р.В., Спектор A.A. 1978) получены поправки к лсевдо-трехмерной модели трещины Перкинса-Ксрна, учитывающие нелокальные эффекты, связанные с конечностью длины трещины и неоднородностью распределений давления и закрепителя по ее длине. В качестве нулевого приближения использовано аналитическое решение плоской задачи о частично закрепленной трещине гидроразрыва при неоднородном распределении по высоте бокового горного давления в массиве.

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

Чнсленпое исследование формы гидравлической трещины в приближении мгновенного осаждения закрепителя проведено в §4. Распределение бокового горного давления в массиве аппроксимировалось кусочно-постоянной зависимостью. В рамках схемы "мгновенного" закрепления трещины, проведепо исследова-

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

В §4 для кусочно-однородного профиля бокового горного давления и при постоянном давлении закачки вязкой степенной жидкости в рамках Р-ЗБ (псевдо-трехмерной) трещипи

£\р(х, - а(К)](1 - С2Г1УЧ - Кг/Щ^) « О

h(x, t) = ho/ cos <p, ip

где С = £)> Co = ho/h(x, i).

найдено автомодельное решение задачи о распространении трещины гидроразрыва

Р(х = о, t) = const, X

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

Одними из основных результатов этой главы является разработанная методика определения высоты и раскрытия трещины

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

В шестой главе строится теория гидравлического разрыва изотропного пороунругого пласта плоскими и дискообразными трещинами в насыщенных жидкостью изотропных пороупругих средах. Внимание исследователей на важность учета пороупругих эффектов в теории гидроразрыва было обращено в работах Т.J.Boon, E.Detourny 1990, 1993. Предполагается, что трещины расклиниваются потоком вязкой фильтрующейся в пласт жидкостью разрыва. Напряженное состояние и деформации пороуп-ругой среды описываются уравнениями Бпо 1941(§1).

Уравнения деформации пороупругого насыщенного жидкостью тела можно разбить на две группы, одна из которых описывает упругую деформацию среды, а другая - движение в ней жидкости. Методами теории аналитических и обобщенных аналитических функций получено комплексное представление плоских (§2) и осеспмметричных (§3) смешанных граничных задач для пороупругой среды. Развиты методы решения смешанных краевых задач для обобщенных аналитических функций, к которым сводятся задачи о пространственных осеспмметричных дискообразных трещинах гидроразрыва в изотропных пороупругих телах. Используя данпый формализм, проинтегрированы уравнения деформации среды, и получено обобщение классической теории гидроразрыва для раскрытия трещин.

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

и функциональную связь раскрытия трещины 1У с давлением жидкости разрыва в трещине р и поровым давлением в среде

ттг, ^ 1 I I1 Г Ъ