Page 26 - 6747

24
Для обчислення фактору стисливості за рівнянням стану (2.53)

необхідно визначити густину  при заданих тисках (р, МПа) і температурі (Т,
м
К).

Густина  із рівняння стану (2.53) обчислюється за методом Ньютона в
м
наступному ітераційному процесі:

1) початкову густину визначають за формулою

3
 ) 0 (  9  10 p /[RT 1 , 1 ( p  7 , 0 )], (2.66)
м п
де приведений тиск обчислюється із виразу

3

p  10 R T , 0 ( 28707  , 0 05559 ), (2.67)
пк пк пк
p  p / p , (2.68)
п пк

Псевдокритичні густину  , температуру Т  і фактор Пітцера  
пк пк
обчислюються за формулами (2.55), (2.57) і (2.60);

2) густина на k-му ітераційному кроці визначається із виразів

3
 (k )  [ 10 p  RTz (k ) 1   (k ) 1  ] /[RT 1 (  A )], (2.69)
м м 1
 (k )   (k ) 1     (k ) , (2.70)
м м м

де z ( k ) 1 обчислюється із рівняння стану (2.53) при густині на

ітераційному кроці (k-1), тобто при  ( k ) 1 , а безмірний комплекс А 1 визначається
м
із рівняння

r S k
l
A   k (  c)1  k / T , (2.71)
1 kl п п
k  l1 0
3) критерії закінчення ітераційного процесу

 (k ) /  ( ) k  10  6 , (2.72)
м м

Якщо критерій (2.72) не виконується, то необхідно продовжити

ітераційний процес, починаючи з пункту 2) алгоритму.

Після визначення фактору стисливості при робочих і стандартних умовах

за формулою (1.22) обчислюється коефіцієнт стисливості.