Промышленность \ Добыча, хранение и транспортировка нефти и газа

Методы повышения эффективности процессов добычи и транспорта газа, страница 17

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

Математическая постановка задачи: по имеющейся матема­тической модели магистрального газопровода (см. пример II.1) и измеренным с помехами параметрам общего расхода газа по магистральному газопроводу и давлениям в точках замера найти оценку относительного коэффициента гидравлического сопротивления.

Из конечно-разностной аппроксимации ^исходной системы дифференциальных уравнений с использованием метода стоха­стической аппроксимации следует алгоритм определения коэф­фициента g, где l = \/il в момент s для 1-й нитки газопровода.



k—1

i=l   n=l


N,


fe-l


\MPk.*-\^Pkti


,S-1 — Pk-l,


(II.3)

где

 = 1 /j] ^'


N

 = Qo,s-. -


m-i   Nt

 2


T


m


T


 *"1 ~~ Pm'


 Pm_i >s_i      Pm-\ ,s) '


(II.4)


.S—1 =

Qo6i4 — количество газа, перекачиваемого по всей системе; kv — коэффициент расхода по каждой нитке газопровода; L — число ниток;

knl = Ui  /Uo   ,

36


Do — диаметр газопровода, относительный расход которого ра­вен единице. Обычно принимается Do = 700 мм;

х.^ — 1..   в - ' У~^' •
1At   /а '   ™      с   V   Axt   '                                                                                ;

Sk = s\gn(pl х , i — Pi , i) l/^lpf., , ,■—p2 e , I;                                                                                                   (H-5)

ft   w      :

2   2 Qin — количество отборов на & участках   газопровода;

Аде —шаг дискретизации по длине газопровода; At — интервал квантования во времени.

Принципиальная блок-схема алгоритма аналогична принци­пиальной блок-схеме алгоритма оперативной оценки отдельных участков газопровода (см. пример II.1).

Для работы алгоритма необходима следующая информация:
И — число точек замера давления; М — число элементарных
участков; L—интервал дискретизации по времени; D[\ : М,
1 : z] — число ниток газопровод^; д;[1 : М, (1 : z] — диаметр газо­
провода; г —длина элементарного участка; ц[\ : z] — число то­
чек попутных отборов; |[0, 1 : М, \ \ z] — начальное значение па­
раметра; е — допустимая абсолютная погрешность; р[0 : М],
Q[\ :М], Т[{ : М\ z[i:M] — соответственно: давления,, попутные
отборы, абсолютные температуры и коэффициенты сжимаемо­
сти для каждого элементарного участка.                ;; i

Размерности и диапазоны изменения входных и выходных констант программы приводятся ниже.

Константа

ния

Константа

НИЯ

Я

0—2Qp,

е

0,1—0,3

Л!

0—200

0,3—0,5

А/, с

0—7200

р, кгс/см2

20-10*—76.10*

z

2—5

Q, кг/с

50—150

D, ы

0.7-1.4

т5, к

220—350

Ц

0—200

0,8—1,0

Ах, м

10000—50 000

0,2-1,1

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


Длина участка Ллг1=Дх2=Л;са, м 15 000    За эталонный  участок  выбран
Величина интервала по времени  3600                й £08

Число элементарных участков        3
Число временных слоев                  15

Диаметр газопровода DX=D2=   1,4 =D3, м


Скорость звука с, м/с Попутных отборов нет.

 у                 р

 первый. £0=0,8 для всех эле- ментарных участков.


375


37


Результаты измерения расхода и давлений   приведены   в табл. П.З.

Скачать файл