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

5 = 3,6^-l,3|l-^-1'5|o.                                                                . 13)

44


Получим алгоритм   динамического   коэффициента участка в момент s

Is = 2£s_i — ls^2 — Л5_1 ехр {— ys2 sign

ГГ                        1

kN


 для &-го

x


Xsigu(pl-{ ,   l

:                                      (П.14)

где Qo — расход в начале газопровода; Qn — бтбор газа п-м
потребителем; N>—,число мест отборов до &-го элементарного
участка;                                 ■-...                                                                -•■'■-

 — шаг дискретизации по длине МГ; Дг— интервал кванто­вания по времени; с — скорость звука в газе;

As «= 4s-i ехр (— ySs sign"vл в);                                                                                        (II. 15)

лг и /г — постоянные величины,   выбираются   таким   образом, чтобы

mn = ~ внутренний диаметр газопровода;

-^ -1*-> ~ А^ ^ х


X si


Р2.ж_, -


X



X (^-si


+ /t) [.(11.16) 45


Если вид дрейфа выбран правильно, то Ys-^О при s->-oo.

Время сходимости алгоритма зависит от выбора начальных параметров go, ^о, g-i. Для оптимального выбора начальных па­раметров необходимо накопить информацию за число шагов, равное по крайней мере числу задаваемых параметров (в дан­ном случае трем), и затем определить g0, Ао, S-i, которые и принимают за начальные приближения.

Принципиальная блок-схема алгоритма представлена на рис. II.7.

1

Вычисление



Запись Лдовнешнюю

память. ПечатьЛ

(^Начало)

I   ВдодI данныхI

±

Определениеплощади сечениятруды\ расстоя­ниеотначалагазопро­водадок- и. точкиизмере­ниядавления,величины

(l-x)h_________

Нет


I

Нет


Определение Величины

Запись В ОЗУ информации, необходимой, для следующего шага идентификации

I


Рис. П.7. Блок-схема алгоритма оценки эквивалентного коэффициента гидрав­лического сопротивления

Для работы алгоритма необходима следующая информация; Н — число трчек замера давления; - М —• число элементарных участков; At — интервал дискретизации по времени; N — интер­вал текущего сглаживания; D — диаметр газопровода; Ах — интервал дискретизации по координате; ц — число точек попут­ных отборов газа; г — допустимая абсолютная погрешность; т]о[О, 1 : М], ц-\[—1, 1 : М] — начальные значения параметра для М элементарных участков на двух временных слоях (двумер­ный массив); ро[\ : Н], А0[0, 1 : М] — начальные значения дав­лений и параметра; L — число попутных отборов газа до ^-го элементарного   участка;   ps[0 : Щ,   Qs~i[0 : Н] — соответственно

46


давления и попутные отборы; Ts[l: M], zJ[l:M] — абсолютные температуры и коэффициенты сжимаемости газа для каждого элементарного участка.

Описание принципиальной блок-схемы

1.  Ввод необходимых данных ОЗУ.

2.  Определение скорости звука в газе и коэффициентов (Зь,
для первого элементарного участка.

3.  Определение параметра А   для первого   элементарного
участка.

4.  Определение коэффициента уп,$, для первого элементар­
ного участка.

5.  Подготовительные операции, необходимые для определе­
ния параметра |ft(S.

6.  Определение оценки параметра Zk,s для первого элемен­
тарного участка на 5-м шаге итерации.

7.  Определение момента останова счета   (аналогично при­
меру II.1).

При выполнении условий останова — определение X и вы­дача его на печать. В случае невыполнения условий останова — запись в ОЗУ параметров, необходимых для работы алгоритма на следующем временном шаге и для остальных элементарных участков.

8.  Аналогично выполняются пп, 1—7 для остальных элемен­
тарных участков.

9.  При поступлении новой информации аналогично выполня­
ются пп. 1—8.