Разработка и совершенствование технологических решений по повышению эксплуатационных показателей горизонтальных скважин и боковых горизонтальных стволов

Размещено на http://www.allbest.ru/

АВТОРЕФЕРАТ

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

кандидата технических наук

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

Разработка и совершенствование технологических решений по повышению эксплуатационных показателей горизонтальных скважин и боковых горизонтальных стволов

ГИЛАЕВ РУСЛАН ГАНИЕВИЧ

Краснодар - 2009

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

Актуальность проблемы

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

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

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

Цель работы

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

Задачи исследования

1. Экспериментальные и теоретические исследования гидравлических

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

2. Оценка потерь давления на трение и смешивание периферийного притока с основным потоком в горизонтальном стволе скважины.

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

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

Научная новизна

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

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

3. Аналитически установлено, что гидравлические сопротивления в пористых каналах, работающих на приток (добывающие скважины) выше в 1, 5-2, 6 раза, чем для непроницаемых каналов и каналов, работающих на нагнетание (нагнетательные скважины).

Практическая ценность работы

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

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

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

4. Разработаны методические указания для студентов по выполнению практических заданий по дисциплине «Горизонтальное бурение» для нефтегазовых специальностей 130501, 130503, 130504 и 130602.

5. Результаты диссертационной работы в виде методик вошли в научно-технические отчеты «Разработка и совершенствование методик расчета интенсивности искривления ствола скважины и плотности перфорации горизонтальных скважин и боковых стволов» (2006-2008 гг.).

Апробация работы

Материалы диссертационной работы докладывались на 4-й Международной конференции «Освоение и добыча трудноизвлекаемых и высоковязких нефтей (г. Анапа, Краснодарский край, 2004 г.), на научно-техниче-ских советах ДОАО «Нефтетерммаш» (2005 г.), ООО «Кубаньгазпром» (2006 г.), на кафедре прикладной математики Кубанского государственного университета (2008 г.), на совместном заседании кафедр «Нефтегазового промысла» и «Оборудование нефтяных и газовых промыслов» Кубанского государственного технологического университета (2008 г.).

Публикации результатов работы

Основное содержание диссертационной работы отражено в 5 печатных работах, в том числе 2 статьи в изданиях, рекомендованных ВАК РФ.

Структура и объем работы

Диссертационная работа состоит из введения, четырех разделов, основных выводов и рекомендаций, списка использованных источников (97 наименований). Работа изложена на 98 стр. машинописного текста, содержит 14 таблиц и 18 рисунков.

Автор выражает благодарность за ценные консультации и полезные

советы при обсуждении результатов работы: д.т.н. Гарушеву А.Р., д.т.н. Кошелеву А.Т., сотрудникам и аспирантам кафедры НГП КубГТУ.

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

скважина ствол фильтрационный гидродинамический

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

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

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

В 1856 году А. Дарси установил основной закон фильтрации жидкости, а в 1858 году Ж. Дюпюи рассмотрел задачу водопритока к вертикальной скважине (колодцу), заложив своими работами основы подземной гидромеханики. В 1886 году Ф. Форхгеймером решена задача о притоке воды к горизонтальному рву, аналогу горизонтальной скважины, а в 1908 году им же для бесконечного ряда колодцев, лежащих на расстоянии «b» от контура питания, получена формула изменения давления в пласте

, (1)

где Q - суммарный дебит скважины; а - половина расстояния между скважинами; b - расстояние до контура питания; k, h, - соответственно проницаемость, мощность пласта, вязкость жидкости; c - постоянная интегрирования.

Уравнение (1) получено методом наложения стоков и источников для однородного изотропного пласта и удовлетворяет уравнению Лапласа.

Дальнейшее развитие подземной гидравлики и нефтяной подземной гидромеханики связано с именами Л.С. Лейбензона, Н.Н. Павловского, Г.Н. Каменского, М. Маскета, П.Я. Полубариновой-Кочиной, В.Н. Щелкачева, И.А Чарного, Ю.П.Борисова, Б.Б Лапука и др.

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

Развитие теории и практики наклонного и кустового наклонного бурения шло высокими темпами. В 1935 году под руководством инж. Н.С. Тимофеева и К.Ф. Михайлова пробурена наклонная скважина под дно Каспийского моря, а уже в 1947 году под руководством инж. А.М. Григоряна и В.А. Брагина пробурили на Краснокамском месторождении скважину № 240 с двумя горизонтальными ответвлениями длиной по 30 м. В Куйбышевской области в 1957 году при глубине скважины 530 м забурен горизонтальный ствол длиной 148 м.

Примерно с 1948 по 1960 гг. начались специальные исследования по разработке нефтяных месторождений с помощью ГС и многозабойных скважин. Крупные гидродинамические исследования проводились во ВНИИнефти.

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

Приведем некоторые из них. В работах В.Г. Григулецкого и Б.А. Никитина предложена обобщенная формула, построенная на модели Ю.П. Борисова, и проводится ее сравнение с формулами С. Джоши и Ренарда - Дюпюи.

Расчеты выполнены для длины ГС от 100 до 250 м. Показано, что формулы Борисова и Джоши дают примерно одинаковые значения дебитов ГС.

А.И. Ибрагимов и А.А. Некрасов сравнивают рекомендованную ими формулу с формулами В.П. Меркулова, Ю.П. Борисова, С. Джоши и Ренарда. Численные эксперименты и результаты сравнения показывают, что при длине горизонтального ствола м формулы С. Джоши и Ю.П. Борисова дают одни и те же результаты. По формуле В.П. Меркулова получаются несколько заниженные результаты, а Ренарда - очень завышенные. В работе А.И. Ибрагимова и А.А. Некрасова рекомендованы формулы притока для оценки продуктивности ГС, расположенной в середине трехмерного ограниченного пласта, а также для ГС, расположенных в одной плоскости под углом друг к другу и для параллельных стволов в середине пласта. В первом случае формула сравнивалась с указанными выше, а на фактических данных она не оценивалась.

В большой и содержательной работе В.Л. Богданова, Ю.П. Батурина и др. дан анализ результатов бурения и эксплуатации ГС на Федоровском месторождении. Здесь приведен статистический материал по 174 ГС, даны рекомендации по повышению эффективности эксплуатации месторождения ГС и по комплексу геофизических и гидродинамических исследований ГС.

А. Бронзов, Н. Бадовский и др. рассмотрели статистику предыдущей работы по Федоровскому месторождению, построили графики зависимости дебитов ГС от длины горизонтального участка, показали, что для этого месторождения «оптимальная» длина ГС составляет 100-200 м.

В работах С.К. Сохошко, А.П. Телкова, С.И. Грачева и др. рассмотрены вопросы притока к пологой скважине, предложены модели расчета изменения перепада давления с учетом потерь напора между сечениями, внутри которых заключено перфорационное отверстие, потерь напора в самом отверстии, и их взаимное влияние (интерференция) друг на друга. Движение жидкости между начальной и конечной точкой пологой скважины описано уравнением Бернулли. Полученные результаты представлены графиком изменения расхода через перфорационные отверстия и депрессии по стволу скважины для конкретных значений параметров продуктивного пласта, вязкости нефти и др. Однако гидравлические сопротивления здесь не учтены.

В работе В.А. Иктисанова и Д.Г. Яраханова рассмотрен экономический критерий зависимости дебита и стоимости строительства скважины от длины горизонтального ствола. В работе дана ссылка на уравнение притока Батлера - Цупруновича, по которому рассчитывается приток. Расчетные профили притока по длине ствола соответствуют равномерному. Высокая стоимость бурения горизонтального ствола снижает чисто дисконтированный доход (ЧДД), а, следовательно, и длину ГС.

Т.О. Акбулатов и др. рассмотрели всю зону течения жидкости к ГС в виде двух областей: радиального и параллельного течения. Рассмотрен нестационарный приток и граница зоны возмущения «сшивается» между областями приравниванием расхода. Это не совсем верно, так как в точке сшивания должны быть приравнены давления и первые производные давления. Принципиально данная схема пригодна для исследования ГС и определения скин-эффекта.

Одной из причин низкой эффективности горизонтальных скважин Ф.А. Агзамов, Т.О. Акбулатов и др. считают кольматацию стенок скважины твердой фазой буровых растворов. Авторы данной статьи рекомендуют для расчетов свою формулу. Сравнивая формулы притока к галерее с формулами З.С. Алиева и В.В. Шеремета для ГС, они получают отношение дебитов 0, 98. Отсюда следует вывод, что формула З.С. Алиева дает дебит для галереи скважин. Далее сравнивается формула Ю.П. Борисова с единичным дебитом галереи скважин. Отношение дебитов составляет 3, 3, что означает, что формула Ю.П. Борисова дает 3-х кратное завышение дебита галереи. Предложенная авторами формула дает дебит, который составляет 0, 6-0, 9 от дебита галереи.

К.О. Каширина в своих расчетах по дебитам ГС единичной длины дает прямо противоположные ответы. Расчеты дебитов ГС единичной длины дают следующие значения в м³/сут: С. Джоши - 0, 23; К.О. Каширина - 0, 22; Ренард - 0, 21; Ю.П. Борисов - 0, 20; Giger - 0, 20; В.П. Пилатовский - 0, 17; В.Д. Лысенко - 0, 12 и В.С. Евченко - 0, 09.

В работах Р.В. Стельмака со ссылками на С. Економидеса, А. Хилла, С. Арбади, В. Диккена вводится индекс удельной производительности и ошибка «производительности». Индекс производительности изменяется вдоль ствола в соответствии с изменением плотности перфорации, проницаемости пласта и т.д., а ошибка «производительности» представляет собой отношение потерь давления на ГУ к депрессии на забое. Построены кривые производительности ГС для различных депрессий и даны ошибки производительности скважин в зависимости от длины ГУ. Для снижения гидравлических сопротивлений на ГУ скважины в США устанавливают «EQUALIZER», снижающие скорость основного потока и выравнивающие скорость притока по длине горизонтального участка ГС.

Используя статистические данные по ГС Федоровского месторождения, работы А.С. Бронзова, Р.В. Стельмака, З.С. Алиева, А.И. Ибрагимова, нами построена фактическая кривая зависимости дебита ГС от длины ГУ. На рисунке 1 приводится фактическая кривая (1) и диапазон рассеивания расчетных данных по формулам Ю.П. Борисова, В.П. Меркулова, С. Джоши, Ренарда.

Рисунок 1 - Фактическая кривая и диапазон рассеивания расчетных данных.

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

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

, (2)

где Q - дебит ГС, тонн/сут; L - длина горизонтального участка, м.

Среднеквадратичное отклонение составляет всего тонн/сут. Кривая (2) имеет асимптоту при тонн/сут.

Следовательно, для «коротких» ГС с длиной ГУ менее 250 м все рассматриваемые модели дают примерено равные результаты. Это свидетельствует о том, что:

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

модели Ю.П. Борисова, А.И. Чарного и С. Джоши наиболее пригодны для расчета дебита ГС, т.к. результаты расчетов для «коротких» () горизонтальных стволов по ним практически совпадают, а остальные модели являются их модификацией;

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

Анализ работ, выполненных в диссертации, показал, что для совершенствования и повышения эффективности ГС необходимо сосредоточить внимание на следующих направлениях:

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

оценка потерь давления на трение и смешивание периферийного притока с основным потоком в горизонтальной скважине (ГС);

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

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

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

Вопросам гидродинамики и тепломассообмена пористых каналов и труб в последние годы уделяется значительное внимание. Большое место здесь уделяется экспериментальным исследованиям при принудительном притоке (оттоке) жидкости через пористые стенки. Специфика решаемых задач по гидродинамике и тепломассообмену (охлаждение ракетных двигателей, электрокабелей, бортовых приборов и др.) требовала принудительного подвода и отвода больших масс охладителей в пористых трубах. Поэтому основные экспериментальные исследования проводились на установках небольшого размера с высокими скоростями подвода и отвода охлаждающего агента. В работах А.С. Бермана, Р.М. Тейлора, С.В. Юаня, Г.Л. Вайсберга, С.М. Тарга, И.П. Гинзбурга и др. рассмотрены различные постановки и методы решения этих задач с граничными условиями, учитывающими подвод или отвод жидкости (газа) на поверхности. При этом течение жидкости в пористых каналах описывается уравнением Навье - Стокса в приближении теории пограничного слоя с учетом инерционных или вязких составляющих. Показано, что течение жидкости в пористых трубах полностью определяются двумя критериями подобия. Это критерии Рейнольдса для течения внутри канала () и периферийного притока (оттока), который отнесен к диаметру трубы (). Здесь U и V - продольная и поперечная скорости жидкости (м/с); Y - коэффициент кинематической вязкости (м²/с); d - диаметр трубы (м).

Используя вышеизложенное и учитывая гидродинамические особенности горизонтальных скважин, нами определены диапазоны изменения геометрических и динамических критериев по следующим исходным данным: длина горизонтального ствола м; диаметр (ширина) горизонтального участка м и м; дебит скважины 1000 м³/сут; кинематическая вязкость нефти м²/c.

Тогда для добывающих скважин (притока) получаем: максимальная и минимальная продольная скорость м/с, м/с; максимальная и минимальная поперечная скорость притока м/с, м/с; числа Рейнольдса , ; число Эйлера ; отношение длины к диаметру .

По расчетам примерно те же порядки R_e, R_n, L/d и E соответствуют нагнетательным скважинам. Видно, что режим течения жидкости на горизонтальном участке охватывает широкий диапазон, а превосходит размеры экспериментальных установок в упомянутых работах.

Далее в работе дается описание экспериментальной установки, гидродинамически подобной горизонтальным скважинам. Установка принципиально не отличается от установки Ф.А. Шавелева по исследованию гидравлики трубного течения. В качестве рабочего элемента применялись перфорированные полиэтиленовые трубки длиной 40 м и диаметром 4 см с перфорационными отверстиями диаметром 1 и 2 мм. Измерялся расход (см³/с) в начале трубки Q_А и конечный (транзитный) расход Q_Т. Разность составляла расход через 40 отверстий, расположенных через 10 см.

Потери напора для перфорированной трубки определялись по формуле , а для неперфорированной трубки - (где А₁ и А₂ - удельное сопротивление трубок; - длина трубки; Q_А и Q₀ - соответствующие расходы). Отношение , т.е. равно отношению коэффициентов гидравлических сопротивлений. В таблице 1 приведены значения начального и конечного расхода и отношение .

Таблица 1 - Значения расхода и отношения коэффициентов гидравлических сопротивлений

№ n/n	QА	Q0	QТ		№ n/n	QА	Q0	QТ
1	10	8	2	0, 358	9	73	57	16	0, 422
2	17	14	3	0, 401	10	76	57	19	0, 437
3	22	17	5	0, 436	11	78	56	22	0, 454
4	37	25	12	0, 476	12	84	57	27	0, 475
5	48	35	13	0, 448	13	89	58	31	0, 499
6	53	40	13	0, 435	14	96	63	33	0, 487
7	63	48	15	0, 431	15	98	64	34	0, 489
8	74	58	16	0, 421	16	108	71	37	0, 487

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

Для установления зависимости скорости оттока жидкости от диаметра перфорационных отверстий по длине трубки через каждые 5 м были установлены микроманометры. Полученные результаты подтвердили предлагаемые в работах В.С. Генкина, В.В. Дильмана, В.Ф. Климкина и А.Д. Ренкина формулы для расчета скоростей принудительного вдува или оттока для коротких трубок. Для малых перепадов давления зависимость скорости линейна от Р, а для высоких скоростей зависит от .

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

, (3)

изменение давления

, (4)

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

Для добывающих скважин в формулах (3) и (4) достаточно подставить (- ) с учетом того, что .

Гидравлические сопротивления равны:

, (5)

где Р - перепад давления по длине канала; - плотность жидкости; U - средняя скорость течения жидкости в канале.

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

для нагнетательных скважин:

, (6)

для добывающих скважин:

, (7)

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

Рисунок 2 - Изменение коэффициента гидравлических сопротивлений в зависимости от kL

Графики построены для величин k и L в диапазоне 1/м, м. На рисунке 2 для добывающих (кривая 1) и нагнетательных скважин (кривая 2) относятся к ламинарным потокам, а кривые 1 и 2 - к квадратичному течению. Прямая 3 соответствует потоку в непроницаемом канале, а пунктирная прямая 4 соответствует равномерному оттоку и притоку. Видно, что при постоянном k (например, 1/м) увеличение L от 0 до 500 м приводит к увеличению гидравлических сопротивлений для добывающих скважин и трещин гидроразрыва в 1, 5 раза для ламинарного потока и в 2, 6 раз для квадратичного течения. Эти изменения (снижение гидравлических сопротивлений) для нагнетательных скважин составляют 1, 3 и 1, 6 раз соответственно.

Полученные результаты показывают, что

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

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

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

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

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

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

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

Наши исследования посвящены разработке методики расчёта таких профилей, лишенной перечисленных недостатков.

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

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

1 - глубина по вертикали Н; Н_в - глубина зарезки наклонного ствола; А - проектное смещение;

2 - Н, Н_в, А и R - радиус искривления ствола скважины на интервале набора зенитного угла;

3 - А, R и - зенитный угол в конце интервала его набора.

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

На все варианты расчёта составлены программы в виде электронных таблиц (в среде EXCEL).

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

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

1. Точность расчёта и построения фактической трассы скважины недостаточно высокая;

2. Не полностью используется информация, полученная в результате расчёта фактической трассы скважины;

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

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

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

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

Ранее отмечалось, что распределение давления и расхода по длине ГУ горизонтальной скважины носит нелинейный характер. Перепады давления внутри и снаружи в системе «скважина - пласт» приводят к неравномерности притока (оттока) флюида по длине пористого (перфорированного) канала.

Как показали наши экспериментальные исследования на перфорированных трубках при течении жидкости, скорость оттока (расход) на начальных участках выше, чем в конце. Промысловые исследования, выполненные совместно с ООО «Кубаньгазпром» на наклонно-горизонтальной скважине № 102 Кущевская показали, что для добывающих скважин приток неоднороден по длине ГУ. Основная зона притока (1/3 перфорированной длины) приходится на начало скважины. Исследования, проведенные зарубежными специалистами (Norsk Hydro и Baker Oil Tools) на ГС в Северном море с длиной ГУ от 350 до 3500 м показали, что приток не распределяется равномерно по длине ГУ и появляется опасность преждевременного прорыва воды (или газа). Для однородных пластов место прорыва меняется в зависимости от проницаемости по длине ГУ.

В случае однородного по проницаемости пласта равномерность притока по длине ГУ можно получить неравномерным распределением плотности перфорации (пористости трубы) по длине ствола. Здесь под «пористостью» трубы понимается отношение суммарной площади отверстий на единицу длины трубы к ее боковой поверхности. При этом плотность перфорации («пористость») необходимо увеличивать от начала ГУ к концу ствола. Схема расчета приведена на рисунке 3.

Здесь х_k, х_k-1, … х₁ середины интервалов, на которые разбивается вся длина ГУ. Число интервалов k выбирается от 5 до 10 в зависимости от длины труб хвостовика. Длина каждого интервала равна .

В основу расчета плотности перфорации положена известная формула

Рисунок 4 - Схема расчета плотности перфорации

, (8)

где i - номер интервала, i = 1, 2, …, k; x_i и x_i+1 - расстояние от начала ГУ до середины интервала; N_i и N_i+1 - плоскость перфорации на предыдущем и последующем интервалах от конца скважины.

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

Очевидно, конструкторские и технологические соображения по проектированию и изготовлению «пористых» труб вносят коррективы в результаты расчета. Так, с позиций изготовления удобно располагать отверстия в ряду через 60^о, 45^о или 40^о, т.е. по 6, 8 или 9 отверстий в ряду, а диаметры отверстий d = 5; 6; 8; 10; 12; 14; 20 мм. Все расчеты в работе выполнены для длины трубы хвостовика, равной 10 м.

Предлагаемый алгоритм и программа расчета дизайна фильтра-хвостовика применима для случая переменных по длине ГУ фильтрационных характеристик коллектора. В этом случае расчет ведут для средней проницаемости по всей длине ГУ, а затем пересчитывают плотность перфорации из соотношения (где N_i - плотность перфорации по расчету; - средняя проницаемость по всей длине ГУ; k_i - средняя проницаемость на i-том участке).

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

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

ОСНОВНЫЕ ВЫВОДЫ И РЕКОМЕНДАЦИИ

1. Анализ расчетных моделей дебита горизонтальных скважин показал, что:

ни для одной модели нет статистических оценок (среднеквадратическое отклонение, мера адекватности и др.) между расчетными и фактическими данными;

наиболее предпочтительными для небольших длин ГУ ( м) являются модели Ю.П. Борисова, А.И. Чарного и С. Джоши, которые дают примерно равные результаты, а остальные модели являются их модификацией;

недоучет гидравлических сопротивлений по длине ГУ приводит к несоответствию расчетных и фактических дебитов для скважин с протяженностью ГУ 250 м и более;

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

2. Экспериментальные исследования на полиэтиленовых трубках, гидродинамически подобных ГС, показали, что:

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

при небольших скоростях течения жидкости (расходах) внутри перфорированного канала скорость оттока (притока) жидкости пропорциональна перепаду давления внутри и снаружи канала Р, а при больших скоростях (расходах) основного потока скорость оттока (притока) пропорциональна .

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

4. Расчетными методами установлено, что гидравлические сопротивления в пористых каналах, работающих на приток (добывающие ГC) в 1, 5-2, 6 раза выше, чем для непроницаемых каналов и каналов, работающих на отток (нагнетательные ГС).

5. Для обеспечения равномерного притока (оттока) жидкости по длине ГУ для однородных по проницаемости пластов плотность перфорации («пористость») необходимо увеличивать от начала ГУ к концу ствола.

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

7. Разработана методика, алгоритм и программа расчета дизайна «фильтра-хвостовика», которые применимы для однородных и неоднородных по проницаемости пластов. Алгоритм и программа могут быть также использованы для расстановки проволочных фильтров с различной скважностью.

8. По результатам работы разработаны и внедрены в учебный процесс методические указания по дисциплине «Горизонтальное бурение» для студентов нефтегазовых специальностей 130501, 130503, 130504 и 130602.

ЛИТЕРАТУРА

1. Кузнецов М.В., Окс С.В., Гилаев Р.Г. ДОАО МЗ «Нефтетерммаш» - надежный партнер в области поставки нефтепромыслового оборудования // Сб. докладов 4-й Междунар. конференции «Освоение и добыча трудноизвлекаемых и высоковязких нефтей». - Краснодар: Изд-во ЭТВИ, 2004. - С. 445-446.

2. Вартумян Г.Т., Гнездов Г.В., Гилаев Р.Г. Вероятностный подход образования трещин гидроразрыва пласта // Инженер нефтяник, 2009. - № 8. - С. 15-18.

3. Гилаев Р.Г. Гидродинамические особенности горизонтальных скважин и трещин гидроразрыва / Г.Т. Вартумян, Р.Г. Гилаев, А.Т. Кошелев, А.В. Гнездов, С.В.Смык // Строительство нефтяных и газовых скважин на суше и на море. - М.: ОАО «ВНИИОЭНГ», 2009. - № 8 - С. 15-18.

4. Гилаев Р.Г. Гидродинамические сопротивления пористых каналов большой протяженности / Строительство нефтяных и газовых скважин на суше и на море. - М.: ОАО «ВНИИОЭНГ», 2009 - № 8 - С. 18-20.

5. Арутюнов А.А., Гилаев Р.Г. Горизонтальное бурение: метод. указания по выполнению практических заданий по дисциплине для студентов нефтегазовых специальностей 130501, 130503, 130504 и 130602

Размещено на Allbest.ru

...