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

Исследование равновесия в системах газ-жидкость. Теоретические основы и экспериментальные методики. Моделирование химико-технологических процессов. Уч

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
1 Мб
Скачать

Окончание табл. 23

1

2

3

4

5

6

7

 

198,5

72,5

500

764

22

1,032

 

190,5

65

500

764

22

1,036

 

201

75

500

764

22

1,032

 

200,5

75

500

764

22

1,036

 

239

164

295

755

20,5

1,011

 

164,5

76

350

755

20,5

1,016

 

220,5

131

355

755

20,5

1,019

 

214,8

125,5

350

755

20,5

1,007

 

125,5

35,8

350

755

20,5

1,003

 

223

97,8

500

755

20,5

1,026

Водород

202

77,5

500

759

21,0

1,037

201

76,5

500

759

21,0

1,037

 

218,5

93,8

500

759

21,0

1,036

Азот

200,5

76

500

759

21,0

1,037

215,2

90,2

500

759

21,0

1,033

 

225,2

100

500

759

21,0

1,032

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

погрешности [62]

=1

- среднее арифметическое

ср =

- стандартное отклонение

 

=1( ср)

=

( 1)

61

 

= ,

 

 

 

доверительный интервал

 

 

 

 

 

=

100%

 

 

Метрологические

 

 

 

относительная погрешность объема

установки

 

 

 

 

 

 

 

 

 

ср

 

Таблица 24

 

характеристики измерений объема установки

 

 

 

 

 

 

Vсреднее (среднее выборки), дм3

 

 

 

1,02

 

 

 

 

 

 

7

 

Стандартное отклонение

 

 

 

0,01

 

 

 

 

1

 

 

 

 

 

 

 

Количество параллельных определений

 

22

 

 

 

 

Квантиль распределения Стьюдента (двусторонняя поста-

2,08

 

новка задачи, доверительная вероятность Р = 0,95)

 

 

Абсолютная погрешность, дм3

 

 

 

0,00

 

 

 

 

 

 

5

 

Относительная погрешность, %

 

 

 

0,50

 

 

 

 

 

 

 

 

Для оценки предельной погрешности разового измерения объема установки как результата косвенного измерения, рассчитываемого по формуле (24), использовался подход, изложенный в справочном руководстве [63]. Согласно рекомендациям при вычислении предельной погрешности предполагается, что систематические ошибки измерения различных величин суммируются наиболее неблагоприятным образом, то есть взаимно усиливают друг друга, в то время как на практике может иметь место частичная взаимная компенсация погрешностей. Формулы для расчета погрешностей косвенно измеряемой величины F с использованием величин систематических ошибок (погрешностей) результатов непосредственных измерений для наиболее простых и часто встречающихся функциональных зависимостей приведены в табл. 25.

62

Таблица 25

Расчет предельной погрешности косвенного измерения в зависимости от функциональной связи параметров

Функциональ-

 

Абсолютная системати-

 

Относительная сис-

 

ная зависи-

 

ческая погрешность δF

 

тематическая по-

 

мость

 

 

грешность εF

 

 

 

 

 

 

 

 

 

 

 

= 1 + 2

 

 

 

δ 1 + δ 2

 

 

δ 1

+ δ22

 

 

 

 

 

 

 

 

 

 

 

 

 

δ

 

 

δ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 1 2

 

 

 

δ

 

 

 

δ

 

 

1

 

 

 

 

2

 

 

 

 

 

1 + 2

 

 

δ 1

δ2

 

 

 

= 1 2

 

δ 1 2 + δ 2 1

 

 

 

 

 

 

 

 

1

+

 

 

2

 

 

=

 

 

 

δ

 

+

 

δ

 

 

δ 1

+

δ 2

 

1

 

 

 

1

 

1∙ 2 2

 

 

1

 

 

2

 

 

 

2

 

 

δ2

 

 

 

 

 

 

1

δ

2

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

δ−1

 

 

 

δ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Исходя из

уравнения (24) предельная

погрешность разового

 

измерения объема установки складывается

из погрешности

величин

Vводы, Рбаром, Рнач

и Ркон.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Паспортная погрешность средства измерения барометрическо-

го давления ± 2 мм рт. ст..

 

Примем среднее значение атмосферного

метра (см. табл. 21) относительная погрешность перепада давления Р

= (Рнач – Ркон) может быть представлена в виде

нач + кон

где δРнач = δРкон = 1 деление

нач − кон

=(усл. ед.).

давления равным 750 мм рт. ст., тогда относительная погрешность со-

ставляет 0,3 %. В соответствии с табл. 25 и характеристиками мано-

При этом согласно данным табл. 23, средний перепад давления Рнач - Ркон по всем экспериментам составляет 117 делений (усл. ед.).

Соответственно расчетное значение εР составляет 2/117 или 1,71 %.

63

Абсолютная погрешность Vводы составляет для мерной колбы объемом 500 см3 величину 0,60 см3, соответственно относительная погрешность равна 0,60/500 или 0,12 %. Суммируя относительные погрешности, получаем величину предельной погрешности измерения объема установки, равную 2,13 %. Действительная (экспериментальная) погрешность измерения объема установки по данным табл. 24 составляет

0,5 %.

Формулы для расчета случайных погрешностей косвенно определяемой величины приведены в табл. 26.

Таблица 26

Расчет случайной погрешности косвенного измерения в зависимости от функциональной связи параметров

Функциональ-

Абсолютная случайная

Относительная слу-

ная зависи-

чайная погрешность

погрешность F

 

 

мость

 

 

 

 

 

 

 

 

εF

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 1

 

+ 2

 

 

 

1 2

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

 

= 1

 

2

( ) + ( )

 

 

 

 

 

1

+ 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

 

 

 

1 2

 

2

2

 

 

 

 

 

 

 

 

 

 

 

= 1

2

( ) + ( )

2

 

 

 

1

2δ

 

 

 

 

 

 

 

 

12

 

 

δ

 

 

 

 

2

1 2

 

 

1

 

 

 

 

2

 

 

 

 

 

(

 

) + (

 

)

 

 

 

+

 

 

 

 

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δ

 

 

 

δ

 

 

 

 

 

δ

 

 

 

 

 

δ

 

 

= 1

1

 

 

1

 

 

2

 

 

 

1

 

 

2

 

 

 

2

2

1 +

 

2

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

1 +

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.5.2. Исследование равновесия газ-жидкость на примере системы водород – этанол

Для отработки методики в качестве объекта исследования была выбрана система водород – этанол, поскольку физические свойства

64

растворителя, требующиеся для расчетов, подробно отражены в литературе.

1.Барометрическое давление в ходе эксперимента 756 ± 2

ммрт. ст.

2.Рабочая температура 27,5 ± 0,5 °С.

3.Полный объем автоклава 1,027 ± 0,005 дм3.

4.Масса жидкости 0,651 кг (бинарная смесь этанол – вода, содержание этанола, рассчитанное через плотность раствора при ком-

натной температуре (23,7 °С, ρ = 0,813538 г/см3), составляет 93,6 % об. или 90,46 % мас., или 78,76 % мол. (расчет по программе

AlcoDens [64]).

5.

Плотность жидкости при рабочей температуре 810,2 кг/м3

(расчет по программе AlcoDens).

 

 

6.

Свободный от жидкости объем автоклава, занятый газом при

рабочей температуре, 0,223 дм3.

 

 

7.

Верхний предел измерений манометра 1 кгс/см2.

8.

Количество делений шкалы манометра 250 дел.

9.

 

1 735,56

 

Цена деления шкалы манометра:

 

ценаделения =

 

250

= 2,94 ммрт. ст.

10. Начальное давление избыточное 183 деления.

11. Конечное давление избыточное 128,5 деления. 12. Конечное= 756давление+ 183 2,942237полное = 1134,1 ммрт. ст.

кон

13. Давление насыщенного пара растворителя при рабочей температуре 59,6 мм рт. ст. (расчет по уравнению Антуана для бинарной смеси этанол – вода, содержащей 78,76 % мол. этанола). Параметры уравнения Антуана для этанола и воды приведены в табл. 27.

Таблица 27

Параметры уравнения Антуана (Antoine Equation)

Соединение

Интервал

Коэффициенты уравнения Антуана

температур,

 

[9]

 

 

К

А

В

С

Этанол

273 - 351,70

5,37229

1670,409

-40,191

 

 

 

 

 

Вода

273 - 303

5,40221

1838,675

-31,737

 

 

 

 

 

65

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

газа в жидкости (nГ составляет 10-4 nS), соответственно погрешность, возникающая при не учете мольной доли газа в жидкости составляет порядка 0,01 % отн.

14. Абсолютное количество растворенного газа (без учета газа, растворе( нногопри атмосферном) (183давлении128,5):) 2,942237 133,322

Г = нач кон Г 8,314 300,65 = 1,909=10−3моль =

15. Растворимость газа в жидкости (без учета газа, растворенного при атмосферном давлении), то есть количество растворенного

газа, измеренное по перепаду давления Рнач – Ркон:

моль/кг

Г = =

0,651

= 2,932 10−3

Г

1,909∙10−3

 

.

16.По аналогичному алгоритму обсчитываются результаты серии измерений при различных значениях Рнач.

17.Далее часть растворителя массой 181 г была слита из автоклава, после чего была выполнена серия измерений при той же температуре.

18.По графику зависимости моляльной концентрации растворенного водорода от его парциального давления (рис. 23) определяется величина растворимости газа в жидкости при рабочей температуре и барометрическом давлении (система соединена с атмосферой, избыточное давление отсутствует, при этом состав газовой фазы над раствором – рабочий газ и пары растворителя). Напомним, что нулевое значение по ординате соответствует парциальному давлению газа 696,4 мм рт. ст.

Как видим, величины свободного члена и углового коэффициента линейной регрессии, определенные по результатам двух серий экспериментов, практически не зависят от объема газовой фазы в установке (относительная погрешность свободного члена между двумя сериями на уровне 5 %, углового коэффициента - на уровне 6 %).

66

0,006

C*,

 

 

 

 

 

 

моль/кг

 

 

 

 

 

 

0,004

y = 7,9309E-06x - 5,5674E-03

 

 

 

 

 

 

 

0,002

 

R2 = 9,967E-01

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

0

200

400

600

800

1000

1200

1400

-0,002

 

 

 

 

 

 

РН2,

 

 

 

 

 

 

мм рт.ст

-0,004

 

y = 8,4452E-06x - 5,8697E-03

 

 

 

 

 

R2 = 9,918E-01

 

 

 

-0,006

 

 

 

 

 

 

 

Рис. 23. Моляльная концентрация растворенного водорода в зависи-

 

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

 

(масса растворителя 0,651 кг и 0,470 кг; t = 27,5°С;

 

 

Рбаром = 756 мм рт.ст.)

 

 

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

Таблица 28

Метрологические характеристики измерений для системы этанол – водород (t = 27,5 °C, Рбаром = 756 мм рт. ст.)

 

 

 

 

Характеристика

 

mS = 0,651

mS =

 

 

 

 

 

 

 

кг

 

0,470 кг

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

2

 

3

 

 

Величина

достоверности

аппроксимации

 

 

 

 

 

 

R2

 

 

Г

 

 

 

 

 

0,9967

 

0,9918

Г

 

 

 

Б

 

Г

Б

 

 

 

 

 

 

 

 

 

 

 

СmБ

как свободный член линейной

 

-3

5,8697

 

1

=

=

5,5674 10

 

 

 

 

 

 

 

 

0

-3

 

 

регрессии у = a + b x, моль/кг

 

 

 

 

 

 

67

Окончание табл. 28

1

 

2

 

3

 

1/(Кн MS) как угловой коэффициент ли-

7,9309 10

-6

8,4452 1

 

нейной регрессии у = a + b x

 

0-6

 

2

 

4,4290 10

-9

1,5888 1

 

Дисперсия оценки s0

 

 

0-8

 

Количество параллельных определений

10 (+1)*

 

10 (+1)*

 

Квантиль распределения Стьюдента

2,26

 

2,26

 

 

 

 

 

 

 

Стандартное отклонение sa

1,5981 10

-4

2,8217 1

 

0-4

 

Доверительный интервал

 

3,6118 10

-4

6,3771 1

 

(абсолютная погрешность) a

 

0-4

 

Относительная погрешностьεa, %

6,5

 

10,9

 

Стандартное отклонение sb

1,5349 10

-7

2,5595 1

 

 

0-7

 

Доверительный интервал

 

3,4689 10

-7

5,7845 1

 

(абсолютная погрешность) b

 

0-7

 

Относительная погрешностьεb, %

4,4

 

6,8

 

* – десять экспериментальных измерений и

точка «0» для

барометри-

 

ческого давления

 

 

 

 

 

Исходя из методики эксперимента, предельная погрешность

разового измерения CmГ*

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

величин: Т, Pбаром, ms, VГ,

Р = Рнач - Ркон. Величины «плотность рас-

творителя» и «давление насыщенных паров растворителя» брались как табличные значения, то есть могут рассматриваться как const, не имеющие погрешности. Абсолютная погрешность при взвешивании (погрешность массы растворителя) составляла 1 г, что при массе растворителя порядка 500 г составляет 0,2 %. Максимальная погрешность объема газовой фазы при загрузке порядка 50 % объема системы жидкостью может быть принята равной абсолютной погрешности суммарного объема установки, то есть 0,005 дм3, при объеме газовой

равной 0,5 дм3. Соответственно относительная погрешность VГ равна 0,005/0,5 или 1 %. Относительная погрешность измерения температуры 0,5 %. Относительная погрешность измерения барометрического давления равна 2/756 или 0,3 %. Относительная погрешность измере-

ния перепада давления Р при разнице Рнач - Ркон 20 делений (см.

68

табл. 26) составляет 2/20 или 10 %. Таким образом, предельная погрешность величины nГ* составляет 0,2 + 1,0 + 0,5 + 0,3 + 10,0=12,0 %,

а погрешность величины СmГ* равна 12,0 + погрешность массы растворителя или 12,2 %. Как видно из табл. 28, максимальное значение относительной погрешности коэффициентов уравнения линейной рег-

рессии СmГ* = f(Ргаза) составляет 10,9 %, что не превышает предельной погрешности, равной 12,2 %.

2.5.3. Исследование равновесия газ-жидкость на примере системы водород – толуол

Данная система была выбрана в качестве базового объекта для определения надежности предлагаемой методики. В работе использовался электролитический водород марки «Б» (ГОСТ 3022-80 изм. № 1) с содержанием основного вещества не менее 99,95 % об. Характеристики исходного растворителя приведены в табл. 29. Перед экспериментом растворитель дополнительно перегонялся при атмосферном

давлении и выдерживался 24 ч над свежепрокаленным при t = 400 °С гамма-оксидом алюминия. Параметры уравнения Антуана для толуола приведены в табл. 30.

 

Таблица 29

Характеристики толуола марки «чда» (ГОСТ 5789-78)

 

 

 

Наименование показателя

Норма

Массовая доля основного вещества, % мас., не менее

99,50*

 

 

 

 

Массовая доля воды, % мас., не более

0,03

 

Массовая доля общей серы (S), % мас., не более

0,0003

 

Массовая доля нелетучего остатка, % мас., не более

0,001

 

Массовая доля тиофена, % мас., не более

0,0001

 

* – при расчетах принималась за 100 %; молярная масса толуола 92,1384 г/моль

Таблица 30

Параметры уравнения Антуана (Antoine Equation) для толуола

Интервал темпера-

Коэффициенты уравнения Антуана [9]

тур, К

 

 

 

А

В

С

273,13 - 297,89

4,23679

1426,448

-45,957

 

 

 

 

69

Равновесие газ – жидкость для данной системы было исследовано в интервале температур 14 – 59 °C. В качестве примера приведены результаты эксперимента при температуре 46,2 °C (табл. 31).

Таблица 31

Mетрологические характеристики измерений для толуола (Рбаром = 762 мм рт. ст., масса растворителя 537 г)

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

 

Г

 

Г

Б

Г

Б

0,9997

 

 

Cm как свободный член линейной регрессии

3,4316 10-3

 

Б

=

=

у = a + b x, моль/кг

 

 

 

1/(Кн MS) как угловой коэффициент линейной регрес-

5,0161 10-6

сии

 

 

 

 

 

у = a + b x

 

 

 

 

Количество параллельных определений

5 (+1)

Квантиль распределения Стьюдента

2,78

 

 

 

 

2

 

 

2,33125 10-

Дисперсия оценки s0

 

 

10

Стандартное отклонение sa

 

3,5393 10-5

Доверительный интервал

 

9,8392 10-5

(абсолютная погрешность) a

 

 

Относительная погрешностьεa, %

 

2,9

Стандартное отклонение sb

 

3,1346 10-8

Доверительный интервал

 

8,7142 10-8

(абсолютная погрешность) b

 

 

Относительная погрешностьεb, %

 

1,7

Значения константы Генри, рассчитанные по величине свободного члена CmБ для парциального давления1 газа, Г)баром = Рбаром

Рs и по величине углового коэффициента приведены в табл. 32.

70

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]