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


Принципиальная блок-схема алгоритма расчета Кэ между КС представлена на рис. II.5.

Исходная информация для расчета: число участков 7; число ниток МГ—1; число временных слоев — 9; Д£=120 мин — ин-



МассиЛ состояния перемычек

1[1]

Оцвнидание

коэффициента

Лэ для к-го

элементарного

участка


Ночей,

Оценивание Я3 для участий между КС

I


Оценаданае коэффициента

для к-го элементарного участка

Нет

Оценивание коэффициента Лэ оля участка между K~MuQi+i)-M

замерами


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

тервал дискретизации по времени; |о = 7 — начальное значение параметров; 1=0,1—величина,   используемая   при остановке

процесса оценивания; р = 0,60 — плотность транспортируемого газа; ДХ2=10 000 м — интервал дискретизации по длине МГ; <7g2=0 — дисперсия ошибок измерения величины массового рас­хода газа; ар2 = 0 — дисперсия ошибок измерения давления; £ = 0,7 — диаметр газопровода.

42


Результаты расчета приведены в табл. II.5.

Таблица     II.5


Время


1 11,3019


8,7347


7,1315


8,1293


8,5638


8,0113


8,6505


8 8,4634


8.42


Истинное значение |э=8,33.

Сходимость   £э   к   истинному   значению   приводится   на рис. II.6.



Рис. II.6. Кривые схо­димости эквивалент­ного коэффициента гидравлического со­противления


ист



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

Рассматривается i-й элементарный участок в начальный мо­мент времени. Расчет \ для всех остальных участков и других моментов времени проводится аналогично.

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

2.  Расчет £э для М параллельных ниток газопровода:
а) площадь сечения трубопровода

б) расстояние от начала газопровода до kточки измере­ния давления

k хь = V Ахл

в) нахождение £>=(/—х)Ах аналогично 2а, 26 для всех то­чек измерения давления.

3.  Определение уп.

4.  Нахождение ii = £o+{*} согласно алгоритму.

5.  Проверка на «останов счета» — проводится аналогично
примеру II. 1. При выполнении условий останова данные выда-

43


ются на печать. В противном случае — запоминание информа­ции, требуемой для следующих временных слоев.

Пример И.4. Оперативная оценка динамически изменяюще­гося коэффициента гидравлического сопротивления магистраль­ного газопровода Я.       .   .

Рассматривается газопровод, коэффициент гидравлического сопротивления которого изменяется во времени. Изменение ко­эффициента гидравлического сопротивления во времени обус­ловлено образованием и постепенным нарастанием гидратной пробки, загрязнением газопровода песком и другими частица­ми, скоплением влаги и конденсата, а также кажущимися (фиктивными) изменениями за счет неучета некоторого коли­чества газа при вычислении I в случае аварийной утечки. Ско­рость изменения Я различна.

При медленных изменениях коэффициента гидравлического сопротивления имеется возможность «следить» за трендом Я. В этом случае необходимо использовать алгоритм динамиче­ской оценки.

Математическая,, постановка задачи оперативной оценки ди­намического коэффициента . гидравлического сопротивления участка магистрального газопровода формулируется следую­щим образом.

По имеющейся модели газопроводами измеренным значени­ям расхода в начале газопровода л давлениям вдоль трассы магистрального газопровода определить оценку параметра X в момент t, изменяющегося во времени согласно известной функциональной зависимости

где вектор коэффициентов а заранее не известен.

Задача решается по методу динамической стохастической аппроксимации.

На основании использования алгоритма статистической оцен­ки коэффициента

sух

и выражения тренда коэффициента [25]

l(t) = -~t*~Bt~l0,                                                                                        (ini)

где

Л = 2-1,3^1 -±уи510,                                                                (ПЛ2)