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

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

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

где

р = ризб + Ре = 2,942 + 0,101 = 3,043

МПа — абсолютное

давление

газа.

 

 

 

 

 

 

 

 

В качестве второго примера рассчитаем критическое давление и темпера­

туру

газовой смеси

с плотностью

при

нормальных

условиях

ря. см =

= 0,7763

кг/м* для газовой

смеси,

не содержащей азот

и углекислый газ,

а также

содержащей

азот с

объемной концентрацией 1% и углекислый газ

с концентрацией 0,5%.

 

 

 

 

 

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

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

(40) и

(41):

 

 

 

 

 

Ркр =

48,48— 1,8 ря. см = 48,48— 1,8-0,7763 =

 

 

 

= 48,48 — 1,39734 = 47,08 кгс/см2 = 4,617

МПа;

 

 

 

Т„р = 87,5 + 155,24 ря. см = 87,5 + 155,24 0,7763 =

 

 

 

 

=

87,5 + 120,5 =

208,0 К.

 

 

Для

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

держанием азота 1% и углекислого газа 0,5% предварительно определяем

концентрацию

азота и углекислого газа в долях от единицы:

 

 

 

 

 

N2 = l % : 100 =

0,01

и

С02 =

0,5% : 100 =

0,005.

 

 

 

газа

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

определяем по формулам (39)

и

(40):

 

 

 

 

 

 

 

 

 

 

Ркр = 48,51 — 1,808 Рн. с« — 11,826N2 + 30, !68COj =

 

 

 

 

 

=

48,51 — 1.808- 0,7763 — ll ,826-0,01 + 30,168-0,005 =

 

 

 

 

 

 

 

 

=47,14 кгс/см2=4,623 МПа;

 

 

 

 

 

 

 

7VP =

87,5 + 155.24ря. см — 148,35N2 — 88,25СОг =

 

 

 

 

 

= 87,5 + 155,24-0,7763 — 148,35-0,01 -88,25-0,005 = 206,1 К.

 

 

 

 

Расчет коэффициента сжимаемости природного газа. Рассчитаем коэффи­

циент сжимаемости

природного газа по формулам (46)

и

(65) при

следую­

щих

значениях

параметров

газа:

плотность

при нормальных

условиях

рн =

0,76 кг/м3,

избыточное

давление

ризе

=

0,78 МПа,

температура

t =

=

—9,5 °С, барометрическое

давление

ре =

0,1

МПа.

Концентрация

 

азота

N i

составляет

1%, углекислого газа

СОз — 0,5%.

 

 

 

 

 

 

ем

 

Для определения коэффициента сжимаемости предварительно определя­

по формуле

(15)

абсолютное давление

р,

а также

критические и

приве­

денные давления и температуры:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Р = Ризб +

ре =

0,78 +

0,1 =

0,88 МПа;

 

 

 

 

 

 

 

 

T =

i + 273,15 = —9,5 +

273,15 = 263,65 К.

 

 

 

 

Критические

давление

и

температуру

определяем

по

формулам

(48) и

(43):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ркр =

4,757 0,1773рв — 1,160N* +

2,958С02 =

 

 

 

 

 

=

4,757 0,1773

0, 7 6 1,160

0,01 +

2 ,9 5 8 - 0,005 =

 

 

 

 

 

 

 

 

=

4,62 МПа =

47,11

кгс/см2;

 

 

 

 

 

 

 

 

Гкр =

87,5 + 155,24рн — 148,35N2 — 88,25COj =

 

 

 

 

 

= 87,5+ 155,24-0,76 — 148,35-0,01 — 88,25.0,005 =203,56 К.

 

 

 

и

Приведенные

давление

и

температуру

рассчитываем по формулам

(39)

(38):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Р п р — : р /р к р =

0 ,88/ 4,62 =

0, 190;

 

 

 

 

 

 

 

 

 

г пр =

Г /Г к р =

263, 65/ 203,56 =

1, 295.

 

 

 

 

 

Для расчета коэффициента сжимаемости Z по формуле (46) находим

значения коэффициентов ре и Т е из выражений (48) и

(47):

 

 

 

Рс= 0,6714рпр=0,6714-0,190=0,128;

Г с= 0,718927^=0,71892-1,295 =0*931.

и

 

Так

как

значения

рс <

1,3,

а

Т с <

1.09,

значение

F

получаем

по

фор­

муле (54):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F = 1 — 75'10-5рс

{ро*-а[2 — e-20d.09-ic)] +

 

 

 

 

 

 

 

4-1756(1,09 — Гс)М1.69 — рс2) } =

1 — 0,00075 0,128 X

 

 

 

 

 

X { 0.1281,3 [2 — 2,7~го<,'0#-°>931>] +

1756(1,09 — 0,931)4 X

 

 

 

 

 

 

 

 

X (1,69—0,1282)}= 1,000

 

 

 

 

 

 

 

 

 

Значение 9 |

находим по формуле

(52):

 

 

 

 

 

 

 

 

 

 

 

 

0, =

(

1’61353^,Г<:1~~.2'?1323— +

3,30378) /0.9312 =

3,3566.

 

 

Значение 02 находим по формуле

(53):

 

 

 

 

 

 

 

 

 

 

 

 

 

02 =

(■4,57657!r^_+_26_j8 27_

__ I3)3185)j

/ (т-с01) =

 

 

 

 

= ^4,57657/ О^П2 4- 26,5827

_

13 3I85 j , { (0931.3,3566) =

6,65244.

 

 

Далее по формуле (51) получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

Bo ^ 4 5 O _ - e ,0 2 2

Q 2 _

5 0 p c F

)

/ 0 1 =

^ —450 — 3,3566-6,652442 , х

 

 

X 6,65244 — 50-0,128-1) /3.3566 =

 

67,87/3,3566 =

20,2211.

 

 

 

По

формуле

(50)

определяем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В 2 -V

В

0 +

V

W

 

20,2211+

+^

У207 .22112=

+ V5.01353 =

3,5129.

 

Значение В \

находим по формуле

(49):

 

 

 

 

 

 

 

 

 

 

В, =

(300 — 0 ,0 22) / (90.) =

(300 — 3,3566-6,652442) / (9-3,3566) =5,0135.

 

Зная

коэффициенты Во,

B t,

В2, ©г и 0 2, по формуле

(46)

определяем

коэффициент

сжимаемости

природного

газа

Z

в

рабочих

условиях:

 

 

7

_

(1 4 - 1,32-10-3/Гс3-25)2

_

 

_

 

 

 

(1 4- 1,32-10-3/0,9313-25)2

v

 

 

 

В ,/В 2 — В2 4 -0 2 /3

с

 

 

5,0135 / 3,513 — 3,513 4- 6,65244 / 3

 

 

 

 

 

 

 

X 0,128 =

0,9747.

 

 

 

 

 

 

 

 

(64)

Рассчитываем коэффициент

сжимаемости

1

по

приближенной

формуле

при

абсолютном

давлении

 

р =

0,88

 

МПа

и температуре

газа

Т =

= 263,6 К:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z =

1 — 1,57р Ц Т — 198) =

1 — 1,57• 0,88/ (263,6 — 198) =

0,979.

 

ния

Сравнивая значения Z , определенные по точной формуле (44),

и значе­

Z ,

полученные по приближенной

формуле (62), видим, что они отли­

чаются друг от друга

не более чем на

0,5%,

а

это вполне

допустимо

при

выполнении большинства практических расчетов.

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

веденного

давления рпр> и

восстанавливается перпендикуляр до

пересечения

с кривой,

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

заданной приведенной температуре

Т в р .

Из точки пересечения перпендикуляра с кривой приведенной темпера­

туры

проводится горизонтальная прямая до пересечения с осью ординат, на

которой отсчитывается

искомое значение коэффициента сжимаемости Z.

Если

заданное значение

приведенной температуры находится между двумя

42

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

мация этих значений.

 

 

 

 

 

 

 

 

 

 

 

при р а р =

3,41

В качестве примера найдем коэффициент сжимаемости Z

и Г„р =

1,71.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отложим на горизонтальной оси значение приведенного давления рир —

= 3,41 (рис. 4) и восстановим

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

чения с

кривой приведенной

температуры

Тпр =

1,71.

Так как

на рис. 4

имеются только кривые

Гпр =

 

1,7 и 1,8, то

значение Гпр =

1,71

находится

методом аппроксимации.

 

перпендикуляра,

соответствующего

р а р =

3,41,

Из

точки

пересечения

с кривой Т Пр — 1,71, проведем

горизонтальную прямую

до пересечения с вер­

тикальной осью Z, на которой и отсчитываем искомое значение Z

= 0,86.

Значение 2, определенное по графику, приведенному

на рис. 4, по исход­

ным данным

предыдущего

примера рпР = 0,190 и

Тпр =

1,295,

равно

0,98.

По данным расчета 2 =

0,979.

 

 

 

 

 

 

 

 

 

 

 

 

 

Расчет динамической вязкости природного газа. Определим динамиче­

скую вязкость природного

газа

следующего

состава,

%: метан — 92, этан —

4, пропан — 2,

бутан— 1,

пентан— 1;

при

температуре

20 °С,

избыточном

давлении газа

риаб = 3,0 МПа

и барометрическом давлении р в = 0*1 МПа.

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

20°С (293К).

Для метана

j.ii =

1,13-10—5 Па-с;

 

для

 

этана

р2 — 0,95Х

Х10-5 Па-с;

для пропана

р3 =

0,82■ 10~' Па-с;

для

бутана

р4 = 0.75Х

Х10-5 Па-с; для пентана

ps =

0,70-10~5 Па-с.

 

 

 

 

 

 

 

 

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

природного газа: для метана

 

Гкр i =

190,66 К,

Ркр i =

4,889

МПа;

для

этана Ткр г =

305,46 К,

ркр s =

5,145МПа;

для пропана

 

7’кр 3 =

369,90 К,

ркр з =

4,483

МПа; для

бутана

 

Гкр 4 = 425,05 К,

 

ркР 4 =

4,003

МПа;

для

пентана Гкр s= 469,50 К,

Ркр s= 3,554 МПа.

 

 

 

 

 

 

 

 

 

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

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

(35) и

(37):

 

 

 

 

 

 

Гкр. СМ— (1/100) (92 ■190,66 + 4• 305,46 +

2 ■369,90 +

 

 

 

 

 

+1 • 425,05 +

1 - 469,50)

= 204,0 К;

 

 

 

 

 

 

Ркр. ем = (1/100) (92-4,889+ 4-5,145+ 2-4,483 +

+1 -4,003 + 1 -3,554) = 4,71 МПа.

и

Приведенные

температуру

и

давление

 

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

(38)

(39):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Гпр. см =

(t +

273,15) /204,4 =

293,15 / 204,4 =

1,44;

 

 

 

 

 

рпр. см ” Р / Ркр =

3,1 /4,71 =

0,66,

 

 

 

где абсолютное значение

давления

р =

рИзв'+ ре =

3 + 0,1 = 3,1

МПа. По

полученным приведенным

 

температуре

и давлению смеси

газов

найдем из

рис. 5, б поправочный

множитель

Сц

при

Гпр. см =

1,44

и рпр. си =

0,66,

который равен

1,02.

 

 

смеси

Цсм.н

при

рабочей

температуре

 

=

Определяем

вязкость

 

t =

20°С и нормальном давлении по формуле

(69);

 

 

 

 

 

 

Пси. Н =

(10-5/ 100) (92-1,13 +

4-0,95 +

2-0,82 +

 

 

 

+ 1-0,75+ 1-0,70)

=

(1,04 + 0,04 +

0,02 +

0,01 + 0,01)-10-5 =

 

 

 

 

 

 

 

=

1Д1-10-5 Па-с.

 

 

 

 

 

 

 

Вычисляем вязкость

смеси р в рабочих условиях по. формуле

(68):

 

 

р =

Ср. Рем. н =

1,02-1,11

10-5 =

1,13.10-5 Па-С.

 

 

 

При необходимости' выражения динамической вязкости в килограмм-си­

лах— секундах

на квадратный

метр воспользуемся

выражением

(66):

 

 

 

1,13-10-5

П а-с/9,80665 =

1,15-10-»

кгс-с/м*.

 

 

43

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

условиях при содержании метана

92% (м- == 1,13- 10—в Па с) практически

равна вязкости метана (цметана =

1,13-10—® Па-с) при нормальном давле­

нии и рабочей температуре, в связи с чем в дальнейшем при практических расчетах в диапазоне температур от —30 до + 100вС и давлений от 0 до 10 МПа вязкость смеси природного газа можно принимать равной вязкости

метана при нормальном

давлении и рабочей

температуре, т. е. р, = цСм. н.

Расчет влажности газа. В качестве примера определим по рис. 6 абсо­

лютную влажность газа

при температуре его

точки росы + 5 °С и абсолют­

ном давлении

0.5 МПа.

Для

этого из точки

на оси абсцисс со значением

температуры

точки росы

+ 5 °С

восстановим

перпендикуляр до пересечения

с кривой, соответствующей давлению 0.5 МПа, и из точки пересечения про­ ведем горизонтальную прямую до пересечения с осью ординат. При этом искомое значение абсолютной влажности газа равно 1,4 г/ма, Из рис. 6 так­ же можно определить массу воды, выделяющуюся из каждого кубометра на­

сыщенного влагой газа

при снижении

температуры от t, до <2-

t\ = +10 вС

до

Так,

при снижении

температуры

стенок трубопровода от

U =

—10 °С абсолютная влажность, определенная по рис. 6,

изменяется

от

10 до 2 г/м3, т. е. из каждого 1 м3

газа выделится 8 г воды.

 

1.5. ПОГРЕШНОСТИ ИЗМЕРЕНИЯ РАСХОДА ГАЗА

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

Предельная погрешность измерения сухого газа по методу переменного перепада давления с использованием дифманомет- ров-расходомеров, у которых предельная погрешность б<з<4Р) за­ дана по шкале перепада давления, может быть определена выражением

«<?<*.) = У Ч + Ч + 8’е + (1/4) (8»н + Ъ\р + Ц + VT +pz). (95)

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

bQivrP) = у ъ \ + 5? + 8* е +VYM + (1/4) (8рн + 82р + 82г + Щ , (96)

где 6<з(4Р) и бQ (V iji) — предельные погрешности измерения рас­ хода сухого газа при использовании дифманометров-расходоме- оов со шкалой перепада давления и расхода соответственно, %;

8а, 8е,

8Др, ..., бz — предельные погрешности

измерения а, е,

Ap,'...t

Z и других величин, входящих в формулы (17)—(21),

.%; бке — предельная погрешность от влияния

изменения числа

Рейнольдса в рабочем диапазоне расходов,.%.

 

44

flредельйая погфШйбсть определения коэ ффщиентя расхода Диафрагмы ба рассчитывается по формуле

К - у ь 1 + Ча + Чо + Чт + Чп+ + S- + U

(97)

где баи— предельная погрешность определения исходного коэф­ фициента расхода диафрагмы, принимаемая равной 0,3%, т. е.

и=0,3%; б ad— предельная погрешность измерения внутрен­ него диаметра отверстия диафрагмы d, %; 6 »D — предельная погрешность измерения внутреннего диаметра измерительного трубопровода D, %; 6„ш — предельная погрешность влияния шероховатости внутренней поверхности трубопровода перед диафрагмой, %; бan— предельная погрешность, вызванная при­ туплением входной кромки диафрагмы, %; бах,— погрешность от влияния местных гидравлических сопротивлений перед диа­ фрагмой, %; бае— предельная погрешность от влияния экс­ центриситета между осью диафрагмы и осью измерительного трубопровода, %; бан — предельная погрешность от влияния уступа от стыковки труб внутри прямого участка трубопровода перед диафрагмой, %.

При внутреннем диаметре трубопровода перед сужающим устройством больше 300 мм значениями погрешностей бащ и оап можно пренебречь и принять бащ “ б ап= 0. При диаметрах трубопровода меньше 300 мм погрешности б1Ш и dan опреде­ ляются по методике, приведенной в Правилах [11].

Значения погрешностей 6ad и б«о, вызванных отклонением от расчетных значений диаметра сужающего устройства d и внутреннего диаметра трубопровода D для диафрагм, сопел и сопел Вентури, равны 0,025% при 0,4 и 0,05% при т < 0,4. Для сопел Вентури б«/=0,05% для всех допустимых значений т . Значения погрешности 6 « D = 0 , 1 5 % д л я всех типов сужаю­ щих устройств в диапазоне допустимого значения т. При числе Рейнольдса R e^lO 6 в рабочей зоне изменения расхода газа коэффициент коррекции на изменение числа Рейнольдса /Ске= 1, а погрешность от изменения числа Рейнольдса бяе=0.

Погрешность

от изменения числа

Рейнольдса

в рабочем

диапазоне расходов

 

 

 

 

 

бне

(1 — Кце) 6(1.

 

 

(98)

где Кяе — коэффициент

коррекции на

изменение

числа Рей­

нольдса в рабочем диапазоне расходов;

б^ — погрешность

определения динамической вязкости газа, 6(1= 5%.

 

Коэффициент Кце может быть определен выражением

/Сне =

[С + В(106/ Remin)0,7S] / +

5 ),

(99)

где

 

 

 

_____

С = (0,5959 + 0,0312т1'05 — 0,184 т 4) 1 т 2;

 

В =

0,0029 т 1-2s / У 1 — т 2;

 

 

45

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

Погрешность бав от влияния эксцентриситета между осью сужающего устройства и осью измерительного трубопровода перед сужающим устройством определяется следующим обра­ зом. Если эксцентриситет

е < 0 ,00057) /(0,1+2,3 т 2),

(100)

то погрешностью бае можно пренебречь, т. е. можно принять, что 6 ае = 0. Если значение эксцентриситета е находится в пре­ делах

0,00057)/ (0,1 + 2,3m2) <

0,0057)/ (0,1 + 2 ,3 т2),

(101)

то погрешность бае не превышает 0,3% и можно принять, что

*««*=0,3%.

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

Н = (h/D) • 100 < 0,3%,

(102)

то погрешностью 6 „н можно пренебречь и принять 6ан=0. Если относительная высота выступа Н между двумя участками со­ стыкованных труб больше 0,3%, но соответствует условиям

д < ° . 2 ( с т ш , <) » « < * > %

.(103)

то погрешность б»ы=0,2%. Здесь I — расстояние от переднего торца сужающего устройства до выступа, h — высота уступа.

Погрешность баь от влияния местных сопротивлений до и после диафрагмы определяется по табл. 5—10. Подставив зна­ чения погрешностей bad, 8ад, 8аш, §<,«, °ае, §«я, и baL в формулу (97), получим значение предельной погрешности 6 а, определения коэффициента расхода а сужающего устройства.

Погрешность определения коэффициента расширения е оп­ ределяется по формуле

S£ = [ . ( l - e ) / e ] | / ^ + S2p + 82+ S 0,

(104)

где е — коэффициент расширения; б*, бдр, 6р — предельные по­ грешности определения показателя адиабаты к, перепада дав­ ления Ар и давления р, %; 6о=2Др/р для т ^ 0 ,5 6 и 60=4Ар/р для 0,56^т^0,64 .

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

0,5(Qnp/ Q )5/4^O»

(195)

где Qnp — предельное значение расхода газа по шкале расхода

46

ДифманоМетрафйсХодомера; Q— текущее значение рйсхбдй М- за в поверяемой точке; S Y T p ~ класс точности дифманометра-

расходомера.

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

ности определен

по

шкале

перепада давления,

погрешность

измерения бдр определяется формулой

 

 

 

Ыр =

0,5 (Др„Р/Ар) S*p,

(106)

где

6 др — предельная

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

ния,

%; АрП — предельное

значение шкалы дифманометра по

шкале перепада

давления;

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

в поверяемой точке;

— класс точности дифманометра.

Погрешности

регистрирующих дифманометров

рассчитыва­

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

ЧР= 0.25 ^ L 5 д р ) 2 + 0,256ПК + 0,0012ДтдР (107)

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

Ь2у~Р = 0,25 ( ~ S y s p j + 0,256пп + 0,0012Ат2^ . <108)

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

&УТр = 0,25 l ^ ~ SvTpJ + 0,256пк .+ 0,0012Дх2^, (109)

где бпк и бпп— приведенные погрешности корневого и пропор­ ционального планиметров соответственно, %; Дтдр и Дт/дJ — абсолютные погрешности хода картограмм дифманометров с классом точности по перепаду давления и расходу соответ­ ственно, мин.

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

бР = 0,5 (рпр / p)Sp, (НО)

где 6р — предельная погрешность измерения давления газа, %; Рпр— предельное значение шкалы манометра; р — текущее зна­

чение измеряемого

давления; 5р — класс

точности манометра.

Для регистрирующих манометров

 

 

Sp2 = 0,25

( - ^

)2 + 0,2561 +

0,0012ДтР2,

(111)

47

Где бдп— погрешность планиметра,

%;

Атр — абсолютная по­

грешность хода диаграммы манометра, мин.

определяется

Погрешность отбора 6о

перепада

давления

по формуле

 

 

 

 

бо = 2Др/р

для

0,56.

(112)

Для показывающих термометров

 

 

 

бг 2 = 0,5[Nt /

(273,15 +

*)]S t,

(113)

где Ьт— предельная погрешность измерения температуры газа,

%;

Nt — диапазон шкалы

измерений термометра; St — класс

точности термометра, %.

 

 

 

 

Для регистрирующих термометров

 

 

 

6г2 = 0,25 (-273j 5^

t

St J + 0,256m..

+ 0,0012Ат*2,

(114)

где

Ат* — абсолютная

погрешность хода

диаграммы

термо­

метра, мин.

 

 

 

 

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

рассчитывается по формуле

 

6рн = 50Арн/ рн,

(115)

где Дрн — максимальная абсолютная

погрешность определения

рн, равная половине разряда последней значащей цифры в таб­

личном

значении рн;

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

при нормальных

условиях.

 

 

Погрешность определения коэффициента сжимаемости Z

природных газов в диапазоне температур от

—25 до 80 °С и

давлений от 0 до 8 МПа рассчитывается по формуле

бz =

У (1 - Z ) W

+ 16бг2 + 46^ + 0,046*, + 0,003б2Оа +

 

 

5 2 та6л'

(116)

где Z — коэффициент сжимаемости газа; 6N2 и бсо2— погреш­ ности определения молярной концентрации азота и углекислого газа в природном газе, %; 2 табл— погрешность определения табличных значений коэффициента сжимаемости, %. Для прак­ тических расчетов можно принять 6гтабл=0,25%.

Подставив полученные значения составляющих всех погреш­ ностей в формулу (95) или (96), получим значения предельных погрешностей измерения расхода газа бв(/д£> или 6<з(дР).

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

В то же время в практике учета газа в газовой промышлен­ ности производится не разовое измерение мгновенного расхода

48

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

обработки нескольких суточных диаграмм уменьшается в yt раз и определяется формулой

б = б; / У i + AI'CHCT,

(117)

где б — результирующая погрешность определения объема газа за несколько суток по данным обработки i суточных диаграмм, %; б,-— погрешность измерения и обработки одной суточной диаграммы, %; i — число обработанных диаграмм; А/Сист— си­ стематическая составляющая погрешности измерения и обра­ ботки одной суточной диаграммы, %.

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

Аг’сист = (1/3) бг.

(118)

Подставив значение А^ист из (118) в

(117), получаем

6 = б</уГ+ (1/3)бг.

(119)

В качестве примера рассчитаем по формуле (119) резуль­ тирующие погрешности определения объема газа бзо и бзб5 (за месяц и за год) по результатам обработки 30 и 365 суточных диаграмм соответственно, снятых с дифманометра-расходомера класса 1,5 (б*=1,5%);

бзо =

1,5 / У 30 +

(1/3) • 1,5 =

0,27% + 0,5% = 0,77%,;

бзо5 =

1.5/У365 +

(1/3) -1,5 =

0,08% + 0,5%,= 0,58%.

Из рассмотренного примера видно, что при обработке боль­ шого числа суточных диаграмм дифманометра-расходомера класса 1,5 результирующая погрешность определения объема газа за отчетный период— 1 мес— не превышает 0,77%, а за год — 0,58%; это подтверждается практическими данными экс­ плуатации пунктов учета газа.

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

4 Зак. 1626

Глава 2

СУЖАЮЩИЕ УСТРОЙСТВА И ИЗМЕРИТЕЛЬНЫЕ ТРУБОПРОВОДЫ

2.1. СПОСОБЫ ОТБОРА ДАВЛЕНИЙ ОТ СУЖАЮЩИХ УСТРОЙСТВ

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

В СССР распространены два способа отбора давлений — угловой и фланцевый. При угловом способе отбора давление отбирается непосредственно у диафрагмы с помощью угловых отверстий (рис. 10, б) или кольцевых камер (рис. 10, в). При фланцевом отборе (рис. 10, е ) давления p i и р2 отбираются через отверстия во фланцах, находящихся на равном расстоя­ нии l\ = h = 25 мм или на расстоянии l\ = h= Q k)D от соответ­ ствующей плоскости диафрагмы. В США применяется метод суженной струи (рис. 10, г), при котором давление р\ отбира­ ется через отверстия, расположенные в сечении АА (рис. 10, а), а давление р2— в сечении ББ. При этом давление р^ от­ бирается на расстоянии /i= (0,5-Ь2)£> от передней плоскости диафрагмы, а р2 — на расстоянии 12= (0,2—0,8 )D от задней плоскости диафрагмы.

Радиальный способ отбора (рис. 10, д) весьма близок ме­ тоду суженной струи. В этом случае давление Pi отбирается на расстоянии l\—D, а давление р2 — на расстоянии /2=0,5Z> от соответствующей плоскости диафрагмы.

Влияние способа отбора давления может быть учтено по

формуле

а = а„. УК, ,

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

На рис. 11 приведена зависимость изменения коэффициента от модуля т при переходе от углового отбора к другим способам. Как видно из рисунка, при модуле диафрагмы т ^ 0 , 5 перепады давления при угловом способе отбора практически равны перепадам при отборах по методу суженной струи, ра­ диальном и фланцевом, что позволяет говорить о равенстве

60