Расчет многофазных смесей

Курсовая  работа

На тему «Расчет многофазных смесей»

Содержание

С.

Введение 4

1.Конструктивная схема и принцип работы двухвинтового насоса 6

2. Математическая  модель двухвинтового мультифазного насоса 7

3. Работа  насоса на многофазной смеси 10

4. Силы, действующие на винты  13

5. Изменение    величины    утечек    из-за  изменения    формы

радиального зазора 14

6. Определение  давлений в камерах винтов 15

7. Характер    распределения     давления    вдоль    винта     при

перекачивании жидкости  16

8. Характер   распределения    давления    вдоль     винта     при

перекачивании    смеси жидкости и газа   16

9 Расчётная  часть 17

Список  литературы   27

Введение

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

Мультифазная технология, обладает рядом важных достоинств, среди которых можно выделить:

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

-Продление срока эффективной эксплуатации истощающихся месторождений.

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

-Снижение негативного воздействия на окружающую среду за счёт  эффективного использования попутного газа и отсутствие факела.

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

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

Способностью  перекачивать многофазные жидкости обладают несколько типов насосов, а именно,  одновинтовые, двухвинтовые, поршневые и гелико-осевые (helico-axial). Наибольшее распространение среди этих типов насосов получили двухвинтовые насосы.

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

-  Кинематическая вязкость перекачиваемой жидкости от 10^-6 до 10^-3 м²/с (1÷135°Е).

- Содержание попутного газа до 90%  по объему при условиях на входе в насос.

- Содержание попутного газа до 100% (газовые пробки) в течение 5 мин.

- Содержание механических примесей в перекачиваемой среде – 0,2% максимальный размер частиц до 1 мм.

- Содержание сероводорода в газе до 2%.

- Давление на входе в насос может достигать 25 бар (2,5МПа).

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

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

1 Конструктивная схема и принцип действия двухвинтового насоса

Конструктивно двухвинтовой насос представляет собой  два винта вставленных в обойму (рис. 1,2) [1,2]. На каждом винте имеется  два участка винтовой нарезки  направленные навстречу друг другу. Жидкость через всасывающий трубопровод  подводится в корпус насоса, где  разделяется на два потока, захватывается  винтами и попадает в камеры, образованные нарезками винтов и обоймой и  переносится в этих камерах в  осевом направлении к центру, где  находится  полость нагнетания, 

Рисунок 1 - Конструктивная схема мультифазного двухвинтового насоса

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

Рис.2. Гидравлическая часть двухвинтового мультифазного насоса

Рисунок 2а - Винты двухвинтового насоса и направление движения потока жидкости

2 Математическая модель двухвинтового мультифазного насоса

Двухвинтовой  насос относится к насосам  объёмного типа. Подача насосов объёмного  типа определяется геометрией рабочей  камеры и частотой совершения рабочих  циклов.

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

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

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

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

Q_насос= Q_геом - 4Q_l1

Коэффициент 4 соответствует общему количеству винтовых участков нарезки на винтах, по два на каждом винте (рисунок 2).

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

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

В работе G.Vetter и M.Wincek 1993г. [3] экспериментальным путем было установлено, что через диаметральные зазоры проходит около 80% утечек, остальная часть проходит по боковым зазорам между винтами. Поэтому в формулу для определения количества утечек через зазоры введен поправочный коэффициент, учитывающий это. Соответственно полное количество утечек можно определить, умножив выражение для определения утечек через диаметральные зазоры на коэффициент 1/f, где f=0.8.

Работу двухвинтового насоса можно представить как работу многоступенчатого насоса, где каждая камера, образованная нарезками винтов и обоймой, представляет собой своеобразную ступень. Общее дифференциальное давление (разница между давлением на выходе и входе в насос), создаваемое насосом, складывается из давлений, создаваемых каждой ступенью в отдельности. При перекачивании только жидкости количество жидкости перетекающей между камерами одинаково Q_l1=Q_l2=…=Q_ln, поэтому разность давлений p₁=p₂=…=p_n=(P_вых-P_вход)/n  в каждой  камере  будет  одной  и  той  же  (n – число замкнутых камер на каждом участке винтовой нарезки).

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

- местные потери при входе жидкости в зазор;

- потери на трение по длине зазора.

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

Для величины утечек жидкости через зазор по наружному  диаметру винтов получим следующее  выражение:

,

где ρ-плотность жидкости; ζ-коэффициент потерь при входе жидкости в зазор; λ-коэффициент потерь на трение зависит от рода жидкости и от режима течения жидкости в зазоре; d_k-гидравлический диаметр d_k=2 δ; δ-величина радиального зазора между винтами и обоймой; l-ширина ленточки нарезки винта; f_зазор-площадь зазора, .

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

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

          Для определения режима течения  используется критерий Рейнольдса. Принимается, что при числе Рейнольдса меньшем 1200 режим течения – ламинарный, при числе Рейнольдса большем 1200 – турбулентный.

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

Для ламинарного  режима коэффициент потерь вычисляется  по формуле .

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

.

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

Коэффициент местных потерь при входе в  зазор ζ принимаем равным  1,5.

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

3 Работа  двухвинтового насоса на многофазной  смеси

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

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

Рисунок 5- Схематичное представление движения двухфазной жидкости в камерах винта

Распределение давления вдоль винта в этом случае будет иметь параболический характер.

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

Рисунок 6 - Распределение давлений в камерах винта при перекачивании жидкости с большим содержанием газа

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

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

  - расход смеси газа и жидкости,

где  Q_ж-  количество жидкости, поступающей в насос,

Q_г.вх- количество газа на входе.

Объёмный  к.п.д. насоса равен-  .

где Q_геом- геометрическая подача насоса.

- количество газа, поступающего в насос с учетом давления на входе.

  - газосодержание в перекачиваемой смеси.

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

где  p_iaбс- абсолютное давление в i-ой камере винта.

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

4 Силы, действующие  на винты

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

Рис.4. Распределение  давления жидкости вдоль винта

В работе [4] предложена формула для определения  радиальной силы, действующая на каждый винт:

где t - шаг нарезки винта; D₂- наружный диаметр нарезки винта; Δρ-дифференциальное давление насоса; ξ-отношение внутреннего диаметра нарезки к наружному ;

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

Рисунок 7 - Форма диаметрального зазора при деформации винтов под действием давления

5 Изменение  величины утечек из-за изменения формы радиального зазора

 Величина утечек через эксцентрический зазор определяется известным из гидравлики соотношением :

,

где e - величина эксцентриситета, Q₁- расход жидкости через концентрический зазор.

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

6 Определение  давлений в камерах винтов

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

где  Q_{l i} - количество утечек через зазор в i-ой камере рисунке 3; Q_газi - количество газа в i-ой камере; E_ж - модуль объёмной упругости жидкости; V_t -объём жидкости в i-ой камере.

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

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

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

,

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

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

Количество  жидкости в i-ой  камере можно представить, как баланс утечек поступающих из (i+1)- ой камеры и выходящих из данной камеры в (i-1) - ю, а также уменьшением объёма газа (рисунка 3, 5).

7 Характер распределения давления вдоль винта при перекачивании жидкости

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

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

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

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

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

Список  литературы

1. Рязанцев  В.М. Роторно-вращательные насосы  с циклоидальными зацеплениями. М.:Машиностроение, 2005г. – 346 c.

2. Bornemann pumps. Multiphase Pumps@Systems

3. Vetter, G and Wincek, M., 1993 “Performance Prediction of Twin Screw Pumps for Two-Phase Gas/Liquid Flow”, Pumping Machinery -1993, FED-154, ASME, pp. 331-340.

4. Allan J. Prang, Paul Cooper. Enhanced multiphase flow predictions in twin screw pumps. Proceedings of the twenty first intrenational pumps user symposium, 2004, p.p. 69-76.

5. ОАО  «Ливгидромаш», Каталог насосов для нефтедобывающей и нефтеперерабатывающей отраслей, 2006г.,  28 - 32стр.