Авторефераты по всем темам  >>  Авторефераты по земле

На правах рукописи

УДК 622.276.1/.4 Доманюк

Федор Николаевич РАЗРАБОТКА АНАЛИТИЧЕСКИХ МЕТОДОВ ПРОГНОЗИРОВАНИЯ ПРОИЗВОДИТЕЛЬНОСТИ ГОРИЗОНТАЛЬНЫХ И СЛОЖНОПРОФИЛЬНЫХ СКВАЖИН

Специальность 25.00.17 - Разработка и эксплуатация нефтяных и газовых месторождений

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата технических наук

Москва - 2012

Работа выполнена в Российском государственном университете нефти и газа имени И.М. Губкина

Научный консультант: доктор технических наук, профессор Золотухин Анатолий Борисович

Официальные оппоненты: доктор технических наук Крянев Дмитрий Юрьевич кандидат технических наук Ибрагимов Ильдар Ильясович

Ведущая организация: Институт проблем нефти и газа РАН ОАО

Защита диссертации состоится л13 марта 2012 г. в ауд. 731 в 15 час.

00 мин. на заседании диссертационного совета Д.212.200.08 при Российском государственном университете нефти и газа имени И.М. Губкина по адресу:

119991, Москва, Ленинский проспект, 65.

Автореферат размещен на интернет-сайте РГУ нефти и газа имени И.М. Губкина www.gubkin.ru л7 февраля 2012 г. и направлен для размещения в сети Интернет на сайте Министерства образования и науки РФ по адресу referat_vak@mon.gov.ru л7 февраля 2012 г.

С диссертацией можно ознакомиться в библиотеке РГУ нефти и газа имени И.М. Губкина.

Автореферат разослан: л7 февраля 2012 г.

Ученый секретарь диссертационного совета, д.т.н., проф. Сомов Б.Е.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

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

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

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

Задачи исследования Для достижения поставленной цели в работе решены следующие задачи:

1. Исследование особенностей стационарного притока жидкости к скважинам с прямолинейной и изогнутой в вертикальной плоскости траекторией ствола.

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

3. Решение задачи прогнозирования производительности скважины с синусоидальным профилем ствола в круговом пласте.

4. Анализ влияния параметров пласта и траектории ствола на продуктивность скважины с синусоидальным профилем и оценка области эффективного использования таких скважин.

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

6. Решение задачи определения оптимальной траектории и длины ствола протяженной скважины в анизотропном пласте.

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

Научная новизна 1. Поставлена и решена задача о притоке жидкости к скважине с синусоидальным профилем ствола.

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

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

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

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

Защищаемые положения 1. Методы расчета продуктивности скважин с прямолинейным профилем ствола при произвольном их положении относительно кровли и подошвы пласта.

2. Точные и приближенные решения задачи притока жидкости к скважине с синусоидальным профилем ствола.

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

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

Апробация работы Основные результаты исследования представлены на следующих конференциях и семинарах:

1. Научно-практическая молодежная конференция Новые технологии в газовой отрасли: опыт и преемственность, Москва, Газпром ВНИИГАЗ, 67 октября 2010 г.

2. Российская нефтегазовая техническая конференция SPE Передовой опыт и инновационные технологии при разработке зрелых месторождений и освоении новых регионов, Москва, 26-28 октября 2010 г.

3. Научно-практическая конференция Математическое моделирование и компьютерные технологии в процессах разработки месторождений, добычи и переработки нефти, Уфа, 26-28 апреля 2011 г.

4. Первый российский нефтяной конгресс, Москва, 14-16 марта 2011г.

5. Научно-технический семинар Актуальные вопросы проектирования разработки месторождений углеводородов, Москва, Газпром ВНИИГАЗ, 26 мая 2011 г.

6. Научные семинары кафедры Разработки и эксплуатации нефтяных месторождений РГУ нефти и газа имени И.М. Губкина, 2011 г.

Публикации По теме диссертационной работы опубликовано 6 печатных работ, из них 3 статьи в изданиях, включенных в Перечень российских рецензируемых научных журналов ВАК РФ.

Объем работы Диссертационная работа состоит из введения, 4 глав и основных выводов по работе. Содержание работы изложено на 125 страницах машинописного текста, включает 46 рисунков, 7 таблиц и список литературы из 1наименований.

Благодарности Автор выражает огромную благодарность своему научному руководителю, доктору технических наук Золотухину Анатолию Борисовичу, а также Андрющенко Н.В., Хруленко А.А. и Булаевой К.В. за ценные советы, консультации и поддержку в процессе выполнения работы.

СОДЕРЖАНИЕ РАБОТЫ

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

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

Вопросам стационарного притока жидкости и газа к горизонтальным и наклонным скважинам посвящены работы З.С. Алиева, К.С. Басниева, В.В.

Бондаренко, Ю.П. Борисова, С.Н. Бузинова, С.Н. Закирова, А.Б. Золотухина, Л.С. Лейбензона, В.Д. Лысенко, В.П. Меркулова, В.П. Пилатовского, А.М.

Пирвердяна, П.Я. Полубариновой-Кочиной, Б.Е. Сомова, А.П. Телкова, В.П.

Табакова, И.А. Чарного, В.А. Черных, В.И. Щурова, D.K. Babu, H. Cinco-Ley, F.M. Giger, A.S. Odeh, J.Besson, R.M. Butler, M.J. Economides, K. Furui, S.D.

Joshi, J. Lu, D. Michelevicius, G. Renard, J.M.Dupuy, A.C. Van der Vlis и других исследователей.

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

Рассматривается модель кругового пласта, вскрытого скважиной с прямолинейным профилем ствола. Решение задачи определения производительности скважины при линейном законе фильтрации сводится к решению уравнения Лапласа в трехмерной постановке с соответствующими граничными условиями. Основная идея метода заключается в том, что скважина моделируется в виде прямолинейной цепочки из N сфер с одинаковыми радиусами rc, ориентированных вдоль заданной траектории, при этом поверхности соседних сфер соприкасаются друг с другом. Впервые такой способ моделирования скважины был предложен в работах D. Michelevicius и А.Б.

Золотухина.

Условие непротекания на кровле и подошве пласта задается с использованием метода отображения источников-стоков. Для этого каждая сфера зеркально отражается относительно кровли и подошвы пласта бесконечное число раз.

Рис.1. Моделирование кровли и подошвы пласта методом отображения стоков Из общей теории фильтрации известно, что распределение давления вблизи единичной сферической скважины записывается в виде:

qB, (1) P C 2 2 4 k x x0 y y0 z z где q - дебит сферы; k - проницаемость; - вязкость жидкости; В - объемный коэффициент; С - константа, определяемая из граничных условий;

x0, y0, z0 - координаты центра сферы;

С использованием принципа суперпозиции получено решение для определения давления в любой точке пласта М (x,y,z) при работе прямолинейной цепочки сфер, произвольно ориентированной относительно кровли и подошвы пласта:

2 2 N x xi y yi z 2hj zi B, PМ C q i 4k i1 j 2 2 x xi y yi z 2hj (h 2i) zi (2) где qi - дебит i-ой сферы; i - расстояние от центра i-ой сферы до срединной плоскости пласта; h - толщина пласта; xi, yi, zi - координаты центра iой сферы; j - порядковый номер сферы-отображения со стороны кровли (берется со знаком л+) и подошвы пласта (берется со знаком л-).

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

N N 4Nk Pk Pз sij , (3) qi B s1 i1 j где sij - коэффициенты, зависящие от траектории скважины в пространстве:

sij 2 2 xs xi rc2 z 2hj zi xs xi rc2 z 2hj h 2i zi .

Rk z 2hj zi Rk z 2hj h 2i zi (4) Анализ уравнения (3) показал, что для случая прямолинейной конфигурации траектории скважины бесконечные суммы с незначительной потерей точности могут быть аппроксимированы простыми аналитическими зависимостями. Уравнение производительности скважины с прямолинейным профилем, вскрывшей анизотропный пласт круговой формы и залегающей под углом к вертикали, получено в следующем виде:

Pk Pз 2 kHh , (5) Q B 2hsin ln eRk h ln 2Lsin L e 1 2 / h 1 r c где L - длина скважины; kH - проницаемость в горизонтальном направлении; - угол наклона скважины к вертикали; = ((h - Lcos)/4h)0,5 - параметр, учитывающий неполноту вскрытия пласта; = (kH/kV)0,5 - коэффициент анизотропии; kV - проницаемость в вертикальном направлении; e - основание натурального логарифма; Т = arctg(tg/) - угол наклона к вертикали в эквивалентной изотропной среде; LТ = Lcos/cosТ - длина скважины в эквивалентной изотропной среде. Уравнение (5) справедливо при выполнении условий: 1,4Rk > Lsin; L > h.

Сопоставление точного решения (3) с приближенным аналитическим уравнением притока жидкости к прямолинейной скважине (5) свидетельствует о приемлемой точности последнего: максимальная относительная погрешность приближенного решения не превышает 3%.

Отметим, что полученное решение для определения производительности прямолинейной скважины справедливо при любом положении ствола относительно непроницаемой кровли и подошвы пласта. Если = 0, то уравнение (5) переходит в решение для совершенной по степени вскрытия наклонной скважины. В том случае, когда = 0,5, уравнение преобразуется в решение для ГС. Если = 0 и L = h уравнение (5) преобразуется в известное решение Дюпюи для вертикальной скважины.

Решение (5) позволяет провести количественный анализ влияния угла наклона скважины прямолинейного профиля к вертикали на ее производительность в пластах с различной вертикальной анизотропией. На рис. 2 представлена зависимость отношения дебитов скважины прямолинейного профиля и горизонтальной скважины от угла наклона прямолинейного ствола к вертикали. Когда = 90, ствол скважины занимает строго горизонтальное положение и рассматриваемое отношение равно 1. Уменьшение зенитного угла в изотропной среде приводит к незначительному уменьшению производительности скважины с прямолинейным профилем по отношению к горизонтальному положению ствола. Ситуация меняется в пластах с ярко выраженной вертикальной анизотропией: уменьшение приводит к увеличению дебита. Так в пласте с kH/kV = 50 производительность совершенной по степени вскрытия наклонной скважины почти на 10% превышает дебит ГС той же длины. Увеличение длины ствола приводит к ослаблению эффекта прироста дебита скважины прямолинейного профиля при ее повороте (уменьшении ) в вертикальной плоскости. Это объясняется тем, что диапазон значений угла поворота при таком соотношении длины скважины и толщины пласта составляет всего несколько градусов.

Рис.2. Зависимость отношения дебитов наклонной и горизонтальной скважин от угла наклона к вертикали : L/h = Во второй главе подробно рассматриваются вопросы, связанные с притоком жидкости к скважинам с искривленным в вертикальной плоскости профилем ствола, приводятся точные и приближенные решения для расчета производительности таких скважин.

Из общей теории фильтрации жидкости и газа к ГС известно, что наиболее привлекательными кандидатами для горизонтального бурения являются относительно тонкие пласты с хорошей вертикальной проницаемостью. В пластах с высоким отношением kH/kV применение ГС становится малоэффективным в силу очень высоких фильтрационных сопротивлений в области конвергенции линий тока вблизи ствола скважины. В настоящее время разработаны технологии бурения, которые обеспечивают искусственное искривление траектории ствола скважины. Многократное искривление горизонтального ствола в плоскости, перпендикулярной напластованию пород, снижает влияние низкой вертикальной проницаемости за счет уменьшения области, в которой течение флюида происходит вдоль вертикальной координаты. Для обозначения таких скважин будем использовать термин Ускважины с синусоидальным профилем стволаФ. Опыт применения скважин с синусоидальным профилем на месторождениях шельфа Аляски свидетельствует об их высокой эффективности при разработке пластов с низкой проницаемостью в вертикальном направлении.

Анализ литературных источников показывает, что проблема прогнозирования производительности сложнопрофильных скважин слабо изучена. В настоящее время известны всего три аналитические модели, описывающие приток к таким скважинам. В моделях R. Kamcom et al. и В.П. Табакова рассматривается пласт с нулевой вертикальной проницаемостью, вскрытый синусоидальным стволом с упрощенной траекторий ствола (рис.3).

Рис. 3. Схематизация траектории синусоидальной скважины:

1 - траектория сложнопрофильной скважины; 2 - упрощенная траектория Допущение об отсутствии потока в вертикальном направлении ограничивает область применения этих моделей сильнослоистыми пластами с очень низкой проницаемостью в вертикальном направлении. В работах В.А.

Черных получено точное решение задачи о притоке жидкости к синусоидальной скважине, справедливое, однако, только в изотропных средах.

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

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

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

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

Pk Pз 2 kh , (6) Q B eRk h 2hsin ln ln 2Lс sin Lс e 1 rc n где Lс n - длина синусоидального ствола; n - количество i iпрямолинейных участков; - длина прямолинейного участка; - угол наклона прямолинейного участка к вертикали.

Сопоставление уравнения (6) с точным решением системы уравнений (3) показало, что максимальная погрешность приближенного решения в широком диапазоне исходных параметров не превышает 3%.

Решение (6) справедливо в изотропной среде. Для учета вертикальной анизотропии был использован известный по работам М. Маскета способ масштабного преобразования координат. В вертикально-анизотропной среде уравнение (6) преобразуется к виду:

Pk Pз 2 kHh , (7) Q B eRk h 2hsin ln ln 2Lc sin Lc e 1 rc где Т - угол наклона прямолинейного участка к вертикали в эквивалентной изотропной среде; LсТ = Lсcos/cosТ - длина синусоидальной скважины в эквивалентной изотропной среде. Решение (7) справедливо при выполнении условий: 1,4Rk > Lsin; L > h.

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

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

Значительное влияние на производительность скважины с синусоидальным профилем оказывает толщина пласта. Увеличение толщины приводит к росту параметра . Эффект увеличения во многом объясняется расширением диапазона изменения угла : чем больше толщина пласта, тем меньше может быть угол наклона прямолинейного участка к вертикали.

Рис. 4. Влияние угла наклона и вертикальной анизотропии на гидродинамическую эффективность двухцикловой синусоидальной скважины, вскрывающей пласт на всю толщину Полученные выше решения (6)-(7) справедливы для упрощенной траектории синусоидальной скважины. С целью количественной оценки эффекта искривления скважины в работе проведена серия численных экспериментов для участков синусоидальной скважины различной кривизны. Моделирование криволинейных участков проводилось дугами окружностей равного радиуса, рис. 5. Степень кривизны траектории в модели регулировалась изменением радиусов сопрягаемых окружностей: чем меньше радиус окружностей, тем больше длина искривленного участка скважины, заключенного между точками 1 и 3, и тем выше степень кривизны траектории в сравнении с ее прямолинейной аппроксимацией. Определение производительности синусоидальных скважин с криволинейной и упрощенной прямолинейной траекториями ствола выполнено в программном комплексе Matlab с использованием решения (3).

Рис. 5. Моделирование траектории скважины с синусоидальным профилем Анализ результатов исследования показывает, что изменение кривизны траектории между точками 1 и 3 незначительно влияет на производительность сложнопрофильной скважины и в большинстве случаев отклонением траектории ствола от прямолинейного приближения можно пренебречь.

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

Проблеме моделирования продуктивности ГС с учетом гидравлических потерь давления по стволу, открытому притоку, посвящены труды З.С. Алиева, В.А. Черных, B.J. Dikken, R.A. Novy, K. Seines, K. Aziz, A.D. Hill, V.R.

Penmatcha, E. Ozkan и других исследователей.

Моделирование производительности ГС с учетом изменения профиля давления вдоль рабочего участка ствола основывается на взаимосвязи процессов течения флюида в пористой среде и стволе скважины, и сводится к решению соответствующей системы дифференциальных уравнений. В настоящее время широкое распространение получили численные методы решения рассматриваемой задачи. Однако численное решение исходной системы дифференциальных уравнений требует проведения многократных итераций - в большинстве случаев это приводит к неприемлемому росту времени машинного расчета. Существующие аналитические решения задачи течения жидкости в системе лпласт - горизонтальный ствол получены для модели ГС бесконечной длины и в большинстве случаев неприменимы для практических расчетов (B.J.Dikken, R.A. Novy). В связи с этим в работе предлагается полуэмпирический подход к решению рассматриваемой задачи, в основу которого положена обработка результатов экспериментов, проведенных на численной модели системы лпласт - горизонтальный ствол.

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

С использованием разработанной численной модели была проведена обширная серия численных экспериментов. В процессе расчета в широком диапазоне варьировались проницаемость пласта, прилагаемая депрессия Р, длина горизонтального участка, вязкость и плотность флюида, внутренний диаметр обсадной колонны участка скважины, открытого притоку. При каждом фиксированном значении параметра рассчитывались: истинный перепад давления по горизонтальному стволу Рз и его аппроксимация Рз*, производительность скважины с учетом (Q) и без учета (Qin) потерь давления в стволе. Для приближенной оценки величины перепада давления в горизонтальном стволе использовался параметр Рз*, рассчитанный по формуле Дарси-Вейсбаха для труб с непроницаемыми стенками, в которой в качестве расхода бралась половина величины дебита скважины, определенного без учета потерь давления. Впервые такой подход к оценке гидравлических потерь в открытом стволе был предложен в работах Hill et al.

Анализ полученных результатов показал, что снижение производительности ГС вследствие гидравлических потерь в стволе является функцией безразмерного параметра Рз /Р (рис. 6). Обработка результатов численных экспериментов позволила установить однозначную полиномиальную связь между безразмерными параметрами Q/Qin, Pз/P и Pз*/P.

Рис. 6. Зависимость относительного снижения дебита ГС от безразмерного перепада давления в горизонтальном стволе Предлагаемый алгоритм расчета производительности ГС состоит из следующих шагов:

1. Для заданной величины депрессии P рассчитывается дебит ГС Qin без учета скважинной гидравлики.

2. По уравнению Дарси-Вейсбаха определяется значение параметра Pз*:

2QinL Pз* , (8) d где - коэффициент гидравлического сопротивления; - плотность флюида; L - длина скважины, d - внутренний диаметр скважины.

3. Рассчитывается величина отношения потерь давления в стволе к величине депрессии на пласт:

Pз P 0,000314ln 100Pз* / Р 0,0111ln 100P* / Р . (9) з 0,11342ln 100Pз* / Р 0,2859ln 100P* / Р 0,276 з 4. Вычисляется производительность скважины с учетом потерь давления на горизонтальном участке:

Q Qin 0,14 Р / Р 0,6435 Р / Р 0,9987 . (10) зз Расхождение между результатами расчетов по предлагаемой методике и численными решениями задачи в рассматриваемом диапазоне исходных параметров не превышает 2,5%.

Проведенные исследования также показали, что для предварительной оценки влияния потерь давления в стволе на производительность ГС может быть использован параметр Pз*/P. Установлено, что если Pз*/P < 0,07, снижение производительности по отношению к модели, не учитывающей гидравлические потери, составляет менее 5%.

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

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

На рисунке 7 приведена зависимость параметра , представляющего собой отношение дебитов синусоидальной и горизонтальной скважин равной длины, от h и kH/kV. Из рис. 7 следует вывод о том, что наиболее привлекательными кандидатами для бурения синусоидальных скважин являются пласты значительной мощности с высоким контрастом проницаемостей в горизонтальном и вертикальном направлениях. Например, применение синусоидальных скважин в пластах с h > 40 м и kH/kV > 60 (Rk/L = 1,5) даст более чем 20% преимущество в дебите таких скважин по отношению к ГС той же длины.

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

Рис. 7. Эффективность применения скважин с синусоидальным профилем ствола: Rk = 1500 м, L = 1000 м, rc = 0,1 м Предлагаемый подход к решению данной задачи заключается в расчете чистого дисконтированного дохода (NPV) как функции пускового дебита и динамики его изменения во времени, стоимости строительства скважины, налоговых отчислений и эксплуатационных затрат. Задача сводится к нахождению таких параметров траектории и длины ствола скважины, которые обеспечивают максимум функции NPV:

(11) NPV , L,n max где n - количество наклонных участков упрощенной синусоидальной траектории.

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

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ 1. Получены приближенные аналитические решения задачи притока жидкости к скважине с прямолинейной траекторией ствола, произвольно ориентированной относительно кровли и подошвы пласта. Показано, что в вертикально-анизотропной среде уменьшение угла наклона скважины прямолинейного профиля к вертикали может приводить к значительному приросту дебита.

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

3. Разработана аналитическая методика оценки продуктивности протяженных нефтяных скважин с учетом потерь давления в горизонтальном стволе. Выявлены критерии значимости влияния гидравлических потерь в открытом притоку горизонтальном стволе на производительность нефтяной ГС. Показано, что если отношение перепада давления по длине рабочего участка скважины к величине депрессии на пласт Pз*/P составляет менее 0,07, снижение производительности ГС вследствие гидравлических потерь давления на трение не превышает 5%.

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

Список опубликованных работ по теме диссертации 1. Доманюк Ф.Н., Золотухин А.Б. Определение дебита скважины с прямолинейным профилем в вертикально-анизотропном пласте // Нефтяное хозяйство, 2011, №5, с.92-95.

2. Доманюк Ф.Н. Стационарный приток жидкости к скважине с волнообразным профилем // Нефтепромысловое дело, 2011, №7, с. 21-26.

3. Доманюк Ф.Н. Моделирование продуктивности нефтяных скважин со сложной траекторией горизонтального ствола // Труды РГУ нефти и газа имени И.М. Губкина, 2011, №3, с.37-47.

4. Доманюк Ф.Н. Исследование стационарного притока жидкости к скважинам сложного профиля. Тезисы докладов 2-ой научно-практической молодежной конференции Новые технологии в газовой отрасли: опыт и преемственность, Москва, Газпром ВНИИГАЗ, 2010, с.14.

5. Доманюк Ф.Н. Производительность горизонтальной скважины со сложным профилем в анизотропном пласте // Тезисы докладов IV научно-практической конференции УМатематическое моделирование и компьютерные технологии в процессах разработки месторождений, добычи и переработки нефтиФ. - М.: ЗАО Издательство Нефтяное хозяйство, 2011. - с.30.

6. Доманюк Ф.Н. Оценка эффективности применения змеевидных скважин при разработке залежей УВ // Материалы 2-го Научнотехнического семинара УАктуальные вопросы проектирования разработки месторождений углеводородовФ. - М.: ООО ВНИИГАЗ, 2011.

- с. 22.

Соискатель: Доманюк Ф.Н.

e-mail: fedor.domanyuk@gmail.com Подписано в печать 31 января 2012 г.

Объем 1,2 п.л.

Тираж 100 экз.

Заказ № Отпечатано в Центре оперативной полиграфии ООО Ол Би Принт Москва, Ленинский пр-т, д.    Авторефераты по всем темам  >>  Авторефераты по земле