Многофазный поток в скважинах
..pdfМоделирование двухфазного потока и расчет градиента давления будет успешным только в том случае, если каждую составляющую градиента давления рассматривать как функцию, зависящую от режима потока. Гравитационная составляющая напря мую зависит от плотности р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 трения |
^ I’L 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) |
где 5С и 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