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

книги / Приборы и средства учета природного газа и конденсата

..pdf
Скачиваний:
1
Добавлен:
12.11.2023
Размер:
9.96 Mб
Скачать

числом Рейнольдса, представляющим собой отношение сил инерции к силам вязкости потока, т. е.

 

 

Re =

HcpDp / |i,

 

 

(73)

где «Ср— средняя скорость

потока,

м/с; р — плотность

газа

в рабочих условиях, кг/м3; D —диаметр трубопровода, м; р —

динамическая вязкость среды в рабочих условиях, Па-с.

среды

В

формуле (73)

единица динамической

вязкости

(газа)

р приводится

в СИ. Так как

в ряде

источников

дина­

мическая вязкость дается в МКГСС, то для перевода единиц вязкости из одной системы единиц в другую следует пользо­ ваться соотношением (66).

Формулу (73) можно переписать, если ввести в нее кине­

матическую вязкость газа v= p/p:

 

Re = «СрD / v.

(74)

Если среднюю скорость газа выразить через объемный рас­ ход в рабочих условиях Qaб, а площадь трубопровода S через л£)2/4, то формула (74) будет иметь вид

 

R e = l,273Q0e/ (Dv).

(75)

Учитывая, что v=p/p, получаем

 

л „

1273<2обР

!273QHPB

 

К е =

Dp

= ----- Dp-------’

(/Ь )

где Qos и QH— объемный расход газа в рабочих и нормальных

условиях, м3/с;

р и рн— плотность газа

в рабочих и нормаль­

ных условиях,

кг/м3; D — внутренний

диаметр трубопровода,

мм; р — динамическая вязкость газа, Па-с.

Если объемный расход газа в рабочих условиях Qoe выра­ зить в кубических метрах в час, динамическую вязкость газа —

в килограмм-силах — секундах на квадратный метр,

а внутрен­

ний диаметр трубопровода D — в миллиметрах, то

выражение

(76) для определения числа Рейнольдса принимает вид

Re = 0,0361 Q06p / (Dp).

(77)

При известном расходе газа QH и плотности газа рн в нор­ мальных условиях число Рейнольдса определяется формулой

Re =

0,0361 QHp„/ (£>|х),

(78)

где QB— расход газа

в нормальных условиях, м3/ч;

р« — плот­

ность газа в нормальных условиях, кг/м3; D — внутренний диа­

метр трубопровода, мм; р — динамическая вязкость

газа в ра­

бочих условиях, кгс-с/м2.

 

Если динамическую вязкость газа р в рабочих условиях вы­

разить в паскаль-секундах, формула (77) принимает вид

 

Re = 0,354 Q„p„/ (Dp).

(79)

В качестве примера рассчитаем число Рейнольдса при следующих дан-

31

НыХ. Расход газа в нормальных

условиях

= 100 000

м3/ч, плотность

газа

в

нормальных

условиях

рн =

0,7172

кг/м3,

диаметр

трубопровода

D =

=

400 мм,

динамическая

вязкость

в

рабочих

условиях

pi = 1,13 X

Х10~* кгс-с/м2.

 

 

 

 

 

 

(78)

 

 

 

 

 

Число Рейнольдса определим по формуле

 

 

 

 

 

Re =

0,0361

Q H P H

=

0,0361

100000-0,7172

5,728 -10s.

 

 

£)ц

400-1,13-10-*

 

 

Проверим

значение Re по

формуле

(79)

при

выражении

динамической

вязкости р, в паскаль-секундах, связанной соотношением 1,13-10—6 кгс-с/м2х

X9,806 =

1,108-10-* Па-с.

 

 

Число Рейнольдса, определяемое по формуле (79), равно

 

Re =

0,354

QHPH

100 000-0,7172

5,728-10е.

Оц

°.354 400-1,108-10 -5

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

а = Qp/ QHK-

(80)

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

Поскольку коэффициент расхода зависит от геометрии су­ жающего устройства, то в практических расчетах под реаль­ ным расходом Qp понимают расход через сужающее устрой­ ство, имеющее заданную геометрию. При этом коэффициент расхода теряет свой первоначальный смысл, определяемый фор­ мулой (80), и становится коэффициентом пропорциональности между расходами Q06 и QH, определяемым одной из формул (12), (17) —(19), и правой частью этих формул. В дальнейшем под коэффициентом расхода а будем понимать коэффициент пропорциональности в указанных формулах расхода, который, в частности, для сопел, может быть больше 1 (рис. 7).

Коэффициент расхода а для диафрагм может изменяться от

ат т=0,6

при т = 0,05 до

а т ах = 0>8 при т = 0,7 (рис. 7, а), соп­

л а — от

amin = 0,98 при

т = 0,05

до ат ах=1,2 при т=0,65

(рис. 7, б).

 

(пропорциональности) в ос­

Значения коэффициента расхода

новном определяются экспериментальным путем. Коэффициент расхода а увеличивается с увеличением модуля т и уменьша­ ется с его уменьшением.

32

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

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

Так как вязкость газа изменяется несущественно, то основ­ ное влияние в формуле (73) на число Рейнольдса оказывает скорость потока г>ср или его объемный расход Q0 e, а также плотность газа р. Зависимость коэффициента расхода а от ■числа Рейнольдса при различных значениях модуля т приве­ дена на рис. 7.

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

щее число Рейнольдса Remin

было больше

граничного

числа

.Рейнольдса Rerp.‘

 

 

 

 

Remin == QmlnPmln /

(^flmax)

^ RCrp,

(81)

где .Qmin— минимальный измеряемый

расход

газа через

сужа-

3 Зак. 1626

33

ющее устройство; pmin— минимально возможная плотность про­ ходящего газа; ртах — максимально возможная динамическая вязкость газа.

Как видно из рис. 7, а, минимально допустимое граничное число Рейнольдса для диафрагм составляет 30000 при /п=0,05 и 1000 000 при т = 0,7. Граничное число Рейнольдса для сопел практически не зависит от модуля т и составляет 200 000 (рис. 7, б).

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

Сказанное справедливо только для трубопроводов с глад­ кой внутренней поверхностью, а для диафрагм необходимо еще и условие остроты входной кромки отверстия. При шерохова­ тости трубопровода и недостаточной остроте входной кромки диафрагмы должны быть введены поправки к значениям а. Коэффициент расхода при гладкой внутренней поверхности тру­ бопровода и острой входной кромке, постоянный в области Rerp<Re<Remax, называется исходным коэффициентом рас­ хода аи, который является функцией модуля сужающего уст­ ройства т или произведение та для различных типов сужаю­ щих устройств. Зависимости исходного коэффициента расхода сси сужающих устройств от модуля т и произведения та при­ ведены на рис. 8.

Значения максимального числа Рейнольдса Remax определя­ ются формулой

Remax = QmaxPmax / (Пр,щ1п)>

(82)

где Qmax— максимальный измеряемый объемный

расход газа,

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

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

<J = OHKHKIUI

(83)

где а — фактический коэффициент расхода сужающего устрой­ ства с учетом влияния неостроты входной кромки и шерохова­ тости трубопровода; К:1 и Кш— поправочные множители на неостроту входной кромки диафрагм и на шероховатость трубо­ провода.

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

34

Рис. 8. Зависимость исходного коэффициента расхода от m и та.

а, в —диафрагма; 6, г — сопло и сопло Вентури.

трубопровода меньше 300 мм. При диаметре трубопровода больше 300 мм коэффициент шероховатости можно принять равным 1, при диаметре трубопровода 100 мм К т = 1,014.

Коэффициент расхода диафрагмы увеличивается и при не­ достаточной остроте входной кромки. У диафрагмы с диамет­ ром отверстия меньше 150 мм кромка не должна иметь замет­ ного невооруженным глазом притупления. Такие условия соот-

3*

35

гп=6,64 +

ветствуют радиусу притупления кромок диафрагмы -т-0,07 мм.

Если входная кромка не отвечает этим требованиям, необ­ ходимо вводить поправочный коэффициент К& значение кото­ рого можно определять по формуле

/( „ = 1 + 0 ,4 /4 ,

(84)

где d — диаметр отверстия диафрагмы, мм.

Расчеты показывают, что влияние неостроты входной кром­ ки сказывается на коэффициенте расхода диафрагмы в трубо­ проводе диаметром до 300 мм. Для трубопровода диаметром свыше 300 мм влияние притупления кромок весьма мало и с достаточной для практических расчетов точностью коэффи­ циент Ки можно принять равным 1. Изложенное показывает, что для трубопроводов диаметром более 300 мм коэффициент расхода диафрагмы можно принять равным исходному коэф­ фициенту расхода аи, т. е. а = а и.

Для выполнения весьма точных расчетов значения коэффи­ циентов Кн и Кш даны в таблицах, приведенных в Правилах{11].

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

определен по одной из нижеприведенных формул

в

диапа­

зоне чисел Рейнольдса Remin< R e< 1 0 8.

 

 

 

Для углового способа отбора давлений

 

 

«у = ---- = L =

Г 0,5959

+ 0,0312 т 1-05 — 0,1840 т< +

У

1 — m2

L

/

10» \ 0.75-]

 

 

 

 

 

 

(85)

 

+ 0,0029 т 1-25 (-RS-J

J .

 

Для фланцевого отбора давлений

 

 

 

аФ= ■

[ 0,5959 +

0,0312 m1-05 — 0,1840 т* +

+ 0,0029 m1.25

106 \ 0,75

2,286 т г

0,8562 m'-5

1

(86)

(-£■ )

 

 

 

D

J*

 

 

 

 

 

 

 

Сравнивая формулы (85) и (86), можно формулу (86) для определения коэффициента расхода при фланцевом отборе дав­ лений выразить через коэффициент расхода при угловом спо­

собе отбора давлений,

определяемый по формуле (85), т. е.

Иф

,

2,286

т 7

0,8562 т 1'5

(87)

ССу Т

£ ( 1 _

т 2)

D

 

 

 

 

 

где ссу и ссф — коэффициенты

расхода диафрагмы

при угловом

и фланцевом отборах

давлений

соответственно;

т — модуль

диафрагмы; D — внутренний диаметр измерительного трубопро­

вода, мм; Re — фактическое

число Рейнольдса,

соответствую­

щее действительному расходу газа в измеряемой точке; Remin—

36

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

При использовании формул (85), (86) для определения коэффициента расхода а минимальное число Рейнольдса Remin выбирается в зависимости от модуля т:

0,05 <

т <

0,20 — Rem,„ =

7 -10<;

(88)

0,20 <

т <

0,59 — Remi„ =

6 ■105;

0,59

m < 0,64 — Remm = 8 • 10s.

 

Для расчета коэффициента расхода

а по формулам

(85),

(86) значение модуля диафрагмы т предварительно опреде­ ляют по номограммам, приведенным в приложении 32 Правил (1 1 ], в зависимости от вычисленного по формуле (26) коэффи­ циента С и заданного предельного перепада давления ДрП на диафрагме. Число Рейнольдса определяют по одной из фор­ мул (74)—(79).

При известных значениях модуля т и диаметра D трубо­ провода значение диаметра отверстия диафрагмы d может быть определено по нижеприведенной формуле, полученной из фор­

мулы (2 2 ):

 

d = Dim.

(89)

В качестве примера рассчитаем значения коэффициентов расхода диа­ фрагм ау и аф по формулам (85) и (87) при одинаковых модулях т = 0,5, минимальных числах Рейнольдса Remi„ = 105 и внутренних диаметрах изме­ рительных трубопроводов D = 720 мм для углового и фланцевого отборов давлений.

Для углового способа отбора давлений

 

а у = —

■ 1

Г[0,5959 + 0,0312 ■0,51-05 — 0,1840 • 0,54 +

 

У1 _ о,52

[

 

 

 

 

+ 0.0029-0,б' 2», /-1 ^1 ')°'751 =

(0,5959 +

0,0312-0,4830 —

 

I 1°5 }

\

У1 -0 ,2 5

 

 

 

-0,1840-0,0625 +

0,0029-0,42-10°-75) = 1,1547(0,5959 +

0,0150— 0,0115 +

 

+ 0,0012-5,62) =

1,1547-0,6061 = 0,6998.

 

Для фланцевого отбора давлений

 

 

 

а* =

ау •+ 2,286 от2

0,8562 т 1-»

= 0,6998 +

2,286- 0,52

 

D(1 — /и2)

 

D

 

720(1 — 0,5=)

0,8562-0,51-® =

0,6998 4

2,286-0,25

0,8562-0,3536

= 0,6998 +

720

 

720-0,75

720

 

 

+

0,5715

— З . 3-028- =

0,6998 +

0,0011 — 0,0004 =

0,7005.

 

540

720

 

 

 

 

 

Как видно из полученных результатов, коэффициент расхода при флан­ цевом отборе давлений несколько больше (на 0,035%), чем коэффициент расхода при угловом способе отбора давлений.

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

37

ьр/р,

btfp,

Рис. 9. Номограммы для определения коэффициента изменения плотности газа е в зависи­ мости от Др/pi, т и у.

а — большие значения Др/pf, 6 — малые значения Др/Pi.

устройство и является функцией относительного перапада дав­

лений Лр/рь модуля т и

показателя

адиабаты

природного

газа к.

е можно

найти по

приведенным

Поправочный множитель

на рис. 9 номограммам и известным значениям Др/pi, т и я следующим образом. Откладывают на верхней горизонтальной

оси значение Ap/pi

и из этой точки опускают перпендикуляр

к соответствующей

прямой показателя адиабаты я измеряемого

газа. Из точки пересечения значений AP/P I и я проводят гори­ зонтальную линию до пересечения с графиком значения модуля т выбранного сужающего устройства (сопла, диафрагмы). Из полученной точки опускают перпендикуляр на нижнюю горизон­ тальную ось е, пересечение с которой дает искомое значение е. При расчете сужающего устройства определяют минимальные и максимальные значения emtn и ет ах, соответствующие Дртш/ / Ридах и APmax/ApiminГрадуировка расходомера производится при среднем значении еср, вычисленном для наиболее вероят­ ного расхода, равного 2/з от максимального.

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

Исм = (1/ 100) {ацц гЪ*” (90)

38

где fli, Й2. •••> —содержание отдельных компонентов в смеси, %; «1, «2, .... и,-— показатели адиабаты для компонентов смеси.

Значения показателей адиабаты основных компонентов при­ родного газа и воздуха приведены в табл. 2 .

Значение поправочного множителя на расширение е может

быть определено и по приведенным ниже формулам:

 

е =

1 -

(0,41 + 0,35т2)Ар / (put);

 

(91)

и = 1,29 +

7,17-10-6 [2575 + (346,23 — Г)2]ри

(92)

где Ар — перепад давления на

сужающем устройстве,

МПа;

Pi — абсолютное давление газа

перед сужающим

устройством,

МПа; и — показатель

адиабаты

природного газа;

m — модуль

сужающего устройства; Т — абсолютная температура газа, К. При практических расчетах перепад давления Ар на диа­ фрагме обычно выражается в килограмм-силах на квадратный

метр, а давление — в килограмм-силах на

квадратный

санти­

метр.

(91) имеет вид

 

 

В этом случае формула

 

 

е = 1 — (0,41 +

0,35 т 2)Ар-10~4 /

(рце),

(93)

где Ар — перепад давления на диафрагме,

кгс/м2; рi — давле­

ние газа перед диафрагмой, кгс/см2.

 

 

Если в формуле (92) для определения показателя адиабаты х давление газа pi также выразить в килограмм-силах на квад­ ратный сантиметр, то ее можно переписать в следующем виде:

х = 1,29+0,704-10-6[2575+(346,23—Г)2^!,

(94)

где pi —давление газа перед диафрагмой, кгс/см2.

Для обеспечения достаточной точности измерения расхода газа значение в не должно существенно изменяться в рабочем диапазоне расхода — от 0,3 Qmax до Qmax-

Диапазон изменения Ае (%) определяется выражением

Ае =

(вшах — 8mln) ■100 / вшах.

Допустимое значение Ае не должно превышать 1%, при ко­

тором систематическая

погрешность измерения расхода газа

в диапазоне от 0,3 Qmax до Qmax не более 0,5%.

Так как поправочный множитель е в основном зависит от отношения Ap/pi и ш, то максимальный перепад давления Арщах на сужающем устройстве, соответствующий Qmax, не должен превышать значения, определяемого по номограммам, показанным на рис. 9. Следует отметить, что влияние перепада давления Ар на систематическую погрешность измерения рас­ хода газа особенно сильно сказывается при небольших давле­ ниях (до 20 кгс/см2).

С тем чтобы погрешность измерения расхода газа от изме­ нения поправочного множителя е не превышала 1 % в рабочем диапазоне изменения расхода, целесообразно при рабочих дав­

лениях в трубопроводе до 4 кгс/см2 иметь максимальный пере­

пад

давления Дртах не более 0,15 кгс/см2, при давлениях до

12

кгс/см2 — 0,4 кгс/см2, а при давлениях до 30 кгс/см2

1,1

кгс/см2.

При давлениях свыше 30 кгс/см2 и максимальном перепаде давления не более 0,63 кгс/см2 изменение поправочного мно­ жителя е практически не влияет на изменение систематической погрешности при измерении расхода газа.

В тех случаях, когда необходима высокая точность измере­ ния расхода газа, необходимо автоматическое введение попра­ вок в показания расходомера на изменение поправочного мно­ жителя по формулам (91) или (93) с учетом текущих значений Ар и pi и расчетных формул (92) или (94) для определения показателя адиабаты х.

1.4. ПРИМЕРЫ РАСЧЕТА ОСНОВНЫХ ПАРАМЕТРОВ ПРИРОДНОГО ГАЗА

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

Расчет плотности газа при нормальных условиях. Определим плотность газа при нормальных условиях при следующих данных по составу и объем­ ному содержанию его компонентов, %: метан — 93, этан — 4, пропан— 1, бутан — 0,6, пентан — 0,2, углекислый газ— 0,2, азот— 1.

По табл. 2 определим плотность входящих в смесь газовых компонентов

при нормальных условиях, кг/ма: метан — 0,6681,

этан— 1,2600, пропан —

1,8659,

бутан — 2,4947, пентан — 3,1633, углекислый газ— 1,9767, азот —

1,1889.

 

 

 

 

Плотность газовой смеси определяем по формуле (29):

 

+

рн. с„ =

(1/100) (93-0,6681 +4-1,2600 +

1-1,8659 +

кг/м3.

0,6-2,4947 +

0,2-3,1633 + 0,2-1,9767 + 1 • 1,1889) = 0,735

Расчет критических и приведенных температур

и давлений.

Определим

приведенные температуру и давление смеси природного газа следующего со­

става,

%:

метан — 92,

этан — 4,

пропан — 2, пентан — 1;

при

температуре

Т = 293 К

и избыточном давлении

 

рИзб = 2,942 МПа. Барометрическое да­

вление р е — 0,101325 МПа.

 

 

 

 

 

 

 

 

 

Из табл. 2 возьмем критические температуры и давления компонентов

природного

газа: для метана 7’к р

i =

190,66 К, ркр i = 4,889 МПа; для

этана

7’к р г =

305,46 К,

рнр 2 = 5,145

МПа;

для

пропана

Тир з = 369,90 К,

Р к р з =

4,483 МПа; для бутана

7’кр4 = 425,05

К, ркр4 =

4,003

МПа;

для

пентана 7’кр 5 =

469,50 К, pKP s =

3,554 МПа.

 

 

 

 

 

По этим критическим температурам и давлениям вычисляем критические

температуры и давления смеси по формулам (36) и

(37):

 

 

 

 

 

7’к р .

см =

(1/Ю0) (92-190,66 +

4-305,46 +

2-369,90 +

 

 

 

 

 

+ 1 • 425,05 +

1 - 469,50) =

204 К;

 

 

 

 

 

Ркр. см =

(1/100) (92-4,889 +

4-5,145 +

2-4,483 +

 

 

 

 

 

+

1 -4,003 +

1-3,554) = 4,87 МПа.

 

 

 

Приведенные температуры и давления выбранной смеси газов для задан­

ного избыточного давления смеси

р„зв =

2,942 МПа и

температуры

газа

Т — 293 К

определяем по формулам

(38)

и (39):

 

 

 

 

 

 

 

ГЛр. см = 7 7 7’кр. ои =

293/204 =

1,44;,

 

 

 

 

 

Рп р. см = р / Р кр . см = 3,043/4*87 =

0,625,

 

 

40