Бесплатный автореферат и диссертация по наукам о земле на тему
Разработка методов прогнозирования технологических режимов магистральных газопроводов в условиях неопределенности спроса на газ
ВАК РФ 25.00.19, Строительство и эксплуатация нефтегазоводов, баз и хранилищ

Автореферат диссертации по теме "Разработка методов прогнозирования технологических режимов магистральных газопроводов в условиях неопределенности спроса на газ"

На правах рукописи

004603350

Казак Константин Александрович

РАЗРАБОТКА МЕТОДОВ ПРОГНОЗИРОВАНИЯ ТЕХНОЛОГИЧЕСКИХ РЕЖИМОВ МАГИСТРАЛЬНЫХ ГАЗОПРОВОДОВ В УСЛОВИЯХ НЕОПРЕДЕЛЕННОСТИ СПРОСА НА ГАЗ

Специальность: 25.00.19 - «Строительство и эксплуатация нефтегазопроводов, баз и хранилищ»

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата технических наук

- з июн 2910

Москва-2010

004603350

Работа выполнена в Обществе с ограниченной ответственностью «Научно-исследовательский институт природных газов и газовых технологий - Газпром ВНИИГАЗ»

Научный руководитель:

доктор технических наук, профессор Лесных Валерий Витальевич

Официальные оппоненты: доктор технических наук

Черний Владимир Петрович

доктор технических наук, профессор Марон Вениамин Исаакович

Ведущая организация:

Московский филиал ОАО «Гипроспецгаз»

Защита состоится «09» и(£)к$ 2010 года в 13 часов 30 минут на заседании диссертационного совета Д 511.001.02 при Обществе с ограниченной ответственностью «Научно-исследовательский институт природных газов и газовых технологий - Газпром ВНИИГАЗ» (ООО «Газпром ВНИИГАЗ») по адресу: 142717, Московская область, Ленинский район, поселок Развилка

С диссертацией можно ознакомиться в библиотеке ООО «Газпром ВНИИГАЗ»

Автореферат разослан « 30 у> ЗЬрё/Я 2010 года

Ученый секретарь диссертационного совета, кандидат технических наук

И.Н. Курганова

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность работы. В настоящее время в связи со строительством магистральных газопроводов рассчитанных на высокое давление, в том числе морских газопроводов («Голубой поток», «Северный поток», «Южный поток») возрастает необходимость разработки более точных технологических моделей магистрального транспорта газа. Это относится и к прогнозированию граничных условий при расчете режимных показателей газопровода (температура, давление и расход газа). Так, например, при относительном отклонении давления в начальном узле на 1% значение расхода в конечном узле может иметь отклонение до 5%. То есть, незначительная ошибка прогноза в результате газодинамического расчета может увеличиться на порядок. Таким образом, должны быть разработаны методы идентификации и адаптации значений давления, температуры и расхода газа в граничных узлах моделируемого участка газопровода, что позволило бы получать прогнозируемые режимные показатели в пределах точности измерительных приборов.

Основной задачей системы магистрального транспорта газа является надежность обеспечения поставок газа потребителям. Следует отметить, что потребители природного газа, как правило, привязаны к узлам газотранспортной системы (ГТС). В зависимости от объемов потребления распределение газа в узлах ГТС может осуществляться непосредственно на промышленные объекты или через газораспределительные станции. При формировании планов развития и реконструкции ГТС необходимо исходить из достоверных прогнозов спроса на газ в этих узлах, чтобы наиболее точно определять динамику развития мощностей газопроводов. В кризисных условиях это особенно важно для минимизации неиспользуемых резервов ГТС и обеспечения надежности поставок газа потребителям в соответствии с их требованиями. В связи с этим должны быть разработаны методы прогнозирования потребности в газе, позволяющие адаптировать развитие системы магистрального транспорта газа к текущей ситуации в газовой промышленности.

Таким образом, в силу указанных выше причин возникает необходимость разработки новых, а также развития и адаптации существующих моделей и методов прогнозирования в задачах функционирования и развития систем магистрального транспорта газа.

Цель работы состоит в развитии методов краткосрочного и перспективного прогнозирования в граничных узлах газопроводов их технологических параметров для моделирования функционирования и развития систем магистрального транспорта газа.

Основные задачи исследования:

• разработка моделей краткосрочного прогнозирования диспетчерских данных в условиях неизменяющегося технологического режима работы ГТС, при ограничениях на время оценки;

• анализ результатов прогноза и разработка на его основе вероятностной методики оценки работоспособности датчиков и соответствующей телеметрии;

• разработка адаптивных моделей краткосрочного прогнозирования диспетчерских данных, зависящих от управляющих параметров, в условиях изменяющегося технологического режима работы ГТС;

• обоснование моделей долгосрочного прогнозирования спроса на газ с учетом различных экспертных оценок для потребителей, имеющих контракты на поставку газа;

• анализ применимости центральной предельной теоремы для определения различных уровней спроса на газ (пессимистичный, наиболее вероятный, оптимистичный) в задачах долгосрочного прогнозирования спроса при наличии больших выборок потребителей с применением международной классификации существующих и перспективных потребителей по категориям спроса;

• анализ подходов и разработка методики экспертной многопараметрической характеристики проектов поставок газа для численной оценки приоритетов реализации проектов в динамике.

Научная новизна. В диссертационной работе впервые проведены комплексные исследования по разработке обобщенных подходов к задачам краткосрочного и долгосрочного прогнозирования характеристик узловых точек газотранспортных систем с применением современного математического аппарата. Наряду с традиционными методами прогнозирования (регрессионный анализ; модели авторегрессионного проинтегрированного скользящего среднего; статистические, теоретико-множественные и теоретико-вероятностные подходы) используются нейросетевые алгоритмы и методы нечеткой логики. При решении задач долгосрочного прогнозирования впервые применена международная классификация потребителей по категориям спроса в зависимости от оценок вероятности реализации различных проектов использования газа. Построена универсальная экспертная многопараметрическая характеристика проектов поставок газа для формирования приоритетов их реализации в виде набора нечетких переменных как составной части экспертной системы.

Практическая значимость. Алгоритмы краткосрочного прогнозирования могут быть использованы для целей оперативного управления работой газотранспортных систем и согласования технологических характеристик в узловых точках магистральных газопроводов, а также в составе задачи по газодинамическому моделированию магистрального транспорта газа. Статистический анализ краткосрочных прогнозов позволяет выявлять датчики и соответствующие им станции телеметрии, функционирующие с большой вероятностью некорректно и требующие проверки работоспособности.

Модели долгосрочного прогнозирования позволяют определять вероятностную динамику перспективного спроса на газ на возможных объектах его потребления сроком на десятки лет. Результаты долгосрочного прогнозирования спроса на газ предназначены для формирования задач развития и реконструкции ГТС на перспективу. Использование разработанных методиче-

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

Результаты диссертационного исследования использованы при разработке Генеральной схемы развития газовой промышленности Боливариан-ской Республики Венесуэлы, а также при разработке математических моделей в составе программного комплекса по газодинамическому моделированию магистрального транспорта газа.

Защищаемые положения:

• методика краткосрочного прогнозирования показаний датчиков телеметрических станций ГТС (диспетчерских данных), состоящая из: корректирующей фильтрации данных; построения модели регрессии, представляемой в виде линейного тренда и периодической составляющей; корректировки прогноза динамически пересчитываемыми моделями авторегрессии;

• методика применения нейросетевого прогнозирования с использованием усовершенствованной реализации алгоритма обратного распространения ошибки для краткосрочного прогнозирования показаний датчиков, зависящих от неявных управляющих параметров, в период изменения технологических режимов работы ГТС;

• вероятностная методика оценки корректности функционирования датчиков телеметрических станций на основе анализа статистики погрешности регрессионных прогнозов с последующей рекомендацией к проверке работоспособности;

• методика долгосрочного прогнозирования спроса на газ по отдельным потребителям с применением регрессионных моделей, эмпирически оптимизированным выявлением периодической составляющей, построением гладкой корректирующей функции в виде кубического сплайна для опорных значений, с возможностью учесть экспертные предположения о будущих значениях и с возможностью взвешивать известные значения;

• методики построения различных уровней спроса на газ (пессимистичный, наиболее вероятный, оптимистичный) в задачах долгосрочного прогнозирования спроса при наличии больших выборок потребителей с применением международной классификации существующих и перспективных потребителей по категориям спроса;

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

Апробация исследований. Материалы исследований по теме диссертационной работы докладывались и обсуждались на секциях Ученого совета ООО «Газпром ВНИИГАЗ», на III Международной научно-технической конференции «Газотранспортные системы: настоящее и будущее (GTS-2009)»

(Москва, 2009 г.) и на III Всероссийской научно-технической конференции «Безопасность критичных инфраструктур и территорий» (Екатеринбург, 2009 г.).

Публикации. Материалы диссертации опубликованы в 12-ти работах, из них 10 работ опубликовано в журналах, включенных в «Перечень...» ВАК Минобрнауки РФ, в том числе опубликована одна монография.

Структура и объем диссертации. Диссертация состоит из введения, 4-х глав, выводов и списка литературы из 59 наименований (2 на английском языке). Работа изложена на 157 страницах машинописного текста, содержит 52 рисунка и 6 таблиц.

КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ

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

В первой главе дано описание ГТС как объекта прогнозирования, представлен обзор подходов к прогнозированию и методов построения прогнозов. Подобные математические модели развиваются уже давно. Это, прежде всего, классические методы математической статистики: анализа и прогнозирования временных рядов (Бокс Дж., Дженкинс Г.); а также более современные: нейросетевые модели, в том числе и прогнозирующие (Гор-бань А.Н, Weigend A.S.), и методы нечеткой логики (Заде JI.A., Рыжов А.П.). Применительно к газовой отрасли: экономико-математическая модель динамического моделирования прогнозного спроса на газ на основании эластичности спроса (Кисленко H.A., Кудрявцев A.A.), модели нейросетевого прогнозирования среднесрочного газопотребления (Агаев Н.Б.); а также задачи моделирования, прогнозирования и оптимизации характеристик функционирования и развития трубопроводов, которые широко представлены в работах Галиулина З.Т., Леонтьева Е.В, Сухарева М.Г., Ставровского Е.П., Яковлева Е.И., Казака A.C., Черняева В.Д., Сощенко А.Е., Сарданашвилли С.А. В этих работах дано теоретическое описание и показаны возможности практического применения моделей регрессии и авторегрессии, а также экспертного анализа для рассматриваемого класса задач. Однако в силу вышеуказанных причин существующие подходы прогнозирования требуют дальнейшего развития и возникает необходимость разработки и адаптации моделей и методов прогнозирования в задачах функционирования и развития систем магистрального транспорта газа.

В главе формулируется общая постановка задачи прогнозирования как построения некоторой прогнозирующей функции. Описываются различные принципы отбора исходных экспертных и статистических входных данных. Приводятся методы последующей корректировки отобранных данных: сгла-

живание данных, корреляционный анализ, редукция. Рассматриваются традиционные для решения задач экстраполяции (прогнозирования) временных рядов модели и методы математической статистики. Даются общие рекомендации по выявлению составляющих (тренд, сезонная составляющая и случайный шум) и приведению временного ряда к стационарному путем дискретного дифференцирования ряда.

Сформулирована модель множественной регрессии в обобщенном виде и описаны основные методы нахождения неизвестных коэффициентов модели: метод наименьших квадратов, метод взвешенных наименьших квадратов с рекомендацией по выбору весовых коэффициентов (для задач прогнозирования) в виде неубывающей по времени последовательности, метод максимального правдоподобия. Приводятся модели авторегрессии, скользящего среднего и их вариации (АРСС и АРПСС), модель спектрального анализа. Даны рекомендации по применению данных моделей в задачах прогнозирования применительно к объектам ГТС.

Показано, что для ряда задач прогнозирования (например, с большим количеством входных параметров или с плохо выраженной аналитической зависимостью между входными и выходными данными) целесообразно применять нейронные сети. Дано описание природы нейронных сетей, сформулирована математическая модель нейронной сети: искусственный нейрон, весовые коэффициенты и смещения, функция активации, архитектура и процесс функционирования нейросети. Поставлена задача обучения нейронной сети в общем виде, приводятся основные группы алгоритмов обучения нейронной сети как задачи многоэкстремальной невыпуклой оптимизации.

Одним из наиболее распространенных алгоритмов обучения многослойных нейронных сетей является алгоритм обратного распространения ошибки. Это итеративный градиентный алгоритм обучения, который минимизирует функционал ошибки Е = ||у — уЛ Ц2 для реального выхода нейросети у и желаемого ук по обучающей паре векторов (хк,ук). Особенностью данного алгоритма является метод вычисления градиента функционала Е как многомерной функции от весов и смещений нейросети, которое происходит от выходов нейронной сети к ее входам в направлении, обратном прямому распространению сигналов в обычном режиме функционирования. Данный процесс был реализован в виде усовершенствованного векторно-параметрического алгоритма (см. рисунок 1). Здесь п - количество слоев нейросети; /.(•) - активационная вектор-функция некоторого /-го слоя нейросети, составленная из всех функций активации нейронов данного слоя (подразумевается, что производная от вектор-функции берется покомпонентно); XV, - матрица весовых коэффициентов г'-го слоя, в строках которой стоят наборы весов соответствующих нейронов; Ь, - вектор смещения г-го слоя, составленный из значений смещений всех нейронов слоя; хкф - входной вектор (набор входных сигналов) /-го слоя, получаемый при подаче на вход всей сети вектора хк.

Рисунок 1. Блок-схема алгоритма вычисления градиента функционала ошибки

Отметим, что введенные в данном алгоритме параметры связаны соотношениями, описывающими аналитически функционирование нейронной сети, которые могут быть использованы при ее практической реализации: х(1) = х,

усо=/;-К-ХСО+Ь,), _/=],...,«, (1)

»0+1) = у (Л,

У = У(и).

После того, как посчитан градиент ошибки, дальнейший процесс обучения становится стандартной задачей оптимизации, решаемой с помощью различных вариаций метода градиентного спуска. При этом могут использоваться всевозможные алгоритмы динамического изменения длины шага вдоль антиградиента и случайного «выбивания» из локального минимума.

Предложены возможные методы решения проблем обучения нейронной сети:

- варьирование количества нейронов сети с целью их минимизации при условии сохранения точности обучения;

- разбиение обучающих примеров на обучающее и контролирующее подмножества;

- обучение нейронной сети по сглаженным данным, с постепенным уменьшением степени сглаженности.

Даны общие рекомендации по преобразованию входных-выходных параметров для более эффективного использования области наибольшей чувствительности функции активации, а также повышения их различимости путем равномерного перераспределения данных значений по некоторому интервалу.

В главе также показано, что повысить эффективность экспертных систем в задачах прогнозирования можно за счет использования нечетких или лингвистических входных переменных, влияющих на прогнозируемую величину. Также нечеткая логика может использоваться для представления знаний эксперта с помощью правил для лингвистических переменных.

Во второй главе рассматриваются методы краткосрочного прогнозирования технологических параметров внешнего воздействия на ГТС. Под технологическими параметрами внешнего воздействия на газотранспортные системы понимается динамика значений давления, температуры и массового расхода газа в узлах входа и выхода рассматриваемой ГТС. Следует отметить, что расчеты режимов работы газотранспортных систем должны обеспечивать адекватное моделирование различных технологических ситуаций, которые включают в себя, в том числе: эксплуатацию газопроводов в штатных режимах; изменение объемов поставок газа поставщиками; изменение объемов потребления газа; изменение технологических параметров на входах и выходах ГТС.

Прогнозное моделирование динамики давления, температуры и массового расхода газа предложено осуществлять на основе ретроспективных данных газотранспортной диспетчерской службы в четыре этапа:

1) сглаживание (фильтрация) замеренных данных;

2) построение регрессионных моделей для прогнозирования трендов сглаженных временных рядов;

3) корректировка получаемых прогнозов моделями авторегрессии;

4) прогнозирование технологических характеристик при изменении управляющих воздействий методами нейросетевого анализа.

Сглаживание ретроспективных данных осуществляется для тех точек исходного ряда х(/), которые не удовлетворяют определенному условию гладкости. Для рассматриваемой задачи в работе определяется условие гладкости следующим образом:

1) строится новый временной ряд

г(/) = х(! +1) - 2х(1) + х{1 -1), (2)

который является дискретным аналогом второй производной по времени / от исходного ряда х(1);

2) из множества всех значений ряда отбрасываются 5% (от общего количества членов ряда - И) минимальных значений и столько же мак-

симальных значений. Для оставшегося множества находятся его минимальный и максимальный элементы - Ь . и I „ соответственно. То есть,

Ш111 Шал '

если Ъ - отсортированный по возрастанию массив значений {г(/)}, то

1^=2(1 + 1^.0.05]), 4„=г(ЛМЛЛ0.05]), (3)

где [о] - операция округления числа до ближайшего целого. Отметим, что величины 1т(п и ¿тах являются аналогами квантилей уровней 0.05 и 0.95 соответственно для множества значений {г(/)}; 3) далее по найденным значениям ЬтЫ и £шах определяются точки исходного ряда х(1), для которых соответствующие точки ряда г(/) удовлетворяют условию

^ 2(0 < ¿тах. (4)

Данные точки ряда х(1) и только они удовлетворяют условию гладкости.

Точки ряда х(г), неудовлетворяющие условию гладкости - неравенству (4), заменяются сглаженными значениями. Для этого предложено использовать модифицированный метод экспоненциально взвешенного сглаживания с некоторым окном локального усреднения. Как показали практические расчеты и анализ их результатов, наиболее эффективно использовать окно ширины 7 с центром непосредственно в точке сглаживания:

¿е-'-хО + О

т = -• (5)

IX

(—3

На рисунке 2 приведен пример реализации предложенного метода частичного сглаживания для диспетчерских данных.

2.99

2.94 -1-------.-----1

0 24 48 72 96 120 144 168 192 время(час)

♦ Исходный рад О Точки коррекции -Сглаженный рад

Рисунок 2. Сглаживание показаний давления на датчике Р-1

Показано, что краткосрочное прогнозирование динамики технологических характеристик газотранспортной системы может быть основано на использовании методов регрессионного анализа. При этом в качестве исходных временных рядов используются результаты частичного сглаживания замеренных данных, по которым с помощью методов обработки статистических данных могут быть построены различные регрессионные модели (уравнения). В работе получены уравнения множественной регрессии для экстраполяции ретроспективных данных на последующие временные шаги, а также для построения моделей авторегрессии, позволяющих корректировать краткосрочные прогнозы режимных характеристик работы ГТС с целью повышения их точности.

В результате проведенных предварительных расчетов и их анализа выявлено, что в качестве уравнений регрессии для оценки краткосрочных трендов динамики технологических характеристик ГТС целесообразно использовать семейство функций вида:

к

ДО = а0 + «,/+£( «2, С0!3Д' + яшД/), (6)

где - неизвестные коэффициенты регрессионной модели, кото-

рые могут быть найдены методом наименьших квадратов. Выражение под знаком суммирования представляет частичную сумму ряда Фурье, состоящую из К гармоник. Последовательность коэффициентов {Д}*,, зависящая

от предполагаемого периода исходного временного ряда, определяет, какие именно гармоники войдут в частичную сумму ряда Фурье. Натуральный параметр К и вещественная последовательность {Д}*, подбираются эмпирически с применением регрессионного анализа.

Проведенные исследования показали, что для аппроксимации ретроспективных технологических характеристик эксплуатации ГТС функцией вида /(г) достаточно первых 16-ти гармоник ряда Фурье (К = 16). При этом предполагаемый период динамики рассматриваемых технологических характеристик ГТС можно считать приближенно равным 48 часов. Соответствующие коэффициенты Д (/ = 1,...,16) определяются следующим образом:

л - о

Используя последовательность коэффициентов {Д}", были вычислены значения коэффициентов а0,...,аи и построено уравнение регрессии с первыми 16-ю гармониками ряда Фурье. Было показано, что для практического использования регрессионной модели в частичной сумме ряда Фурье построенного уравнения регрессии достаточно отставить четыре наиболее значимые гармоники. В качестве степени значимости некоторой гармоники предложено рассматривать ее коррелированность с исходными данными. Отобранным

четырем гармоникам ряда Фурье соответствует подпоследовательность коэффициентов |Д|4 , для которой пересчитываются коэффициенты результирующей регрессионной модели а0,...,а,. Полученные подобным образом уравнения регрессии с четырьмя гармониками используются для прогнозирования значений расхода и термобарических характеристик газотранспортной системы в последующие моменты времени.

Уравнение регрессии для датчика давления Р-1 ./(О = 2.959 + 0.00004-1 +

+ 0.0008-о«! —Л + О.ОООООб-вт (—+ 148 ) 148 )

+ 0.002-СО8^^ + О.ООЗ-мп^/^ + ^

+ О.ОООЗ-сов^-^-/^ + 0.0004-яп^^ +

+ 0.00004-сов^^ + 0.0003- эт^^.

Показана возможность улучшения прогноза динамики давления, температуры и расхода газа, полученного при помощи регрессионных моделей, за счет применения моделей авторегрессии - динамически пересчитываемыми моделями АРПСС( 1,1,0). Для этого проводится последовательная корректировка каждого прогнозного значения, вычисленного с помощью некоторого регрессионного уравнения /(/), от начала прогнозного периода и далее на каждом временном шаге прогнозирования. При этом корректировка значений уравнения регрессии на каждом временном шаге определяется по р предыдущим отклонениям регрессионной модели /(/) от прогнозируемого временного ряда и(0:

КО = "(0-/(0- (9)

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

Предложенный в работе метод корректировки прогноза обеспечивает получение более точных статистических оценок прогнозируемых показаний датчиков телеметрических станций по сравнению с теми, которые получаются с помощью экстраполяции уравнениями регрессии. Поскольку речь идет о статистических оценках, то не исключается, что в некоторых случаях прогнозные оценки из уравнения регрессии могут оказаться ближе к фактическим показаниям, чем их корректировка с помощью моделей авторегрессии (см. рисунок 3).

2.98

2.97

" 'А ♦ / \

» к \

♦ ♦

2.95

2.94

О

24

48

72

96

120

144

168

время(мае»

♦ Ретроспектныщешлчення — Регрессионная модель

- Начало прогноза

Корректировка авторегресспей

Рисунок 3. Динамика сглаженных ретроспективных и полученных на основании уравнения регрессии (8) прогнозных значений давления на датчике Р-1. Средняя погрешность прогноза: регрессия - 0.2%, авторегрессия - 0.13%

Исследования показали, что задача прогнозирования массового расхода газа вполне корректна только для датчиков, работающих в одном технологическом режиме. Переключение между рабочими режимами вызывает существенные скачки в показаниях практически всех датчиков, которые не всегда могут быть спрогнозированы при помощи моделей регрессии и авторегрессии. В этом случае было предложено использовать показания датчиков расхода в качестве управляющих параметров для нейросетевого прогнозирования значений на датчиках давления и температуры. При этом независимо от изменений технологических режимов показания температурных датчиков имеют незначительные и достаточно гладкие колебания, и для них можно построить вполне качественный прогноз, используя только модели регрессии и авторегрессии. Таким образом, для нейросетевого прогнозирования были выбраны датчики давления газа, для которых были характерны существенные скачки в ретроспективных данных, связанные с изменениями рабочего режима.

В работе рассмотрено два варианта построения прогноза: прогнозировались непосредственно показание датчика в следующий момент времени и изменение показания датчика на следующем временном шаге.

Была построена нейронная сеть, которая прогнозирует показания датчика давления на один шаг (час) вперед. В качестве входных параметров нейронной сети использовались ретроспективные (или уже спрогнозированные) показания датчика давления, а также ретроспективные и перспективные

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

Для прогнозирования предложено рассматривать «конусовидные» трехслойные персептроны. Общее правило определения количества нейронов по слоям (ЛГ - количество рассматриваемых датчиков расхода газа):

1-ый слой: «25^+24 (количество входов),

2-ой слой: «0.5-(25^+24),

3-ий слой: 1.

Обучение нейронной сети осуществляется при помощи усовершенствованной реализации алгоритма обратного распространения ошибки. В работе обосновано использование метода динамического изменения длины шага вдоль антиградиента в зависимости от скорости обучения нейросети как для каждого из обучающих примеров в отдельности, так и для всего обучающего множества. Для ускорения обучения или недопущения попадания в локальный минимум функционала ошибки нейронной сети применяется метод случайного смещения. Пример нейросетевого прогнозирования приведен на рисунке 4.

Для повышения качества получаемого прогноза были предложены следующие методики «дообучения»:

1) сеть обучается обычным образом на обучающем наборе, состоящем из двух подмножеств. Первое подмножество состоит приблизительно из 10% от исходных (основных) обучающих примеров, отобранных случайным образом, на которых нейросеть обучалась на этапе основного обучения. Данное подмножество обучающих примеров не позволит нейронной сети «забыть» основные закономерности поведения прогнозируемого временного ряда, которым нейросеть уже была обучена. Второе подмножество обучающих примеров сформировано по актуальным ретроспективным показаниям (почасовые показания за последние трое суток) выбранного для прогнозирования датчика, рассматриваемым в качестве целевых значений выхода нейросети.

2) набор обучающих примеров формируется только на основании актуальных ретроспективных показаний выбранного для прогнозирования датчика, рассматриваемых как целевые значения. При этом значение базовой длины шага обучения нейросети выбирается достаточно малым: р < 0.05. Небольшое значение шага обучения р способствует предотвращению возникновения проблемы переобучения нейронной сети на последних трех сутках.

4.8

4.6

3.2 1 ...

3 ---.-.-.---.-.-.-.-.-.-------.--

0 1 2 3 4 5 в 7 8 9 10 11 12 13 14 15 16

время I сутки I

■"Показания датчика -Предварительный прогноз Окоичателънийпрогно*

Рисунок 4. Нейросетевое прогнозирование показаний датчика давления Р-2 с

ежесуточным «дообучением» нейросети. Погрешность прогноза: предварительный - 3.52%, окончательный - 0.93%

В третьей главе рассматриваются различные аспекты долгосрочного прогнозирования спроса на газ. Известно, что прогноз спроса на газ является вероятностной величиной. В мировой практике принято определять динамику потребности газа с заданным уровнем достоверности. При этом применяют дискретные модели распределения вероятностей реализации проекта и сдвига его во времени.

По международной классификации существует четыре основные категории существующих и перспективных потребителей газа:

1) Базовый спрос - спрос существующих потребителей газа, который увеличивается в соответствии с ростом внутреннего валового продукта (ВВП);

2) Законтрактованный спрос представляет собой увеличение потребления газа уже подключенными клиентами или новыми потребителями газа, при условии наличия подписанных контрактов;

3) Возможный спрос - к категории возможных потребителей газа относятся те клиенты, которые находятся на стадии заключения договоров о поставках природного газа;

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

Дифференциация потребителей по категориям позволяет статистическими (метод Монте-Карло) или теоретико-вероятностными методами оценить пессимистичный (кривая р,0), оптимистичный (кривая р90) и наиболее

вероятный (кривая Р5„) уровни перспективной потребности в газе по годам (прогнозные уровни спроса на газ).

Кривая Р10, описывающая пессимистичный сценарий развития спроса на газ, соответствует квантилям уровня 0.1 в вероятностных пространствах суммарных объемов выполнения контрактных (проектно-контрактных) обязательств по годам перед всеми рассматриваемыми потребителями с учетом задержек поставок газа. То есть в каждый год ниже этой кривой лежат 10% возможных (вероятностных) значений суммарных объемов поставок газа всем потребителям.

Кривая Р50 описывает наиболее вероятный (реалистичный) сценарий развития спроса на газ и соответствует квантилям уровня 0.5 - медианам. Ниже этой кривой, а также, соответственно, и выше лежат 50% возможных значений суммарных объемов поставок газа рассматриваемым потребителям за каждый год.

Кривая Р90 описывает оптимистичный сценарий. В каждый год выше этой кривой лежат только 10% возможных значений суммарных объемов поставок газа (квантили уровня 0.9). То есть вероятность того, что спрос окажется выше указанной кривой, является сравнительно небольшой.

При построении данных кривых предполагается, что поставки газа различным потребителям являются попарно независимыми событиями. Также все вероятностные пространства возможных объемов поставок газа одному потребителю в разные годы рассматриваются как независимые. То есть, каждое предыдущее значение прогнозной кривой не влияет на вычисление последующего значения кривой. Таким образом, для построения прогнозных кривых Р10, Р50, Р,0 для произвольной группы потребителей предложено за каждый год прогнозирования некоторым образом вычислять или оценивать функцию распределения вероятностей суммарного объема поставки газа данной группе потребителей на основании вероятностных пространств выполнения контрактных обязательств перед каждым потребителем в отдельности.

Показано, что на основе метода Монте-Карло возможна оценка распределения вероятностей появления всех возможных исходов (суммарных объемов поставок газа потребителям) за каждый год прогнозирования, а, следовательно, построение кривых Р10, Р50, Р,0 для данной группы потребителей.

Исследования показали, что для многочисленной группы потребителей распределения вероятностей объемов поставок газа каждому потребителю в отдельности таковы, что можно предположить выполнимость условий центральной предельной теоремы. То есть, распределение суммарного объема поставки газа группе потребителей в каждый год / будет достаточно близко к распределению по нормальному закону . При этом данная степень

близости будет прямо пропорционально зависеть от количества потребителей в рассматриваемой группе. Для одного потребителя с вероятностным пространством объемов поставки газа V в некоторый год прогнозирования ма-

тематическое ожидание М и дисперсия О вычисляются следующим образом:

Л/ = £уР(У), ¿> = ^ур (у)-М2. (Ю)

Таким образом, вычислив параметры д и ст,2, можно оценить прогнозные значения Рг(/), как / -квантили нормального распределения /V(д,ст,2):

рг(0 « Л, Р»(0«Л. (И)

где гг - табличное значение, у -квантиль стандартного нормального распределения N(0,1).

Сравним результаты прогнозирования спроса на газ, полученные двумя различными подходами (см. рисунок 5).

Расчеты показали, что применение обоих подходов прогнозирования дает схожие результаты. Более того, близость полученных результатов прогноза, а также принципиальное различие и алгоритмическая независимость используемых подходов, гарантируют качество данных результатов и корректность подходов прогнозирования. Однако заметим, что для многочисленных групп потребителей теоретико-вероятностный подход при построении прогноза спроса на газ является менее ресурсоемким и требует значительно меньших вычислительных мощностей. В силу этого применение вероятностного подхода для подобных задач является более предпочтительным.

Были предложены алгоритмы, направленные на решение задачи корректировки прогнозных спросов на газ по отдельным клиентам, так, чтобы получаемые в итоге результаты прогнозов удовлетворяли условию аддитивности.

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

Для проектов поставок газа потребителям предложен характеризующий их набор параметров, по которым формируются базовые приоритеты реализации (в порядке значимости): (1) категория потребителя; (2) актуальность проекта; (3) прогрессивность сектора экономики; (4) прогрессивность региона', (5) запрашиваемые объемы газа. По данному набору параметров были определены соответствующие нечеткие переменные (семейства нечетких переменных), параметризованные по году прогнозирования /: (от,, ^{проект,/)), (а2,р2(проект)), {а3, ^(сектор,/)), {аА,р4 (регион^)), (а5, ¿¡¡(проект)). Данные переменные были использованы в экспертной системе для формирования базовых приоритетов реализации проектов поставок газа и их дальнейшей корректировки. Непосредственно приоритет реализации проектов определяется как некоторый функционал принадлежности нечеткого множества «Финальная реализация проекта», зависящий от данных нечетких переменных:

0, если (¡.(проект, /) = 0;

а, =1,

1, если ^(проект,!) = 1; (12)

> а.-^(проект, /), в противном случае; 4 ,

. '=1

где е! - заданная точность вычисления (дискретизация) функции принадлежности р., зависящая от необходимой степени детализации элементов нечеткого множества, описывающего семантику нечеткой переменной а;.

Введенные нечеткие переменные позволяют описать алгоритм формирования базовых приоритетов реализации (семейство функций принадлежности проект,/)) и выбора реализуемых проектов поставок газа потребителям в каждый год прогнозирования I. Пусть В(0 - подмножество проектов поставок газа, которые в год прогнозирования I относятся к категории «базовые» (были реализованы ранее). Алгоритм определения значений функционала р(проект,1) и последовательного построения множеств ВД/) на конец каждого года (, пробегающего значения от /0 - года начала периода прогнозирования до Т - года окончания периода показан на рисунке 6. Здесь индекс у отражает зависимость данных множеств от спрогнозированного суммарного спроса на газ на национальном уровне Р,(?). В результате функционирования алгоритма также строится кривая Ьг(0, которая определяет

уровень минимального приоритета для реализуемых в год / проектов поставок газа. То есть, для любого проекта е Вг(/) выполняется неравенство:

[¡(проект, () > ЬДс) • (13)

(Начало алгоритма

Конец алгоритма 1

В,('о) := {"роекты, "базовые" изначально}

Ф,(<0) := 0 Формируем Пг(/„). Определяем Ог00)

Ф(0 - подмножество проектов, которые в год 1 являются финально нереализованными;

П(/) — подмножество проектов, которые являются предстоящими в год

£1(0 = (в(0иФ«))\п(0

Для /:=<0 до Т

Для проектов из Вг(/), ФД/), ПД/) вычисляем //,(проект, <)

Для проектов с &г(1) вычисляем вточке (проект,1)

ДВДО := 0

Я := РДО

.(0

проект € В. (/)

Для всех проектов вычисляем [¡(проект, I)

Выберем ргеС1г(1), т.ч. М(рг, I) = шах {//(проект,/)}

проект е £1^,(1) 1

лвг(о = ЛВ,(/)и{рг},

пдо = п,(0 \{рг},

я = и- УргО),

Ргш = рг

Ц(0 := 1

М'> := Мрп^.0,

ВД0 := В,(0 и ДВД0

- Плановый объем | поставки газа по проекту в год 1

0.5 < с < 1 - некоторый фиксированный | коэффициент допустимости превышения ¡спрогнозированного суммарного спроса ¡на газ Р^/). предлагается брать с = 0.5

¡ДФ^О - появившиеся в год прогнозирования / !финально нереализованные проекты;

/\ДП7(0 - проекты, переставшие быть ' ¡предстоящими в год 1

ВД/ + 1)

ФДГ+1) = ФД0>~>ДФД/ + 1),

ПД/ + 1)

пдг+1) = ПД0\ЛФД' + 1)иДПД<+1)

Рисунок 6. Блок-схема алгоритма формирования базовых приоритетов реализации и выбора реализуемых проектов поставок газа потребителям

На основании сформированных приоритетов реализации проект, I) и кривой L(í) можно определить множества реализованных («базовых») проектов B(í). Таким образом, корректируя семейство функций принадлежности ¿¡(проект, t) и кривую L(f), эксперт может модифицировать изначально сформированные множества В(Г) и получить тем самым уточненный результат прогноза реализации проектов поставок газа потребителям.

В работе также рассмотрена задача прогнозирования спроса на основании ретроспективных значений. При этом при формировании прогноза учитывается, что некоторые значения в плане поставок газа потребителю уже известны (например, могут быть заданы основные предполагаемые вехи в развитии потребителя) и полученный прогноз должен гладко проходить через известные значения. Условие гладкости возникает из свойств любого естественного процесса, которым являются спрос на газ и газопотребление, не иметь существенных скачков и резких изломов.

Опишем математическую постановку данной задачи. Пусть Р - множество ретроспективных значений объемов поставок газа, П - множество известных перспективных (плановых) значений. Необходимо построить прогнозирующую функцию

F(t) = R(t) + Q( t), (14)

где R(t) - функция множественной регрессии, построенная по множеству точек РиП, Q{t) - корректирующая функция, которая является гладкой и для любой точки (/,у)еП обеспечивает равенство

F(t) = v. (15)

Функция множественной регрессии R(t) образуется суммированием линейного тренда RL(t) и гармонической составляющей R"(i):

RL(t) = c0+ct-t, R"(t) = a-cos{2Ktlp) + b%m(2ntlp). (16) Анализ показал, что функции множественной регрессии подобного вида составляют вполне оптимальный класс регрессионных функций с точки зрения долгосрочного прогнозирования.

Исходя из эмпирических соображений, период гармонической составляющей р должен удовлетворять следующим условиям: период должен быть кратен году и не может быть существенно больше глубины всех имеющихся данных (ретроспективных и перспективных), в противном случае его определение будет некорректной задачей. Отсюда вытекает, что

/,6N, 2<Р<У2 (max{/}-min{/}). (17)

Из практических ограничений на существующие ретроспективные данные и формируемый прогноз, параметр р по порядку не превосходит 102. Следовательно, вполне корректно определять р методом полного перебора.

В работе предложен ряд модификаций построения прогнозирующей функции, основанных, прежде всего, на экспертных заключениях и предположениях. Одна из модификаций заключается в выделении во множестве известных перспективных значений объемов поставок газа П некоторого под-

множества, которое не накладывает ограничения на прогнозирующую функцию в виде равенства (15). Выделенное подмножество представляет экспертные предположения относительно возможных перспективных объемов поставок газа и может содержать за один и тот же год несколько возможных объемов поставок газа - несколько различных экспертных предположений, например, выдвинутых разными экспертами. Следующая модификация построения прогнозирующей функции заключается в определении для каждой точки множества РиП весового коэффициента ю, > 0, характеризующего степень значимости рассматриваемой точки для построения прогноза. При этом неизвестные параметры прогнозирующей модели будут определяться методом взвешенных наименьших квадратов. Весовые коэффициенты для множества перспективных значений объемов поставок газа П определяются экспертно и отражают степень качества соответствующего экспертного предположения или заключения - элемента множества П. Как правило, в силу своей актуальности, точки множества П имеют вес ц > 1. В свою очередь точки множества Р имеют вес <у, < 1. Учитывая специфику рассматриваемой задачи, в частности, относительно небольшую возможную глубину ретроспективных данных, весовые коэффициенты со, множества точек Р предлагается определять по следующей формуле:

Отметим, что данная последовательность весовых коэффициентов является возрастающей и выпуклой вверх, то есть с увеличением срока давности значимость (весовые коэффициенты) более новых ретроспективных данных убывает медленнее, чем значимость старых данных.

На рисунке 7 представлен результат применения предложенного алгоритма с использованием описанных модификаций построения прогнозирующей функции. ___________ ___________

(18)

50

20

25

30

35

40

45

10

15

5

0

1980 1985 1990 1995 2000 2005

2010 2015 2020 2025 2030 ♦ Линейный тренд -♦-К 1но"Лественнля регрессия -♦-Прогно-шрующля '¡^нкцня

-•-Ретроспекгпвные'Шпчения • Ичвестные перспектнвные-шлчення О Возможные перспективные чнлчення

Рисунок 7. Построение прогнозирующих функций /?'(/), Я(!) и

В четвертой главе рассмотрены вопросы практического применения разработанных методов прогнозирования на примере газотранспортной системы Боливарианской Республики Венесуэлы. Общая протяженность ГТС Венесуэлы составляет 3695 км и включает в себя пять подсистем: Anaco -Jose - Puerto La Cruz, Anaco - Puerto Ordaz, Anaco - Caracas - Barquisimeto, La Toscana - San Vicente, Ulé - Amuay. В настоящее время в эксплуатации находится одна линейная компрессорная станция Altagracia (РСА) общей установленной мощностью 20 МВт. Замеры в узловых точках ГТС проводятся более чем на 2000 датчиках давления, температуры и расхода газа.

Разработанные подходы по прогнозированию динамики давления, температуры и массового расхода газа позволили провести вероятностную оценку работоспособности соответствующих датчиков и телеметрической аппаратуры. Для этого были выявлены датчики, для которых модели регрессии и авторегрессии дали слишком хороший прогноз (средняя погрешность прогноза равна или очень близка к 0% - характерен для датчиков, значения которых не менялись во времени за исключением, возможно, единичных выбросов) или слишком плохой прогноз (значение средней погрешности прогноза достаточно велико). Выявлено, что подобные прогнозы могут быть обусловлены следующими причинами: неисправность датчика или телеметрии, несоответствие технических параметров оборудования обрабатываемой информации (например, измеряемые значения значительно выходят за рамки рабочего диапазона датчика), в некоторые периоды времени замеры не производились, данные по телеметрии некорректно передавались на диспетчерский пункт.

Был проведен анализ реальных данных для 1052 датчиков давления, температуры и расхода газа за полугодовой период их функционирования. Для каждого датчика в автоматическом режиме с помощью разработанной компьютерной программы был построен ретроспективный прогноз на основании регрессионной модели и далее вычислено значение средней погрешности прогноза. Перед построением прогноза заведомо неверные показания датчиков экстраполировались предыдущими адекватными значениями. Границы адекватности показаний датчиков выбирались исходя из технологических ограничений на эксплуатацию газопровода:

- показания датчиков давления должны быть не меньше 1.37 МПа;

- показания датчиков температуры не меньше 283.148 К; показания датчиков расхода строго больше 0 кг/с.

Первоначально из дальнейшего анализа были исключены датчики с нулевой погрешностью построенного прогноза (слишком хороший прогноз). Для каждого из трех типов датчиков (давление, температура, расход) на основании полученной статистики погрешности прогноза была построена соответствующая функция плотности распределения погрешности. Далее был проведен экспертно-статистический анализ плотностей распределений погрешностей прогнозов и в результате по каждому типу датчика определено допустимое значение погрешности прогноза. Прогнозы показаний датчиков с

погрешностью, превосходящей допустимое значение, считаются слишком плохими прогнозами, и соответствующие датчики, а также связанную с ними телеметрию рекомендовали к дополнительной аппаратурной проверке корректности функционирования. Аналогично были рекомендованы для дополнительной проверки работоспособности датчики, по которым был построен слишком хороший прогноз (см. таблица 1). На рисунках 8 и 9 приведен пример выявления некорректно функционирующих датчиков давления.

70 •

1.1 ■ || 1

датчики давления

Рисунок 8. Статистика погрешности прогноза для датчиков давления

в 0.1

20 25 30 35 40 погрешность прогноза <Л и I

Рисунок 9. Плотность распределения погрешности прогноза для датчиков давления

Анализ результатов применения разработанной вероятностной методики оценки корректности функционирования датчиков телеметрических станций показал, что для получения наиболее достоверных оценок необходимо применять данную методику в периоды работы ГТС в одном технологическом режиме. Как отмечалось выше, в случае наличия изменений технологических режимов модели регрессии и авторегрессии могут дать не вполне точный прогноз. Наиболее чувствительны к смене режимов работы ГТС датчики расхода газа. Тем не менее, для всех проанализированных датчиков ГТС Венесуэлы были получены вполне приемлемые результаты по оценке работоспособности (см. таблицу 1).

Таблица 1. Сводные данные по оценке работоспособности датчиков

Датчики давления Датчики температуры Датчики расхода

Общее количество датчиков 675 236 141

Количество исключенных из анализа датчиков 95 37 62

Количество проанализированных датчиков 580 199 79

Допустимое значение погрешности прогноза 12% 1.5% 25%

Количество датчиков рекомендованных к проверке 28 7 19

Доля датчиков рекомендованных к проверке 4.8% 3.5% 24%

В работе подробно рассмотрен пример практического применения разработанных методов оценки прогнозного спроса на газ в условиях неопределенности в рамках формирования программы развития и реконструкции газотранспортной системы Республики Венесуэлы.

На рисунке 10 представлена перспективная схема газотранспортной системы Венесуэлы, на которой отмечены пять основных потоковых узлов: Guiria, Puerto Ordaz, Barcelona, Morichal, Maracaibo. К 2030 году в данных узлах ГТС должны обеспечиваться поставками газа наиболее крупные потребители. Динамика максимальной потребности в газе по основным потоковым

-•"Cruiru -»-Puerto Ordaz -♦-Barcelona -«-Morich.il -'-Maracaibo

Рисунок 11. Динамика максимальной потребности в газе по узлам ГТС Венесуэлы

Полученные результаты прогнозов спроса на газ в узлах ГТС были использованы для последующего расчета перспективных потоков газа и объемов развития и реконструкции магистральных газопроводов. На основании данных по текущей производительности и динамики потоков по участкам ГТС Венесуэлы в период до 2030 года были оценены необходимые вводы мощностей по годам. В частности, на суше необходимо сооружение 3980 км газопроводов различных диаметров и 8 КС с 32 ГПА общей мощностью 502 МВт. Суммарные объемы ввода морских газопроводов оцениваются в 1330 км.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ РАБОТЫ

1. Разработана регрессионная модель, описывающая динамику поведения технологических параметров ГТС, и ее корректировка моделями авторегрессии для прогноза на 24 часа.

2. Обосновано применение нейронных сетей для прогнозирования технологических параметров ГТС в случае наличия изменений технологических режимов работы ГТС.

3. Разработана вероятностная методика оценки работоспособности датчиков телеметрических станций ГТС.

4. Выполнено построение прогнозных уровней спроса на газ (пессимистичный, наиболее вероятный, оптимистичный) на основе эффективного применения международной классификации существующих и перспективных потребителей газа.

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

6. Разработан алгоритм экспертной многопараметрической характеристики проектов поставок газа для определения приоритетов их реализации, согла-

сованных с прогнозным уровнем спроса на газ по стране в целом, на основании методов нечеткой логики.

7. Разработанные методы прогнозирования технологических параметров магистральных газопроводов успешно применены для повышения эффективности функционирования и устойчивого развития газотранспортной системы Республики Венесуэлы в условиях неопределенности спроса на газ.

Основные публикации по теме диссертации:

1. Самсонов P.O. Анализ перспектив развития экспорта российского природного газа / P.O. Самсонов, H.A. Кисленко, К.А. Казак, A.C. Казак // Нефть, газ и бизнес. - 2007. - № 1-2. - С. 9-14.

2. Казак A.C. Применения экспертных систем для формирования прогнозов спроса на газ / A.C. Казак, H.A. Кисленко, К.А. Казак, В.Н. Башкин // Нефть, газ и бизнес. - 2007. - № 10. - С. 67-72.

3. Казак К.А. Корректировка диспетчерских данных методом частичного сглаживания / К.А. Казак, В.В. Кнотько, H.A. Кисленко, A.C. Казак // Системы управления и информационные технологии. - 2007. - № 4.1(30). -С. 152-155.

4. Казак К.А. Регрессионный анализ технологических характеристик ГТС / К.А. Казак, В.В. Кнотько, H.A. Кисленко, A.C. Казак // Системы управления и информационные технологии. - 2007. - № 4.1(30). - С. 155-159.

5. Казак К. А. Вероятностная методика оценки работоспособности датчиков/ К.А. Казак, В.В. Кнотько, H.A. Кисленко, A.C. Казак // Системы управления и информационные технологии. - 2007. - № 4.2(30). - С. 249-251.

6. Казак К.А. Нейросетевое прогнозирование параметров ГТС / К.А. Казак, В.В. Кнотько, H.A. Кисленко, A.C. Казак // Системы управления и информационные технологии. - 2008. - № 1 (31). - С. 56-60.

7. Казак К.А. Стохастические модели долгосрочного прогнозирования спроса на газ при наличии больших выборок потребителей / К.А. Казак, A.C. Казак, H.A. Кисленко, И.С. Гитман, КБ. Кудрявцев // Системы управления и информационные технологии. - 2008. - № 1.3(31). - С. 363-366.

8. Казак К.А. Формирование приоритетов реализации проектов поставок газа потребителям на основании методов нечеткой логики / К.А. Казак, И.С. Гитман, A.C. Казак, H.A. Кисленко // Системы управления и информационные технологии. - 2008. - № 2.2(32). - С. 256-261.

9. Казак К.А. Системный анализ потребления газа // Наука и техника в газовой промышленности. - 2010. - № 2. - С. 19-23.

10. Казак К.А. Оценка адекватности вариантов развития газовой отрасли // Наука и техника в газовой промышленности. - 2010. - № 2. - С. 105-110.

11. Чионов A.M. Система теплогидравлического моделирования режимов работы магистральных газопроводов при различных способах их прокладки / A.M. Чионов, И.О. Стурейко, К.А. Казак // Тезисы докладов III Международной научно-технической конференции «Газотранспортные системы: настоящее и будущее (GTS-2009)» - М.: Газпром ВНИИГАЗ, 2009.-С. 72.

12. Samsonov R.O. Master plan methodology for gas industry development/ R.O. Samsonov, A.S. Kazak, V.N. Bashkin, K.A. Kazak et al. - Moscow: Scientific World, 2007. - 304 p.

Подписано к печати 29.04.2010 г.

Заказ № 5845 Объем 1,75 уч.-изд.л. Ф-т 60x84/16. Тираж 100 экз.

Отпечатано в ООО «Газпром ВНИИГАЗ», 142717, Московская область, Ленинский р-н, п. Развилка, ООО «Газпром ВНИИГАЗ»

Содержание диссертации, кандидата технических наук, Казак, Константин Александрович

ВВЕДЕНИЕ.

ГЛАВА 1. ОБЗОР ПОДХОДОВ И МЕТОДОВ

ПОСТРОЕНИЯ ПРОГНОЗОВ.

1.1. Описание объектов исследования.

1.2. Общая постановка задачи прогнозирования.

1.3. Анализ, отбор и корректировка входных параметров для функции прогнозирования.

1.4. Экстраполяция (прогнозирование) статистическими методами.

1.5. Нейросетевое построение прогнозирующей функции.

1.6. Экспертные системы и элементы нечеткой логики в прогнозировании.

ГЛАВА 2. КРАТКОСРОЧНОЕ ПРОГНОЗИРОВАНИЕ ТЕХНОЛОГИЧЕСКИХ ПАРАМЕТРОВ

ВНЕШНЕГО ВОЗДЕЙСТВИЯ НА ГТС.

2.1. Постановка задачи.

2.2. Корректировка диспетчерских данных методом частичного сглаживания.

2.3. Регрессионный анализ технологических характеристик ГТС.

2.4. Корректировка прогноза с помощью моделей авторегрессии.

2.5. Нейросетевое прогнозирование технологических характеристик газового потока.

ГЛАВА 3. ДОЛГОСРОЧНОЕ ПРОГНОЗИРОВАНИЕ

СПРОСА НА ГАЗ.

3.1. Постановка задачи.

3.2. Долгосрочное прогнозирование спроса на газ при наличии больших выборок потребителей.

3.3. Формирование приоритетов реализации проектов поставок газа потребителям на основании методов нечеткой логики.

3.4. Долгосрочное прогнозирование спроса на газ для «базовых» потребителей.

ГЛАВА 4. ПРАКТИЧЕСКОЕ ПРИМЕНЕНИЕ

РАЗРАБОТАННЫХ МЕТОДОВ ПРОГНОЗИРОВАНИЯ ТЕХНОЛОГИЧЕСКИХ ХАРАКТЕРИСТИК В УЗЛАХ ГАЗОТРАНСПОРТНОЙ СИСТЕМЫ.

4.1. Структура ГТС Республики Венесуэлы.

4.2. Вероятностная методика оценки работоспособности датчиков

ГТС Венесуэлы на основе анализа погрешности прогноза.

4.3. Анализ неопределенности прогноза спроса на газ для оценки объемов развития и реконструкции ГТС Венесуэлы.

ВЫВОДЫ.

Введение Диссертация по наукам о земле, на тему "Разработка методов прогнозирования технологических режимов магистральных газопроводов в условиях неопределенности спроса на газ"

В настоящее время в связи со строительством магистральных газопроводов рассчитанных на высокое давление, в том числе морских газопроводов («Голубой поток», «Северный поток», «Южный поток») возрастает необходимость разработки более точных технологических моделей магистрального транспорта газа. Это относится и к прогнозированию граничных условий при расчете режимных показателей газопровода (температура, давление и расход газа). Так, например, при относительном отклонении давления в начальном узле на 1% значение расхода в конечном узле может иметь отклонение до 5%. То есть, незначительная ошибка прогноза в результате газодинамического расчета может увеличиться на порядок. Таким образом, должны быть разработаны методы идентификации и адаптации значений давления, температуры и расхода газа в граничных узлах моделируемого участка газопровода, что позволило бы получать прогнозируемые режимные показатели в пределах точности измерительных приборов.

Основной задачей системы магистрального транспорта газа является надежность обеспечения поставок газа потребителям. Следует отметить, что потребители природного газа, как правило, привязаны к узлам газотранспортной системы (ГТС). В зависимости от объемов потребления распределение газа в узлах ГТС может осуществляться непосредственно на промышленные объекты или через газораспределительные станции. При формировании планов развития и реконструкции ГТС необходимо исходить из достоверных прогнозов спроса на газ в этих узлах, чтобы наиболее точно определять динамику развития мощностей газопроводов. В кризисных условиях это особенно важно для минимизации неиспользуемых резервов ГТС и обеспечения надежности поставок газа потребителям в соответствии с их требованиями. В связи с этим должны быть разработаны методы прогнозирования потребности в газе, позволяющие адаптировать развитие системы магистрального транспорта газа к текущей ситуации в газовой промышленности.

Таким образом, в силу указанных выше причин возникает необходимость разработки новых, а также развития и адаптации существующих моделей и методов прогнозирования в задачах функционирования и развития систем магистрального транспорта газа.

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

В диссертационной работе были исследованы следующие основные задачи:

• разработка моделей краткосрочного прогнозирования диспетчерских данных в условиях неизменяющегося технологического режима работы ГТС, при ограничениях на время оценки;

• анализ результатов прогноза и разработка на его основе вероятностной методики оценки работоспособности датчиков и соответствующей телеметрии; разработка адаптивных моделей краткосрочного прогнозирования диспетчерских данных, зависящих от управляющих параметров, в условиях изменяющегося технологического режима работы ГТС;

• обоснование моделей долгосрочного прогнозирования спроса на газ с учетом различных экспертных оценок для потребителей, имеющих контракты на поставку газа; анализ применимости центральной предельной теоремы для определения различных уровней спроса на газ (пессимистичный, наиболее вероятный, оптимистичный) в задачах долгосрочного прогнозирования спроса при наличии больших выборок потребителей с применением международной классификации существующих и перспективных потребителей по категориям спроса;

• анализ подходов и разработка методики экспертной многопараметрической характеристики проектов поставок газа для численной оценки приоритетов реализации проектов в динамике.

В работе впервые проведены комплексные исследования по разработке обобщенных подходов к задачам краткосрочного и долгосрочного прогнозирования характеристик узловых точек газотранспортных систем с применением современного математического аппарата. Наряду с традиционными методами прогнозирования (регрессионный анализ; модели авторегрессионного проинтегрированного скользящего среднего; статистические, теоретико-множественные и теоретико-вероятностные подходы) используются нейросетевые алгоритмы и методы нечеткой логики. При решении задач долгосрочного прогнозирования впервые применена международная классификация потребителей по категориям спроса в зависимости от оценок вероятности реализации различных проектов использования газа. Построена универсальная экспертная многопараметрическая характеристика проектов поставок газа для формирования приоритетов их реализации в виде набора нечетких переменных как составной части экспертной системы.

Разработанные алгоритмы краткосрочного прогнозирования могут быть использованы для целей оперативного управления работой газотранспортных систем и согласования технологических характеристик в узловых точках магистральных газопроводов, а также в составе задачи по газодинамическому моделированию магистрального транспорта газа. Статистический анализ краткосрочных прогнозов позволяет выявлять датчики и соответствующие им станции телеметрии, функционирующие с большой вероятностью некорректно и требующие проверки работоспособности.

Представленные в работе модели долгосрочного прогнозирования позволяют определять вероятностную динамику перспективного спроса на газ на возможных объектах его потребления сроком на десятки лет. Результаты долгосрочного прогнозирования спроса на газ предназначены для формирования задач развития и реконструкции ГТС на перспективу. Использование разработанных методических подходов и алгоритмов позволяет формировать пессимистичные, реалистичные и оптимистичные варианты спроса на газ с тем, чтобы определить возможный диапазон количественных оценок стратегии развития магистрального транспорта газа, а также ранжировать имеющиеся проекты поставок газа по стратегическому приоритету реализации.

Результаты диссертационного исследования были использованы при разработке Генеральной схемы развития газовой промышленности Боливарианской Республики Венесуэлы, а также при разработке математических моделей в составе программного комплекса по газодинамическому моделированию магистрального транспорта газа.

В соответствии с основными задачами исследования, представленными выше, в диссертационной работе защищаются следующие положения:

• методика краткосрочного прогнозирования показаний датчиков телеметрических станций ГТС (диспетчерских данных), состоящая из: корректирующей фильтрации данных; построения модели регрессии, представляемой в виде линейного тренда и периодической составляющей; корректировки прогноза динамически пересчитываемыми моделями авторегрессии;

• методика применения нейросетевого прогнозирования с использованием усовершенствованной реализации алгоритма обратного распространения ошибки для краткосрочного прогнозирования показаний датчиков, зависящих от неявных управляющих параметров, в период изменения технологических режимов работы ГТС;

• вероятностная методика оценки корректности функционирования датчиков телеметрических станций на основе анализа статистики погрешности регрессионных прогнозов с последующей рекомендацией к проверке работоспособности;

• методика долгосрочного прогнозирования спроса на газ по отдельным потребителям с применением регрессионных моделей, эмпирически оптимизированным выявлением периодической составляющей, построением гладкой корректирующей функции в виде кубического сплайна для опорных значений, с возможностью учесть экспертные предположения о будущих значениях и с возможностью взвешивать известные значения;

• методики построения различных уровней спроса на газ (пессимистичный, наиболее вероятный, оптимистичный) в задачах долгосрочного прогнозирования спроса при наличии больших выборок потребителей с применением международной классификации существующих и перспективных потребителей по категориям спроса;

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

Материалы диссертации опубликованы в 22-х работах, из них 18 работ опубликовано в журналах, включенных в Перечень ведущих рецензируемых научных журналов и изданий ВАК Минобрнауки РФ, в том числе опубликована одна монография. Также материалы исследований по теме диссертационной работы докладывались и обсуждались на секциях Ученого совета ООО «Газпром ВНИИГАЗ» и различных конференциях:

III Международная научно-техническая конференция

Газотранспортные системы: настоящее и будущее (GTS-2009)» (Москва, 2009 г.);

III Всероссийская научно-техническая конференция «Безопасность критичных инфраструктур и территорий» (Екатеринбург, 2009 г.).

Диссертационная работа состоит из введения, четырех глав, выводов и списка литературы.

Заключение Диссертация по теме "Строительство и эксплуатация нефтегазоводов, баз и хранилищ", Казак, Константин Александрович

выводы

В ходе выполнения диссертационной работы получены следующие основные результаты:

1. Разработана регрессионная модель, описывающая динамику поведения технологических параметров ГТС, и ее корректировка моделями авторегрессии для прогноза на 24 часа.

2. Обосновано применение нейронных сетей для прогнозирования технологических параметров ГТС в случае наличия изменений технологических режимов работы ГТС.

3. Разработана вероятностная методика оценки работоспособности датчиков телеметрических станций ГТС.

4. Выполнено построение прогнозных уровней спроса на газ (пессимистичный, наиболее вероятный, оптимистичный) на основе эффективного применения международной классификации существующих и перспективных потребителей газа.

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

6. Разработан алгоритм экспертной многопараметрической характеристики проектов поставок газа для определения приоритетов их реализации, согласованных с прогнозным уровнем спроса на газ по стране в целом, на основании методов нечеткой логики.

7. Разработанные методы прогнозирования технологических параметров магистральных газопроводов успешно применены для повышения эффективности функционирования и устойчивого развития газотранспортной системы Республики Венесуэлы в условиях неопределенности спроса на газ.

Библиография Диссертация по наукам о земле, кандидата технических наук, Казак, Константин Александрович, Москва

1. Ширяев А.Н. Вероятность. В 2-х кн. М.: МЦНМО, 2007. - 552+416 с.

2. Бокс Дж., Дженкинс Г. Анализ временных рядов. Прогноз и управление. Вып. 1,2.-М.: Мир, 1974.-405+197 с.

3. Крамер Г. Математические методы статистики. М.: Мир, 1975. - 648 с.

4. Тутубалин В.Н. Теория вероятностей и случайных процессов. М.: Изд-во МГУ, 1992.-400 с.

5. Ван дер Варден Б.Л. Математическая статистика. М.: Изд-во иностранной литературы, 1960. - 435 с.

6. Андерсон Т. Статистический анализ временных рядов. М.: Мир, 1976. -757 с.

7. Дженкинс Г., Ватте Д. Спектральный анализ и его приложения. -М.: Мир. Вып. 1, 1971. Вып. 2, 1972.-318+288 с.

8. Эконометрика. Учебник / Под ред. Елисеевой И.И. М.: Финансы и статистика, 2006. — 576 с.

9. Тутубалин В.Н. Границы применимости (вероятностно-статистические методы и их возможности). М.: Знание, 1977. - 64 с.

10. Тюрин Ю.Н., Макаров А.А. Статистический анализ данных на компьютере. М.: ИНФРА-М, 1998. - 528 с.

11. Тихонов А.Н., Уфимцев М.В. Статистическая обработка результатов эксперимента. М.: Изд-во МГУ, 1988. - 172 с.

12. Чураков Е.П. Математические методы обработки экспериментальных данных в экономике. Учебное пособие. М.: Финансы и статистика, 2004. -240 с.

13. Большев Л.Н., Смирнов Н.В. Таблицы математической статистики. -М.: Наука, 1983.-416 с.

14. Архипов Г.И., Садовничий В. А., Чубариков В.Н. Лекции по математическому анализу. М.: Высшая школа, 1999. - 695 с.

15. Гельфанд И.М. Лекции по линейной алгебре. М.: Добросвет, МЦНМО, 1998.-320 с.

16. Бахвалов Н.С., Жидков Н.П., Кобельков Г.М. Численные методы. -М.: БИНОМ. Лаборатория знаний, 2008. 636 с.

17. Соболь И.М. Численные методы Монте-Карло. М.: Наука, 1973. - 311 с.

18. Заде JT.А. Понятие лингвистической переменной и его применение к принятию приближенных решений. — М.: Мир, 1976. — 165 с.

19. Рыжов А.П. Элементы теории нечетких множеств и измерения нечеткости. М.: Диалог-МГУ, 1998. - 82 с.

20. Саати Т. Принятие решений. Метод анализа иерархий. М.: Радио и связь, 1993.-320 с.

21. Саати Т., Керне К. Аналитическое планирование. Организация систем. -М.: Радио и связь, 1991.-224 с.

22. Уоссермен Ф. Нейрокомпьютерная техника: Теория и практика. -М.: Мир, 1992.- 184 с.

23. Горбань А.Н., Россиев Д.А. Нейронные сети на персональном компьютере. Новосибирск: Наука, 1996. - 276 с.

24. Weigend A.S., Huberman В.А., Rumelhart D.E. Predicting the Future: A Connectionist Approach // International Journal of Neural Systems. 1990. Vol. l.-P. 193-209.

25. Круглов B.B., Борисов B.B. Искусственные нейронные сети. Теория и практика. М.: Горячая линия-Телеком, 2001. - 382 с.

26. Крисилов В.А., Чумичкин К.В. Ускорение обучения нейронных сетей за счет адаптивного упрощения обучающей выборки // Одесский национальный политехнический университет. Институт компьютерных систем.

27. Крисилов В.А., Кондратюк А.В. Преобразование входных данных нейросети с целью улучшения их различимости // Одесский национальный политехнический университет. Институт компьютерных систем.

28. Агаев Н.Б. Нейросетевые модели объема газопотребления в распределительных системах газоснабжения // Газовая промышленность. -2006. №2.-С. 37-40.

29. Сухарев М.Г. Методы прогнозирования. Учебное пособие. М.: РГУ нефти и газа, 2009. - 208 с.

30. Березина И.В., Ретинский B.C. Оперативное управление системами газоснабжения. М.: Недра, 1985. - 192 с.

31. Галиуллин З.Т., Леонтьев Е.В. Интенсификация магистрального транспорта газа. М.: Недра, 1991. - 272 с.

32. Сухарев М.Г., Ставровский Е.Р. Оптимизация систем транспорта газа. -М.: Недра, 1975.-277 с.

33. Сухарев М.Г., Горлов Д.В. Новая методология моделирования единой системы газоснабжения // Доклады РАН. Энергетика. 2002. № 4.

34. Казак А.С., Седов В.И., Орехова И.В., Яковлев Е.И. Оперативный контроль магистральных газопроводов. М.: Недра, 1989. - 289 с.

35. Черняев В.Д., Яковлев Е.И., Казак А.С., Сощенко А.Е. Трубопроводный транспорт углеводородного сырья. — М.: ВНИИОЭНГ, 1991. 344 с.

36. Сарданашвили С.А. Расчетные методы и алгоритмы (трубопроводный транспорт газа). М.: Нефть и газ, 2005. - 576 с.

37. Кисленко Н.А., Кудрявцев А.А., Постников В.Д., Гришко В.В., Анохина О.С. Анализ эффективного спроса на газ в рамках программ газификации // Газовая промышленность. 2004. № 8.

38. Самсонов P.O., Кисленко Н.А., Казак К.А., Казак А.С. Анализ перспектив развития экспорта российского природного газа // Нефть, газ и бизнес. -2007. № 1-2.-С. 9-14.

39. Самсонов P.O., Кисленко Н.А., Казак К.А., Казак А.С. Формирование краткосрочных и долгосрочных прогнозов спроса на газ в России // Нефть, газ и бизнес. 2007. № 3. - С. 20-26.

40. Казак А.С., Кисленко Н.А., Казак К.А., Башкин В.Н. Применения экспертных систем для формирования прогнозов спроса на газ // Нефть, газ и бизнес. 2007. № 10. - С. 67-72.

41. Казак К.А., Казак А.С. Разработка метода расчета неустановившихся режимов транспорта газа по участку трубопровода при возникновении утечки // Системы управления и информационные технологии. 2007. № 2.2(28). - С. 237-240.

42. Казак К.А., Казак А.С. Моделирование нестационарных режимов газопроводов с граничными условиями по давлению и температуре в начале участка и массовому расходу на конце // Системы управления и информационные технологии. 2007. № 3(29). - С. 8-10.

43. Казак К.А., Казак А.С., Башкин В.Н. Оценка адекватности вариантов развития газовой отрасли // Системы управления и информационные технологии. 2007. № 3.1(29). - С. 146-149.

44. Казак К.А., Казак А.С. Идентификация мест утечки в линейной части магистральных газопроводов на основе неустановившихся изотермических моделей // Системы управления и информационные технологии. 2007. № 3.2(29). - С. 250-254.

45. Казак А.С., Казак К.А. Идентификация местоположения утечек на основе неизотермических моделей течения газа // Системы управления и информационные технологии. 2007. № 3.2(29). — С. 254-258.

46. Samsonov R.O., Kazak A.S., Bashkin V.N., Kazak К.А. et al. Master plan methodology for gas industry development. Moscow: Scientific World, 2007. -304 p.

47. Казак K.A., Кнотько B.B., Кисленко H.A., Казак А.С. Корректировка диспетчерских данных методом частичного сглаживания // Системы управления и информационные технологии. 2007. №4.1(30). - С. 152155.

48. Казак К.А., Кнотько В.В., Кисленко Н.А., Казак А.С. Регрессионный анализ технологических характеристик ГТС // Системы управления и информационные технологии. 2007. № 4.1(30). - С. 155-159.

49. Казак К.А., Кнотько В.В., Кисленко Н.А., Казак А.С. Применение авторегрессионной модели для улучшения качества прогноза // Системы управления и информационные технологии. — 2007. № 4.2(30). — С. 247249.

50. Казак К.А., Кнотько В.В., Кисленко Н.А., Казак А.С. Вероятностная методика оценки работоспособности датчиков // Системы управления и информационные технологии. 2007. № 4.2(30). - С. 249-251.

51. Казак К.А., Кнотько В.В., Кисленко Н.А., Казак А.С. Нейросетевое прогнозирование параметров ГТС // Системы управления и информационные технологии. 2008. № 1(31). - С. 56-60.

52. Гитман И.С., Казак К.А., Кисленко Н.А., Казак А.С., Кудрявцев И.Б. Долгосрочное прогнозирование спроса на газ при наличии больших выборок потребителей // Системы управления и информационные технологии. 2008. № 1.3(31). - С. 349-352.

53. Казак К.А., Гитман И.С., Казак А.С., Кисленко Н.А. Формирование приоритетов реализации проектов поставок газа потребителям на основании методов нечеткой логики // Системы управления и информационные технологии. 2008. № 2.2(32). - С. 256-261.

54. Казак К.А. Системный анализ потребления газа // Наука и техника в газовой промышленности. 2010. № 2. - С. 19-23.

55. Казак К.А. Оценка адекватности вариантов развития газовой отрасли // Наука и техника в газовой промышленности. — 2010. № 2. — С. 105-110.