Принципиальная блок-схема алгоритма расчета Кэ между КС представлена на рис. 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)
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.