1470
.pdfАрр =10"3 - (187,4 —0,154 • р,) |
------ 1 |
(8.33) |
|
Ч Р ат |
J |
Величина р, называется фиктивной потому, что легкие углево дороды, азот, углекислый газ при стандартных условиях не могут быть жидкостью.
Плотности индивидуальных углеводородов pj для подстановки в формулу (8.32) берутся из табл. 8.1.
Величины фиктивных плотностей метана и этана в общем слу чае зависят от компонентного состава жидкой фазы и вычисляются по формулам
Pc, —(l 0,1-Хс.)
8,24 + 2,14 -М,
(
МсУ*
Рс = !°3 0,242 + ________ I (8.34) 125,27 + 2,04-М С3+;
где Мс - молярная масса жидкой фазы без метана и не углеводо
родных компонентов; Мс |
- величина М с |
без этана. |
|
Более просто величины рС1 и рС2 |
находятся по методике Сиб- |
||
нии нп |
|
|
|
Р С| |
= 0,46 • р, - |
1 3 ,8 ; |
|
рС2 = 0,30 -рС) +248,0. |
(8.35) |
Однако в этом случае для нахождения р, по формуле (8.32) тре буется решать трансцендентное уравнение.
Учитывая незначительное содержание метана и этана в жидкой фазе при интересующих нас давлении и температуре, с достаточной для инженерных расчетов точностью можно принять рС1 =459,3 кг/м3, рС2=519,5 к г / м 3.
Плотность условного компонента С6+ в жидкой фазе вычисляет
ся по формуле Крего |
|
|
|
|
103-М С |
|
|
____________ ч+_____ |
(8.36) |
||
Рсб+ 47,409 +0,953-Мс |
|||
|
|||
* |
4+ |
|
10. Б-762 |
289 |
Плотность газовой фазы при условиях фазового перехода
Р |
273 |
|
Рг —Рну ‘ |
(8.37) |
|
Т ’ |
||
Z C * Р) |
где рнуплотность газа при нормальных условиях
Р<у=мг/ 22,4;
zcкоэффициент сжимаемости газа
0,0241 ■Рпр
zn =1 — (8.38)
1 -1 ,68-Т11|)+ 0,78 -Т ’„+0,0107-Т:пр
пр
Рпр, Тпрприведенные соответственно давление и температура
Рпр |
(8.39) |
Рм Д ф - критические давление и температуры i-ro |
компонента |
(табл.8.1). |
|
Для условного компонента С6+ критические параметры вычис
ляют по корреляционным зависимостям |
|
Р .* . = 7 ,7 7 -9 ,5 -1 0 -’ - Т ^ ; |
|
Т , ^ = 19,25 ■(in MCfc )2 + 44,06 - In Мс__ - 7 0 . |
(8.40) |
Динамическая вязкость жидкой фазы приближенно может быть найдена по формуле Кондела-Монроэ
|
(8.41) |
\ i=l |
у |
где Pj - динамическая вязкость i-ro компонента (табл. 8.1).
Для фиктивного компонента Сб+ величину динамической вязко сти можно вычислить по формуле СибНИИ НП
|
10 |
0 .0 1 0 5 -М г |
- 2 ,8 |
^ = 1 7 , 1 - М ^ |
Сб+ |
(8.42) |
290
Для жидкой фазы, содержащей растворенный газ, можно также пользоваться формулой
8,6 Ю-3 |
КР |
|
(1 .0 3 -р ж/ |
при |
Рж < 782— |
р .)3-75 |
м |
|
М-s |
|
(8.43) |
4,69-10"5 |
^-700 КГ |
|
[(1,03 - р, / |
при |
рж> 782— |
р .)7-35 |
м |
а для жидкостей, разгазированных при атмосферном давлении
1,12-10'
М*о = |
9,09 ’ |
(8.44) |
(1 .0 3 -р ,о/ |
р .) |
|
где р„ - плотность воды.
Вформулах (8.42)-(8.44) динамическая вязкость рассчитывается
вмПас.
Учитывая приближенный характер расчета р по формулам (8.43), (8.44), можно рекомендовать для нахождения динамической вязкос ти при температуре Т подставлять в них плотность жидкой фазы при этой же температуре.
Динамическая вязкость газовой фазы при нормальных условиях находится по формуле
i X |
i - y i - л /н |
|
Цг = ^ |
----------------• |
(8-45) |
Ъ г ф л ,
i= l
Вязкость условного компонента С6+ в газовой фазе можно при нять равной 0,0057 мПа-с.
Давление насыщения жидкой фазы при заданной температу ре находится методом последовательных приближений из урав
нения |
|
^ K f •х1=1. |
(8.46) |
291
§ 8.4. Расчет параметров нефтегазовых и газоконденсатных смесей по результатам разгазирования проб
Нередко динамическая вязкость и плотность де газированной нефти (или стабильного газового конденсата) заданы. В этом случае свойства нестабильных жидкостей рассчитывают в зависи мости от количества растворенного газа и с учетом его плотности.
Количество растворенного газа характеризуется величиной га зового фактора Г. Это отношение объема газа, выделившегося из нестабильной жидкости при стандартных условиях, к объему жид кой фазы, приведенное к нормальным условиям.
Расчетным путем величина газового фактора при нормальных условиях находится как
_ |
Мг ■р |
• е |
(8.47) |
г __ |
I г Ж |
с т |
|
|
|
|
М . - Р н у - О - О ’
где £„ - мольная доля газовой фазы , найденная по формуле (8.17) при стандартных условиях.
Величина газового фактора связана с давлением насыщения Р5 формулой
<8‘48)
где хр, ур - эмпирические коэффициенты, определяемые по резуль татам разгазирования проб нестабильной жидкости.
При Ур=1, коэффициент Хр=Кр коэффициенту растворимости газа в углеводородной жидкости при Т=293К. Приближенная оценка ве личины Кр в нм3/(м 3-МПа) для нефти может быть сделана по формуле
К„ * 115,4 • 101г • Ар61 ■ех р (-7 ,76 • Дн) , |
(8.49) |
где Дн , Дг —относительные ( по воде ) плотности нефти (при 293 К) и растворенного газа.
Более точно, хотя и более сложно, коэффициент растворимости
газа в нефти при 293К в нм3/(м 3 |
МПа) может быть вычислен следу |
|
ющим образом |
|
|
Кр = 0,0093 + 0,934 |
■t, + 3,634 • 10~5 • t, • t2 , |
(8.50) |
где t, , t2 —расчетные коэффициенты:
t, = 0,172 - 0,030z2 +1,006z, + 0,624 • 10 '3 • z, ■z2;
292
t 2 = 0 , 1 7 9 + 0 , 999z3- 0 , 0 2 4 z 4 + 0 , 6 5 4 • 10 "3 z 3z 4 ;
z, = 2,11 2 - 0,345y, + 0,678y3 + 0,04 ly,y3 + 0,904 ■1 O '2 yf - 0,023y2;
z 2 = 2 , 5 3 2 + 0 , 3 6 0 y 2 - 0 , 7 2 6 - 1 0 - 2 y i + 0 , 0 3 5 y j y 2 - 0 , O i l y 2 + 0 , 1 0 5 -10"3y f ; z 3 = 5 ,0 2 7 + 0 , 0 6 2 y 3 - 0 , 3 2 6 y 4 + 0 , 0 3 7 y 3y 4 + 0 , 2 0 5 - l ( T 2y 3 + 0 , 7 8 8 - 1 0 ' 2у 2 ; z 4 = 1 1 ,2 6 3 - 0 , 15 2 y 4 + 0 , 14 6 y , + 0 , 2 1 5 y , y 4 - 0 , 12 l y 2 - 0 , 0 6 9 y ? ;
y, = 1 3 2 , 6 5 8 - 1 0 4 , 6 5 0 x N - 1 5 3 , 1 0 - A 2 + 2 3 2 , 5 0 x 2 ;
у2 = -22,685 + 36619,20 • Лг -1 1 2 ,690xN - 93620,0 ■ArxN + 617,90x2 ; y 3 = 3 5 , 4 0 1 - 2 8 2 , 1 4 0 X N - 5 4 , 7 0 X C + 3 2 1 , 1 0 X N X C + 5 8 7 , 5 0 X 2 + 2 4 , 3 3 0 X 2 ;
y 4 = 1 4 0 , 1 5 5 + 3 6 , 5 3 0 x c - 4 3 , 0 1 O x 2 - 1 7 2 , 7 0 - A 2 ;
xN , xc - объемная доля в растворенном газе соответственно азота и метана.
При известном газовом факторе давление насыщения газонасы щенной нефти при температуре 293 К
P s * » = P „ + - ^ - |
( 8 . 5 2 ) |
К Р
С увеличением температуры величина давления насыщения растет
P S = P S293 + M T - 2 9 3 ) , |
(8-53) |
где Ks - температурная поправка давления насыщения , для нефтей (МПа/град)
К5 = 10^- 26,2 + ^ з - ( 1 9 ,4 - Д н -1 2 ) |
(8.54) |
Рзт |
|
Для газовых конденсатов величина давления насыщения рас считывается по формулам (3.23) , (3.24).
Плотность нестабильных жидкостей может быть рассчитана по формуле
+ Рну Г
(8.55)
В к
где р - плотность стабильной жидкости; Вк — объемный коэффициент, показывающий во сколько раз
увеличивается объем стабильной жидкости при растворении газа в ней. Для нефтей
Вк =(1 + 19,4-д ;-1В • Др49 • Г)” , |
(8,56) |
|
а для нестабильного газового конденсата |
|
|
— **” ®В *Р |
’ *Р » |
(8.57) |
|
293
где а , , в , , с, - расчетные коэффициенты , полученные при обработ ке номограммы Стендинга
а в =1(Г2 -[7 4 ,8 -7 ,99-(Д г -1 )];
ав = 10"4 • [8,65 + 2 ,6 2 '(Д г —1)]; |
(8.58) |
св = 10'3 -[2,37+ 0,35-(Аг -1 )];
Аг - относительная плотность газа по воздуху, Аг = р ну /1,205. Величина Вк вычисляется по формуле (8.56) при 293К. Кинематическая вязкость газонасыщенной нефти У8(мПа-с) при
температуре Т может быть рассчитана через вязкость дегазирован ной нефти при той же температуре v (мПа-с) по формуле
v s =v-B^ |
(8.59) |
где Z —расчетный коэффициент
Z = -0,375 ■(in v)1’358 • (In Вк)"°’33 • Г0’418, |
(8.60) |
где v —относительная (по воде) кинематическая вязкость дегазиро ванной нефти при расчетной температуре.
Динамическая вязкость газовой фазы составляет (мПа с)
Цг =0,005 l[l + рс1 (1,1 —0 .25Рст)] • [о, 037 + Tnp (l —0,104Тпр)] -
(8.61)
. 1 + ___ ---------
30(Т„р -1 )
где Кц - расчетный коэффициент, равный
к м= 0,0051 [1 + р„ (1,1 - 0 ,2 5 р „)]; рст - плотность газа при стандартных условиях.
§ 8.5. Особенности технологического расчета изотермической перекачки нестабильных жидкостей
При транспортировке нестабильных жидкостей в каждой точке трубопровода должно поддерживаться давление, пре-
294
вышающие давление насыщения Ps при температуре перекачки. Поэтому уравнение баланса напоров в данном случае имеет вид
|
f |
P |
^ |
Nr |
+ Н 2S + п • H g — 1,02 ■i ■L + Az + N 3 • |
Нкп + ----- |
,(8.62) |
|
„ 'б |
PsP -g |
j |
|
VrsP |
p ° |
где pSp —плотность нестабильной жидкости при расчетной температуре; Ней, is “ соответственно расчетный напор одной станции и гидрав лический уклон при перекачке нестабильной жидкости.
Отсюда необходимое расчетное число насосных станций
п = 1,028, ■L+ AZ + N .- Q H .- H » ) |
H - N ,- H „ |
(8.63) |
н CTS |
н crS |
|
Нетрудно видеть, что если опустить индекс «s», то формула (8.63) ничем не отличается от аналогичной формулы (5.32) для случая пе рекачки дегазированных нефтей. Поэтому все рассуждения относи тельно округления расчетного числа насосных станций, приведен ные в главе 5, справедливы и для случая перекачки нестабильных жидкостей. Однако необходимо иметь в виду, что из-за невозмож ности закачки нестабильных жидкостей в атмосферные резервуары, число эксплуатационных участков N 3 не может быть больше одного.
При давлении, превышающем давление насыщения, нестабиль ные жидкости ведут себя как капельные. Поэтому расчет величены is и HctS может быть выполнена по формулам (5.18) и (5.30) соответ ственно, но при подстановке в них параметров нестабильной жид кости: расхода Qs и кинематической вязкости vs.
Расчетный часовой расход нестабильной жидкости Qs находится по формуле
Q s = |
o s |
(8.64) |
|
24-Np -ps ’ |
|||
|
|||
где Gs - массовый годовой расход нестабильной жидкости |
|
||
G s = G r |
(1 + рг Г /р 2И); |
(8.65) |
р293 ~ плотность стабильной жидкости при 293 К.
Поскольку при транспортировке нестабильных жидкостей под пор основных насосов увеличен на Ps/Pspg ,то рабочее давление на выходе головной насосной станции
Pps — PsP ‘ g ' ( m MHS ' ^ M H S + ^ 2s ) + ^ s • |
( 8 .6 6 ) |
295
Так как величина Ps может достигать 2...3 МПа то число основ ных насосов mMHS может быть меньше 3-х.
§ 8.6. Гидравлический расчет трубопроводов, транспортирующих нефтяные эмульсии
Из подавляющего большинства скважин из влекается обводненная нефть. При этом образуются высокодиспер сные стойкие эмульсии, которые ведут себя как однородные жидко сти, и грубодисперсные неустойчивые нефтяные эмульсии.
Гидравлический расчет трубопроводов для перекачки высоко дисперсных эмульсий выполняются так же, как и для однофазных жидкостей, но с использованием физических свойств стойких не фтяных эмульсий.
Для грубодисперсных неустойчивых нефтяных эмульсий гидрав лический расчет трубопроводов ведется с учетом эффекта гашения тур булентных пульсаций дисперсионной среды каплями дисперсной фазы.
Исходными данными для гидравлического расчета трубопрово дов, транспортирующих нефтяные эмульсии, являются его внут ренний диаметр d и длина L, разность геодезических отметок конца и начала трубопровода Дг, плотности нефти р и пластовой воды рв, динамическая вязкость нефти рж и пластовой воды |iB,межфазное натяжение онв, объемные расходы нефти QHи пластовой воды QB.
Последовательность гидравлического расчета трубопроводов для перекачки нефтяных эмульсий следующая.
Объемная доля воды в эмульсии
Рв = |
Q , |
(8.67) |
|
QH +Q„
Критическое расходное содержание воды в эмульсии, при кото ром происходит инверсия фаз
|
(Нж/УвГ5 |
р..ч> |
(8.68) |
|
■ -ф ж Л ь )0,5 |
Если Рв< РВ1ф , то тип эмульсии - «вода в нефти», в которой нефть —дисперсионная среда (в дальнейшем индекс «с»), вода —дис персная фаза (индекс «ф») и объемная доля дисперсной фазы рф=рв.
Если же Рв>Рвкр, то тип эмульсии — «нефть в воде», в которой дисперсионная среда —вода, дисперсная фаза —нефть, а объемная
296
доля дисперсной фазы рф=1-рв. Плотность эмульсии
Р э = Р с - 0 - Р ф ) + Р ф - Р ф - |
(8.69) |
Динамическая вязкость эмульсии
Р з = Р с -0 -Р ф Г 2’5 |
(8.70) |
Средняя скорость эмульсии
m , 4 - w - + Q .) |
(8.71) |
’Ttd2
Число Вебера
We = |
(8.72) |
Первое приближение величины среднего объемно-поверхност ного диаметра капель эмульсии (без учета эффекта гашения турбу лентности)
d, = l,4 -d - We0,6 |
(8.73) |
Дополнительное напряжение сдвига эмульсии при рф>0,524
то =(0,195-рф- 0 ,1 0 2 ) - ^ - . |
(8.74) |
“ l |
|
Если рф < 0,524, то t0=0. |
|
Параметр Ильюшина |
|
И = - ^ Й - . |
(8.75) |
Число Рейнольдса при течении эмульсии |
|
Re. = со -d-p, |
(8.76) |
i . иЛ |
|
^ | 1 + б |
|
Наличие или отсутствие эффекта гашения турбулентности опре деляется по параметру Медведева
297
Md= Re3- W e1,2 • рф /p 3. |
(8.77) |
Если Md < 0,46, то эмульсия является плотной и эффект гаше ния турбулентности отсутствует (у, = 0), если же Md > 0,46, то эмуль сия неустойчива и эффект гашения турбулентности имеет место (у,=1). В последнем случае необходимо вычислить уточненное значение объемно-поверхностного диаметра капель эмульсии
d ,.= |
0.4 ’ |
(8.78) |
[ ( l —Рф ) • ( l —0,863 • Рф • M°,ls)] |
|
|
где М, - расчетный параметр |
|
|
м ,= |
м 2 Х |
(8.79) |
d-p c -CT4 ’ |
||
|
с ~нв |
|
При найденной величине d,. уточняются значения т0, И и Re3 по формулам (8.74)...(8.76).
Коэффициент гидравлического сопротивления при течении
эмульсии |
|
|
|
|
64 |
|
при |
Re3 < 2320 |
|
Re, |
|
|||
|
|
(8.80) |
||
X = |
|
|
||
0,3164 |
|
при 2320 < Re3 < 105 |
||
(1 + 1.125-у,-рф)-Ке°-3 |
||||
|
|
|||
Потери давления в трубопроводе при движении эмульсии |
||||
^ |
L-C02 |
|
^ |
|
АР = рэ |
X -------- + g • Az |
(8.81) |
||
|
3 2d |
6 |
J |
§ 8.7. Гидравлический расчет трубопроводов для перекачки газоэмульсионных смесей
Газоэмульсионной называется трехфазная смесь нефти, газа и воды. Она образуется на участках промысловой сети сбора между устьями скважин и сепарационными установками.
298