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

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

Is = 2£_s_i — ls^₂ — Л₅_1 ехр {— y_s2 sign

ГГ                        1

kN

 для &-го

x

Xsigu(pl-_{ ,   l

:                                      (П.14)

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

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

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

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

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

-^ -¹*-> ~ ^А^ ^ ^х

X si

Р2._ж_, -

X

X (^-si

+ /t) [.(11.16) 45

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

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

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

1

Вычисление

(^Начало)

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

±

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

(l-x)h_________

Нет

I

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

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

I

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

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

46