Численное исследование восходящего нефтегазового потока в вертикальной скважине с установкой электроцентробежных насосов тема автореферата и диссертации по механике, 01.02.05 ВАК РФ

Бородин, Станислав Леонидович АВТОР
кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ
Тюмень МЕСТО ЗАЩИТЫ
2012 ГОД ЗАЩИТЫ
   
01.02.05 КОД ВАК РФ
Диссертация по механике на тему «Численное исследование восходящего нефтегазового потока в вертикальной скважине с установкой электроцентробежных насосов»
 
Автореферат диссертации на тему "Численное исследование восходящего нефтегазового потока в вертикальной скважине с установкой электроцентробежных насосов"

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

БОРОДИН Станислав Леонидович

ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ ВОСХОДЯЩЕГО НЕФТЕГАЗОВОГО ПОТОКА В ВЕРТИКАЛЬНОЙ СКВАЖИНЕ С УСТАНОВКОЙ ЭЛЕКТРОЦЕНТРОБЕЖНЫХ НАСОСОВ

01.02.05 - механика жидкости, газа и плазмы

АВТОРЕФЕРАТ

диссертации на соискание учёной степени кандидата физико-математических наук

3 МАЗ Ш

Тюмень-2012

005018800

Работа выполнена в Тюменском филиале Федерального государственного бюджетного учреждения науки Института теоретической и прикладной механики им. С.А. Христиановича Сибирского отделения Российской академии наук.

Научный руководитель: кандидат физико-математических наук,

доцент Мусакаев Наиль Габсалямович

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

профессор

Кислицын Анатолий Александрович

(Тюменский государственный университет, заведующий кафедрой микро- и нанотехнологий)

кандидат технических наук Черемисин Николай Алексеевич

(ООО «Тюменский нефтяной научный центр», старший эксперт)

Ведущая организация: Институт механики Уфимского научного

центра РАН

Защита состоится 17 мая 2012 года в 11— на заседании

диссертационного совета ДМ 212.274.09 в Тюменском

государственном университете по адресу: 625003, г. Тюмень, ул. Перекопская, д. 15а, ауд. 410.

С диссертацией можно ознакомиться в библиотеке Тюменского государственного университета по адресу: 625003, г. Тюмень, ул. Семакова, д. 10.

Автореферат разослан «^6» апреля 2012 года.

И.о. учёного секретаря диссертационного совета доктор физ.-мат. наук

Татосов А. В.

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

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

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

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

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

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

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

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

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

Практическая и теоретическая ценность. В диссертационную работу вошли результаты исследований, выполненных в соответствии с планами фундаментальных исследований РАН, а также в рамках государственного контракта №02.516.11.0004 по ФЦП «Исследования и разработки по приоритетным направлениям развития научно-технологического комплекса России на 2007-2012 годы», государственного контракта №119-ДОН по Тюменской областной целевой научно-технической программе и подержанных грантом РФФИ № 08-01-90300-Вьет_а, грантами Президента РФ для государственной поддержки ведущих научных школ РФ (руководитель научной школы - академик РАН Р.И. Нигматулин), фантом Губернатора Тюменской области на реализацию проекта по фундаментальным и прикладным научным исследованиям.

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

На защиту выносятся:

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

2. Результаты численных экспериментов по изучению влияния определяющих параметров на структуру потока в подъёмной колонне

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

Апробация работы. Результаты, приведённые в диссертации, докладывались и обсуждались на следующих конференциях:

• Всероссийская научно-техническая конференция «Нефть и газ Западной Сибири», Тюмень, 2007;

• VI и VII Международная научная школа-конференция «Актуальные вопросы теплофизики и физической гидрогазодинамики», Украина, Алушта, 2008, 2009;

• International Workshop on «Thermal Hydrodynamics of Multiphase Flows and Applications», Hanoi-Vietnam, 2009;

• VII Международная конференция «Химия нефти и газа», Томск, 2009;

• Российская конференция «Многофазные системы: природа, человек, общество, технологии», посвящённая 70-летию академика Р.И. Нигматулина, Уфа, 2010;

• XV Международная конференция по методам аэрофизических исследований, Новосибирск, 2010;

• X Всероссийский съезд по фундаментальным проблемам теоретической и прикладной механики, Нижний Новгород, 2011.

Результаты работы докладывались автором на семинарах ТюмФ ИТПМ СО РАН под руководством профессора А.А. Губайдуллина.

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

Объём и структура работы. Диссертация состоит из введения, трёх глав, заключения, списка литературы и двух приложений. Работа изложена на 94 страницах, включает 27 рисунков и 1 таблицу. Список литературы содержит 99 наименований.

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

Автор выражает глубокую признательность научному руководителю Наилю Габсалямовичу Мусакаеву за понимание, терпение и поддержку в науке и жизни.

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

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

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

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

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

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

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

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

ф

"Г = аг

<1г

"

т,с, +от с +-

(и)

А,—т,{1щ 1(М) + 1{Ш)), А —т1ш),

£ -А:

1 +

' -Р(м)

, (Р-Рт) ,

В' = ~

р,т =р.ехр(-Т,/Т).

Здесь р. Тир- давление, температура и плотность двухфазной смеси; ^ - сила трения между двухфазным потоком и стенками канала, отнесённая к единице длины канала; g - ускорение свободного падения; mi (/ = /,§)- массовый расход /-ой фазы; р,° (/ = /,§)-истинная плотность /-ой фазы; /(У) (/ = Л/, ¿) - удельная теплота перехода у'-ой компоненты из жидкой фазы в газовую; ккл (/ = I, g; _/ = М, V) - массовая концентрация в /-ой фазе >ой компоненты; С{1) -постоянная Генри; и = М,Ь)~ газовая «постоянная» >ой

компоненты в газовой фазе; ()„- интенсивность отвода тепла, отнесенная к единице длины канала; с, (/ = /,&)- удельная теплоёмкость /-ой фазы; р, и Т. - эмпирические аппроксимационные параметры. То, к какой фазе относится соответствующий параметр, указывается нижним индексом / или I- жидкая фаза, g- газовая фаза.

Для замыкания вышеприведенной системы уравнений необходимо задание интенсивности теплоотвода ()„, скачка давления и температуры газожидкостного потока в результате работы УЭЦН и, при необходимости, роста твёрдых отложений на внутренних стенках подъёмной колонны скважины.

Примем, что радиальные температурные поля в системе «скважина - окружающая порода» являются квазистационарными. Тогда для интенсивности отвода тепла, отнесённой к единице длины канала, можно записать следующее выражение:

к , , ^ с}„ = 2л(т-ты) —— + —,

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

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

При сделанных допущениях уравнения, описывающие внешнюю тепловую задачу в мёрзлых породах, начальные и граничные условия могут быть записаны в виде:

Здесь Та, т и Тех! (2) - температуры окружающей скважину мёрзлой породы и протаявшей породы, соответственно; -геотермальная температура; г - время; Хвкп и Хан <2> - коэффициенты температуропроводности окружающей скважину мёрзлой породы и

^,0=^(2)=^. (, = 0)

протаявшей породы, соответственно; г - радиальная координата; ка,т и Лет»о- коэффициенты теплопроводности окружающей скважину мёрзлой породы и протаявшей породы, соответственно; в- радиус протаявшей зоны; Тех1(и2)- температура фазового перехода в мёрзлых породах (температура таяния льда); Рехцц - плотность мёрзлой породы; ¿(12) - удельная теплота протаивания мёрзлой породы.

Для расчёта скачков давления (Ар) и температуры (Д7) газожидкостного потока в результате работы УЭЦН получены следующие выражения:

Лр = А0+А,-т + А2-т2,

где А0, А,, Л2- коэффициенты полинома, рассчитываемые из характеристики центробежного насоса, т - массовый расход смеси, N— мощность погружного двигателя УЭЦН.

При необходимости в модели можно учесть рост парафиновых отложений на внутренних стенках скважины. В работе приведён вывод соотношения для интенсивности такого роста (Шагапов В.Ш., Мусакаев Н.Г., 1998):

где 6Х- толщина твёрдых отложений; Те- температура начала кристаллизации тяжёлого компонента; Д- коэффициент теплопередачи между твёрдыми отложениями и стенкой скважины; Т„- температура внутренней стенки подъёмной колонны скважины; р,°- истинная плотность твёрдой фазы; Ь, - удельная теплота парафинообразования; £)- коэффициент диффузии; 5Л- число Шервуда; Я - радиус поперечного сечения канала; кКИ)е - равновесная концентрация тяжёлой компоненты; Та- температура на поверхности твёрдой фазы.

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

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

_ Ы-ЛР пг/р

/и, с, +тяс1

Для стабильной работы УЭЦН необходимо, чтобы на глубине подвески электронасосного агрегата структура потока была пузырьковой, в этой связи изучено влияние различных параметров на структуру потока в подъёмной колонне скважины (рис. 1).

2000 ¡500 1000 500 0

___ а)

2000 1500 1000 500 0

i б)

——>■

0 100 200 300 т, т/сут

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

Области параметров, расположенные ниже рассчитанных кривых, соответствуют пузырьковой структуре потока.

Из рис. 1а видно, что при увеличении газового фактора интервал, в котором существует пузырьковая структура потока, сужается, это приводит к необходимости устанавливать электронасосный агрегат на большей глубине. После достижения газовым фактором некоторого значения пузырьковая структура потока престаёт существовать по всей высоте скважины и использование УЭЦН в подобных обстоятельствах без дополнительного оборудования (газосепаратор) приведёт к неисправностям электронасосного агрегата. Как видно из рис. 16 и1в, увеличение дебита и снижение давления на забое скважины приводят к сужению интервала, в котором существует пузырьковая структура потока. Поэтому при установке насоса необходимо учитывать зависимость глубины его погружения от массового расхода насоса.

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

2000 1500 1000 500

0.

к Нсфп, Нола-ч........ V V Пок|- а) \! .....\ жишП.л

\ ! \ га-} N. \ •........

•^Очу 1

| ! гч

273 283 293 303 Те 323 Т„]

В)

п ^ N ; Поко-б)

1Ьч Г1 состоянии | \ А га}

к|м[вскЦии ! 1 : : : 1 .......

| | Ч|

• ■ ! 1 ! ;

273 283 293 303

2000 1500 1000 500

°273 283 293 303 Тс 323 Г., К

Рис. 2. Профили температуры внутренней стенки подъёмной колонны скважины Тк для различных флюидов в межтрубном пространстве (а), различного поведения газа (б) и

различной толщины теплоизолирующего слоя <й (в)

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

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

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

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

Рис. 3. Профили температуры внутренней стенки подъёмной колонны скважины Т„ для различной глубины погружения электронасосного агрегата (а) и массового расхода насоса (б). Линии 1, 2, 3 и 4 (а) соответствуют глубине подвески Нр = 500, 1000, 1500 и 2000 м; Мр = 50 т/сут. Линии 1, 2 и 3 (б) соответствуют расходу Мр = 50, 75 и 125 т/сут; Нр = 2000 м.

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

В связи с вышеизложенным, проведена серия расчётов движения фронта таяния мёрзлых пород окружающих скважину с УЭЦН.

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

Рис. 4. Изменение во времени радиуса протаивания окружающих скважину мёрзлых пород. Зависимости построены для глубины 100 метров. Линии 1, 2 и 3 (а) соответствуют случаям, когда в межтрубном пространстве скважины находится нефть, газ и вода. Линии (б) соответствуют случаям: 1 -отсутствия теплоизоляции; 2 - АЛ = 5 мм, в качестве материала теплоизоляции в расчетах взят пенополиуретан с коэффициентом теплопроводности 0,0067 Вт/(м К); 3 -ДА = 10 мм, пенополиуретан; 4 - ДА = 2 мм, Thermo Coat с коэффициентом теплопроводности 0,001 Вт/(мК).

Из рис. 46 видно, что при увеличении теплоизолирующего слоя или при использовании материала с меньшим коэффициентом теплопроводности радиус протаявшей зоны уменьшается. Как видно из рис. 4а, если в межтрубном пространстве находится нефть (коэффициент теплопроводности в расчётах равен 0,13 Вт/(м-К)), то радиус протаивания максимален, если вода (0,68 Вт/(м К)) -минимален, газ же с наименьшим коэффициентом теплопроводности занимает промежуточную позицию. Т.е. увеличение теплопроводности вещества в межтрубном пространстве скважины не ведет однозначным образом к повышению темпов роста радиуса протаявшей зоны. Для изучения этого интересного факта в работе была проведена серия расчётов с различными коэффициентами теплопроводности вещества А™ в межтрубном пространстве (рис. 5а).

0,7 0,5 0,3 ОД

А • ! а)

1 1 • V •

• • • 1 1 1 1 1 1 1 1 1 1

• 2 1 1 4} 1 1 |1 3

т-----;— -г*-!^!---1-----

0,001 0,01 0,1 А^т, Вт/(м-К)

300 290 280 270

£

б)

¡3

0,001 0,01 од л_вт/о^)

Рис. 5. Значения радиуса протаивания мёрзлых пород (а) и температуры внутренней стенки подъёмной колонны скважины (б) при различном коэффициенте теплопроводности вещества в межтрубном пространстве. Значения рассчитаны для глубины 100 метров. Вертикальные линии 1, 2, 3 и 4 - коэффициенты теплопроводности нефти, газа, воды и газа в состоянии термогравитационной конвекции, соответственно.

Из рис. 5а видно, что увеличение Лпп до определенного значения приводит к большим значениям зоны протаивания, затем же происходит обратное - радиус в имеет с ростом Л™, все меньшие значения. Причиной максимума в(Атп), вероятно, является «конкуренция» двух факторов в величине теплоотдачи от скважины в окружающие породы, а именно, коэффициента теплопередачи р через систему труб скважины и разницы температур ДТ в стволе скважины и в окружающих породах. С увеличением А™ коэффициент теплопередачи р растет, однако величина АТ снижается вследствие больших потерь тепловой энергии скважинного потока до начала зоны многолетних мёрзлых пород. Таким образом, для уменьшения радиуса протаивания можно использовать в межтрубном пространстве вещество и с большим коэффициентом теплопроводности, однако для предотвращения выпадения парафина в стволе скважины необходимо поддержание температуры стенки подъёмной колонны скважины выше температуры начала кристаллизации парафина. Это достигается, если в межгрубном пространстве находится флюид с малым коэффициентом теплопроводности (рис. 56).

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

Рис. 6 Изменение со временем радиуса протаивания окружающих скважину мёрзлых пород. Зависимости построены для глубины 100 метров. Линии 1, 2 и 3 (а) соответствуют Мр = 50, 75 и 100 т/сут; Нр- 2000 м. Линии 1 и 2 (б) соответствуют Нр = 2000 и 1000 м; Мр - 50 т/сут.

Из рис. 6а видно, что увеличение массового расхода УЭЦН ведет к росту зоны протаивания. Данное обстоятельство можно объяснить тем, что увеличение массового расхода смеси приводит к повышению скорости движения фаз, и поток, двигаясь к многолетнемёрзлым породам, успевает отдать меньшее количество тепла в окружающие горные породы, т.е. подходит более тёплый. В результате большему значению Мр соответствует большая разница температур АТ в зоне многолетнемёрзлых пород и, следовательно, больший нагрев окружающих пород. Результаты расчетов, представленные на рис. 66 показывают, что подвеска насоса на большей глубине приводит к незначительному повышению темпов роста радиуса протаявшей зоны.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ

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

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

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

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

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

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

ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ

1. Мусакаев Н.Г., Бородин С.Л. Теоретическое исследование особенностей двухфазного течения в оснащённой электроцентробежным насосом скважине // Вестник Нижегородского университета им. Н.И. Лобачевского. 2011. № 4 (2). С. 502-504.

2. Губайдуллин A.A., Мусакаев Н.Г., Бородин С.Л. Математическая модель восходящего газожидкостного потока в вертикальной скважине // Вестник Тюменского государственного университета. 2010. № 6. С.68-75.

3. Губайдуллин А.А., Мусакаев Н.Г., Бородин C.J1. Компьютерное моделирование процессов в оснащенной УЭЦН нефтегазовой скважине // Известия вузов. Нефть и газ. 2010. № 5. С.59-65.

4. Мусакаев Н.Г., РоманюкС.Н., Бородин C.J1. Теоретическое исследование закономерностей увеличения радиуса протаивания многолетних мёрзлых пород вокруг скважины с установкой электроцентробежных насосов // Вестник Кыргызского национального университета имени Жусупа Баласагына. 2011. Серия 3. Естественно-технические науки. С. 102-106.

5. Musakaev N.G., Borodin S.L. Computational investigation of preventive methods against formation of paraffin deposits in the vertical wells equipped with an installation of electric centrifugal pumps // Proceedings of 15th International Conference on the Methods of Aerophysical Research: Abstracts, Novosibirsk, Russia, November 01-06, 2010. Parti. Pp. 190-191.

6. Shagapov V.Sh., Borodin S.L, Gubaidullin A.A., Duong Ngoc Hai, Musakaev N.G. The Theoretical Research of an Upward Two-Phase Flow with Phase's Changes in a Vertical Well // Vietnam Journal of Mechanics. 2010. Vol. 32, No. 4. Pp. 211-221.

7. Н.Г. Мусакаев, С.Л. Бородин. Численное исследование процессов, происходящих в вертикальной скважине с установкой центробежных электронасосов // Тезисы докладов Российской конференции «Многофазные системы: природа, человек, общество, технологии», посвящённой 70-летию академика Р.И. Нигматулина. Уфа: Изд-во Нефтегазовое дело, 2010. С. 156.

8. Бородин С.Л., Мусакаев Н.Г. Программа «PVTWell», свидетельство о государственной регистрации программы для ЭВМ № 2010616663. 2010.

9. Мусакаев Н.Г., Бородин С.Л. Численное исследование восходящего газонефтяного потока в добывающей скважине, оснащенной УЭЦН // Материалы VII международной конференции «Химия нефти и газа», Томск, 21-26 сентября 2009. Томск: Изд-во Института оптики атмосферы СО РАН, 2009. С. 335-338.

10. Губайдуллин А.А., Мусакаев Н.Г., Бородин С.Л., Зыонг НгокХай. Численное исследование восходящего двухфазного течения с фазовыми переходами в вертикальной скважине // Сборник научных статей «Современная наука: исследования, идеи, результаты, технологии». 2009. № 1. С. 69-71, специальный выпуск по материалам VII научной конференции «Актуальные вопросы теплофизики и физической гидрогазодинамики», Алушта, Украина, 21-25 сентября 2009.

11. Shagapov V.Sh., Borodin S.L., Gubaidullin А.А., Duong Ngoc Hai, Musakaev N.G. Mathematical Modeling of Upward Flow of a Liquid-Gas

Mixture in a Vertical Well // Proceedings of International Workshop on «Thermal Hydrodynamics of Multiphase Flows and Applications», May 5-6 2009, Hanoi-Vietnam, Pp. 161-172.

12. Мусакаев Н.Г., Губайдуллин A.A., Бородин С.Л., Зыонг НгокХай. Математическое моделирование восходящего двухфазного потока при наличии фазовых превращений и внешнего теплообмена // Материалы VI международной научной школы-конференции «Актуальные вопросы теплофизики и физической гидрогазодинамики», Алушта, Украина, 22-28 сентября 2008. Вып. 6, 4.1. С. 52-55.

13. Бородин С.Л., Мусакаев Н.Г. Математическое моделирование восходящего газожидкостного потока в вертикальной скважине // Материалы Всероссийской научно-технической конференции «Нефть и газ Западной Сибири», Тюмень, 15-17 октября 2007. Т.1. С.105-106.

Подписано в печать 10.04.2012. Формат 60x90 1/16. Усл. печ. л. 1 Тираж 100 экз. Заказ № 1649.

ЗАО «Полиграфмаркет» 625026, г. Тюмень, ул. Мельникайте, 106 оф. 106. Тел./факс: (3452) 40-38-50, 22-96-66. E-mail: poligraf@tyumen.ru. http://www.poligrafmarket.ru

 
Текст научной работы диссертации и автореферата по механике, кандидата физико-математических наук, Бородин, Станислав Леонидович, Тюмень

61 12-1/972

ТЮМЕНСКИЙ ФИЛИАЛ ФЕДЕРАЛЬНОГО ГОСУДАРСТВЕННОГО БЮДЖЕТНОГО УЧРЕЖДЕНИЯ НАУКИ ИНСТИТУТА ТЕОРЕТИЧЕСКОЙ И ПРИКЛАДНОЙ МЕХАНИКИ им. С.А. ХРИСТИАНОВИЧА СИБИРСКОГО ОТДЕЛЕНИЯ РОССИЙСКОЙ АКАДЕМИИ НАУК

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

Бородин Станислав Леонидович

ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ ВОСХОДЯЩЕГО НЕФТЕГАЗОВОГО ПОТОКА В ВЕРТИКАЛЬНОЙ СКВАЖИНЕ С УСТАНОВКОЙ ЭЛЕКТРОЦЕНТРОБЕЖНЫХ НАСОСОВ

01.02.05 - Механика жидкости, газа и плазмы

Диссертация на соискание учёной степени кандидата физико-математических наук

Научный руководитель: канд. физ.-мат. наук, доцент Н.Г. Мусакаев

Тюмень - 2012

ОГЛАВЛЕНИЕ

ОГЛАВЛЕНИЕ..........................................................................................................2

ВВЕДЕНИЕ.................................................................................................................3

ОБОЗНАЧЕНИЯ.......................................................................................................8

ГЛАВА 1. ОБЗОР ЛИТЕРАТУРЫ......................................................................12

1.1. Многофазный поток в скважине......................................................................12

1.2. Эксплуатация скважин установками погружных электроцентробежных насосов........................................................................................................................22

1.3. Отложения парафинов.......................................................................................24

1.4. Заключение.........................................................................................................26

ГЛАВА 2. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ.....................................................27

2.1. Схема скважины.................................................................................................27

2.2. Математическая модель восходящего течения двухфазной смеси в вертикальном канале.................................................................................................31

2.3. Теплообмен скважинного потока с окружающей породой...........................45

2.4. Изменение параметров потока вследствие эксплуатации скважины установкой электроцентробежных насосов...........................................................49

2.5. Отложения парафинов.......................................................................................52

2.6. Постановка начальных и граничных условий.................................................56

2.7. Заключение.........................................................................................................57

ГЛАВА 3. ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ...................................................58

3.1. Алгоритм расчётов.............................................................................................58

3.2. Апробация...........................................................................................................62

3.3. Численные эксперименты.................................................................................65

3.4. Заключение.........................................................................................................79

ЗАКЛЮЧЕНИЕ.......................................................................................................81

ЛИТЕРАТУРА.........................................................................................................83

Приложение 1............................................................................................................93

Приложение 2............................................................................................................94

ВВЕДЕНИЕ

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

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

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

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

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

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

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

Практическая и теоретическая ценность. В диссертационную работу вошли результаты исследований, выполненных в соответствии с планами фундаментальных исследований РАН, а также в рамках государственного контракта №02.516.11.0004 по ФЦП «Исследования и разработки по приоритетным направлениям развития научно-технологического комплекса России на 2007-2012 годы», государственного контракта №119-ДОН по Тюменской областной целевой научно-технической программе и подержанных грантом РФФИ № 08-01-90300-Вьет_а, грантами Президента РФ для государственной поддержки ведущих научных школ РФ (руководитель научной школы - академик РАН Р.И. Нигматулин), грантом Губернатора Тюменской области на реализацию проекта по фундаментальным и прикладным научным исследованиям.

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

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

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

Объём и структура работы. Диссертация состоит из введения, трёх глав, заключения, списка литературы и двух приложений. Работа изложена на 94 страницах, включает 27 рисунков и 1 таблицу. Список литературы содержит 99 наименований.

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

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

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

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

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

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

• Всероссийская научно-техническая конференция «Нефть и газ Западной Сибири», Тюмень, 2007;

• VI Международная научная школа-конференция «Актуальные вопросы теплофизики и физической гидрогазодинамики», Украина, Алушта, 2008;

• VII Международная научная школа-конференция «Актуальные вопросы теплофизики и физической гидрогазодинамики», Украина, Алушта, 2009;

• International Workshop on «Thermal Hydrodynamics of Multiphase Flows and Applications», Hanoi-Vietnam, 2009;

• VII Международная конференция «Химия нефти и газа», Томск, 2009;

• Российская конференция «Многофазные системы: природа, человек, общество, технологии», посвященная 70-летию академика Р.И. Нигматулина, Уфа, 2010;

• XV Международная конференция по методам аэрофизических исследований, Новосибирск, 2010;

• X Всероссийский съезд по фундаментальным проблемам теоретической и прикладной механики, Нижний Новгород, 2011.

Результаты работы докладывались автором на семинарах Тюменского филиала Федерального государственного бюджетного учреждения науки Института теоретической и прикладной механики им. С.А. Христиановича СО РАН под руководством профессора A.A. Губайдуллина.

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

ОБОЗНАЧЕНИЯ

То, к какой фазе относится соответствующий параметр, указывается нижним индексом / или g;l- жидкая фаза, g - газовая фаза.

То, к какой компоненте относится соответствующий параметр, указывается нижним индексом Н, М или L в скобках; (Н) - тяжёлая компонента, (М) - средняя компонента, (L) - лёгкая компонента.

с - удельная теплоёмкость двухфазной смеси (м2/(с2-К));

О 0

Ci (i = /, g) - удельная теплоёмкость z'-ой фазы (м /(с -К));

О 9

са - удельная теплоёмкость материала погружного электродвигателя (м /(с -К));

2, 2

сеХ( ~ удельная теплоёмкость окружающей скважину породы (м /(с -К));

9 О

Cext(i) ~ удельная теплоёмкость мёрзлой породы (м /(с -К)); сех, (2) - удельная теплоёмкость протаявшей породы (м2/(с2-К)); D - коэффициент диффузии (м2/с); Ei (i = I, g)- удельная полная энергия z'-ой фазы (м2/с2);

Eig - интенсивность обмена энергией между жидкой и газовой фазами в единице объёма смеси и в единицу времени (кг/(м-с3));

Fig - составляющая межфазной силы, зависящая от скоростной

неравновесности между жидкой и газовой фазами, отнесённая к единице 2 2

объёма (кг/(м -с )); Fr - число Фруда;

Fw - сила трения между двухфазным потоком и стенками канала, отнесённая к

2 2

единице длины канала (кг/(м -с )); g - ускорение свободного падения

(м/с );

G(l) ~ постоянная Генри (кг/(м-с2)); GF- газовый фактор (м3/т); Н- глубина скважины (м);

Нр - глубина подвески электронасосного агрегата (м);

Нимп - глубина многолетних мёрзлых пород (м);

2 2

ii (i = l, g) - удельная энтальпия z'-ой фазы (м /с ); У/g - интенсивность газовыделения (кг/(м3-с));

JigiL) - интенсивность выделения растворённого в жидкости газа (кг/(м3-с)); J\g{M) ~ интенсивность испарения жидкой фазы (кг/(м3-с)); J is - интенсивность твёрдых отложений (кг/(м3-с)); js - интенсивность образования парафиновых отложений (кг/(м2-с)); К- коэффициент продуктивности скважины (м-с);

ki(j) (i = l> £> j = M,L)- массовая концентрация в z'-ой фазе у'-ой компоненты; кцщ ~ текущая концентрация тяжёлой компоненты; кцще - равновесная концентрация тяжёлой компоненты; L - длина окружности поперечного сечения канала (м);

Lt (z = I, g) - длина части окружности поперечного сечения канала, соприкасающейся с z'-ой фазой (м);

1( л (у = М,Ь) - удельная теплота перехода у'-ой компоненты из жидкой фазы в 2 2

газовую (м /с );

- удельная теплота парафинообразования (м2/с2);

^(1,2) ~ удельная теплота протаивания мёрзлой породы (м2/с2); Ьр - длина электронасосного агрегата (м);

Ьц - длина погружного двигателя установки электроцентробежных насосов (м); т - массовый расход смеси (кг/с);

тр - расход погружного насоса (массовый расход двухфазной смеси при откачке погружным электроцентробежным насосом) (кг/с); т(р> - массовый расход лёгкой компоненты (кг/с);

Щь)Р ~ массовый расход лёгкой компоненты при откачке погружным

электроцентробежным насосом (кг/с);

т, (г = /,£-)- массовый расход /-ой фазы (кг/с);

т1{]) (г = I, g; у = М, Ь) - массовый расход у'-ой компоненты, которая содержится в г-ой фазе (кг/с);

Ма - масса погружного двигателя (кг); Ы- мощность погружного двигателя (кг-м2/с3); АЫ - число Нуссельта; р - давление двухфазной смеси (кг/(м-с2));

и = М,Ь) - парциальное давление у'-ой компоненты в газовой фазе (кг/(м-с2));

рь - давление в пласте (кг/(м-с )); р^ давление на забое скважины (кг/(м-с2));

р1р - давление на забое скважины при откачке погружным электроцентробежным насосом (кг/(м-с2));

Ар - скачок давления двухфазной смеси вследствие работы электронасосного агрегата (кг/(м-с ));

р* - эмпирический аппроксимационный параметр; Рг - число Прандтля;

Qw - интенсивность отвода тепла, отнесенная к единице длины канала (кг-м/с3); д, (г = /, - поток тепла от ¿-ой фазы через стенки канала, отнесённый к единице длины окружности поперечного сечения канала (кг/с3); ца - интенсивность теплообмена между двухфазным потоком и поверхностью

о

твёрдой фазы (кг/с );

Я - радиус поперечного сечения канала (м);

- радиус поперечного сечения канала (м); Ям - радиус скважины (м);

^нкт - радиус насосно-компрессорных труб (м); /?эк - радиус труб эксплуатационной колонны (м); Яс - радиус скважины (м);

- радиус погружного двигателя (м);

- газовая «постоянная» газовой

фазы (м2/(с2-К));

/?(/) (у = М, Ь) - газовая «постоянная» у'-ой компоненты в газовой фазе (м /(с2-К));

R( л° (j = M,L) - газовая постоянная /-ой компоненты в газовой фазе, если

2 2

рассматриватьу-ую компоненту, как идеальный газ (м /(с -К)); Re - число Рейнольдса;

Rig - межфазная сила из-за сил трения, давления, сцепления между жидкой и газовой фазами, отнесённая к единице объёма (кг/(м2-с2));

Rgi - межфазная сила из-за сил трения, давления, сцепления между газовой и жидкой фазами, отнесённая к единице объёма (кг/(м2-с2)); г - радиальная координата (м); S - площадь поперечного сечения канала (м2);

St (i - /, g) - площадь поперечного сечения канала, занимаемого z'-ой фазой (м2); Sh - число Шервуда; t - время (с);

Т- температура двухфазной смеси (К);

Ti - температура внутренней стенки канала либо твёрдой фазы на внутренней стенке канала (К);

Td - температура погружного двигателя (К);

AT - скачок температуры двухфазной смеси вследствие работы электронасосного агрегата (К);

Те - температура начала кристаллизации тяжёлого компонента (К); Text - температура окружающей скважину породы (К); Text (i) - температура окружающей скважину мёрзлой породы (К); Text (2) - температура окружающей скважину протаявшей породы (К); Text(\,i) ~ температура фазового перехода в мёрзлых породах (К); Tf - температура окружающей скважину породы на забое (К); Tgeo ~ геотермальная температура (К);

Th - температура окружающей скважину породы на устье (К); TN - температура на внешней границе скважины (К);

Тг - температура в интервале от внутренней стенки канала, по которому движется двухфазная смесь, до внешней границы скважины (К); Ts - температура в слое твёрдых отложений (К);

Tw - температура внутренней стенки подъёмной колонны скважины (К); Т„ - температура на поверхности твёрдой фазы (К); Т* - эмпирический аппроксимационный параметр;

2 2

и¿ (i = I, g) - удельная внутренняя энергия /-ой фазы (м /с ); Vi (i = /, g) - скорость z'-ой фазы (м/с);

v¡g - скорость массы, претерпевающей фазовый переход из жидкости в газ и находящейся в газовой фазе (м/с);

Vg¡ - скорость массы, претерпевающей фазовый переход из газа в жидкость и находящейся в жидкой фазе (м/с);

W- общая приведённая скорость двухфазной смеси (м/с); W¡ - приведённая скорость жидкой фазы (м/с); Wg - приведённая скорость газовой фазы (м/с); z - вертикальная координата (м);

z(y) (у =M,L)~ коэффициент сжимаемостиу'-ой компоненты в газовой фазе;

и

а - объёмное газосодержание;

/3 - коэффициент теплопередачи через систему труб скважины

(кг/(с -К));

рс1 - коэффициент теплопередачи между погружным двигателем и двухфазным потоком (кг/(с3-К));

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