Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Многофазный поток в скважинах

..pdf
Скачиваний:
135
Добавлен:
15.11.2022
Размер:
17.52 Mб
Скачать

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

д р \

Р т рд,

(4.275)

dL /)

гравит.

 

где

 

 

Р т р = P L H L + р д { 1 - H i ) .

(4.276)

Значение объемного содержания жидкости H i определяют по уравнению (4.273). Составляющей градиента давления по ускорению в пузырьковом режиме потока

можно пренебречь, то есть

I)/ уск.

 

(4.277)

 

 

 

Составляющая градиента давления по трению в затрубном пространстве рассчи­

тывается из уравнения

 

 

 

(*Е)

4 / '

V l

(4.278)

V .^ 'V трения

 

 

 

где скорость смеси vm равна

 

 

(4.279)

V-m = V S L

+ V s g .

 

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

PTPVmdh

NReTP =

Р т р

Вязкость смеси дтр вычисляют в соответствии с объемным содержанием фаз:

Р т р = P L ^ L + Р д { 1 — А/,).

Здесь объемное содержание жидкости без учета эффекта проскальзывания A i ляется следующим образом:

V S L

Аъ

V S L + VS g

(4.280)

(4.281)

вычис­

(4.282)

Объединяя уравнения (4.274), (4.275), (4.277) и (4.278), получаем следующее выражение для общего градиента давления:

dp

= рт рд +

4 / '

(4.283)

dL

-т -рт р

 

dh

 

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

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

В случае, когда толщина пленки является постоянной, модель строится на осно­ ве упрощенной модели Фернандеса и др. [39] для вертикального восходящего потока в трубе. Впоследствии Каэтано и др. расширили область применения модели, учитывая характеристики пробкового режима потока в затрубном пространстве. В том случае, когда пузырек Тейлора состоит из одной шапки, а толщина пленки варьируется, мо­ дель опирается на метод МакКуиллана и Уэлли [44], описывающий процесс перехода из пузырькового в пробковый режим потока. Ранее уже упоминалось, что длина труб добывающих скважин настолько велика, что на всем протяжении скважины можно счи­ тать пробковый режим потока установившимся. Такой подход может привести лишь к незначительным ошибкам. Поэтому мы не будем рассматривать модель развивающе­ гося пробкового режима потока.

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

С к о р о с т ь п о с т у п а т е л ь н о г о п е р е м е щ е н и я п у з ы р ь к а Т е й ­ л о р а . Каэтано и др. пришли к выводу, что скорость подъема пузырька Тейлора можно

рассчитать, применив концепцию эквипериферического диаметра:

 

 

V TB = 1,2VM + 0,345 у/g(dc + dt).

(4.284)

Б а л а н с м а с с

в п р о б к е и п л е н к е . Баланс масс в жидкой фазе между

сечениями А -А и В-13 (рис. 4.25) выражается соотношением

 

 

(V TB - V L L S )HL LS = (VT B + V L T B )HL T B ,

(4.285)

где VI T B РеальнаЯ

скорость жидкости в пленке, окружающей

пузырек Тейлора.

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

Здесь Н ьтв — объемное

содержание жидкости по всему сечению трубы

с учетом

толщины пленки.

 

 

Аналогично рассчитывается баланс масс для газовой фазы:

 

(V TB vgLs)( 1 - H L LS) = (VTB + V9T B )( 1 - H L T B )-

(4.286)

О б л а с т ь п р о б к и

ж и д к о с т и . Предполагается, что газ и жидкость в рай­

оне пробен жидкости взаимодействуют так же, как в пузырьковом потоке. Опираясь на метод Барни и Браунера [50], Каэтано и др. заключили, что объемное содержание жидкости в обдасти пробки можно рассчитать так же, как и объемное содержание жидкости на гРанице между пробковым и пузырьковым режимами потока. Принима­ ется, что значение объемного содержания жидкости на границе пузырькового и проб-

vsl — приведенная скорость жидкости vSlJ - приведенная скорость газа

VLLS реальная скорость жидкости в пробке жидкости v,)Ls ~ реальная скорость газа в пробке жидкости

VLTH ~ реальная скорость жидкости в пленке жидкости

V,JTH реальная скорость газа в пузырьке Тейлора vTH — переходная скорость пузырька Тейлора

HLLS — объемное содержание жидкости в пробке жидкости Lls - длина пробки жидкости

Lu, - длина пленки жидкости Lsv - длина блока пробки

Рис. 4.25. Окончательно установившийся пробковый поток [67]

кового режимов потока постоянно. Поэтому объемное содержание жидкости в проб­ ке также постоянно и совпадает со значением в переходных условиях. Таким обра­ зом, объемное содержание жидкости в области пробки в концентрическом затрубном пространстве H u s имеет значение 0,80, а в эксцентрическом затрубном простран­ стве — 0,85.

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

Поскольку в любой точке пробки жидкости общий объемный дебит постоянен,

баланс масс для пробки жидкости описывается выражением

 

( V S L + vsg)A = v u s H n s A + V g b s il - H L L S ) A .

(4.287)

Здесь А представляет собой площадь сечения потока в затрубном пространстве:

A = ^ - d f ) .

(4.288)

Объединяя уравнения (4.177) и (4.287), получаем следующее выражение для реальной скорости жидкости в пробке:

 

 

 

1/4

 

 

 

(,PL ~

Pg)9CL

(4.289)

V L L S =

( V S L + V S g )

- 1 , 5 3

H L L 1/2(1 - H LLs).

 

 

 

PI

 

О б щ и й

б а л а н с

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

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

 

L LS

 

L L F

(4.290)

V S L = V L L S H L L S Lsu

-

VL T B H LTB L su

где длина блока пробки L su равна

 

 

 

 

 

Lsu = L

L F

+

S

(4.291)

 

 

L L -

 

Объединяя уравнения (4.290) и (4.291), получаем:

 

L L S _

VSL + V L T B H L T B

(4.292)

Lsu

V L L S H L L S + V L T B H L T B

 

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

L L S _

VSg + VgTBC L - H L T B )

(4 293^

L s u

V 9 L S { 1 - H u s ) + V 3T B { 1 - H L T B )

 

О б л а с т ь п л е н к и

ж и д к о с т и . Можно считать, что пленка жидкости, окру­

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

£ (1 - С м ) / С м

VL

(4.294)

V L T B =

Ц С м 4 p L

1/3

 

А

С к

,g{PL-Pg)PL.

Значения индексов С к и С м зависят от режима потока в зоне пленки. Число Рейнольдса ДЛЯ ЛЛенки жидкости берут равным

N n . , TB =

,4.295)

Для случая ламинарного потока (N ^ LTB < 1000) индексы Ск и См были определены У°ЛЛИсом [18] аналитическим путем, они равны 0,9086 и 1/3 соответственно. Для

турбулентного потока (N^eLTB > 1000) Фернандес и др. [39] предложили использовать значения 0,0682 и 2/3. Эти же значения Брод [41] получил экспериментальным путем.

Если пузырек Тейлора является полностью развитым, то толщина пленки име­ ет постоянное предельное значение и объемное содержание жидкости в зоне пленки можно определить по геометрии потока в зависимости от толщины пленки 5:

 

4S(dc -

<5)

(4.296)

 

d2c ~ d j

 

 

С р е д н е е о б ъ е м н о е

с о д е р ж а н и е

ж и д к о с т и . Среднее объемное со­

держание жидкости в блоке пробки H is u рассчитывается по уравнению

 

H L S U = (

H L L S + (1 - у— ) H L T B -

(4.297)

\ L s u J

L su

 

Отношение L i s / L s u определяют из (4.292) или (4.293), а значение H I T в

— по урав­

нению (4.296) для полностью развитого пузырька Тейлора. Толщину пленки 5 рассчи­ тывают итерационным методом.

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

Гравитационная составляющая градиента давления в блоке пробки зависит от плот­ ности проскальзывания и равна

(4М8>

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

PL S = PL H U S + рд(1 - H LLs).

(4.299)

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

= PL —HLLL(уьтв + V T B ){V L T B + V L L S )-

(4.300)

LSU

 

Когда пузырек Тейлора полностью развит, Н ы в — это среднее объемное содержание жидкости, а У ы в — средняя скорость пленки.

Составляющая градиента давления по трению равна

(4.301)

\ d L J трения

^ IL S {VS‘ + V S L ) , { I ^

 

Значение коэффициента трения Фаннинга / ' определяют методом, описанным в главе 2. Соответствующее ему число Рейнольдса для блока пробки рассчитывают по уравне­ нию (4.280), в котором значение ртр заменяется на значение p i s . Плотность в блоке пробки рассчитывается по уравнению (4.299).

Объединяя уравнения (4.298)-(4.301), получаем следующее выражение для вычис­ ления общего градиента давления:

d p \

2f

 

 

М ) о бш.

P L s { v s g + VS L ) 2 L L S + P b S g L b S +

 

 

dh

 

 

 

+ P L H L T B (;V L T B + V T B ) { V L T B + V L L S )

1

(4.302)

 

 

L s u

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

Г~>

А

А'

с

i ' ^ ’г

^ к Газ и жидкость

Периметр внешней

Толщина внешней пленки

Пленки со стороны

жидкости в любой ее точке, <5(.

обсадной трубы^ S,.

Однородная смесь газа

 

Периметр Внешней пленки

и жидкости

Толщина внутренней пленки

на границе со смесью

газа и жидкости Srl

жидкости в любой ее точке, 6,

Периметр внутренне

пленки на границе Со сме с ы о Р а з р е з А -А газа и >ки^кости, Sh

Периметр внутренней пленки со сТоРоПы стволовОй трубы^ S,

Рис. 4.26. Кольцевой режим потока в концентрическом затрубном пространстве [67]

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

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

dр \

Sc

S a

n

(4.303)

хг)

+ Т с ~л---- TCiH

P L 9 — 0,

d L ) об.пл.

A c f

A c J

 

 

где (dp/dL)cj — общий градиент давления для внешней пленки жидкости, тс и та — силы касательного напряжения на стенке обсадной трубы и на границе между плен­ кой жидкости и смесью газа и жидкости. S c и Sci — периметры внешней пленки со стороны обсадной трубы и внешней пленки на границе со смесью газа и жидкости соответственно. A cj — общая площадь сечения внешней пленки жидкости.

Аналогично выписывается условие сохранения линейных импульсов для внутрен­ ней пленки жидкости:

dp

,

S t

Sti

(4.304)

dL

+ п ~аГ,

” Ti; л 7 '>W =

 

Все параметры данного уравнения имеют те же значения, что и в предыдущем уравне­ нии.

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

ф \

S a

Su

,

n

(4.305)

+ т

а ^

- г „ —

+ р С9 =

°,

dL ) внутр.

 

 

 

 

 

где (dp/dL) внутр. — общий градиент давления для смеси внутри затрубного простран­ ства, рс — плотность смеси внутри затрубного пространства, А с — площадь затрубного пространства, занятая смесью.

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

/ ф \

dp

d p \

\ ^ /

dL

(4.306)

внутр.

об.пл.

С учетом данного уравнения равновесия, а также уравнений сохранения линейных импульсов (4.303)-(4.305), мы приходим к следующим выражениям для суммарного импульса:

~ тс - г £ -

+

A CJ

+

TCil T '

+

ТЧ 1 Г

~ (pL ~ Рс^д = 0

(4.307)

Acf

 

 

A c

 

A c

 

 

и

 

 

 

 

 

 

 

 

 

+ TliX i

+

TCi%

+

TliJ 3

~ {pL ~ pc)9 = °-

(4-308)

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

Ac/ 7T<*c(dc

(4.309)

и

А С = \ [d2 c- d j - 45c(dc - 6С) - 46t(dt + <5,)],

(4.311)

где и St — толщина внешней пленки и толщина внутренней пленки соответствен­ но.

Смоченные периметры внешней и внутренней пленок (как со стороны контакта со стенкой трубы, так и со стороны контакта с другой фазой) равны:

Sc — 7Гб?с,

Sci — тг(dc 25с), Со II

и

Su = тг(dt + 26t).

(4.312)

(4.313)

(4.314)

(4.315)

Основываясь на понятии гидравлического диаметра и используя уравнения (4.309)—(4.315), можно получить соотношения для гидравлического диаметра внеш­ ней пленки жидкости, внутренней пленки жидкости и внутренней части затрубного пространства:

dcfh ~ 4<5с ( l

d c ) ’

II

+ 1—т

и

 

 

[ 4 - 4 - 46c(dc -

6C) - 4St(dt + $t)]

(dc —25C) +

(dt + 2St)

(4.316)

(4.317)

(4.318)

Сила касательного напряжения на стенки равна

т = / 'р у ,

(4.319)

где т — сила Касательного напряжения на стенке, р — плотность той фазы, которая смачивает стен*<у? v — реальная средняя скорость фазы. При этом коэффициент трения Фаннинга / ' определяют, пользуясь формулой Блазиуса:

/ ' = C N £е,

(4.320)

где п принимает значение — 1 для ламинарного потока и п = -0 ,2 5 для турбулентного потока. Коэффициент С рассчитывают с учетом геометрии затрубного пространства (см. главу 2). Г[ри в ы ч и с л е н и и ЧИсла Рейнольдса необходимо использовать значения диаметров dcf}t или dt/д в зависимости от того, где рассчитывается сила касательного напряжения.

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

V 2

 

Гг = f'iPc- f ,

(4 .3 2 1 )

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

Pc = PL H LC + рд( 1 - Ньс)

(4.322)

и

 

VC = v s c j ^ -

(4.323)

Значения H LC и VSC рассчитывают по уравнениям (4.332)

и (4.338), значение А с —

по уравнению (4.311), А — площадь поперечного сечения затрубного пространства. Каэтано и др. модифицировали корреляцию Уоллиса для коэффициента межфазно­

го

трения,

которая используется для

тонких пленок

или для

случаев интенсивно­

го

захвата

газом капелек жидкости

в вертикальном

потоке

затрубного простран­

ства (уравнения (4.220) и (4.221)). В результате были получены следующие выражения для межфазных коэффициентов трения:

f'i

= f'sc

( l

+ 3 0 o Q

(4 .3 2 4 )

И

 

 

 

 

f'i

= f'sc

( l

+ 3 0 o | ) ,

(4.325)

где f'Sc — коэффициент трения, полученный по уравнению (4.320) для числа Рейнольд­ са, вычисленного по приведенной скорости смеси, плотности смеси и гидравлическому диаметру внутренней части затрубного пространства (уравнение (4.318)).

С учетом формулы для касательного напряжения и геометрических параметров системы (уравнения (4.309)-(4.325)) приходим к следующему уравнению сохранения линейного импульса для внешней пленки жидкости и смеси внутри затрубного про­ странства:

,/

Vcf

 

+

f'sc{ 1+ mtc)Pcf

 

 

 

- J c P L

2 Sc(dc -

6c)

 

 

 

 

[ ( 4

-

2SC)

4 ( 4

-

2<5C)

 

 

+

 

 

 

 

+

 

 

 

 

 

 

4 ( 4 -

4 )

[d2c - d 'f - 4 4 ( 4

-

4

) -

4 4 ( 4

+ 4 ) ]

+ fs c [ 1 + з о о | ) ^

 

 

4 ( 4

-

2 4 )

 

 

 

 

 

 

 

 

 

 

 

[d2c - d 2 -

4

4

( 4

-

4 ) -

4 4 ( 4 + 4 ) ]

- ( P L - p c ) g = 0- (4.326)

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

Jtf

f'iPL

+

f'sc

( 1 + 3 0 0 ^ -)

pc~§-

2

Stidk + St)

 

dt J

^

2

 

I (<k + 26t)

 

4(dt + 2SC)

 

6t(dt -

+

 

 

+

 

St)

[dl - d \ - 46c(dc - Sc) - 4St(dt + 5t)]

+

f'sc ( 1 + 300“k ) PC ~Y

 

4{dc - 25c)

 

 

 

 

 

[dc ~ d%—4Sc(dc — Sc) — 4St(dt + St)]

 

 

 

 

 

~ { p L ~ P c )g = 0. (4.327)

У р а в н е н и я н е п р е р ы в н о с т и

фа з . Непрерывность жидкой фазы, пред­

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

VSL = v t f H b t f + v cf H Lcf + V S L F E ,

(4.328)

где HLft реальное объемное содержание жидкости во внутренней пленке, Н и/ — реальное объемное содержание жидкости во внешней пленке и F # — долевой объем жидкости, захваченный газовым ядром.

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

VSg = Vgc{ 1 -

(4.329)

где vgc — реальная скорость газа в затрубном пространстве, Н

— общее содержание

жидкости в затрубном пространстве.

 

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

ний (4.309) и (4.288) получаем:

 

H Lcf = П 1-!)

(4.330)

<4(1- к 2)

'

где К — соотношение диаметров труб dt/dc, образующих затрубное пространство. Аналогично рассчитывается реальное объемное содержание жидкости во внутрен­

ней пленке:

1+!

 

Н и / = ~4 г К

(4.331)

dc

(1 - К 2)

 

Реальное объемное содержание жидкости внутри затрубного пространства, заполнен­ ного однородной смесью газа и жидкости, равно

тт_ VS L FE

LC V S L FE + vsg