Бесплатный автореферат и диссертация по наукам о земле на тему

Гидродинамика двухфазных смесей в процессах бурения нефтяных и газовых скважин

ВАК РФ 25.00.15, Технология бурения и освоения скважин

Автореферат диссертации по теме "Гидродинамика двухфазных смесей в процессах бурения нефтяных и газовых скважин"

На правах рукописи УДК 622.24

ИСАЕВ ВАЛЕРИЙ ИВАНОВИЧ

ГИДРОДИНАМИКА ДВУХФАЗНЫХ СМЕСЕЙ В ПРОЦЕССАХ БУРЕНИЯ НЕФТЯНЫХ И ГАЗОВЫХ СКВАЖИН

Специальность 25.00.15 - Технология бурения и освоение скважин

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

- 3 ДЕК 2009

Москва-2009

003487056

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

Научный консультант доктор технических наук,

профессор Леонов Е.Г.

Официальные оппоненты:

доктор технических наук, главный научный сотрудник Шеберстов Е.В.

доктор технических наук, профессор Войтенко В.С.

доктор технических наук, профессор Гусман А.М.

Ведущая организация: Буровая компания ОАО

«Газпром» ООО «Бургаз»

Защита состоится 08 декабря 2009 г. в 15.00 часов в ауд. 731 на заседании Диссертационного Совета Д.212.200.08 при Российском государственном университете нефти и газа имени И.М. Губкина по адресу: 119991, ГСП-1, Москва, В-296, Ленинский проспект, д.65.

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

Автореферат разослан " 06 " _2009 г.

Ученый секретарь диссертационного совета доктор технических наук, профессор

Б.Е. Сомов

Общая характеристика работы

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

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

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

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

Цели диссертационной работы

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

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

- совершенствование одномерных моделей течения двухфазных смесей в элементах циркуляционной системы скважины;

- экспериментальное исследование процесса истечения газовых струй через слой НЖ и ННЖ;

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

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

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

2. Обобщены законы гидростатики ньютоновских, неньютоновских жидкостей и многофазных смесей из них.

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

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

5. Дано обоснование перехода от ламинарного течения к турбулентному при движении вязкопластической жидкости (ВПЖ) в трубах.

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

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

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

Достоверность полученных результатов

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

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

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

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

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

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

Результаты работы докладывались, демонстрировались и одобрены на всесоюзных, всероссийских и международных съездах, конференциях и симпозиумах. Автор докладывал основные разделы диссертации на научных семинарах: в институте механики МГУ им. М.В. Ломоносова; в институте проблем нефти и газа РАН; на кафедре бурения нефтяных и газовых скважин и на кафедре нефтегазовой и подземной гидромеханики РГУ нефти и газа имени И.М. Губкина, "Актуальные проблемы нефтегазовой гидромеханики и разработки месторождений углеводородов" РГУ нефти и газа имени И.М. Губкина.

Публикации

Основные материалы диссертации опубликованы в 43 печатных работах; из них 19 статей вышли в журналах, рекомендованных ВАК для публикации материалов по докторским диссертациям. Два учебника выпущены издательством "Недра"; один учебник, четыре учебных пособия и статьи напечатаны в других издательствах. Всего автором опубликовано 126 работ.

Основные защищаемые положения

1. Модель двухфазной гидродинамики процессов бурения.

2. Обобщение закона гидростатики для ньютоновских и неньютоновских растворов и их многофазных смесей.

3. Методика управления скважиной при бурении на депрессии.

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

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

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

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

8. Формула для расчёта коэффициента гидравлических сопротивлений при турбулентном течении вязкопластического раствора в трубах.

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

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

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

Благодарности

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

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

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

Содержание работы

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

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

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

Бакланова и Р.И. Шищенко, Д. Верслюиса, Т.Ф. Мура, Г.Д. Уайльда, А.П. Крылова, A.A. Арманда, Т. Поэттмана, П. Карпентера и многих других ученых учитывались зависимости для истинного содержания фаз.

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

Возрастающее внимание к двухфазным средам способствовало изданию монографий и учебников теоретического и прикладного характера для нефтегазовой направленности как в России, так и за рубежом следующих авторов: В.Г. Багдасаров (1947); В.А. Архангельский (1958); С.С. Кутателад-зе, М.А. Стырикович (1958); К.В. Виноградов (1964); М.А. Гейман и В.И. Мусинов (1965); Г.Ф. Агаев (1966); Д.Ф, Файзулдаев (1966); А.О. Межлумов и Н.С. Макурин (1967); С. Л. Coy (1967); В.А. Мамаев, Г.Э. Одишария, Н.И. Семёнов и A.A. Точигин (1969); Г.Б. Уоллис (1969); М.С. Винарский и Н.М. Гончаренко (1969); Д.Ф. Файзуллаев, P.C. Гурбанов и Я.М. Расизаде (1970); Дж. Хьюитт и Н. Холл-Тейлор (1970); В.А. Амиян и Н.П. Васильева (1972); Г.В. Циклаури, B.C. Данилин (1973); А.О. Межлумов (1976); Д. Баттерворс и Г. Хьюитт (1977); Р.И. Нигматулин (1978, 1987); В.А. Мамаев, Г.Э. Одишария, О.В. Клапчук и др. (1978); Д.Ф. Файзуллаев, А.И. Умаров и A.A. Шакиров (1980); В.Н. Николаевский (1984); А.Х. Мирзаджанзаде А.Х. и В.М. Ентов (1985); Е.Г. Леонов и В.И. Исаев (1987); H.A. Гукасов (1988); В.Д Малеванский и Е.В. Шеберстов (1990); А.И. Булатов, А.Г. Аветисов (1993-1996); А.И. Гриценко, О.В. Клапчук и Ю.А. Харченко (1994); А.А.Точигин и Г.Э. Одишария (1998); В.И. Ямпольский. (1999); Дж.П. Брил и X. Мукерджи (1999); Ю.М. Басарыгин, А.И. Булатов, Ю.М. Проселков (2000); В.А. Сахаров и М.А. Мохов (2004); Л.Н.Полянин и В.П. Дробков (2004) и др.

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

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

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

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

а) Схема ЦС: 1 - кольцевое пространство (КП); 2 - бурильные трубы (БТ); 3 - утяжелённые бурильные трубы (УБТ); 4 - забойный двигатель; 5 - долото; 6 - замок; 7 - обсадная колонна; 8 - открытый (необсаженный) ствол; 9 - перекрытый слабый пласт;

10 - горная порода под башмаком последней спущенной обсадной колонны;

11 - вскрываемый пласт.

б) Распределение давлений в элементах ЦС (I - гидростатическое, II - при циркуляции в КП, III - при циркуляции в БК): 1-2, 3-4, 4-7 - за БТ; 2-3 - за замками; 7-8 - за УБТ; 8-9 - за двигателем; 9-12 - в долоте; 12-13 - в забойном двигателе; 13-14 - в УБТ; 14-15,16-17-внутри БТ; 15-16-взамках.

Значения давлений: 5, 10 - пластовые давления рпл1, рга2; 6, 11 - давления гидроразрыва (поглощения) ррь рр2 в горной породе и нижнем вскрываемом пласте; 18 -забойное гидростатическое давление; 9 - забойное давление при циркуляции (промывке); 19 - гидростатическое давление в КП под башмаком обсадной колонны; 20 - давление в КП при промывке под башмаком; 21 - давление в стояке; 1 - давление в КП на устье.

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

распределения параметров: 1) расходов фаз; 2) давлений; 3) плотностей; 4) напряжений; 5) концентраций фаз; 6) температур; 7) геометрических размеров элементов ЦС (длина, диаметр и расположение в пространстве, глубина расположения, радиус и толщина пластов); 8) характеристик компрессоров и насосов, цементировочных агрегатов и смесительных машин (подача, давление); 9) прочностных характеристик элементов системы; 10) характеристик подъемного механизма буровой установки (скорость и ускорение при спускоподъемных операциях); 11) характеристик забойных двигателей (перепад давления при различных расходах промывочной смеси); 12) гранулометрического состава выносимого из скважины шлама.

Распределения п.п. 1-6 связаны друг с другом общими уравнениями гидродинамики в области распределений п.п. 7-12, существующих при бурении. Описание гидродинамических процессов бурения сводится к нахождению соотношений, связывающих распределения п.п. 1-12.

На рис. 2 приведен перечень основных процессов 1.1-1.3 и 2.1-2.5 и связанных с ними задач 1.1.1-1.3.4 и 2.1.1-2.5.1, которые приходится рассматривать при бурении. Для них необходимо изучать распределения п.п. 1-12 при установившихся и неустановившихся течениях в элементах системы скважина - пласт. При решении конкретной задачи находят одно или больше распределений п.п. 1-12 так, чтобы они не противоречили остальным. Например, распределение давлений в подземной части ЦС, которое часто приходится находить при осуществлении гидромеханического процесса бурения с промывкой жидкостью или смесями.

На рис. 16 построено искомое распределение (эпюра) давлений в ЦС некоторой вертикальной скважины при бурении при заданной компоновке бурильной колонны с учетом условий: давление в стояке рст не превышает допустимого давления насоса рдоп, т.е. выполняется соответствие распределениям п.п. 2 и 9; давление в необсаженных частях скважины выше давлений в проявляющих пластах рПЛ1, рПЛ2 и не превышает давлений поглощения или гидроразрыва рр1, рр2: рПЛ1 < Р < рРь Рш.2 < Р < РР2, т.е. выполняется соответствие распределениям п. 2 и 9; расходы жидкости (^щ в КП и на забое (2за5 обеспечивают вынос шлама и являются одними из значений распределения п. 1; разность распределений давлений в трубах р,р и КП ры, удовлетворяет условиям прочности труб рпр-. ¡Ртр-РшК рПр; т.е. соответствует п.п. 2 и 9.

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

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

основном, приводятся постановки и решения задач, полученные автором (см. пп. 1.1,1.3 и 2.3 на рис. 2).

1. Установившиеся процессы

11. Бурение с промывкой жидкостью

1.2. Бур «ии ее про дли ком гамм

1.1.1. 1.1.2.

Выносная Выбор

способность насоса

потока

1.13. 1.1.4.

Буренке на Гидромонито-

равновесном рное бурение и

давления очистка за&оя

13, Буренке с промывкой многофазными смесям»

1.21, Выносная способность газа

1.2.1 Выбор

компрессора

1.3.1 Выносная способность смеси 13.2. Бурение иаравно- новесном далении

1.3.3. Выбор насосов м компрессоров 1.3.4. Бурение на депрессии

2. Неустановившиеся процессы

2.1. Пуск 2.2 С пт скоп о дь- 23. Цементи- 2.4. Осложнения 2.5. Вскрытие в

насосов ёмные операции рование при бурении испытание пластов

2.1.1 Время запуска насосов

22.1. Спуск колони с открытым концом

22.2 Спуск КОЛОНН с закрыть™ концом

23.1 Выбор режимов цемеятнро ваикя

241.

Борьба е гио-водонефте-проянл типмн

2.4Л.

Глушение открытых фонтанов

z:

2.4.3. Борьба с поглощтиямп

IZ.

2S.1.

Определение

параметров

пластов

2.4.4.

Кавернообразо-ванне я сужение ствола скважины

X.

2 A.S. Смятие н прихват колонн

Рис. 2. Перечень основных гидроаэродинамических процессов при бурении

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

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

(2)

(3)

альных уравнений по живому сечению канала площадью Б и замыкающих функций:

уравнений сохранения массы для каждой фазы

|-(5фР1)+^-(8фр1У1) = 112+:1т; ¿(8(1-ф)р2)+^-(8(1чр)ргу2) = :21+12т; (1)

о! ог о1 ох 4 '

уравнения движения ^(8(ФР1У1 + (1 - Ф)Р2У2 ))+^(фр^,2 + (1 - Ф)Р2У22 )) = = Б(<рр, + (1- ф)рг)дсоза-8^ + ^-[фр^2 +(1 - ф)р2У2];

уравнения энергии д д

-(8(фр,Е1 + (1 - Ф)Р2Е2)) + —^(«рр&у, + (1 - Ф)Р2Е2У2)) =

= -—+ (1 - Ф) у2)р) + 3(фр,у1 + (1 - ф)р2У2)ёсоз а + ртш +

где ф и (1-ф) - концентрации первой и второй фаз; р1 и р2 - плотности фаз; V] и \'г - скорости фаз; I\г= - - интенсивность фазовых переходов, ,Гт= 1(гп-+-12т, От«, - заданные суммарные интенсивности потоков массы и потока тепла фаз из внешнего пространства и за счёт конвекции; Е[= в1+ У]2/2 и Е2 = е2+ у22/2 - удельные энергии фаз; е^ Ь.-р/р! и е2= Ь2-р/р2 - внутренние энергии фаз; ^ и Ъ2 - энтальпии фаз; у_ - периметр поперечного сечения канала; ос - угол отклонения канала от вертикали; г - вертикальная координата; I - время;

уравнения состояния для обеих фаз с одинаковым давлением в фазах р=р(рьТ,),р = р(р2,Т,); (4)

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

Ф = ф(Р, Бг, Ые, р, ц, б), (5)

где р = Ql/(Ql+ СЫ -_расходное содержание первой фазы; СЬ и 0>2 -объёмные расходы фаз; р - отношение плотности первой фазы к плотности второй; ц - приведенная по воде вязкость жидкой фазы; г - относительная шероховатость труб; Рг, Яе, 'У/е - безразмерные числа Фруда, Рейнольдса и Вебера смеси, выраженные через параметры входящие в уравнения;

теоретической или эмпирической зависимости для коэффициента гидравлического сопротивления смеси

^т = Рг, Ле, \Уе, р, 11, е). (6)

Величины Ьь Ъ2, (3, Qw находятся с привлечением термодинамических законов и опытных соотношений.

Система (1)- (6) содержит 8 уравнений с восемью неизвестными ф, рь Р2, V], \'2, р, Т и Хт. Следовательно, система замкнута.

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

р = р (рь т, У1,...,У„), р = р (р2,Т,хь...,хп); (7)

уравнение фазовых концентраций Ф^,Г|1,...,Т1П ,кь...,кп)= 0> (8)

где W - мольная доля газовой фазы; т|; - мольная доля ¡-го компонента в смеси; к; - коэффициент распределения для ¡-го компонента;

выражения для мольных концентраций компонентов у\ и X;: = х^х^ЛькО, (1=1,2,...,п); (9)

выражения для летучестей £ компонентов газовой и жидкой фаз: ^ = ^ = Г2г(р,Т), А, =0, (¡=1,2,...,п); (10)

выражения для определения расходного газосодержания р Р = (3 (р,Т,Уь—>Уп ,х,,...,хп); (11)

соотношения для динамических коэффициентов вязкости фаз: Ц1= (р,Т,уь...,уп), = Цг (р,Т,хь...,х„); (12)

соотношения для энтальпий фаз: Ь1 =Ь1(р,Т,уь...,уп), Ь2 = Ь2 (р,Т,хь...,хп); (13)

выражения для поверхностного натяжения на межфазной границе ст= ст(р, Т, уь...,уп, хь...,х„). (14)

Вид соотношений (7)-(14) известен и определяется законами термодинамики. Так как компонентный состав смеси задается, добавилось (1+5п) переменных (\У, х;, у, к^п^) и (7+5п) уравнений. Шесть из этих уравнений служат для определения недостающих величин ¡3, 11ь Ь2, Ць Цг, ст. Таким образом, для многокомпонентной двухфазной смеси система состоит из (15+5п) уравнений с (15+5п) неизвестными и также является замкнутой.

Система уравнений допускает аналитическое решение при обоснованных предположениях. В основном её решения получают численными методами. Для получения решения системы для задач, указанных на рис. 2, и построения, например, эпюры распределения давления, приведённой в правой части рис. 1, следует задать вид функций (4)- (5) и (7)-(14), а также начальные и граничные условия, диктуемые каждой задачей.

Система основных одномерных уравнений, которая используется для решения установившихся течений в трубах и КП в данной работе, имеет более простой вид и является следствием системы (1)- (3):

(Ю! т . сЮ2 т

1Г 12+ 11ш;"121+12

Лг т г Б^г^р, (1-ф[)р2

-.2

+ Х^д 8Б3

/

в

2

_с1_ Аг

,2

ь2 + |-|С2

Ф1Р1 (1-ф])Р2

(15)

где рт = ф!р1+(1"ф1)р2 - плотность смеси; в] = р^р^^ и С2= ргО-ф^Б -массовые расходы фаз, Ош = 02 - массовый расход смеси, равный сумме массовых расходов смеси; Ьь Ь2 - энтальпии первой и второй фаз; <3™ -поступивший в канал приток тепла за счёт конвекции и 1(5 = хк(Т~Т0)~ из

внешнего пространства вместе с обеими фазами; То= Т„+ Г(г- т^) - температура пород в зависимости от глубины г, определённая по геотермическому градиенту Г с учётом температуры Тп нейтрального слоя на глубине т^,.

К уравнениям (15) также добавляются уравнения (4)- (5) и (7)-(14), в которых учитываются физико-химические свойства фаз.

Для систем газ-жидкость уравнения (5) и (6) экспериментально определены в различных научно-исследовательских организациях и рядом авторов. На основе этих работ в диссертации получено обобщение, удобное для составления алгоритмов счёта и последующих расчётов на ЭВМ.

Вид соотношений (5) и (6) зависит от структуры (режимов) течения. Обобщенная запись формул для истинного газосодержания ср для потоков в вертикальных трубах по данным ООО «Газпром ВНИИГАЗ» имеет вид:

А • Р • ^ (Бг ) ПРИ Ут Уа ~ пузырьковый или снарядный режим;

ф _ . А • Р ■ ^ (Ргт ) + {2 (Р) при у„ < ут < V, - кольцевой;

А ■ Р • ^ (Ргт) - £3 (Р) ПрИ уг < ут < уст — дисперсно-кольцевой; ^ при ут > уст =5 м/с - дисперсный,

где А = 0,5 + 0,3-ехр[0,067-(1-р)], = 1-ехр^,4^Ргт/Бг*);

(16)

т=

(1-А)..

-2-0-й)

\0,25

g•a

оф(-7,5-^1=р); Г3©=(1+А-2-р)-е?ф(-7,5-^);

_0,8б-уг-ехр[9'(1-р)]_ ' ~~ 3,3-[1+0,0027-(р-1)]'

р,-р8; " 3,3-[1+0,0027-(Ц-1}] - "ехр[9 (1-р)]' Бг* = 4 • [1 - ехр(0,1 • ц)] - 3 • [1 - ехр(1 - 0,05 • Д)].

Коэффициент сопротивления Хт газожидкостной двухфазной смеси представлен в виде

1т = Ц11ега,б).у(Ргт,\Уеп,р,р;,р), (17)

где А,(Яеш,е) - коэффициент сопротивления однофазного потока; V)/ = Хт/Х -приведенный коэффициент, характеризующий отклонение коэффициента сопротивления смеси от аналогичного для однофазного потока.

По данным ООО «Газпром ВНИИГАЗ» приведенный коэффициент сопротивления в зависимости от режима течения записан в следующем виде

= Е{^(р,р,Р где Е = 1 + 0,03-р; Г4 -

при

ХИзМ] при

при V,

v < v < v уа — "т ^ усг>

(18)

1-Р + 0,03-ехр

-1350 (1-Р)3

и =6,35-щ>илЩ-®Щ, к=ут-((р1 -р8)/(ё-а))'

\0,25

Для расчёта коэффициента сопротивления X для однофазного потока использована формула Черчилля, которая справедлива во всем диапазоне чисел Re и является удобной для составления алгоритмов программ.

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

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

Для вычисления Jm = Jim + J2m используются соотношения вида

, _Рщ,2(2)-Р2. m A'.(L-H)'

(19)

л' = Рго12 -р32 = А, -Um +B1-Um2;Um =Gm/p„,

я-к-И-Т^ Ь-Н

где Аь В] - фильтрационные коэффициенты сопротивления пласта; р^, Тст, рст -давление, температура, плотность при стандартных условиях; Ь - общая длина скважины; Н - длина обсаженной части скважины; к, Ь - проницаемость, толщина пласта; 11К - радиус контура питания; ц - приведенный средний коэффициент вязкости; Z - средний коэффициент сверхсжимаемости.

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

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

Выведено обобщённое уравнение гидростатики многофазных флюидов, из которого следуют частные случаи гидростатик двухфазных смесей, однофазных НЖ и ННЖ, наиболее распространённых в практике строительства и эксплуатации скважин. Осреднённое по живому сечению канала уравнение гидростатики изотермического устойчивого или предельного равновесия многофазных НЖ и ННЖ в вертикальных и наклонных каналах в поле силы тяжести получено как предельный случай уравнения движения многофазной смеси при стремлении скоростей фаз к нулю: <1р 4т„ <1р 4(т0+т,.яп)

¿ = Рсм8-соз« + ^ или ¿ = (20)

где рсм = - плотность смеси; pi - концентрация и плотность ¡-ой фазы (1 = 1...п); тСг= Хо+Ткап- касательное напряжение на стенке канала с учётом напряжения сдвига то и капиллярного напряжения ткш; ё(/) - переменный гидравлический диаметр канала, зависящий от координаты /, совпадающей с осью канала; а - угол наклона канала к вертикали г. В случае непостоянной

кривизны скважины угол а - заданная функция координаты г. Для вертикального канала (а = 0) координаты I и z совпадают.

Вторые слагаемые в правых частях (20) равны нулю для смесей НЖ и ННЖ, не обладающих динамическим напряжением сдвига и при отсутствии капиллярного напряжения. Не в капиллярных каналах это слагаемое не всегда равно нулю, и, следовательно, запись гидростатического уравнения в форме (21) отличается от общепринятой записи. Оно также не равно нулю при наличии действия поверхностных сил различной природы, например, сил поверхностного напряжения (натяжения) между стенкой канала, жидкостью и свободной поверхностью. В силу малости второго слагаемого для относительно больших диаметров труб (не капилляров), величина этого слагаемого незначительна, а в безграничной жидкости равна нулю даже при конечных и достаточно больших тсг различной природы (Ит(т - /d) = 0).

Следует отметить, что при приложении градиента давления к флюидам возникают напряжения и, непосредственно перед их сдвигом, в зависимости от направления действия градиента dp/d/, напряжения достигают своего максимума (|тст/| = т0) при предельном равновесии. Такие флюиды широко используются в нефтепромысловой практике, особенно, при строительстве скважин.

Из (20) следует общеизвестное основное уравнение гидростатики для однофазного флюида в поле силы тяжести, не обладающего динамическим напряжением сдвига (1=1, хт= х0= 0, cpi = cpi = 1, pj = р). Таким образом, для безграничного канала или величины x^/d равной нулю имеем основное уравнение гидростатики, которое обычно приводят в курсах гидродинамики dp/dz = pg.

В гидростатике истинное содержание ср можно представить в виде ф = X/(X+Pi/p2(l-x))- Здесь х = V(l+Tl) - массовое содержание первой фазы; г| = а-ро/рг - массовый коэффициент аэрации; а = Q0/Q2 - расходный коэффициент аэрации; р0 и Q0 - плотность газа и расход (подача компрессоров) при атмосферных условиях, рг и Q2 — плотность жидкости и расход (подача насосов). Величины %, г| и а (при отсутствии растворимости фаз), в отличие от ф, не зависят от давления и их удобно использовать. Следует отметить, что в выражения для концентраций входят динамические переменные. В гидростатических условиях на момент запуска или остановки насосов, когда отсутствуют (стабильная пена) или пока не включились механизмы оседания или всплытия, эти оценки массового содержания и коэффициента аэрации для вычисления плотности справедливы. Для учёта растворимости одной фазы в другой можно использовать закон Генри % = k-р, где % - часть массовой концентрации, перешедшая из одного состояния фазы в другую.

Для интегрирования уравнения (20) и получения основного уравнения гидростатики двухфазной жидкости принято, что термодинамические

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

Pi = a¡ + b¡-p, (22)

где a¡, b¡ - опытные коэффициенты (i = 1,2).

В частном случае, когда первая фаза (например, воздух) подчиняется уравнению состояния реального газа, а вторая (например, вода) - несжимаемая, то ai=0, bi=l/( ZR- Т) и a2=const, b2=0, где R - газовая постоянная, Т -осреднённая температура. Если обе фазы представлены слабосжимаемыми жидкостями, то: а;=рю(1- Pipío), b¡= pi0Pi, где pi - коэффициент сжимаемости; Pío, Pío - постоянные значения. Анализ уравнений состояния жидкостей, газов и твёрдых веществ, встречающихся в бурении, показал, что значения а, и Ъ, лежат в пределах: a¡ = 0-гЗ-103 кг/м3, b¡= 6-10"8-й,3-10"5 кг/(м3/Па).

С помощью выражения (22) получено аналитическое решение (20) в зависимости от знака дискриминанта А = 4Ai Ci - ВД которое при а = const имеет вид

А = кр(Ь, - Ь2); В = Ъ2Х+ Ь,(1-у>кР(а, - а2); С = а2/+ а, (1-х); А, = Ь,Ь2+ Б-крф] - Ь2); В! = а!Ь2+ а2Ь1+Б-В; С, = а,а2+В-С; Б = 4т0 /(а^-сова); Ь - вертикальное расстояние между заданными горизонтальными плоскостями.

Выражения

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

(23) можно проводить с помощью номограмм для определённого флюида.

При кр = 0 получаем гидростатику смеси газа и несжимаемой жидкости или смесь газа со шламом (в момент остановки). Для таких смесей гц = 0 и Ь2 =0. и дискриминант всегда отрицателен Д = -[а2Ь]+В'Ь1(1-х)]2 < 0. Дискриминант может оказаться положительным, так как его знак зависит от величины и знака динамического напряжения сдвига т0.

Если фазы несжимаемые, например, смесь промывочной жидкости со шламом или две несмешивающиеся жидкости, то Ь1 = Ь2 = 0 и дискриминант

h = fi(p)-fi(p0) + f2(PH2(Po) при Д > 0; Ь = fi(p)-fi(po) + f3(p)-f3(Po) при А < 0; h = fi(p)-f1(Po) + f4(P)-f4(Po) при Л = 0

(23)

равен нулю А = 0. Дискриминант равен также нулю, когда и а] = 0 и а2 = 0, то есть флюид является двухкомпонентной однофазной смесью двух газов.

В диссертации приведена классификация гидростатик для разных видов термодинамических уравнений фаз с нулевыми или ненулевыми коэффициентами и Ь; (1 = 1, 2). Из классификации выделены уравнения гидростатик, наиболее востребованные для расчётов в бурении и эксплуатации скважин.

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

Известно из опыта, что однофазная НЖ в зависимости от лиофильной или диофильной поверхности капилляра занимает разные положения равновесия. Радиус капилляра г связан с радиусом мениска го и краевым углом а, а поверхностное натяжение а на стенке капилляра выражается через касательное напряжение ткап и капиллярное поднятие Ь. При равновесии жидкости в капилляре градиент давления равен нулю. Используя зависимость а = ткап • Ь и нулевой градиент давления в (20), получена формула для высоты капиллярного поднятия жидкости, в том числе с * 0, с учётом гидростатического давления газовой фазы р0яЬ

го(Р-РО)'Е + 2'гО'

При То - 0 формула (24) переходит в формулу Жюрена. Формула (24) также легко обобщается и на большее количество фаз.

Для горизонтального капилляра в основном уравнении гидростатики (20) отсутствуют гравитационные силы при То ■*■ 0 или т0 = 0. Из (20) при т0 = 0 также следует, что пока не преодолены капиллярные силы, то есть не создан градиент давления, превосходящий значения правой части уравнения, движения жидкости в капилляре не будет. Вода не обладает динамическим напряжением сдвига (т0 = 0), поэтому предельный градиент сдвига, возникающий при страгивании воды в пористых средах, необходим лишь для преодоления капиллярного сопротивления.

Законы гидростатики применены к задаче цементирования скважин стабильными пеноцементными растворами с использованием решений (23).

Цементирование скважин пеноцементными растворами

Расчёт цементирования скважин пеноцементными растворами позволяет оценить основные параметры его режима при постоянной и

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

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

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

а) без наличия в верхней части кольцевого пространства скважины столба из "чистого" (негазированного) цементного раствора;

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

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

Рис. 3. Режимы работы насосов цементировочных агрегатов и компрессоров СД 9/101 во времени: 1, 2 - расчётная и рекомендуемая ступенчатая подачи насосов ЦА; 3 - подача компрессоров; 4 - степень аэрации. Шаг аппроксимации по подаче насосов выбран согласно паспортной характеристике насосов 9Т.

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

Основные уравнения установившегося движения в каналах ЦС следуют из общей системы уравнений (1)- (3), если в ней принять ф =1. Для полноты изложения даются методики расчёта как ламинарных, так и турбулентных течений НЖ и ННЖ с применением известных решений для ряда технологических задач бурения нефтяных и газовых скважин.

В частности, автором показано, что формула для коэффициента гидравлических сопротивлений X при ламинарном течении ВПЖ в круглых каналах

Х = —— Ле8р

(25)

в результате предельного перехода при т0 0 переходит в решение для НЖ

64 Бе Яе8р

64 Яе

1п8 1-8'

так как

->1.

т0-+0

Бе

'„-о - — ^

Для ламинарного течения ВПЖ в круговой щели известное решение относительно расхода представлено в безразмерном виде

8е = Не/Ке = Зр/(1-1,5р + 0,5р3)) (26)

где Р = 2то./(Н9р/йг) > 0; Бе = т0.л г Н2/(|С>|т1) - число Сен-Венана для течения в круговой щели; Н - раскрытие щели; г - текущий радиус щели; г| -динамический коэффициент вязкости ВПЖ. График зависимости (26) показан на рис. 4.

"0-1 V и к» зе

Рис. 4. Кривая зависимости р = Р(Бе) для круговой щели

Следует заметить, что если Р > 0, то < 0. Этот случай соответствует притоку,если р<0и(2>0-поглощению.

Для расчёта перепада давлений с использованием (26) следует вычислить число Бе = т0Н/(2г|у) при г = гк, из рис. 4 найти I р I и определить

г 2т т

перепад Др = рс-рж по формуле = у1п5 , где 8 = — ; гк - радиус

Н|Р| гс

контура; гс - радиус скважины. Последняя формула преобразована к виду

формулы Дарси-Вейсбаха и получен коэффициент гидравлических

сопротивлений при течении в круговой щели ВПЖ, который переходит в

коэффициент для течения НЖ при т0 —» 0

XsJSeJn» vHp (27)

Re |ß| 1 - 5 r|

Для развитого турбулентного течения ВПЖ в трубах получен коэффициент гидравлических сопротивлений, который переходит в коэффициент гидравлических сопротивлений по Альтшулю для НЖ при т0 -> 0 (Не -> 0):

^«(^ЕЮГ <28)

Известно, что переходный режим для НЖ наступает при критическом числе Рейнольдса ReKp = 2100. Основой для определения чисел ReKp для ВПЖ и степенных жидкостей (СЖ) является гипотеза, что при приближении свойств ННЖ к НЖ, кривые гидростатических сопротивлений приближаются к кривым гидравлических сопротивлений для НЖ. При очень больших числах Re жидкости с любыми физическими свойствами имеют одну асимптотическую область автомодельности гидравлических сопротивлений в гидравлически гладких трубах. Зависимость Блазиуса А = 0,316Re'0,25 на графике Никурадзе (рис. 5) пересекает кривую ламинарного режима А, = 64/Re при числе Re =1187 < Re*,, = 2100, а точка пересечения является началом отклонения течения от ламинарного режима. X

2

J0"' 8

б

4

3,048 2

1,28

10"J

,0J 2Ы 4 6 S104 2 4 6 8,0s 2 V4 6 Rc

Рис. 5. Зависимости для коэффициента сопротивления в гидравлически гладких трубах: 1 - по Пуазейлю; 2, 2' - переходный режим для НЖ и ВПЖ; 3, 4 - для турбулентного режима вязкой жидкости (Не = 0); 4, 5 - для критических режимов; б - по Букингему; 4, 7 - для турбулентного режима ВПЖ в гладких трубах в зависимости от параметра Не

В дальнейшем для ВПЖ в качестве критической кривой принята линия 5 (см. рис. 5), соединяющая две предельные области, общие для НЖ и ВПЖ. Эта линия соединяет две точки: первая с координатами ReKp = 2100 и

А„ф = 0,03048 соответствует началу переходного режима вязких жидкостей и является предельной для начала переходных режимов течений ВПЖ при приближении их свойств к вязким; вторая с координатами ReKp = 3,7-105 и Я,,ф = 0,0128 определяет начало автомодельного режима течения жидкостей. Уравнение линии 5 имеет вид

Л 0,11

~ Re0,168'

Для ВПЖ при ламинарном течении зависимость коэффициента X от чисел Re известна - это формула (24). Подставляя в неё Re = ReKp, получаем формулу, связывающую Хкр, ReKp и параметры Не или Se л 64 Не

(30)

При Rejp > 3,7-105 можно принять л,кр = 0,0128.

Исключением },кр из уравнений (30) и (29), и подстановкой ReKp в зависимость, которая связывает числа Re, Не и Р при течении ВПЖ в трубе Не _ 8(3

Re 1-W*

Зн 3

получена система уравнений для определения критических чисел Rek-p 0,11 _ 64 Не Не _ 1

Re«« ReKT Re^SP i_W (31)

В результате численного расчета системы (31) найдена зависимость ReKp = /(Не), которая хорошо соответствует формуле, полученной Е.М. Соловьёвым при обработке теоретического и опытного материала по течению ВПЖ, включая буровые растворы, Rekp = 2100 + 7,3(Не)0,58.

Подставляя в (29) Хщ, = 0,0128 и заменяя / = Se/8p при больших значениях He/Re её аппроксимацией / = 0,125He/Re, имеем

Re,, = 25^Не-

Переходя к размерным величинам, получаем формулу Б.С. Филатова для критической скорости % = 25ф0/р при автомодельном турбулентном

течении в гладких трубах. При Не < 4-104 эта формула становится неточной.

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

Проводя аналогичные рассуждения, в работе получены выражения для критических чисел при течении в трубах (8 = 0) степенной жидкости (СЖ).

Подобные формулы для ВПЖ и СЖ можно получить и для кольцевого канала при 6^0, однако в настоящее время недостаточно экспериментальных данных, чтобы подтвердить их. В этом случае в первом приближении можно пользоваться для определения чисел ReKp формулой Е.М. Соловьёва, подставляя в неё значение гидравлического диаметра dr = dc - dH.

В шестой главе рассмотрены установившиеся изотермические течения двухфазных смесей в элементах ЦС.

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

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

При установившемся ламинарном или турбулентном течении и при отсутствии фазовых переходов система уравнений (15) для газожидкостной смеси записана в виде:

уравнение движения

|=е(фР1+(1-<?)р2)±^(фр1У12 + (1-ф)рЛ2)-[фр1у1^+(1-ф)рЛ^-]; (32)

где знак плюс берётся для восходящего потока и минус - для нисходящего;

уравнения сохранения массы S<pp,v, = QjP, = m, = const, S(1-9)p2Vj =Q2p2 = m2 = const; (33)

термодинамические уравнения состояния p = zRTp,g; p2 = const, (34)

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

уравнение концентрации для первой фазы <P = q>(ß,Fr); (35)

уравнение для коэффициента гидравлических сопротивлений смеси ^ = Xc(cp,ß,Fr). (36)

Вводя безразмерные переменные p = p/(2lTp2g), 4 = z/(zRT), преобразуем уравнение движения (32) к виду

1 h П2 , 1 ]__iL^fj__П2 1 dtp

dp Р;~ 2 [фр (1 -Ф>J ЖГ [l-ф ф2р/ di; (37)

zRT фр2

где k2 =Q2/(gdS2)- число Фруда; т| = Q0p0/(Q2p2) = ар0/р2 — массовый коэффициент аэрации; Q0, ро - объёмный расход и плотность газа при нормальных условиях (То, ро); а = Q0/Q2 - расходный коэффициент аэрации.

Уравнение (37) можно также разрешить и относительно производной от истинного газосодержания

?.„к

1-ф(1-р)± йф I

фр (1-ф)

zRT фр2 J d£

dE, d J 1 гi2

■=к

zRT ^ 1 — ф ф2р

(38)

Дифференциальные уравнения (37) и (38) в работе использованы для вычисления распределения давления или газосодержания ф вдоль канала.

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

Ф = 0,81Р при Ш - ¡3) ,1,93 при 93, (39)

где Р = <3) /(СЬ + СЫ = У(л + р) - расходное газосодержание. Полученное ими решение сведено к соотношению N = N1 + ^, (40)

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

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

Ф = —РдляО,21<Рг<5; Ф = 1,24р даяРг>5. (41)

■7Рг - 0,45

Коэффициент Х0 для нисходящего потока несколько выше, чем для восходящего и принимается постоянным и равным в среднем 0,06.

Используя полученные значения для ф, Бг, Р и опуская часть инерционного члена в (41), уравнение движения записано в виде

а? '

, где к2=^> п' = тМк-0,45)

[.Л+Р Р и получено его аналитическое решение

1-1 2

1 , (Ы)рг-1П' 1 ,

-г-1П"--- +-5-1П

(1-1)2 а-ад-^1 (1+0

а+ррг+ьу

(42)

(1+1)р-м?У

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

Чтобы удобнее пользоваться формулой (42) для нахождения р на входе в трубы одного типоразмера, если известно давление рт на выходе, она сведена к уравнению относительно числа М, для нахождения которого построена номограмма (рис. 6) N = М + 0,51я [М - 0,217а|+0,5(1 - а)2^ |М+0,217а(1 -а)|,

где К=М0+0,%|Н-0>217а)+0,5(1-а)21§|М0+0,217а(1-а)|-^(1-1)24; (43)

М0=-^(1-а); М = -|-т(1-а); а = 21/(1+1). При а = 0 эта формула переходит в формулу (40).

Рис. 6. Графики для определения чисел М: а - для случаев: 1 - 0 < N < 5 (за ВТ, УБТ, турбобуром, в долоте), 2 - N < 5 в турбобуре), 3 - а > 1, Моо < 0, -1,5 < N < 0 (внутри БТ, УБТ), кривые 1 и 2 построены в положительных координатах, кривая 3 - в отрицательных; б - верхнее семейство кривых: а = О (за УБТ, БТ и турбобуром), а < 1, Моо > 0, (в УБТ и БТ); нижнее семейство кривых: а< 1, Моо < 0 (в УБТ и БТ) (см. рис. 1)

Последовательность расчёта давлений с помощью (43) следующая. По исходным данным (L, db d2, рт и т.д.) определяют числа r(', М0, а, N и знак числа Моо= М0-0,217а. По соответствующей кривой рис. 6 находят число М, затем рассчитывают искомое давление р = 2,ЗМт|'/(1-а), или в размерном виде p = p-(zRTp2g).

Случай М0-0,217сс = 0 соответствует значению градиента dp/d^ = о . При этом течение неустойчивое и при движении нисходящего потока давление по длине участка трубы данного диаметра не меняется (р = const).

На основе общего подхода автором получены методы расчёта и для других элементов ЦС с использованием номограммы (рис. 6). Представленная модель служит для расчётов распределения давления при изотермическом движении двухфазных смесей в элементах ЦС.

Одномерное установившееся неизотермическое течение двухфазных смесей.

В общем случае не удаётся для таких задач найти аналитическое решение. Для удобства получения численного решения система уравнений (1)-(3) разрешена относительно градиентов давления и температуры в виде системы обыкновенных дифференциальных уравнений

¿От йг

£Г

¿г ст.д

аг"

л,-

(44)

1+0,

1 ар

1 50т

О

ш0р

(1-ф)р2р„+Ф(1-р)2р1

^-т ,

-----2 ;ф2=Рш82-^Ф1;Фз=От<Р3(Р.Т.От);

82(рр8+а-р)р,) 20

1

Рр„+(1-Р)Р1

( 2> V (

к ) +(1-р)р, ч

1 2

V т ф ) сТ ЭТ ф

По уравнениям (44) при граничных условиях р = ро и Т = То, при г = го с использованием формулы для притока в виде (45) получены уточнённые распределения давлений по стволу скважин месторождения Карачаганак (рис. 7).

50 р,МПя

1000

Рис. 7. Сравнение результатов расчётов распределения давления в скважинах №109 и №121 месторождения Карачаганак: 1, 3 - по разработанной модели; 2, 4 - по методике Ю.П. Коротаева и др.; ® - по формуле Г.А. Адамова; * - по методике ВолгоУралНИПИга-за; о - промысловые замеры.

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

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

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

В работе используется формула для скорости страгивания частички в виде ус=КЛ/4-(рч-рж)-ач-1/(3-Сигрж) ( где (1ч - диаметр частицы; рж -плотность жидкости; 1 = в(ц.т-эта+соБа) - эффективное ускорение свободного падения, К - коэффициент пропорциональности; коэффициент трения между частицей, стенкой и другими частицами, определяемый экспериментально.

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

Рис. 8. Течение в плоской щели: d, - эквивалентный диаметр частицы шлама, Н - величина зазора, W(r) - профиль скорости течения НЖ или СЖ.

П cpUMKIKutt

сло|

I грммвтяий

СЛФ0

Рис. 9. Течение в плоской щели: ёч - эквивалентный диаметр частицы шлама, Н - величина зазора, Но - размер ядра, - профиль скорости течения ВПЖ

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

В седьмой главе рассмотрено неустановившееся течение двухфазных смесей.

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

В систему одномерных уравнений входят уравнения сохранения массы, движения и энергии в безразмерных переменных:

Э(1~ф(1-р)) 80т _т .

= X

91 54

Ж.Т [Л 8Е, дЕ, _]'

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

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

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

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

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

Для интегрирования системы ■'<■».■ (44) при известном массовом расходе газа От задаются начальные условия в виде р(Н) = р, Т(Н) = Т , где Н - глубина скважины, а р*, Т - давление и температура на забое. В случае аварийного фонтанирования скважины величины р , Т являются неизвестными и нужно задать ещё два условия. Первое является индикаторной кривой пласта рш2-(рТ = А'От+В-От2, где р,ш - известное пластовое давление; А, В - фильтрационные коэффициенты, полученные при исследовании скважин на стационарном режиме. При аварийном фонтанировании на устье скважины (г = 0) может возникнуть критический режим истечения, откуда следует второе условие: Эр/Эг-»оо, сТ/&—ко при г—>0.

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

Течение в скважине принято стационарным, а тепловое поле в породе нестационарным. В момент времени 1 = 0 температура Тс определяется геотермическим градиентом. Из системы уравнений находятся функции р(г) и Т(г) при заданных значениях р\ Т* и Ога,. Считая, что в течение некоторого промежутка времени Т(г) остаётся постоянной, производим расчет внешнего поля температуры и находим новое значение Тс(г). Расчет ведется до изменения Тс(г) на заданную величину ДТ. По значениям Тс(г) и вычисленным р, Т и От находим новые зависимости р(г), Т(г) в скважине и продолжаем расчет внешнего поля температуры. Затем снова находим р(г) и Т(г) в стволе и т. д.

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

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

Во всех случаях распределение температуры в породе описывается уравнением теплопроводности в безразмерной форме

а

9 8, | 1 59; Эг2 г а

(47)

где 6i = Т/Г; 1 = т/т0; г = с,/г0; = а^о/г02; 1 = 1 относится к талой породе, { = 2-к мерзлой, 1 = 3 - к немерзлой; ai - температуропроводность; г0 - радиус скважины; т0- характерное время; Т - температура в стволе; —текущий радиус; т - время.

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

-—2- = а2(92 -1), r = l; в2 = ©} , t = 0; 6/J,t¡) = 0; 02 ->0¡ , г со. йг

Здесь t,,, = t/то соответствует моменту начала протаивания мерзлой породы, 6=ТКД - начальному распределению температуры в мерзлой породе.

На втором этапе, когда в породе две зоны из талых и мерзлых пород, уравнение (47) дополняется условиями

s(tj) = l, 92(r,t¡) = e0(r),r>S(t); = 0^(0,-1), г = 1, б, = 02 = 0, г = S(t);

от

, Э9, , 592 dS п „

= г = S(t); 02 6¡ , г -> оо.

Зг йг dt

В немёрзлой зоне (ниже зоны мерзлоты) при i = 3 имеем ^ = а3(03-1), r = l, 03=0°, t = 0; 03->eO)r-)-oo, а, = к • ib / X* Дt = Х*Тт0 /(r2qpL )(i = 1,2,3); 0О = Т0 / Т; 0° = Т° / Т; S(t) = ?(х) /г0. Здесь к - коэффициент теплопередачи, - коэффициент теплопроводности, q - удельная теплота фазового перехода воды в лёд, Pl - массовое содержание льда в породе, Т0© - распределение температуры в породе при т < тм, T°(z) - начальное распределение температуры в немерзлой породе, Е,(т) - подвижная граница фазового перехода (радиус протаивания).

Вместо условий на бесконечности вводится радиус теплового влияния R(t),

на котором выполняются условия 92 9), 03 -> 6°, 362 / йг + Э93 / Sr = 0, г = R(t).

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

3. Для выполнения вычислений создан пакет программ. В работе приведены результаты расчета для типовой скважины п-ова Ямал, где толщина слоя многолетнемерзлой породы h = 250 м. Мольный состав добываемого многокомпонентного флюида, %: СН4 - 96,37; С2Нб - 2,89; С3Н5 - 0,05; С4Ню - 0,03: С5 - 0,01; С02 - 0,22; N2 - 0,43.

Расчеты свободного дебита сделаны для двух вариантов.

В первом - фонтан возник непосредственно при пуске скважины в эксплуатацию (без учета протаивания). Получено: свободный дебит Qо = 4,22'106 м3/сут, давление и температура на устье скважины р(0) = 0,84 МПа, Т(0) = 247 К. Температура на стенке скважины соответствовала геотермической. Во втором варианте фонтан возник через 352 дня после начала эксплуатации. Получено: свободный дебит Q0 = 4,4-106 м3/сут; р(0) = 1,03 МПа; Т(0) = 257 К.

Для обоих случаев в работе приведены расчётные зависимости изменения давления p(z) и температуры T(z) по глубине скважины. Численные расчёты изменения радиуса протаивания c,(t) от времени приведены в работе для глубин z = 200, 100 и 0 м. Изменение внешнего температурного поля вследствие прогрева и протаивания породы почти не влияет на величину дебита, но заметно сказывается на устьевом давлении и температуре. Из расчетов следует, что на всех глубинах граница t(t) монотонно растет с увеличением времени, причем особенно резко в начальный период

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

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

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

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

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

За начало отсчёта (t=0) процесса глушения скважины принят момент начала подачи задавочной жидкости плотностью рр с постоянным расходом Q3 по наклонной скважине в ствол фонтанирующей скважины с дебитом газа Qo, который известен или может быть найден.

При глушении в пласте и в призабойной зоне скважины движется газ, а в трубном и (или) кольцевом пространстве - газожидкостная смесь. Жидкость, продвигаясь к устью, увеличивает забойное давление р3, уменьшая приток газа. То есть, в стволе скважины имеет место неустановившееся течение двухфазной жидкости. Глушение произойдёт, когда р3 станет равным или больше пластового. Считается, что задавочная жидкость -несжимаемая, газ подчиняется закону (34), фазовых переходов нет, фазы движутся с равными скоростями (ф = |3), движение изотермическое и турбулентное (X = const), инерцией пренебрегается.

При сделанных предположениях система (46) сильно упрощается и получается, как частный случай, система, состоящая из уравнений сохранения массы для газа и задавочной жидкости, уравнений движения газожидкостной смеси и уравнения индикаторной кривой пласта: dpP , dpPv _ Q Зр3(1-Р) | Эр3(1-Р)у о .

St Sz ' 5t dz '

U = g[pP + P3(l-P)].[l + |^.vI];; (48)

p3=p(L,t) = VP™-aQ0(t)-bQ2(t) .

В качестве начальных условий (t = 0) взяты постоянные распределения давлений в пласте и скважине (дебит фонтана Q0 = const), поскольку до начала глушения течение газа в системе скважина - пласт установившееся.

Граничные условия: устьевое давление и расход задавочной жидкости постоянны при I > 0.

В.Д. Малеванский, Е.В. Шеберстов и Е.Г. Леонов использовали систему (43) для определения параметров глушения.

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

Например, на рис. 10 показаны два типичных графика при Ь/с1с =1600. &

4

3

2

1

24 6 8 811 24 6 8 БЬ

Рис. 10. Зависимости у/Ёт от БЬ при рш/ру = 2 и различных безразмерных переменных От для глушения газовых фонтанов. Другие графики в альбоме также выполнены в указанных координатах для различных отношений пластового давления рщ, к забойному давлению р30. На рисунке рг - плотность воды.

Каждый график представляет собой набор кривых для заданных параметров. Чтобы использовать альбом, по исходным данным рассчитывают забойное давление р30 для движения чистого газа, находят параметры рпл/ру, Рпл''(РзёГ), рпл/рзо- Выбирают нужную кривую глушения фонтана. Задавшись определённым объёмом V задавочной жидкости и вычислив число БЬ = У/Ускв, с помощью соответствующей кривой находят значение л/Ёг и производительность насосов по формуле = .

Неустановившееся течение газожидкостной смеси при бурении на депрессии.

Рассмотрено движение раствора в бурящейся скважине с непрерывным поступлением с забоя флюида с заданным расходом <3ПЛ на некоторой глубине для Кончуренского ПХГ.

При бурении с заданной депрессией на пласт требуется определить давление рСТШ1К в стояке и давление ркп в КП перед штуцером для поддержания необходимой депрессии на пласт. Информация об изменении давления в стояке даёт возможность управлять безаварийным процессом бурения.

В диссертации использована изотермическая постановка задачи (32)- (36). Для скважины Кончуринского ПХГ численное решение представлено в виде изменения давления в стояке, которого следует придерживаться

В», бар

Лр устьеиЗб

и

/

Рис. 11. Пример поддержания заданной депрессии при нестационарном течении газированного раствора регулированием штуцера: Др^щр = 6,84 бар - требуемое давление для поддержания депрессии без притока; Ар1ЯТТр = 10 бар - требуемое давление на штуцере с притоком перед герметизацией; Дршггаз= 36 бар - требуемое давление на штуцере (газ на устье); tm!n = 24 мин - минимальное время, через которое газ появляется на устье; 1 мин -интервал времени регулирования.

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

Исследование течения «нулевой подачи».

При переходе с промывки на бурение с продувкой воздухом необходимо рассчитать производительность компрессоров, обеспечивающих вынос жидкости из скважины. При поступлении пластовых вод в скважину производительность компрессоров для выноса жидкости соизмерима с производительностью компрессоров для выноса шлама. До выноса жидкости из скважины существует режим течения - «нулевая подача», при котором скорость жидкой фазы w2 = 0. Поэтому автором были использованы известные опытные данные из литературы, дополненные исследованиями газосодержания ф на лабораторной установке при числах Fr < 0,1.

Установка состояла из стеклянной круглой вертикальной трубы с внутренним диаметром D = 79 мм, к нижнему концу которой был присоединен газопровод. Расход газа Qo измеряли расходомером. Измеряли начальную высоту столба чистой жидкости h и высоту газожидкостного столба Н при установившемся движении газа с заданным числом Фруда Fr„=Q2/(S2gD), где S=ttD2 / 4. За верхнюю границу высоты Н принимали поверхность раздела между газожидкостным столбом и образовавшейся пеной. Вычисляли величину водосодержания 1 - ф = h/H.

Опытные данные аппроксимированы формулами, коэффициенты которых получены методом наименьших квадратов:

1-0,0638^ т,^™ , 1-0,0461757 1_ф =-:-__ при Бг< 130 и 1-(р = -— приРг> 130.

1 + 2,489%/^

1 + 4,477^

Задача распределения давления в режиме «пулевая подача». Для данной изотермической задачи уравнение движения следует из (46), если принять скорость жидкости \у2 = 0 и пренебречь инерционным членом

(49)

^ = 1-Ф(1-р) + ^—,где ф = ф(Рг), р = 1, Яс = Х0 = соп81. гк 2 Фр

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

ф = 2,553к/(2,489к+р) дляБг2 130; ф = 2,553к/(2,489к+р) дляРг< 130 и найден интеграл уравнения (49) с граничными условиями р'=р при \ = 0

и»

х!п

--?1п А + Вр + ср^| + -----

С 2С 2Ср + В-УА

2А24К

{Т

xarctg

2Ср + В + ч/д р В

¡А + Вр + Ср2

1Б7а

!п

2Ср + В-\/1

2Ср + В + Тд

II В — 2 АС

А + Вр + Ср2 +-

- ^ 1 2А л/Д

2Ср + В , —9==— + аЬ

М

2 С

1п А + Вр + С?2

В 2Ср + В

"сЖГъ м

_ при Д>0;

при А< 0;

(50)

4 =

2С

В _ — + Р 2С Н

+ 1п

В _ — +р

2С н

а

+ С

_ ь р--

2С

2 С

В _

+ 21п

В _ — +р

2С н

при Д = 0.

Приведены примеры расчётов давления в обводнённой газовой скважине № 81 Ставропольского месторождения, работавшей при «нулевом режиме».

Для расчетов давления в скважине по формулам (50) предложен простой метод. Результаты расчетов давления столба газированной нефти по предложенной модели удовлетворительно совпадают с промысловыми данными Ю.Г. Валишина для чисел Рг < 0,1. Расчеты по предлагаемой методике также достаточно хорошо совпадают с данными А.Ф. Мусаева.

Истечение струи газа через слой жидкости.

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

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

Схема последней установки, на которой можно проводить компьютерную обработку параметров изучаемого явления, представлена на рис. 12. Бак 1 (350x400x300 мм) из прозрачного органического стекла заполняется водой до некоторого уровня. Сменные насадки устанавливаются на держателе 2, положение которого может изменяться по отношению к поверхности воды как по высоте, так и по углу наклона струи, вытекающей из насадка. Воздух подаётся к насадку от компрессора 3. Для сглаживания пульсаций газа установлен ресивер 4, давление в котором регулируется краном 5 и измеряется образцовым манометром. Давление газа в подводящей трубке к насадку измеряется с помощью датчика 6 типа МРХ фирмы «Motorola». Для измерения расхода газа создан малоинерционный расходомер 7 тахометрического типа с электронным преобразователем

Рис. 12. Схема экспериментальной установки для изучения истечения струи газа через слой жидкости: 1 - бак; 2 - держатель; 3 - компрессор; 4 - ресивер; 5 - кран; 6 - датчик давления;

7 - электронный расходомер;

8 - электронная линейка;

9 - измерительная система PowerLab;

10 - компьютер; 11 - источник света; 12 - волногаситель

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

Параметры бугра (высота Н и ширина D в разных сечениях) измерялись по теневому отражению с помощью сконструированного датчика координат - электронной линейки 8 (рис. 12). Сигналы от датчика давления 6, расходомера 7 и линейки 8 поступали на измерительную систему PowerLab/8sp фирмы «ADInstruments» и далее на компьютер 10. Запись и обработка сигналов осуществлялась при помощи программы Chart, входящей в комплект системы.

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

Проводилась видеосъёмка образования и динамики изменения геометрических размеров бугра на поверхности жидкости при истечении воздуха из щелей с размерами 3,5x0,14 мм, 5,5x0,14 мм, 10x0,17 мм и отверстий с диаметрами 0,71 мм, 1 мм и 1,3 мм при массовых расходах воздуха от 0,01 до 0,73-10"3 кг/с.

сигнала на базе анемометра AM 402.

На держателе 2 (рис. 12) закрепляли насадок и, перемещая держатель по вертикали, устанавливали глубину погружения Нс под поверхность воды. Компрессором 3 нагнетали воздух в ресивер 4 и краном 5 создавали определённое давление в нём. При истечении воздуха из насадка на поверхности воды возникал газожидкостный бугор. Давление измерялось датчиком 6, а расход воздуха - электронным расходомером 7. Сигналы от датчиков поступали на измерительную систему Ро\уегЬаЬ и на компьютер. На экране, установленном на передней стенке бака 1, наблюдали тень от бугра. С помощью линейки 8 получали запись изменения напряжения II и его среднее значение в зависимости от высоты бугра Не за средний промежуток времени. По тарировочной таблице определялись геометрические параметры бугра.

Струя газа, распространяясь в жидкости, приводит в движение пограничные слои жидкости и возбуждает в ней циркуляционное течение. Проходя скачок плотности, струя расширяется. Жидкость, захваченная струёй газа, частично растекается по поверхности и частично поднимается с потоком газа над жидкостью, образуя газожидкостной двухфазный бугор. На рис. 13а представлена типичная фотография характерной формы бугра на поверхности: видна геометрическая особенность картины - диаметр основания бугра Р существенно превышает диаметр выходящей струи._

Рис. 13. Фотография и модели истечения: а) - типичная фотография истечения газа через слой жидкости ((2 = 0,0045 м3/с; Ь = 0,2 м; О0 = 0,03 м; Н = 0,27 м; Я = 0,095 м); б) - предложенная автором геометрическая модель течения, основанная на наблюдениях во время опытов; в) - модель истечения Тепляцкого

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

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

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

Безразмерные параметры менялись в пределах: для круглых отверстий 1,5 < Fr < 149, 43,5 < Нс < 318, 6178 <_Re < 36583, 178 < We <1699; для щелей 1,2 < Fr < 169,8, - от 38,5 < Не < 285,7, 6454 < Re < 41633, от 162 <We< 1900.

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

Результаты опытов сравнивались с экспериментальными данными других авторов. Например, наши опытные данные и опыты М.Дж. Фридла имеют коэффициент корреляции г2 = 0,98 в области пересечения опытных данных.

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

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

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

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

сравнение опытных данных, полученных при истечении газа в НЖ (водопроводная вода) с данными при истечении в ННЖ (раствор ПАА).

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

Основные выводы и рекомендации

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

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

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

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

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

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

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

6. Найдено распределение давления в кольцевом пространстве скважины при переходе с промывки на продувку газом, необходимое для расчёта производительности компрессоров, обеспечивающих вынос жидкости из скважины в режиме «нулевая подача». Исследования при режиме «нулевая подача» дополнены новыми экспериментальными данными по барботажу газа через столб жидкости при малых числах Фруда (Бг < 0.01) и получена зависимость между параметрами течения в широком диапазоне чисел Фруда вплоть до начала выноса жидкости из вертикальной скважины.

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

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

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

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

Основные опубликованные работы по теме диссертации

Публикации в изданиях, рекомендованных ВАК

1. Астрахан И.М., Егорушкин С.А., Исаев В.И., Розенберг Г.Д., Слободкина Ф.А. Моделирование аварийного фонтанирования газоконден-сатной скважины в условиях вечной мерзлоты// Изв. РАН. МЖГ, 1994, № 3. С. 101-108.

2. Исаев В.И. Термогидродинамические модели многофазных течений при строительстве и эксплуатации скважин// Управление качеством в нефтегазовом комплексе, 2005, № 1-2. С. 100-103.

3. Исаев В.И. Обобщённая гидростатика однофазных флюидов и многофазных смесей в поле силы тяжести// НТЖ. Строительство нефтяных и газовых скважин на суше и на море, 2006, № 1. С. 26-32.

4. Исаев В.И., Иванников В.Г., Иванников A.B. Оценка параметров газопроявлений в виде грифонов и экологической обстановки окружающей среды// Управление качеством в нефтегазовом комплексе, 2004, № i.e. 5759.

5. Исаев В.И., Иванников В.Г., Иванников A.B. Экспериментальное мрделирование аварийного фонтанирования газовой скважины через слой жидкости// НТЖ. Строительство нефтяных и газовых скважин на суше и на море, 2005, № 10. С. 35-39.

6. Исаев В.И., Иванников В.Г., Иванников A.B. Внедрение струи газа в слой неньютоновской жидкости// НТЖ. Строительство нефтяных и газовых скважин на суше и на море, 2006, № 2. С. 29-32.

7. Исаев В.И., Леонов Е.Г., Райкевич С.И. Гидростатическое давление двухфазных сред при бурении и эксплуатации скважин// НТЖ. Строительство нефтяных и газовых скважин на суше и на море, 2001, № 1. С. 21-22.

8. Исаев В.И., Фёдоров В.В. Определение функции истинного газосодержания в горловине струйного насоса по данным опытов// Управление качеством в нефтегазовом комплексе, 2004, № 3-4. С. 80-84.

9. Исаев В.И., Фёдоров В.В. Оценка параметров течения многофазных жидкостей в процессе освоения скважин// НТЖ. Строительство нефтяных и газовых скважин на суше и на море, 2005, № 9. С. 43-50.

10. Исаев В.И., Шуть К.Ф. Определение реологических характеристик торфосодержащих растворов// НТЖ. Строительство нефтяных и газовых скважин на суше и на море, 2004, № 12. С. 37-46.

11. Исаев В.И., Шуть К.Ф. Определение подачи насосов для обеспечения выноса шлама при ламинарном течении жидкости, подчиняющейся 3-параметрической реологической модели// НТЖ. Строительство нефтяных и газовых скважин на суше и на море, 2005, № 7. С. 15-17.

12. Басниев К.С., Розенберг Г.Д., Исаев В.И. и др. Карачаганакское месторождение: расчёт свободного дебита скважин// Газовая промышленность, 1992, № 8. С. 30-32.

13. Басниев К.С., Розенберг Г.Д., Исаев В.И., Слободкина Ф.А., Кульпина Н.М., Астрахан И.М. Прогнозирование аварийного фонтанирования скважин при освоении нефтегазовых ресурсов Арктики в целях оценки возможных экологических последствий// НТЖ. Геология, геофизика и разработка нефтяных месторождений, 1996, № 12. С. 31-34.

14. Басниев К.С., Розенберг Г.Д., Исаев В.И., Кульпина Н.М. Расчет эксплуатационных и аварийных режимов в вертикальных и горизонтальных газовых и газоконденсатных скважинах// НТЖ. Геология, геофизика и разработка нефтяных месторождений, 1997, № 1. С. 38-41.

15. Брюховецкий О.С., Степанов Ю.С., Исаев В.И. Математическая модель двухфазных смесей при гидротранспорте// Известия вузов. Геология и разведка, 1990, № 1. С. 59.

16. Исаев В.И., Леонов Е.Г. Определение дебита газовых скважин, аварийно фонтанирующих через слой жидкости// Газовая промышленность, 1976, №9. С. 22-24.

17. Исаев В.И., Леонов Е.Г. Расчет давления по глубине обводненной газовой скважины// Газовая промышленность, 1974, № 4. С. 54.

18. Леонов Е.Г., Исаев В.И., Лукьянов И.П. Теория и методика расчёта расхода бурового раствора различной реологии для очистки от шлама ствола наклонно-направленных скважин// НТЖ. Строительство нефтяных и газовых скважин на суше и на море. М.: ОАО "ВНИИОЭНГ", 2006, № 8. С. 24-31.

19. Петров В.А., Леонов Е.Г., Филатов Б.С., Исаев В.И. Исследование газовых фонтанов и графический метод расчета глушения// Газовая промышленность, 1974, № 8. С. 53-54.

Учебники и учебные пособия для высших учебных заведений

20. Леонов Е.Г, Исаев В.И. Гидроаэромеханика в бурении. Учебник для вузов (гриф министерства высшего и среднего специального образования СССР). М.: Недра, 1987. 304 с.

21. Леонов Е.Г, Исаев В.И. Осложнения и аварии при бурении нефтяных и газовых скважин. Учебник для вузов (гриф министерства образования и науки РФ). В 2 частях. 2-е издание, исправленное и дополненное. Ч. I. Гидроаэромеханика в бурении. М.: ООО «Недра», 2006. 413 с.

22. Исаев В.И., Марков O.A. Управление скважиной. Предупреждение и ликвидация газонефтеводопроявлений (гриф УМО вузов РФ по нефтегазовому образованию). Уч. пособие. М.: ФАЗИС, 2006.134 с.

23. Исаев В.И., Марков O.A. Управление скважиной. Предупреждение и ликвидация газонефтеводопроявлений (гриф УМО вузов РФ по нефтегазовому образованию). Издание 2-е, исправленное и дополненное. Уч. пособие. М.: ФАЗИС, 2007. 146 с.

24. Исаев В.И., Митюшин А.И., Сумбатова А.Р. Лабораторный практикум по изучению процессов истечения газовых струй в жидкость. Уч. пособие. М.: РГУ нефти и газа им. И.М. Губкина, 1999. 50 с.

25. Иванников В.Г., Исаев В.И. Лабораторный практикум по технической гидромеханике. Уч. пособие. М.: ГАНГ им. И.М. Губкина, 1996. 111 с.

26. Астрахан И.М., Исаев В.И., Розенберг Г.Д., Кравченко М.Н. Гидродинамические задачи нефтегазопромысловой экологии. Уч. пособие. М.: РГУ нефти и газа им. И.М. Губкина, 1999. 151 с.

Работы, опубликованные в журналах, которые не входят в список ВАК

27. Исаев В.И. Экспериментальная установка и методика исследования дебита газового фонтана по параметрам газоводяного бугра на поверхности жидкости в кратере// В сб. "Нефть и газ", МИНХ и ГП им. И.М. Губкина, 1972. С. 108-109.

28. Исаев В.И. Гидродинамические проблемы бурения горизонтальных скважин// Геология, бурение, разведка и эксплуатация газовых и газоконден-сатных месторождений. М.: ИРЦ Газпром. Специальный выпуск, 1994, № 14, С. 29-42.

29. Астрахан И.М., Исаев В.И., Слободкина Ф.А. Сопряжённая задача теплообмена при аварийном фонтанировании газоконденсатной скважины в условиях вечной мерзлоты// В сб. «Нефтегазовая гидромеханика» М.: ФГУП. Изд-во "Нефть и газ" РГУ им. И.М. Губкина, 2005. С. 116-125

30. Астрахан И.М., Розенберг Г.Д., Исаев В.И. и др. Расчет эксплуатационных и аварийных режимов вертикальных и горизонтальных скважин при различных условиях вскрытия и различной криологической обстановке// НТЖ. Техника и технология добычи нефти, 1995, № 7. С. 8-21.

31. Бондаренко В.В., Бакшутов B.C., Ангелопуло O.K., Исаев В.И., Леонов Е.Г., Овчинников Н.Т. Крепление скважин пеноцементными тампонажными растворами // Серия "Строительство скважин", 1989 вып. 16. 60 с.

32. Басниев К.С., Розенберг Г.Д., Исаев В.И., Райский Ю.Д., Гиммер П.Р., Кульпина Н.М., Астрахан И.М., Кормишин А.К., Колесников А.Ф., Тимофеев А.Н. Термогидродинамическая модель расчета свободного дебита газоконденсатной скважины. Обзорная информация// Серия "Разработка и эксплуатация газоконденсатных месторождений". М.: ВНИИЭГазпром, 1991. 36 с.

33. Леонов Е.Г., Исаев В.И. Расчет потерь давления в циркуляционной системе скважины при бурении с промывкой аэрированной жидкостью ВИНИТИ, 2061-76. Деп. от 09.06.76.

34. Леонов Е.Г., Исаев В.И. Гидравлические расчеты промывки аэрированными жидкостями при бурении. М.: МИНХ и ГП им. И.М. Губкина, 1980. 48 с.

35. Леонов Е.Г., Исаев В.И. Гидравлические расчеты промывки скважин при бурении. М. : МИНХ и ГП им. И.М. Губкина, 1978.48с.

36. Леонов Е.Г., Исаев В.И. Определение критических чисел Рейнольдса для течения неньютоновских жидкостей в круглых трубах и кольцевых каналах. ВИНИТИ, №1152-78. Деп. от 04.04.78.

37. Леонов Е.Г., Исаев В.И. Уравнения гидромеханических задач буровых процессов. ВИНИТИ, №1329-80. Деп. 08.04.80.

38. Леонов Е.Г., Лукьянов И.П., Исаев В.И. Определение расхода жидкости разной реологии для транспорта единичных твердых частиц в наклонных трубах и кольцевых каналах. ВИНИТИ, №1091-В93. Деп. от 23.04.93.

39. Леонов Е.Г., Исаев В.И., Пономарев Ю.Н. Применение ЭВМ для гидравлических расчетов промывки скважин при бурении. М.: МИНХ и ГП им. И.М. Губкина, труды, выпуск 152,1980. С. 137-160.

40. Леонов Е.Г., Исаев В.И., Фишер В.А. Гидравлический расчет циркуляционной системы при бурении скважин. М.: МИНХ и ГП, "ПОМНИ, 1984. 39 с.

41. Петров В.А., Леонов Е.Г., Филатов Б.С., Исаев В.И. Графический метод расчета глушения открытых газовых фонтанов. Техника и технология геологоразведочных работ, организация производства, вып. 13, М.: ВИЭМС, 1974. С. 19-34.

42. Максимов В.М., Розенберг Г.Д., Исаев В.И., Лимар Е.Е. Теоретические основы прогнозирования последствий аварий на газо- и конденсатопро-водах// НТЖ. Защита от коррозии и охрана окружающей среды, 1995, № 5. С.

43. Бондаренко В.В., Бакшутов B.C., Ангелопуло O.K., Исаев В.И. и др. Способ тампонирования скважин вспененными растворами. А/с. N1521859, опубликованное в БИ №42,1989.

8-21.

Соискатель, E-mail: isa@gubkin.ru

В.И. Исаев

Подписано в печать: 28.10.2009

Заказ № 2850 Тираж - 120 экз. Печать трафаретная. Типография «11-й ФОРМАТ» ИНН 7726330900 115230, Москва, Варшавское ш., 36 (499) 788-78-56 www.autoreferat.ru

Содержание диссертации, доктора технических наук, Исаев, Валерий Иванович

Основные обозначения и принятые сокращения.

Введение.

1. Обзор основных моделей движения двухфазных смесей.

2. Основные задачи гидродинамики двухфазных смесей в бурении.

3. Одномерные уравнения движения двухфазных смесей в элементах циркуляционной системы скважин.„.

3.1. Уравнения неустановившегося течения смесей в бурильных трубах и кольцевых пространствах.

3.2. Уравнения установившегося течения смесей в бурильных трубах и кольцевых пространствах.

4. Обобщённая гидростатика однофазных и многофазных ньютоновских и неньютоновских жидкостей.

4.1. Основные уравнения гидростатики в поле силы тяжести.

4.1.1. Гидростатика однофазных жидкостей.

4.1.2. Гидростатика несжимаемой жидкости при тст =0.

4.1.3. Гидростатика однофазной сжимаемой жидкости газа) при тст =0.

4.1.4. Гидростатика слабосжимаемой жидкости при тст =0.

4.1.5. Гидростатика жидкости с динамическим напряженим сдвига (то Ф 0).

4.1.6. Гидростатика двухфазных флюидов.

4.2. Классификация гидростатик.

4.3. Капиллярные явления.

4.4. Применение законов гидростатики в технологиях строительства скважин.

4.4.1. Некоторые стандартные примеры расчётов давления в скважине.

4.4.2. Методика и расчёт одноступенчатого цементирования скважин стабильными пеноцементными растворами.

5. Установившиеся движения однофазных смесей в элементах циркуляционной системы бурящейся скважине.

5.1. .Уравнения одномерного установившегося движения однофазных смесей в элементах циркуляционной системы скважины.

5.2. Ламинарное течение ньютоновской несжимаемой жидкости в элементах циркуляционной системы скважины.

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

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

5.5. Турбулентное течение жидкости в трубах и концентрических кольцевых каналах.

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

6. Установившиеся течения газожидкостных смесей в элементах циркуляционной системы скважины.

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

6.2. Ламинарное и турбулентное восходящие течения газожидкостных вязкопластических смесей в циркуляционной системе скважины (в БТ и КП).

6.3. Движение газожидкостной смеси в местных сопротивлениях циркуляционной системы скважины.

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

6.5. Методика расчёта подачи и давления насосов и компрессоров при бурении с промывкой аэрированной жидкостью.

6.6. Влияние растворимости газа в жидкости на давление смеси в скважине.

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

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

7. Неустановившиеся движения двухфазных жидкостей в циркуляционной системе скважины.

7.1. Одномерные уравнения неустановившихся двухфазных течений смесей в циркуляционной системе скважины.

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

7.3. Неустановившееся течение газожидкостной смеси при бурении на депрессии.

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

8. Барботаж и струйные течения.

8.1. Исследование течения «нулевая подача».

8.2. Истечение струй газа через слой жидкости.

Введение Диссертация по наукам о земле, на тему "Гидродинамика двухфазных смесей в процессах бурения нефтяных и газовых скважин"

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

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

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

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

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

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

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

Цели диссертационной работы

Основными целями работы являются:

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

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

- совершенствование одномерных моделей течения двухфазных смесей в элементах циркуляционной системы скважины;

- экспериментальное исследование процесса истечения газовых струй через слой НЖ и ННЖ;

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

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

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

2. Обобщены законы гидростатики ньютоновских, неньютоновских жидкостей и многофазных смесей из них.

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

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

5. Дано обоснование перехода от ламинарного течения к турбулентному при движении вязкопластической жидкости (ВПЖ) в трубах.

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

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

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

Достоверность полученных результатов

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

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

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

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

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

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

Основные результаты диссертационной работы докладывались, демонстрировались и получили одобрение на 5-ой и 6-ой Всесоюзных конференциях по теплообмену и гидравлическим сопротивлениям при движении двухфазного потока в элементах энергетических машин и аппаратов (Ленинград, 1974,1978); на 6-ом Всесоюзном семинаре по гидравлике промывочных жидкостей и тампо-нажных растворов (Астрахань, 1978); на Международной конференции (21st petroleum conference and exhibition. Hungary, Siofok, 1990); на 2-ом Всесоюзном школе-семинаре "Разработка месторождений нефти и газа: современное состояние, проблемы, перспективы" (Звенигород, 1991); на Всесоюзной конференции "Проблемы развития нефтегазового комплекса страны" (Красный Курган, 1991); на научно-технической конференции "Проблемы экологически безопасных технологий разведки, разработки и эксплуатации глубокопогруженных месторождений со сложным составом пластовой смеси" (Астрахань, 1991); на научных семинарах "Актуальные проблемы нефтегазовой и подземной гидромеханики" (ГАНГ им. И.М. Губкина, 1992, 1993); на международном семинаре "Горизонтальные скважины" в ГАНГ им. И.М. Губкина (Москва, 1993); на научном семинаре по гидромеханике института механики МГУ под руководством проф. Куликовского А.Г. (Москва, 1993); на научных семинарах по подземной гидромеханике под руководством проф. Басниева К.С.; на научном семинаре кафедры бурения нефтяных и газовых скважин под руководством д.т.н. проф. Ангелопуло O.K.; на научных семинарах академического института проблем нефти и газа (ИПНГ); на научно-технической конференции, посвященной 70-летию первого выпуска российских инженеров-нефтяников (Москва, 1994); на 1-ой российской национальной конференции по теплообмену (Красногорск, 1994); на 4-ой международной деловой встрече "Диагностика-94" (Ялта, 1994); на выставке "Научные технические достижения и интеллектуальная собственность высшей школы, организованной Госкомвузом РФ (Москва,1994); на международном симпозиуме "Методы и средства мониторинга состояния окружающей среды МСОС-95" (Санкт-Петербург, 1995); на заседании кафедры "Гидравлика и гидрофизические процессы геологоразведочных и горных работ" Московской государственной геологоразведочной академии; на научно-технической конференции "Геоэкология в. нефтяной и газовой промышленности" (Москва, ГАНГ им. И.М. Губкина, 1995); на конференции МГРИ им. Серго Ордженикидзе академии (Москва, 1996); на международной конференции по исследованию газа (International Gas Research Conference, Cannes, France, 1995); на всероссийской научной конференции " Фундаментальные проблемы нефти и газа" (Москва, 1996); на 2-ой научно-технической конференции "Актуальные проблемы состояния и развития нефтегазового комплекса России" (Москва, ГАНГ им. И.М. Губкина, 1997); на международном симпозиуме по вопросам теории и практики механики сплошных и сыпучих сред (Москва, 1997); на конференции "Проблемы разработки газовых и газоконденсатных месторождений" (Москва, ГАНГ им. И.М. Губкина, 1998); на 2-ом международном симпозиуме "Актуальные проблемы механики сплошных и сыпучих сред" (Москва, 1999); на 3-ей научно-технической конференции, по-свящённой 70-летию РГУ нефти и газа им. И.М. Губкина "Актуальные проблемы состояния и развития нефтегазового комплекса России" (Москва, РГУ,1999); на 4-ой научно-технической конференции, посвящённой 300-летию инженерного образования в России (Москва, РГУ, 2001); на 8-ом Всероссийском съезде по теоретической и прикладной механике (Пермь, 2001); на международной конференции по газовым технологиям (Амстердам (Нидерланды), 2002); на 5-ой научно-технической конференции "Актуальные проблемы состояния и развития нефтегазового комплекса России" (Москва, РГУ им. И.М. Губкина, 2003); на 6-ой научно-технической конференции, посвящённой 75-лет РГУ нефти и газа им. И.М. Губкина (2005 г.); на 9-ом Всероссийском съезде по теоретической и прикладной механике (Н. Новгород, 2006г.); на 7-ой Всероссийской конференции "Актуальные проблемы состояния и развития нефтегазового комплекса России" (Москва, РГУ нефти и газа им. И.М. Губкина 2007); на Всероссийской конференции, посвящённой 20-летию ИПНГ РАН (Москва, ИПНГ, 2007); на международной (SOF-2007) научно-технической конференции (Москва, ВНИИГАЗ, 2007).

Основные защищаемые положения

1. Модель двухфазной гидродинамики процессов бурения.

2. Обобщение закона гидростатики для ньютоновских и неньютоновских растворов и их многофазных смесей.

3. Методика расчёта управления скважиной при бурении на депрессии.

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

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

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

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

8. Формула для расчёта коэффициента гидравлических сопротивлений при турбулентном течении вязкопластического раствора в трубах.

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

Публикации

Основные материалы диссертации опубликованы в 43 печатных работах; из них 19 статей в журналах, рекомендованных ВАК для публикации материалов по докторским диссертациям; два учебника, выпущенных издательством "Недра", пять учебных пособия выпущены в других издательствах. Всего автором опубликовано 126 работ.

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

Диссертация состоит из введения, 8 глав, основных выводов и рекомендаций, библиографии, включающей 348 наименований работ отечественных и зарубежных авторов. Работа изложена на 331 странице, содержит 85 рисунков и 22 таблицы.

Заключение Диссертация по теме "Технология бурения и освоения скважин", Исаев, Валерий Иванович

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

Библиография Диссертация по наукам о земле, доктора технических наук, Исаев, Валерий Иванович, Москва

1. Абрамович Г.Н. Теория турбулентных струй. М.: Наука, 1984.-716 с.

2. Агаев Г.Ф. Движение газожидкостной смеси в скважине. Баку: Азнефтеиз-дат, 1966. 190 с

3. Агаев Г.Е. Исследование движения твердой частицы при роторном способе бурения// Нефть и газ, № 4,1977. С. 32 36.

4. Айзенбуд М.Б., Дильман В.В. Вопросы гидравлики химических реакторов для систем газ жидкость// Хим. пром., № 3, 1961. С. 51-56.

5. Акбулатов Т.О. Вынос частиц шлама из горизонтального ствола скважины // Нефть и газ, № 1, 2000. С. 34-38.

6. Алиев З.С., Власенко А.П., Андреев С.А., Абрамова Е.С. Определение критического дебита в газовых скважинах.// Газ. пром., №2, 1975. С. 27-30.

7. Аликбеков Б.И., Пирвердян A.M., Сафаров Р.А., Чубанов О.В. К расчёту воздушного подъёмника для вязкопластической жидкости// Нефтяное хозяйство, № 8, 1971. С. 44-46.

8. Аргунов П.П. Исследование работы эрлифта и его расчёт. Сб. «Строительное водопонижение, гидромеханика и физика грунтовых вод». Тр. Ин-та оснований и фундаментов. 1953.

9. Арманд А.А. Сопротивление двухфазного потока в горизонтальных трубах. Труды ВТИ, 15(1), 1946. С. 19-34.

10. Арманд А.А. Исследование механизма движения двухфазной смеси в вертикальной трубе. Сб. Гидродинамика и теплообмен в котлах высокого давления//Изв. АН СССР, 1955. С. 21-34.

11. Арманд А.А., Невструева Е.И. Исследование механизма движения^ двухфазной смеси в вертикальной трубе// Изв. ВТИ, №2, 1950i

12. Архангельский В.А. Движение-газированных жидкостей и газожидкостных смесей в вертикальных трубах. Инж. сб., т.4, вып. 2, 1949.

13. Архангельский В.А. Движение газированных нефтей в системе скважина пласт. М.: Изд-во АН СССР, 1958. 92с.

14. Астарита Д., Маруччи Д. Основы гидромеханики неньютоновских жидкостей. М.: Мир, 1978. 310 с.

15. Астрахан И.М., Егорушкин С.В., Исаев В.И., Розенберг Г.Д., Слободкина Ф.А. Моделирование аварийного фонтанирования газоконденсатной скважины в условиях вечной мерзлоты// Изв.РАН. МЖГ, № 3., 1994. С.101-108.

16. Бабалов М.А. Определение дебита горящей газовой струи по её шуму// Газ. пром., №5 , 1973. С. 12-15.

17. Багдасаров В.Г. Теория, расчёт в практике эргазлифта. M.-JL: Гостоп-техиздат. 1947. 234 с.

18. Базаров И.П. Термодинамика. М.: Высшая школа. 1991. 376 с.

19. Баранов В.Я. Учебный практикум по курсу "Физическая и коллоидная химия". М.: МИНХ и ГП, 1985. 94 с.

20. Басниев К.С., Кочина И.Н., Максимов В.М. Подземная гидромеханика. М.: Недра, 1993. 416 с.

21. Басниев К.С., Исаев В.И., Кульпина Н.М., Розенберг Г.Д. Расчёт эксплуатационных и аварийных режимов в вертикальных и горизонтальных газовых и газоконденсатных скважинах. // Геология, геофизика и разработка нефтяных месторождений. № 1, 1997. С. 38-41.

22. Басниев К.С., Розенберг Г.Д., Исаев В.И., Райский Ю.Д., Гимер П.Р., Кульпина Н.М. Карачаганакское месторождение: расчёт свободного дебита* скважин // Газовая пром. №8, 1992. С. 30-32.ч

23. Белей И.В. К вопросу о работе турбобуров на аэрированной жидкости. М.: МНП, ВНИИБТ, тр. Вып. XXX "Пути интесификации буровых работ", 1972 С. 28-31.

24. Берд Р., Стьюарт В., Лайтфут Е. Явления переноса. М.: Химия, 1974. 688 с.

25. Бондаренко В.В., Бакщутов B.C., Ангелопуло O.K., Исаев В.И., Леонов Е.Г., Овчинников Н.Т. Крепление скважин пеноцементными тампонажными раствора. М.: ВНИИОЭНГ, 1989. Вып. 16. 58 с.

26. Бондаренко В.В., Исаев В.И. Методика расчета цементирования глубоких скважин пеноцементными растворами // Нефтяник, №4, 1987. С. 18.

27. Бондаренко В.В., Исаев В.И. Пути повышения качества цементирования глубоких скважин в регионе. ВНИИЭГазпром, труды МИНГ им. И.М. Губкина, Деп. № 1023-ГЗ/25,25.

28. Бондаренко В.В., Ангелопуло O.K., Исаев В.И. и др. Технология цементирования скважин тампонажными пеноцементными растворами. ВНИИЭГазпром, Деп. №1064-ГЗ/88, от 29.06.88.

29. Бондаренко В.В., Ангелопуло O.K., Бакшутов B.C., Исаев В.И. Способ тампонирования скважин вспененными растворами. А/с. №1521859, БИ №42, 1989.

30. Брилл Дж.П., Мукерджи Ч. Многофазный поток в скважинах. Москва-Ижевск: Институт компьютерных исследований, 2006. 384 с.

31. Бронштейн И.Н., Семендяев К.А. Справочник по математике для инженеров и учащихся втузов. М.: Наука, 1986. 544 с.

32. Брюховецкий О.С. Гидравлическая модель потока гидросмеси для условий напорного транспортирования минерального сырья. М.: Геоинформ. 1990. 87 с.

33. Брюховецкий'О.С., Исаев В.И., Степанов Ю.С. Математическая модель, двухфазных смесей при гидротранспорте // Изв. вузов. Геология и разведка. №1,1991. С.28.

34. Бусройд Р. Течение газа со взвешенными частицами. М.: Мир, 1975. 378 с.

35. Бэтчелор Дж. Введение в динамику жидкости. М.: Мир, 1973. 758 с.

36. Васильев А.С., Талачёв B.C., Павлов В. П., Плановский А.Н. Закономерности истечения струй газа в жидкость // Теор. основы хим. технол. Т. 4, №5, 1970. С. 727-735.

37. Верслуис Д. Математическая теория фонтанирования нефтяных скважин // Нефтяное хозяйство. № 6, 1931.

38. Викторов М.М. Методы вычисления физико-химических величин и прикладные расчёты. Д.: Химия, 1977. 360 с.

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

40. Виноградов К.В. Движение газонефтяной смеси в фонтанных скважинах. М.: Недра, 1964.

41. Воларович М.П., Гуткин A.M. Течение пластично-вязкого тела между двумя параллельными плоскими стенками и в кольцевом пространстве между двумя коаксильными цилиндрами // ЖТФ. Т. XVI. Вып. 3.1946. С. 321-328.

42. Вольфсон В.И., Назаров А.А., Савельев А.С., Белей И.В. К вопросу истечения водовоздушной смеси через цилиндрические насадки. М.: МНП, ВНИИБТ, тр. Вып. XXXI "Осн. проект, глуб. скважин", 1973. С. 142-148.

43. IX Всероссийский съезд по теоретической и прикладной механике. Аннотации докладов. Т. II. Нижний Новгород: Изд-во Нижегородского госуниверситета им. Н.И. Лобачевского, 2006.

44. Вулис JI.A., Джаугастин К.Е., Живов В.Г., Ярин Л.П. Распространение; струй вязкой жидкости в среде со скачком плотности // ПМТФ, №3, 1972. С. 115-122.

45. Газодинамические и термодинамические исследования скважин месторождений Ставрополья. Отчёт по теме 25/64. Ставрополь, 1965.

46. Гейман М.А., Мусинов В.И. Турбинное бурение на аэрированной промывочной жидкости. М.: Недра, 1965. 146 с.

47. Гидродинамика и фильтрация однофазных и многофазных потоков. Тр. МИНХ и ГП им. И.М. Губкина, вып. 101. М.: Недра, 1972. 165 с.

48. Гимер П.Р. Расчет вязкости жидких углеводородных смесей. М.: ИРЦ Газпром. № 1336-гз 92 (Библ. указ.) "Деп. науч. раб." №11, 1992г. 55 е.

49. Гиршович Т.А. Турбулентные струи в поперечном потоке. М.: Машиностроение, 1993. 251 с.

50. Голеевский А.А. Вопросы механики струйного движения жидкостей и газов. М.: Машгиз, 1957. 88с.

51. Гоинс У.К., Шеффилд Р. Предотвращение выбросов: Пер. с англ. М.: Недра, 1987. 288 с.

52. Грей Дж.Р., Дарли Г.С.Г. Состав и свойства буровых агентов (промывочных жидкостей): пер. с англ. М.: Недра, 1985. 509 с.

53. Гришанин К.В. Основы динамики русловых потоков. М.: Транспорт, 1990. 320 с.

54. Гриценко А.И., Акопова Г.С., Максимов В.М. Экология, нефть и газ. М.: Наука, 1997. 600 с.

55. Гриценко А.И., Алиев З.С., Ермилов О.М., Ремизов В.В., Зотов Г.А. Руководство по исследованию скважин. М.: Наука, 1995. 523 с.

56. Гриценко А.И., Клапчук О.В., Харченко Ю.А. Гидродинамика газожидкостных смесей в скважинах и трубопроводах. М.: Недра, 1994. 238 с.

57. Гродзовский Г.Л. О движении мелких частиц в газовом потоке. //Труды. ЦАГИ,- вып. 1155, 1976. С. 3-12.

58. Грязнов Г.С. Конструкция газовых скважин в районах многолетнемёрз-лых пород. М.: Недра, 1978. 137 с.

59. Губин В.Е., Губин В.В. Трубопроводный транспорт нефти и нефтепродуктов. М.: Недра, 1982. 296 с.

60. Гукасов Н.А. Гидравлика газожидкостных смесей в бурении и добыче нефти: Справочное пособие. М.: Недра, 1988. 237 с.

61. Гукасов Н.А. Гидродинамика при креплении скважины. М.: Недра, 1976.

62. Гукасов Н.А. Механика жидкости и газа. М.: Недра, 1996. 445 с.

63. Гукасов Н.А., Кучеров Г.Г. Теория и практика добычи газожидкостных смесей. М.: ООО «ИРЦ Газпром», 2005. 307 с.

64. Гурбанов Р.С., Бахтияров С.И. Обтекание полого цилиндрического тела в вертикальной колонне труб неньютоновской жидкостью // Нефть и газ. №1, 1976. С. 51-53.

65. IX (девятый) Всероссийский съезд по теоретической и прикладной механике. Аннотации докладов. Т. I III. (Нижний Новгород, 22 - 28 августа 2006). Нижний Новгород: Изд-во Нижегородского госуниверситета им. Н.И. Лобачевского, 2006. С. 178, с. 194, с. 233.

66. Де Гроот С., Мазур П. Неравновесная термодинамика. М.: Мир, 1964: 456 с.

67. Дейч М.Е., Филиппов Г.А. Двухфазные течения в элементах энергетического оборудования. М.: Энергоатомиздат, 1987. 328 е.

68. Дерягин Б.В., Чураев Н.В., Муллер В.М. Поверхностные силы. М.: Наука, 1985. 398fc.

69. Достижения в области исследования- теплообмена и» гидравлики двухфазных потоков в элементах энергооборудования. Л. г Наука, Ленингр. отд., 1973.292 с.

70. Дубина М.М., Красовицкий Б.А. Теплообмен и механика взаимодействия трубопроводов и скважин с грунтами. Новосибирск: Наука, 1983. 133 с.

71. Дюнин А.К., Борщевский Ю.Т., Яковлев Н.А. Основы механики многокомпонентных потоков. Новосибирск: Редакционно-издательский отдел сибирского отделения АН СССР. 1965. 76 с.

72. Дядькин Ю.Д., Гендлер С.Г., Смирнова Н.Н. Геотермальная теплофизика. С.-Петербург: Наука, 1993. 256 с.

73. Елин Н.Н. Метод расчета коэффициента сопротивления газожидкостного потока. ВИНИТИ, Деп. №3621-В88-1988. 13 с.

74. Есьман Б.И., Габузов Г.Г. Термогидравлические процессы при бурении скважин. М.: Недра, 1991. 216 с.

75. Задара В.М. Применение методов теории подобия к исследованию движения газожидкостных потоков в горных выработках.// Изв. вузов. Горный журнал. №7, 1974. С. 102-105.

76. Задворных В.Н. Реодинамика нелинейно-вязкопластических буровых растворов в кольцевом пространстве глубокой скважины. Дисс. на соиск. уч. ст. к.т.н. М.: 1986.

77. Зайцев А.И., Бытев Д.О. Ударные процессы в дисперсно- пленочных системах. М.: Химия, 1994. 177 с.

78. Заря А.Н. Концентрация гидросмеси // Изв. вузов. Гор. ж. № 9,1984. С. 57- 60.

79. Инструкция по исследованию газовых скважин. М.: Недра, 1961. 110 с.

80. Иванников А.В. Экспериментальное исследование истечения газожидкостной струи через слой жидкости. Дисс. на соискание уч. ст. к.т.н. (научный рук. Исаев В.И.). М.: РГУ им. И.М. Губкина, 2007. 141 с.

81. Иогансен К.В. Спутник буровика. М: Недра, 1986. 294 с.

82. Исаев В.И. Экспериментальная установка и методика исследования дебита газового фонтана по параметрам газоводяного бугра на поверхности жидкости в кратере. / Сб. "Нефть и газ". М.: МИНХ и ГП, 1972. С. 108-109.

83. Исаев В.И. Гидродинамика горизонтального бурения скважин. Доклад на Всероссийской научной конференции "Фундаментальные проблемы нефти и газа", Москва, 22-25 января, 1996.

84. Исаев В.И. Гидродинамические проблемы бурения горизонтальных скважин. Серия: геология, бурение, разработка и эксплуатация газовых и га-зоконденсатных месторождений. Спец. выпуск, №1-4. М.: ИРЦ Газпром, 1994. С. 29-42.

85. Исаев В.И. Термогидродинамические модели многофазных течений при строительстве и эксплуатации скважин// Управление качеством в нефтегазовом комплексе. №1-2, 2005. С. 100-103.

86. Исаев В.И. Обобщённая гидростатика однофазных флюидов и многофазных смесей в поле силы тяжести// НТЖ. Строительство нефтяных и газовых скважин на суше и на море. №4, 2005. С. 26-32

87. Исаев В. И. Некоторые задачи движения» газожидкостных смесей при-бурении скважин. Дисс. на соиск. уч. ст. к.т.н. М.г: МИНХ и ГП им. И:М. Губкина, 1976. 130 с.

88. Исаев В.И. Метод расчета глушения открытого газоконденсатного фонтана. Тезисы докладов научно-технической конференции "Геоэкология в нефтяной и газовой промышленности". Москва, ГАНГ им. И.М. Губкина, 1995. С. 38.

89. Исаев В.И. Определение расхода промывочной жидкости для выноса шлама в наклонных скважинах. Тезисы докладов науч. конференции посвященной 75-лет МГА-МГГА, Москва, 1993. С.67.

90. Исаев В.И., Иванников В.Г., Иванников А.В. Оценка параметров газопроявлений в виде грифонов и экологической обстановки окружающей среды //Управление качеством в нефтегазовом комплексе. №1, 2004. С. 57 59.

91. Исаев В.И., Иванников В.Г., Иванников А.В. Экспериментальное моделирование аварийного фонтанирования газовой скважины через слой жид-кости//НТЖ. Строительство нефтяных и газовых скважин на суше и на море.- М.: ОАО «ВНИИОЭНГ», №10, 2005. С. 35-39.

92. Исаев В.И., Иванников В.Г., Иванников А.В., Шуть К.Ф. Внедрение' струи газа в слой неньютоновской жидкости // НТЖ. Строительство нефтяных и газовых скважин на суше и на море.- М.: ОАО «ВНИИОЭНГ», №2, 2006. С. 29-32.

93. Исаев В.И., Леонов Е.Г. Гидравлический расчет при промывке скважин аэрированной жидкостью/ТВопросы гидродинамики вязкой и вязкопласти-ческой жидкости. Рязань: 1976. С. 46-61.

94. Исаев В.И., Леонов Е.Г. Определение дебита газовых скважин, аварийно-фонтанирующих через слой жидкости// Газовая промышленность. №9, 1976. С. 22-24.

95. Исаев В.И., Леонов Е.Г., Райкевич С.И. Гидростатическое давление двухфазных сред при бурении и эксплуатации скважин// НТЖ. Строительство нефтяных и газовых скважин на суше и на море. №1, 2001. С. 21-22

96. Исаев В.И, Марков О.А. Управление скважиной. Предупреждение и ликвидация газонефтеводопроявлений. Уч. пособ. (Гриф УМО вузов РФ). М.: Фазис, 2006. 135 с.

97. Исаев В.И., Марков О.А. Управление скважиной. Предупреждение и ликвидация газонефтеводопроявлений. Уч. пособ. (Гриф УМО вузов РФ). 2-е изд. испр. и дополненное. М.: Фазис, 2007. 146 с.

98. Исаев В.И., Матвиенко А.В. Метод расчета перепада давления нелиней-но-вязкопластических жидкостей в кольцевом пространстве. ВИНИТИ, Деп. №1092-В93, от 23.04.93.

99. Исаев В.И., Митюшин А.И., Сумбатова А.Р. Лабораторный практикум по изучению процессов истечения газовых струй в жидкость: Учебное пособие. М.: РГУ нефти и газа им. И.М. Губкина. 1999. 49 с.

100. Исаев В.И., Самохвалов С.Ю. Особенности решения задачи о течении нелинейно-вязких жидкостей в кольцевых каналах. ВИНИТИ, Деп. №4631-В89 от 12.07.89

101. Исследование турбулентных течений двухфазных сред. Под ред. С.С. Кутателадзе. Новосибирск, изд. ИТФ СО АН СССР, 1973. 198 с.

102. Калинин А.Г., Никитин Б.А., Солодкий К.М., Повалихин А.С. Профили направленных скважин и компоновки низа бурильных колонн. М.: Недра, 1995.303 с.

103. Камнер А.А., Яхно О.М. Гидромеханика в инженерной практике. Киев: Техника, 1987. 175 с.

104. Карамышев В.Г, Корнилов; Г.Г. О связи параметров; газожидкостных потоков.//Нефтяное хозяйство. № 4, 1999; С. 37-39;

105. Катц Д.Л., Корнелл Д., Кобаяши Р., Поеттман Ф.Х., Вери Дж.А., Еленба-ас Дж.Р., Уайнауг Ч.Ф. Руководство по добыче, транспорту и переработке природного газа. (Перевод с англ. 1959). М.: Недра, 1965. 676 с

106. Колесников Н.А., Шарипов А.У. Гидромеханические способы повышения технико-экономических показателей бурения. М.: ВНИИОЭНГ, 1989. 50 с.

107. Колтюбинг. Технологическое приложение №1 к журналу «Нефть и Капитал». М.: 2001, 40 с.

108. Костерин С.И. Исследование структур течения двухфазных сред в горизонтальных трубах // Изв. АН СССР, ОТН, №7, 1943. С. 37-45.

109. Коротаев Ю.П. Комплексная разведка и разработка газовых месторождений. М.: Недра, 1968. 290 с.

110. Крылов А.П. Расчёт подъёмников для эксплоатации компрессорных и фонтанных скважин // Нефтяное хозяйство. №2, 1934. С. 24 33.

111. Крылов А.П. Потери трения и скольжения при движении жидкости и газа по вертикальным трубам // Нефтяное хозяйство. №6, 1935. С. 35 42.

112. Крылов В.И., Крецул В.В. Особенности технологии промывки горизонтальных скважин // Нефтяное хозяйство, № 6, 2001. С. 24 — 36 и № 7, 2001. С. 20-24.

113. Куропатенко В.Ф. Модель многокомпонентной среды// Доклады Академии наук. Т. 403, №6. С. 761 763.

114. Киршин В.И. Научно-методические и технологические решения по строительству скважин в условиях депрессии с использованием колпобинговых установок. Автореферат диссертации, на соискание уч. ст. к.т.н. М.; 2006.24 с.

115. Кутателадзе С.С., Стырикович М.А. Гидравлика газожидкостных систем. M.-JL: Госэнергоиздат, 1958. 180 с.

116. Ландау JI.Д., Лифшиц Е.М. Теоретическая физика. В' 10 т. Т.VI1. Гидродинамика. 4-е изд., М.: Наука. Гл. ред. физ.-мат. лит., 1988. 736 с.

117. Левин В.Г. Физико-химическая гидродинамика. М.: Физматгиз, 1959. 699 с.

118. Леонов Е.Г. Методика расчета производительности насосов и компрессоров при бурении скважин с промывкой аэрированной жидкостью // РНТС ВНИИОЭНГ. Бурение. № 8, 1973. С. 5-9.

119. Леонов Е.Г., Исаев В.И. Гидравлические расчеты промывки скважин при бурении. М.: МИНХ и ГП им.И.М. ГубкинаД978. 51 с.

120. Леонов Е.Г., Исаев В.И. Гидравлические расчеты промывки аэрированными жидкостями при бурении. М.: МИНГ и ГП, 1980. 49 с.

121. Леонов Е.Г., Исаев В.И. Гидроаэромеханика в бурении. М.: Недра, 1987. 304 с.

122. Леонов Е.Г., Исаев В.И. Методика подбора диаметров насадок гидромониторных долот при роторном и турбинном способе бурения. ВИНИТИ, 3677-77, Деп. от 19.09. 1977.

123. Леонов Е.Г., Исаев В.И. О переходных режимах течения неньютоновских жидкостей в круглых и кольцевых каналах. Тез. докл. 6-ого Всес. семинара по гидравлике промывочных жидкостей и тампонажных растворов. Астрахань, 1978. С. 41-42.

124. Леонов Е.Г., Исаев В.И. Определение критических чисел Рейнольдса'; для, течения неньютоновских жидкостей в круглых трубах и кольцевых каналах. ВИНИТИ; №1152-78, Деп. от 04.04. 1978.

125. Леонов Е.Г., Исаев В.И. Осложнения и аварии при бурении нефтяных и газовых скважин. Учеб. для вузов (Гриф мин. обр. и науки РФ). В 2 частях.- 2-е изд., испр. и доп. Ч. I. Гидроаэромеханика в бурении. Изд-во ООО «Недра», 2006. 413 с.

126. Леонов Е.Г., Исаев В.И. О фильтрации промывочных и тампонажных жидкостей. Тр. МИНХ и ГП им. И.М. Губкина, вып. 162, 1982.

127. Леонов Е.Г., Исаев В.И. Расчет гидродинамического давления при СПО (учебное пособие). МИНХ и ГП им. И.М. Губкина, 1982. 60 с.

128. Леонов Е.Г., Исаев В.И. Расчет давления газированной жидкости на приеме глубинного насоса// Нефтепромысловое дело, Эи, ВНИИОЭНГ, №16, 1973. С. 24

129. Леонов Е.Г., Исаев В.И. Расчет давления по глубине обводненной газовой скважины // Газовая промышленность, №4, 1974. С. 35.

130. Леонов Е.Г., Исаев В.И. Расчет потерь давления в циркуляционной системе скважины при бурении с промывкой аэрированной жидкостью. ВИНИТИ, №2061-76 Деп. от 9 июня 1976.

131. Леонов Е.Г., Исаев В.И. Определение дебита газовых скважин, аварийно фонтанирующих через слой жидкости. Газовая промышленность, №9, 1976. С.38.

132. Леонов Е.Г., Исаев В.И. Течение пластовых флюидов и горных пород. ВИНИТИ, №5608-81, Деп. от 09.12.81.

133. Леонов Е.Г., Исаев В.И. Уравнения гидромеханических задач' буровых процессов. ВИНИТИ, №1329-80, Деп. 08:04.80.

134. Леонов. Е.Г., Исаев^В.И., Самохвалов С.Ю. Оценка влияния шлама на^ давление бурового раствора в затрубном пространстве скважины при промывке. ВИНИТИ, №5019-81, Деп. от 02.11.81.

135. Леонов Е.Г., Исаев В.И., Лукьянов И.П., Жолмагамбетов Т. Выбор режима дегазации в ДВС-3 // Машины и нефтяное оборудование, 5, 1982. С. 28-32.

136. Леонов Е.Г., Исаев В.И., Лукьянов И.П. Определение расхода жидкости различной реологии для транспорта единичных твердых частиц в наклонных трубах и кольцевых каналах. ВИНИТИ, Деп. №1091-93, от 23.04.93.

137. Леонов Е.Г., Исаев В.И., Пономарёв Ю.Н. Применение ЭВМ для гидравлических расчётов промывки скважин при бурении. В кн.: Технология бурения нефтяных и газовых скважин. Тр. МИНХ и ГП им. И.М. Губкина, вып 152, 1980. С. 137-160.

138. Леонов Е.Г., Исаев В.И., Фишер В.А. Гидравлический расчёт циркуляционной системы при бурении скважин. М.: МИНХ и ГП им. И.М. Губкина, 1984. 40 с.

139. Липатов В.И., Шумилов Л.П. О движении твердых частиц в восходящем потоке раствора. Тр. ВНИИБТ, вып.24, 1970. С 74 79.

140. Литвинов В.Г. Движение нелинейно-вязкой жидкости. М.: Наука, 1982. 374 с.

141. Литвинов А.И. Аналитическое определение темпа нагнетания жидкости при глушении открытых газовых фонтанов // Изв. АН Узб.ССР. Сер. тех:, наук. Добыча и транспорт топлива. №5, 1974. С.58-60.

142. Лойцянский Л.Г. Механика жидкости и газа. М.: Наука. 1987. 720 с.

143. Лурье. М.В., Подоба. Н.А. Модификация* теории Кармана для, расчёта; сдвиговой турбулентности // ДАН СССР. № 10, 1984. С. 570-575.

144. Лутошкин Г.С. Исследование влияния вязкости жидкости и поверхностного натяжения системы «жидкость-газ» на работу эргазлифта. Дисс. на со-иск. уч. ст. к.т.н. ВНИИ: 1955.

145. Лутошкин Г.С., Беленький В.Н., Никитина Н.И. Исследование движения газожидкостных смесей по затрубному пространству. Тр. ВНИИ, вып. 13. М.: Гостоптехиздат, 1958. С. 33-41.

146. Лямаев Б.Ф., Небольсин Г.П., Нелюбов В.А. Стационарные и переходные процессы в сложных гидросистемах. Методы расчёта на ЭВМ. Л.: Ма- . шиностроение. Ленингр. отделение, 1978. 192 с.

147. Максимов В.М. Основы гидротермодинамики пластовых систем. М.:. Недра, 1994. 201 с.

148. Маковей Н. Гидравлика бурения. М.: Недра, 1986. 536 с.

149. Малеванский В.Д. Открытые газовые фонтаны и борьба с ними. М.: Гостоптехиздат, 1963. 212 с.

150. Малеванский В.Д., Шеберстов Е.В. Гидродинамические расчёты режимов глушения фонтанов в нефтяных и газовых скважинах. М.: Недра, 1990. 246 с.

151. Мамаев В.А. и др. Методические указания по расчёту трубопроводов,, транспортирующих нестабильный конденсат в двухфазном состоянии. М.: ВНИИГАЗ, 1974. 146 с.

152. Мамаев В.А., Одишария Г.Э., Семенов Н.И., Точигин А.А. Гидродинамика газо-жидкостных смесей в трубах. М.: Недра, 1969.208 с.

153. Мамаев В.А., Одишария Г.Э., Клапчук О.В.,. Точигин А.А., Семёнов Н.И. Движение газожидкостных смесей в трубах. М.: Недра, 1978. 270 с.

154. МедниковЕ.П. Турбулентный перенос и осаждение аэрозолей. М.: Наука, 1981. 174 с.170: Межлумов А.О. Использование'аэрированных жидкостей" при-проводке скважин. М.: Недра, 1976. 231 с.

155. Мехтиев Е.Х. Об определении гидродинамического давления аэрированной промывочной жидкости на забой. Азербайджанское нефтяное хозяйство, № 11, 1966. С. 20 -22.

156. Мирзаджанзаде А.Х. Вопросы гидродинамики вязко-пластичных и вязких жидкостей в нефтедобыче. Баку: Азернефтнешр, 1959. 409с.

157. Мирзаджанзаде А.Х., Ентов В.М. Гидродинамика в бурении. М.: Недра, 1985. 196 с.

158. Мирзаджанзаде А.Х., Ширинзаде С.А. Повышение эффективности и качества бурения глубоких скважин. М.: Недра, 1986. 278 с.

159. Мищенко И.Т. Скважинная добыча нефти. Учебное пособие для вузов. М.: ФГУП Изд-во «Нефть и газ» РГУ нефти и газа им. И.М. Губкина, 2003. 816 с.

160. Мукминов Р.А., Мавлютов М.Р. Расчет циркуляционной системы при бурении с очисткой забоя скважины аэрированной жидкостью. Нефть и газ, №4, 1965. С. 21 -26.

161. Муравьёв И.М., Крылов А.П. Эксплуатация нефтяных месторождений. M.-JI. Гостоптехиздат. 1949.

162. Муравьёв И.М., Репин Н.И. Исследование движения многокомпонентных смесей в скважинах. М.: Недра. 1972.

163. Накоряков В.Е., Покусаев Б.Г., Шрейбер И.Р. Волновая динамика газо- и парожидкостных сред. М.: Энергоатомиздат, 1990. 248 с.

164. Научные основы прогноза фазового поведения пластовых газоконден-сатных систем / А.И. Гриценко, И.А. Гриценко, В.В. Юшкин, Т.Д. Островская. М.: Недра, 1995. 432 с.

165. Нигматулин Р.И. Основы механики гетерогенных сред. М: Наука. 1978. 336 с.

166. Нигматулин Р.И. Динамика многофазных сред. М.: Наука. Т.1. 1987. 464 с.

167. Нигматулин Р.И. Динамика многофазных сред. М.: Наука. Т.2. 1987г. -359 с.

168. Николаев А.Ф., Охрименко Г.И. Водорастворимые полимеры. JL: Химия, 1979. 144 с.

169. Николаевский В.Н. Механика пористых и трещиноватых сред. М.: Недра, 1984. 232с.

170. Николаевский В.Н. Геомеханика и флюидомеханика. М.: Недра, 1996. 447 с.

171. Одишария Г.Э. Разработка и внедрение технологии магистрального транспорта нестабильного конденсата и природного газа при низких температурах. Дис. на соиск. д.т.н. М.: МИНХ и ГП им. И.М.Губкина, 1980.

172. Островский Г.М. Прикладная механика неоднородных сред. СПб, Наука, 2000. 359 с.

173. Пеньков А.И., Никитин Б.А., Филиппов Е.Ф. Методы регламентирования свойств буровых растворов для горизонтальных скважин. Сб. научн. тр. «Вопросы промывки скважин с горизонтальными участками ствола» // ОАО НПО «Бурение». Краснодар: 1998. С. 9 15.

174. Петров В.А. Исследование движения газожидкостных смесей в трубах при бурении скважин и глушении газовых фонтанов. Дисс. на соискание уч. ст. к.т.н. М.: МИНХ и ГП, 1974.

175. Петров В.А., Леонов Е.Г., Филатов Б.С., Исаев В.И. Графический метод расчета глушения открытых газовых фонтанов // ВИЭМС, Сер. XI. Техн. и технол. геол.разв. работ., орг. производства, 1974, вып.13. С. 19-33.

176. Петров В.А., Леонов Е.Г., Филатов Б.С., Исаев В.И. Исследование газовых фонтанов и графический метод расчета глушения // Газовая промышленность, №8,1974. С: 34-36.

177. Пирвердян A.M. К теории воздушного подъёмника // Нефтяное хозяйство, №4, 1951. С. 18-21.

178. Полянин Л.Н., Дробков В.П. Прикладная гидромеханика восходящих газожидкостных потоков. М.: Энергоатомиздат, 2004. 80 с.

179. Потапов А.Г. Сопротивление при турбулентном течении растворов. В сб. "Бурение глубоких разведочных скважин в осложнённых условиях Нижнего Поволжья", вып.27. Мин. нефт. пром., АН СССР, ИГИРГИ, 1976. С. 27-31.

180. Потапов А.Г., Васильченко С.В., Сутырин А.В. Методика по выбору реологических свойств буровых растворов и технологии очистки горизонтальных скважин. ВРД 39 1.8 - 045 - 2001. М.: ООО «ВНИИГАЗ», 2001. 17 с.

181. Правила 28-64 измерения расхода жидкости, газов и паров стандартными диафрагмами и соплами. М.: Изд. стандартов, 1965.

182. Правила безопасности в нефтяной и газовой промышленности (ПБ 08624-03). Серия 08. Выпуск 4. Федеральное государственное унитарное предприятие «Научно-технический центр безопасности в промышленности Гос-гортехнадзора России», 2004. 312 с.

183. Проблемы механики неоднородных сред и методы их решения для интенсификации технологических процессов. М.: АН СССР, сектор механики неоднородных сред, препринт № 19, 1987. 125 с.

184. Протодьяконов И.О., Люблинская И.Е., Рыжков А.Е. Гидродинамика и массообмен в дисперсных системах жидкость-твердое тело. Ленинград: Химия. 1987. 336с.

185. Пыхачёв Г.Б., Исаев Р.Г. Подземная гидравлика. М.: Недра, 1973.

186. Рабинович Н.Р. Инженерные задачи механики сплошной среды в бурении. М.: Недра, 1989. 270 с.

187. В.ИМяченков, В.П.Мальцев, В.П.Майборода и др. Расчеты машиностроительных конструкций методом конечных элементов: Справочник. Под общей редакцией В.И. Мяченкова. М.: Машиностроение, 1989. 520 с.

188. Рахматулин Х.А. Основы гидродинамики взаимопроникающих движений сжимаемых сред// ПММ, т.20, вып.2, 1956. С. 184-195.

189. Рахматулин Х.А., Сагомонян А.Я., Бунимович А.И., Зверев И.Н. Газовая динамика. М.: Высшая школа, 1965. 722 с.

190. Рид Р., Праусниц Дж., Шервуд Т. Свойства газов и жидкостей: Справочное пособие. JL: Химия, 1982. 592 с.

191. Сахаров В.А., Мохов М.А. Гидродинамика газожидкостных смесей в вертикальных трубах и промысловых подъёмниках. М.: ФГУП изд-во «Нефть и газ» РГУ нефти и газа им. И.М. Губкина, 2004. 398 с.

192. Самохвалов С.Ю., Леонов Е.Г. Оценка скорости осаждения шлама в буровом растворе.Сборник научных трудов. М.: МИНГ им. И.М. Губкина, вып. 214. 1987. С. 16-22

193. Седов Л.И. Методы подобия и размерности в механике. М: Наука, 1965. 195 с.

194. Седов Л.И. Механика сплошных сред. Т. 1. М.: Наука, 1983. 528 с. Т. 2. М.: Наука, 1984. 560 с.

195. Синайский Э.Г. Гидродинамика физико-химических процессов. М.: Недра, 1997. 352 с.

196. Слезкин Н.А. Динамика вязкой несжимаемой жидкости. М.: Гостех-издат, 1955. 520 с.

197. Смолдырев А.Е. Трубопроводный транспорт. М.: Недра, 1970. 272 с.

198. Солдатов A.M. Дренирование скважин с использованием аэрированной, жидкости//Нефтепромысловое дело, № 1, 1968. С. 10-11.

199. Соловьев Е.М. Заканчивание скважин. М.: Недра, 1979. 304с.

200. Coy С. Гидродинамика многофазных систем. М.: Мир, 1971.536с.

201. Справочное руководство по проектированию, разработке и эксплуатации' нефтяных месторождений. Под ред. Ш.К. Гиматудинова. М.: Недра, 1983. 455с:

202. Справочник.инженера по. бурению/Под ред. В.И. Мищенко и Н:А. Сидорова: М.: Недра, т. 1, 1973. С. 446-490:

203. Стернин Л.Е, Шрайбер А.А. Многофазные течения газа с частицами. М.: Машиностроение, 1994. 320 с.

204. Талачев B.C., Чепура И.В., Павлов В.П. Исследование течения жидкости, индуцированного струей газа. Теор. основы хим. технологии, т.6, № 2, 1972. С. 219-231.

205. Тарасов В.К., Гусак Л.Н., Валяев Н.И. Гидравлика двухфазных потоков. М.: МИСИ им. В.В.Куйбышева, 1983. 83с.

206. Тарг С.М. Основные задачи теории ламинарных течений. М.,Л.: Гос. изд-во тех.-теор. лит-ры, 1951. 420 с.

207. Телетов С.Г. Диссертация на соиск. уч. ст. д. т. н., МГУ, 1938.

208. Телетов С.Г. Вопросы гидродинамики двухфазных смесей. Вестник МГУ, № 2, 1958г. С. 15-27.

209. Телетов С.Г. Новые исследования по общим направлениям гидродинамики и энергии двухфазных течений. М.: Атомиздат, 1970. 63 с.

210. Телетов С.Г. Об обработке опытных данных по паро- и газожидкостным смесям и методике эксперимента . В сб. «Гидродинамика и теплообмен при кипении в котлах высокого давления». М.: Изд. АН СССР, 1955. С. 46-64.

211. Телетов С.Г. Уравнение гидродинамики двухфазных жидкостей // ДАН; СССР, т.50, 1945. С. 99-102.

212. Теплицкий И.С. Оценка дебита газового фонтана по размерам газо-водяного бугра на поверхности воды в кратере // ДАН УзССР. № 8,1968. С. 10 -11.

213. Теория и прикладные аспекты гидротранспортирования твердых материалов / И.А. Асауленко, Ю.К. Витошкин, В.М. Карасик, С.И. Криль, В.Ф. Очеретько. Киев: Наукова думка, 1981.364с.

214. Теплопередача в двухфазном потоке. / Под ред. Д.Баттерворса и Г.Хьюитта. М.: Энергия, 1980. 326 с.

215. Технология бурения нефтяных и газовых скважин. Уч. для вузов. Под общей редакцией А.И. Спивака. М.: ООО «Недра Бизнесцентр», 2003. 509 с.

216. Точигин А.А., Одишария Г.Э. Прикладная гидродинамика газожидкостных смесей. М.: ВНИИгаз, Иван. гос. энерг. ун-т. 1998. 400 с.

217. Транспортирование водоугольных суспензий: гидродинамика и температурный режим / А.С.Кондратьев, В.М.Овсянников, Е.П. Олофинский, Б.С.Степин, И.А.Чиненков. М.: Недра, 1988. 213 с.

218. Трошин А.К. Ардалион Николаевич Новосильцев. 1816-1878. (Научно-биографическая литература). М.: Наука, 1996. 76с.

219. Уилкинсон У.Л. Неньютоновские жидкости. М.: Мир, 1964. 216 с.

220. Уоллис Г. Одномерные двухфазные течения. М.: Мир, 1972. 440 с.

221. Устименко Б.П. Процессы турбулентного переноса во вращающихся течениях.-Алма-Ата : Наука, 1977. 22с.

222. Файзуллаев Д.Ф. Ламинарное движение многофазных сред в трубопроводах. Ташкент: Фан, 1966.

223. Файзуллаев Д.Ф., Гурбанов Р.С., Расизаде Я.М. Элементы гидравлики смесей. Ташкент: Фан, 1970. 156 с.

224. Файзуллаев Д.Ф., Умаров А.И., Шакиров А.А. Гидродинамика одно- и двухфазных сред и ее практическое приложение. Ташкент: Фан, 1980. 167 с.

225. Фисенко В.В. Критические двухфазные потоки. М.: Атомиздат, 1978.159 с.

226. Фокс Д.А. Гидравлический анализ неустановившегося течения в трубопроводах. М.: Энергоиздат, 1981. 248 с.

227. Фортье А. Механика суспензий. М.: Мир, 1971. 264с.

228. Франкль Ф.И. К теории движения взвешенных наносов // ДАН СССР.' ХСП, №2. 1953.

229. Фукс М.А. Механика,аэрозолей. М.: АН СССР, 1955. 3 51 с.

230. Хабибуллин З.Ф., Филимонов Ю.Ф. Определение давления аэрированной' жидкости на забой при бурении скважин // Нефть и газ. № 7,1967. С. 10 13.

231. Хаппель Дж., Бреннер Г. Гидродинамика при малых числах Рейнольдса. М.: Мир, 1976. 630 с.

232. Химическая гидродинамика: Спр. Пособие / А.М.Кутепов, А.Д.Полянин, З.Д.Запрянов, А.В. Вязьмин, Д.А.Казенин. М.: Бюро Кван-тум, 1996. 336 с.

233. Хинце И.О. Турбулентность. М.: Гос. изд-во физ-мат. лит-ры, 1963. 680с.

234. Холпанов Л.П., Запорожец В.П. Зиберт Г.К., Кащицкий Ю.А. Математическое моделирование нелинейных термогидрогазодинамических процессов в многокомпонентных струйных течениях. М.: Наука. 1998. 320 с.

235. Хьюитт Дж., Холл-Тейлор Н. Кольцевые двухфазные течения. М.: Энергия, 1974. 407 с.

236. Царевич К., Шищенко Р., Бакланов Б. Глинистые растворы в бурении. Баку, Москва: ОНТИ НКТП СССР, Азнефтеиздат, 1935. 329 с.

237. Циклаури Г.В., Данилин B.C. Адиабатные двухфазные течения. М.: Атомиздат, 1973.

238. Чарный И.А. Неустановившееся движение реальной жидкости в трубах. М.: Недра, 1975.

239. Чарный И.А., Евдокимова В.А., Кочина И.Н. Определение свободного дебита газовых скважин // Газовая промышленность. №4, 1963. С. 3-6.

240. Чарный И.А. О продвижении границы изменения агрегатного состояния при охлаждении или нагревании тел // Изв. АН СССР. ОТН. № 2,1948. С. 187-201.

241. Чугаев P.P. Гидравлика. Учебник для вузов. Ленинград: Энергоиздат, 1982. 672 с.ОАО «Газпром», ООО «ВНИИГАЗ», 2000. 190 с.

242. Чисхолм Д. Двухфазные течения в трубопроводах и теплообменниках. М.: Недра, 1986. 204 с.

243. Шеберстов Е.В., Леонов Е.Г. Расчет давления в скважине при бурении с применением аэрированных жидкостей // Нефтяное хозяйство, № 12, 1968. С. 14-17.

244. Шевцов В.Д. Предупреждение газопроявлений и выбросов при бурении глубоких скважин. М.: Недра, 1988. 200 с.

245. Шищенко Р.И., Есьман Б.И., Кондратенко П.И. Гидравлика промывочных жидкостей. М.: Недра, 1976. 294 с.

246. Шоу С.Ф. Теория и практика газлифта. Пер. с англ. M.-JL, Гостоптехиз-дат: 1948. 198 с.

247. Шлихтинг Г. Теория пограничного слоя. М.: Наука, 1969.742с.

248. Шрайбер А.А. Многофазные полидисперсные течения с переменным фракционным составом дискретных включений // Итоги науки и техники. Сер. Комплексные и специальные разделы механики. М.: ВИНИТИ, 1988,380 с.

249. Шульман З.П., Задворных В.Н., Литвинов А.И. Геодинамика нелиней-но-вязкопластичной жидкости в кольцевых каналах с подвижными стенками. Препринт № 45 АН БССР, инст. тепло- и массообмена им. А.В. Лыкова, Минск: 1987. 51 с.

250. Шумилов Л.П. О транспорте шлама по стволу скважины // Нефтяное хозяйство. № 7, 1966. С. 27 30.

251. Щукин А.К. Теплообмен и гидродинамика внутренних потоков в полях массовых сил. М.: Машиностроение, 1980. 240 с.

252. Якимов Ю.Л. Силы, действующие на малое тело в произвольном потоке несжимаемой жидкости, и уравнения движения двухфазной среды // Изв. АН СССР, МЖГ. № 3, 1973. С. 84-92.

253. Ямпольский В.И. Паронагнетание и движение смесей в трубах. Ижевск: Алфавит, 1997. 316 с.

254. Alves I.N. et al. Modeling Annular Flow Behavior for Gas Wells, paper presented at the 1988 Annual Winter Meeting of ASME, Chicago, 27 November 2 December.

255. Ansari A.M. et a\. A Comprehensive Mechanistic Model for Two-phase Flow in Wellbores. SPEPF, May, 1994, 143; Trans., AIME, 297.

256. Asheim H. MONA, An Accurate Two-Phase Well Flow Model Based on Phase Slippage. SPE Production Engineering, May, 1986. Pp. 221 230.

257. Attu A. Une Solution Generate de L'Ecoulement Diphasique Compressible atravers un Elargissement Brusque // European Journal Mech. Eng. M, Vol. 42, N 1, 1997. Pp. 3-8.

258. Aziz K., Govier G.W., Fogarasi M. Pressure drops in Well Production Oil and Gas. J. Cdn.Pei, Tech. (Juli-Septumber, 1972) 11,38.

259. Baxendell P.B., Thomas R. The Calculation of Pressure Gradients in High-Rate Flowing Wells. JPT (October, 1961) 1023; Trans., AIME, 222.

260. Barnea D., Shoham O., Taitel Y.M. Flow Pattern Transition for Vertical Downward Two-Phase Flow. Chem. Eng. Sci. №37, 1982. 741 p.

261. Barnea D. A Unified Model for Predicting Flow-Pattern Transition for the Whole Range of Pipe Inclinations/ Intl. J. Multihase Flow. 1987, 13,1.

262. Baroczy C.J. A systematic correlation for two-phase pressure drop. Chem. Engng Design, 62(44), 1966. Pp. 232-249.

263. Bert R. Bloodworth, George J. Keely Jr., Peter E. Clark. Mud relaxation» measurements help predict hole cleaning ability. Oil & Gas Journal. Vol.90.№22, . 1992. Pp.73-78.

264. Beggs H.D., Brill J.P. A Study of Two-Phase Flow in inclined Pipes. JPT (May, 1973) 607; Trans. AIME. 255.

265. Bobo A. Roy. New aerated mud equipment can control lost circulation: O.G.J., № 18, 1968. Pp. 76 - 77.283: Bobo A. Roy. New air/mud system can boost performance. O.G.J., № 19, 1968. Pp. 114- 118.

266. Bourgoyne A.T et all Applied Drilling Engineering, Richardson, Texas, Society of Petroleum Engineers, 1986. 214 p.

267. Brauner N., Barnea D*- Slag/Churn Transition in Upward Gas-Liquid Flow. Chem. Eng. Sci. №41, 1986. 159 p.

268. Burkhardt J.A. Wellbore pressure surges produced by pipe mouvement. In: Jour. Petr. Tech., № 6, 1965. pp.595-605.

269. Caetano E.F. Upward Vertical Two-Phase Flow Through an Annulus. PhD dissertation, U. of Tulsa, Tulsa, Oklahoma. 1985.

270. Caetano E.F., Shoham O., Brill J.P. Upward Vertical Two-Phase Flow Through an Annulus. Part I. Single-Phase Friction Factor Taylor Bubble-Rise Velocity and Flow-Pattern Prediction. J. Energy Res. Tech. (March, 1992), 114; 1.

271. Caetano E.F., Shoham O., Brill J.P. Upward Vertical Two-Phase Flow Through an Annulus. Part II. Modeling Bubble, Slug and Annular Flow. J. Energy Res. Tech. (March, 1992), 114; 14.

272. Chang H.C., Smith T.M. Experiment of particle deposition in a horizontal sampling line. Amer.Ind.and Hyg.Asoc.J., 1972, 33,№ 11. Pp.722-728.

273. Chierchi G.L., Ciuccii G.M., Sclocchi G. Two-Phase Flow Vertical Row inOil Wells-Prediction of Pressure Drop. JPT (August, 1974) 927; Trans. AIME. 257.

274. Chisholm D., Sutherland L.A. Prediction of pressure changes in pipeline systems during two-phase flow/ Paper № 4,1. Mech. E./I Chem. E. Joint Symp. on Fluid Mechanics and measurements in two-phase systems, 24-25 Septembre 1969, University of Leeds.

275. Dean D.E., Stiel L.I. // AICHE Journal. Vol 11. 1965 526 p.

276. Diener R., Friedel L. Reproductive accuracy of selected void fraction correlations for horizontal and vertical upflow// Forsch Ingenieurwes, 64, N 4, 5, 1998. Pp. 87-97

277. Dring R.P., Caspar J.R., Suo M. Particle Trajectories in Turbine Cascades. Journal of Energy. Vol.3, №3, May-June, 1979. Pp.161-166.

278. Duns H. Jr., Ros N.C.J. Vertical Flow of Gas and Liquid Mixtures in Wells. Proc. Sixth World Pet. Cong. Tokyo, 1963. Pp.451.

279. Fredrickson A.G., Bird R.B. Non-Newtonian Flow in Annuli. Industrial and Engineering Chemistry. Vol.50, №3, 1958. Pp.347-352.

280. Friedl M.J. Bubble plumes and their interaction with the water surface. Dissertation No 12667. Swiss Federal Institute of Technology Zurich (ETH). 1998.

281. Friedl M.J., Fannelop Т.К. Bubble plumes and their interaction with the water surface.// Applied Ocean Research. 22, No2, 2000. Pp. 119-128.

282. Fannelop Т.К. and Sjoen K. Hydrodinamics of Underwater Blowouts. // Norwegian Maritime Research, v.8, №.4, 1980. Pp. 17-34.

283. Fancher G.H.Jr., Brown K.E. Prediction of Pressure Gradients for Multiphase Flow in Tubing. SPEJ (March, 1963) 59; Trans., AIME, 228.

284. Fernandes R.C., Semait N., Dukler A.E. Hydrodynamic Model for Gas-Liquid Slug Flow in Vertical Nubes. AlChE J. №32, 1986. 981 p.

285. Gabold G., Nguyen J-P. Formulaire du foreur. Sixieme edition. Edition Tech-nip, Paris. 1989. 542 p.

286. Grodde R.H. Rheologie Kolloider Suspension insbesondere der Bohrspuelun-gen Erdoel und Kohle. Bd. 13, 1960, №1, pp. 11-18, №2, Pp. 79-83.

287. Guet S., Ooms G/ Fluid mechanical Aspect of the Gas-lift Technique// Annual? Review of Fluid Mechanics. Vol. 38,2006. Palo Alto (Calif.), 2006. Pp. 225-249.

288. Hagedorn A.R., Brown K.E. Experimental Study of Pressure Gradients Occurring During Continuous Two-phase Flow in small-Diameter Vertical Conduits. JPT (April, 1965) 475; Trans. AIME. 228.

289. Hanks R.W. The Laminar-Turbulent Transition for Fluids with a Jield Stress In: A.I.Ch.E.J., vol.9, № 3. May,1963. Pp.306-309.

290. Hanks W.R., Christiansen E.B. The Laminar-Turbulent Transition in Noniso-thermal Flow of Pseudoplastic Fluids in Tubes. A.I.Ch.E.Journal, Vol.8, №4, pp.467-471.

291. Hasan A.R., Kabir C.S. A Study of Multiphase Flow Behavior in Vertical Wells. SPEPE (May, 1988), 263; Trans., AIME, 285.

292. Hasan A.R., Kabir C.S. Predicting Multiphase Flow Behavior in a Deviated Well. SPEPE (November, 1988), 474.

293. Hasan A.R., Kabir C.S. Two-Phase Flow in Vertical and Inclined Annuli. Int. J. Multiphase Flow. №18, 1992. Pp. 279.

294. Hoefele E.O., Brimacombe J.K. Flow Regimes in Submerged Gas Injection// Metallurgical Transactions B, Vol. 10B, 1979. Pp. 631 -648.

295. Hopkin E.A. Factor affecting cutting removal during rotary drilling. J. Petr. Tech., yunie 1967. Pp. 807 814.

296. Joshi S.D. Horizontal Well Technology. PenWellBooks, Tulsa, Oklahoma, 1991.535 p.

297. Kabir C.S., Hasan A.R. Performance of Two-Phase Gas/Liquid Model in Vertical Well. J.Pet. Sci. and Eng. №4, 1990. P. 273.

298. Kesler H.G., Lee B.I. Improve prediction of enthalpy of fractions // Hydrocarbon Processing, Vol.55, №3, 1976. Pp. 153-158.

299. Kobus H. Bemessungsgrundlagen und Anwendungen fir Luftschleier in Wasserbau.- Bielefeld, E.Schmidt Verlag, 1973.

300. Kobus H. On the use of air bubble screens as oil barriers. Fundam. Tools used environ, probl. 16th Congr, S,o Paulo, 1975.

301. Lockart R.W., Martinelli R.C. Proposed correlation of data for isothermal two-phase two-component flow in pipes. Chem. Engng/ Progr., 45(1), 1949. Pp. 39-48.

302. Lorenz H. Die Arbeitsweise und Berechnung des Druckvasserhebers, Z. D: V. D. J., 1909.

303. McQuillan R.W., Whalley P.B. Flow Pattern in Vertical Two-Phase Flow. Intl. J. Multiphase Flow. №11, 1985. P. 161.

304. Martinelli R.C., Boelter L.M.K., Taylor T.H.M., Thomson E.G., Moen R.H. Isothermal pressure drop for two-phase two-component flow in a horizontal pipe. Trans. Amer. Soc. Mech. Engrng, 66 (2), 1944. Pp. 139 51.

305. Martinelli R.C., Lockart R.W., Putnum J.A. Two-phase two-component flow in the viscous region. Nrans. Amer. Soc. Mech. Engrs, 42(4), 1946. Pp. 681-705.

306. Maximov V., Limar E., Isaev V. Hydrodynamic Study of Underwater Gas Blowouts: Theory and Experiment.// Proceedings Int. Gas Res. Conf., paper EPP-07 (CD-ROM), 5-8 Nov., Amsterdam, Netherlands, 2001.

307. Moore T.V. and Wilde H.D. Jr. Exerimental measurement of slippage in flow through vertical pipes. Trans. Amer. Inst. Mining and Met. Engrs (Pet.Div.), 92, 1931. Pp. 296-319.

308. Mukheijee H., Brill J.P. Pressure Drop Correlations for Inclined Two-Phase flow. J. Energy Res. Tech. (December, 1985) 107, 549.

309. Orkiszewski J. Predicting Two-Phase Pressure Drops in Vertical Pipes.//JPT (June, 1987) 829; Trans., AIME, 240.

310. Okrajni Slavomir S., Azar J.J. The Effects of Mud Rheology on Annular Hole Cleaning in Directional Wells. In: SPE Reprint Series. №30. Directional Drillihg. Society of Petroleum Engineers. 1990. Pp. 46-57.

311. Passut C. A., Danner R. P. Correlation of ideal gas enthalpy, heat capacity and entropy// hid. Eng. Chem. Process Das. and Dev, Vol. 11, № 4,1972. Pp. 543-546.

312. Peng D.Y., Robinson D.B. A new two-constant equation of state // Ind. Eng. Chem. Fundamen. Vol 15, 1976. Pp: 59-64.

313. Poettmann F.H., Carpenter P.G The Multiphase Flow of Gas, Oil and Water through Vertical Flow Strings with Application to the Design of Gas-Lift Installations. A.P.I. Drilling and Production Practice, 257, 1952.

314. Poettmann F.H., Bergman W.E. Density of drilling muds reduced b^air injection. World Oil, 1 august 1955. Pp. 41 43, 45.

315. Rhem B. Don't overlook aerateg mud. O.G.J., 2 dec. 1963. Pp. 150 153.

316. Sahajwalla V., Castillejos A.H., Brimacombe J.K. The Spout of Air Jets Upwardly Injected into a Water Bath.// Metallurgical Transactions B, Vol. 2IB, 1990. Pp. 71-80.

317. Sato Y., Sekoguchi K. Liquid velosity distribution in two-phase bubble flow.// Int.J.Multiphase Flow., vol.2, 1976, № 1.

318. Silvester N.D. A mechanistic Model for Two-Phase Vertical Slug Flow in Pipes. J. Energy Res. Tech. December 1987. Pp. 109,206.

319. Stefan MJ. Versuche uber die Verdampfung Sitzungsberichte Akademie der Wissenschaften, Bd. 68, 1873.

320. Taitel Y.M., Barnea D., Dukler A.E. Modeling Flow Pattern Transitions for Steady Upward Gas-Liquid Flow in Vertical Tubes. AlChe J. № 26, 1980. P. 345.

321. Tchen C.M. Mean value and correlation problems connected with the motion of small particles suspended in a turbulent fluid, автореф. диссерт. The Hague, Dlft, 1947.

322. Turner J.M., Ph.D. thesis, Dartmouth College, Hanover, N.H., 1966

323. Turner J.M., Wallis G.B. Rept. NYO-3114-6, 1965.v ь1345. Users Manual for API 14B.SSCSV Sizing Computer Program. Second edition. API (1978) Appendix B, pp. 38-41.

324. Vogel P., Marx C. Berechnung von Blowoutraten in Erdgassonden. (Determination of of Blowout Rates in Gas Wells). Erdoel-Erdgas, .101. Jg., Heft 10, Oktober, 1985. Pp. 311-317.

325. Warriner R.A. and Cassity T.G. Relief-Well Requirements- To Kill a High-Rate Gas Blowout From a Deepwater Reservoir. Journal, of Petroleum Technol-. ogy, December, 1988. Pp. 1602-1608.

326. Zeidler H. Udo. An experimental analysis of the transport of drilled particles. Soc. Petr. Eng., febr. 1972. Pp. 39 48.