УДК 519.833.2 ББК 22.18
ВЫЧИСЛЕНИЕ ХАРАКТЕРИСТИК ПАССАЖИРОПОТОКОВ В ТРАНСПОРТНЫХ СИСТЕМАХ1
Буре В. М.2
(Санкт-Петербургский государственный университет, Санкт-Петербург) Мазалов В. В.3 (Институт прикладных математических исследований Карельского научного центра РАН, Петрозаводск) Плаксина Н. В.4 (Петрозаводский государственный университет, Петрозаводск)
Статья посвящена исследованию статистической модели оценки интенсивности пассажиропотоков между остановочными пунктами. Предполагается, что при этом известны данные о пассажирах, входящих в общественный транспорт и выходящих из него. В статье предлагается методика оценки распределения входящих на каждой остановке пассажиров по последующим остановкам маршрута. Приведены данные численного моделирования на основе натурных экспериментов.
Ключевые слова: теория игр, пассажиропоток, транспортная система, натурный эксперимент, равновесие.
1 Работа выполнена при финансовой поддержке РФФИ (грант № 13-01-91158-ГФЕН), Программы ОМНРАН и Программы стратегического развития ПетрГУ.
2 Владимир Мансурович Буре, доктор технических наук, профессор (vlb310154@gmail.com).
3 Владимир Викторович Мазалов, доктор физико-математических наук, профессор (vmazalov@krc.karelia.ru).
4Нина Владимировна Плаксина, аспирант (plaksina_nv@mail.ru).
Управление большими системами. Выпуск 47 Введение
В настоящее время в связи со строительством новых микрорайонов города, увеличением количества личного транспорта актуальной становится проблема анализа и оценки пассажиропотоков. Под пассажиропотоком будем понимать движение пассажиров в одном направлении маршрута. Пассажиропоток может быть в прямом направлении и в обратном. Пассажиропотоки играют важную роль при планировании маршрутов движения городского пассажирского транспорта, определении оптимальных интервалов движения городского пассажирского транспорта. Особенностью пассажиропотоков является их неравномерность, т.е. изменчивость по времени (по часам, суткам, дням недели, сезонам года). Пассажирообразующая способность отдельного района определяется в зависимости от количества населения, показателей его подвижности и коэффициентов неравномерности перевозок по времени.
Для получения информации о показателях подвижности населения и коэффициентах неравномерности перевозок применяют различные методы обследования. Натурные методы обследования обладают большой точностью, погрешность таких методов составляет примерно 5% [5]. Главный недостаток таких методов - большие затраты денежных и людских ресурсов для их проведения. Кроме того, на обработку результатов необходимо потратить много времени, поэтому результаты наблюдения могут уже неточно отражать реальную информацию о пассажиропотоках. Для определения пассажиропотоков также применяются энтропийные методы [1]. В энтропийной модели пассажиропотоки вычисляются на основе демографических и социально-экономических данных с использованием информации о территориальном расположении района. Главным недостатком энтропийного метода является допущение стационарности пассажиропотоков, т.е. такая модель не учитывает изменчивость пассажиропотоков от времени суток, дней недели и т.д. Поэтому нередко в совокупности с натурными экспериментами, демографическими
данными используют математические методы для получения актуальной информации о пассажиропотоках [2, 4, 6, 7, 11]. Такие методы также целесообразно применять для определения оптимальных (равновесных) решений при планировании перевозок, используя полученную информацию о пассажиропотоках. Значения численных характеристик пассажиропотоков позволяют решать различные прикладные задачи оптимального управления городским транспортом, такие как необходимое количество автобусов, достаточное для перевозки пассажиров и с учетом интересов пассажиров, а также интересов перевозчиков (чтобы рейсы не были убыточными).
1. Модель
В общем случае транспортная модель может быть представлена графом, вершины которого являются остановками, а ребра описывают транспортные коридоры. Выберем какой-нибудь маршрут транспортного средства, который представляет собой последовательность остановок, соединенных ребрами. Рассмотрим случай, когда на маршруте имеется К остановок. Предположим, что есть потоки пассажиров между этими остановками в прямом направлении. Под направлением будем понимать маршрут следования автобуса. Задача состоит в том, чтобы определить доли пассажиров из общего пассажиропотока и направления движения, по которым перемещаются эти доли. Например, для десяти пассажиров, находящихся на остановке, какое-то количество едет до конечной остановки, а какое-то количество выйдет на промежуточных остановках. Для этого проведем серию экспериментов г, г = 1,...,Ж. В качестве одного эксперимента будем рассматривать одну поездку автобуса от начальной остановки до конечной. Цель эксперимента состоит в фиксации информации о количестве вошедших и вышедших из автобуса на каждой остановке пассажиров. Пусть ж\, ж£, ..., жгк - количество пассажиров, вошедших на остановках г = 1,..., К соответственно в эксперименте г, г = 1,..., N. Пусть уГ, У2, УК - количество пассажиров, вышедших соответственно на остановках ] = 1 ,...,К
79
в эксперименте г, г = 1,...,Ж. Очевидно, что хгк = у\ = 0 в
Г, Уг
любом эксперименте г. В модели , уГ - это наблюдаемые ве
личины. Обозначим через у- - количество пассажиров, которые вошли на остановке i и вышли на остановке у. Это ненаблюдаемые величины, которые и нужно оценить. Тогда информацию о перемещениях пассажиров удобно представить в виде таблицы (таблица 1).
Таблица 1. Информация о количестве вошедших и вышедших из автобуса пассажиров
уи = 0 у12 со 1—1 у у1К -1 у1к х 1
0 уг>2 = 0 у со у2К-2 у2 к х2
0 0 0 ук-1,к-1 = 0 ук-1 ,к хК-1
0 0 0 0 0 хк = 0
уг = 0 у22 ^ СО у у к-1 уК
В таблице 1 приведены результаты эксперимента с номером г. В последнем столбце все элементы равны суммам элементов таблицы из соответствующей строки, и аналогично в последней строке все числа равны сумме чисел в соответствующем столбце; у- - количество пассажиров (эти величины в эксперименте не наблюдаются), вошедших в автобус на остановке с номером i и вышедших на остановке с номером у в эксперименте с номером г:
(1) уГ = Е УГ.
г=1
В эксперименте наблюдаются только элементы последнего столбца и последней строки. Будем считать для всех г, что пассажиропотоки уГ; представляют собой независимые (для разных г), одинаково распределенные случайные величины. Их среднее значение представляет собой интенсивность движения пассажиров по маршруту от остановки г до остановки ]. Обозначим через Р; долю пассажиров, вошедших на остановке г и вышедших на остановке ее можно трактовать как вероятность того, что пасса-80
жир, вошедший на остановке г, выйдет на остановке Составим из этих значений таблицу 2.
Таблица 2. Информация о вероятностном распределении пассажиров
Р11 = 0 Р12 Р13 Р1К-1 Р1К 1
0 Р22 = 0 Р23 Р2К-1 Р2К 1
0 0 0 РК-1,К-1 = 0 РК-1,К 1
В таблице 2 сумма элементов в каждой строке равна единице и представляет собой полиномиальное распределение [3] вероятностей выхода пассажиров, вошедших в автобус на остановке, при этом номер остановки совпадает с номером строки.
Рассмотрим эксперимент с номером г, умножим строку с номером г из таблицы 2 на число хГ, тогда естественно ожидать, что получившаяся после умножения строка соответствует строке с номером г из таблицы 1, так как х[рз представляет собой ожидаемое значение случайной величины , причем случайные величины, стоящие в одном столбце таблицы 1, взаимно независимы, следовательно, их отклонения от ожидаемого значения случайны, независимы и при сложении компенсируют друг друга, поэтому приближенно должно выполняться равенство
3-1
(2) Уг ХР3 •
г=1
Если число экспериментов N велико, то получившиеся уравнения (2) будут, как правило, линейно независимы, так как коэффициенты в разных уравнениях окажутся взаимно независимыми случайными величинами, и решить эти уравнения относительно неизвестных вероятностей (здесь г < ], ] = 1, •••, К) можно методом наименьших квадратов (МНК) [2] либо методом минимизации суммы абсолютных величин разностей (ММС) [7].
2. Модель системы для К = 5
Рассмотрим модель, которая включает пять остановок и пассажиропотоки между остановками в прямом направлении. Исследуем пассажиропотоки по маршруту автобуса, который включает все эти остановки. Согласно модели р12 характеризует пассажиропоток от первой до второй остановки, р 13 - пассажиропоток от первой до третьей остановки и т.д. для каждой остановки и каждого пассажиропотока. Для определения характеристик пассажиропотоков было проведено 100 экспериментов, по итогам которых составлены таблицы перемещений пассажиров (таблицы 3, 4).
Таблица 3. Информация о количестве вошедших в автобус пассажиров
№ остановки № эксперимента
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
1 4 6 5 4 4 5 5 4 5 6 4 4 5 6 4
2 3 2 2 3 3 2 2 3 2 2 3 3 3 2 2
3 1 2 2 1 2 2 1 1 2 1 2 1 1 1 1
4 3 2 2 2 2 2 2 2 2 3 2 3 2 2 2
Таблица 4. Информация о количестве вышедших из автобуса пассажиров_
№ остановки № эксперимента
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
2 3 3 3 3 2 3 2 3 3 2 3 3 3 3 4
3 3 5 4 4 3 4 4 4 3 5 3 3 3 4 2
4 3 3 3 1 4 3 2 1 3 2 4 2 3 3 2
5 2 1 1 2 2 1 2 2 2 3 1 3 2 1 1
В таблице 3, 4 представлена информация о полученных пассажиропотоках для первых 15 экспериментов.
Для поиска неизвестных р^, г < ], ] = 1, ...,К воспользуемся методом наименьших квадратов. Сумма квадратов отклоне-82
НИИ
N К 3-1
(3) 8 = ЕЕ(*2 - Е < Рз )2,
з=1 3=2 г=1
где N = 100, К = 5. Тогда
100 5 / 3-1 \ 2 100
(4) 8 = ЕЕ Уг - Е < Рз = Е N - х1Р12)2 +
3 ру1 = ^|(Уз ~г~ 42
з=1 3=2 \ г=1 / з=1
2 | з ^т ^ ^т ^ ^т ^ \ 2
+ (Уз - ХР13 - х2Р2З) + (У4 - х1р14 - х2р24 - ^Р34) +
+ (У5 - Х>15 - Р25 - Х3Р35 - <Р45) с ограничениями Р12 + Р13 + Р14 + Р15 = 1, Р23 + Р24 + Р25 = 1, Р34 + Р35 = 1, Рз ^ 0, для любых г, j.
Используя метод множителеИ Лагранжа [7], сократим число ограничений, сформировав при этом новую функцию 100
(5) 5 = Е [(У2 - х1р12)2 + (У5 - Х1Р13 - Х2Р23)2 +
т=1
+ (У4 - Х>14 - Х2Р24 - Х3Р34)2 +
+ (У5 - х1Р15 - Х2Р25 - Х3Р35 - Х4Р45)2 - А^Р12 + Р13 +
+ Р14 + Р15 - 1) - А2(Р23 + Р24 + Р25 - 1) - А3(р34 + Р35 - 1), с ограничениями р^ ^ 0, при любых г, j.
Задача состоит в том, чтобы найти минимум функции Б относительно неизвестных параметров р^, А^, для г <
; = 1, ..., к, г = 1, 2, з.
Минимума функция 5 достигает при следующих значениях параметров: Р12 = 0, 48; Р13 = 0, 52; рм = 0; Р15 = 0; Р23 = 0,67; Р24 = 0,33; Р25 = 0; Р34 = 1; Р35 = 0; Р45 = 1.
Общее число пассажиров на остановках колеблется в пределах от 9 до 14 человек. Для пассажиропотока из 11 человек средние значения пассажиропотоков распределены следующим образом (таблица 5).
Теперь найдем неизвестные рз, используя метод минимизации суммы абсолютных величин разностей. Сумма абсолютных
Таблица 5. Информация о среднем значении пассажиропотоков, полученная при помощи метода наименьших квадратов
№ остановки 2 3 4 5 Среднее количество пассажиров, находящихся на остановке
1 2,4 2,6 0 0 5
2 0 2,01 0,99 0 3
3 0 0 1 0 1
4 0 0 0 2 2
Среднее количество пассажиров, вышедших на остановке 2,4 4,61 1,99 2
величин разностей $1 для пяти остановочных пунктов и ста экспериментов:
- х>12| +
100 5 .7-1 100
(6) $1 = Е Е Уз - Е хЗ
2=1 7=2 г=1 2=1
+ У - Х1Р13 - х2р2З| + |у4 - Х>14 - х2р24 - Х3Р341 +
+ |Уз - Ж>15 - х2р25 - Х3Р35 - х4р45|] , с ограничениями Р12 + Р13 + Р14 + Р15 = 1, Р23 + Р24 + Р25 = 1, Р34 + Р35 = 1, Р. ^ 0, для любых г, ]. Введем переменные
(7) Я? =
Уз - Е.- Х2Р., если У.2 - х2Рг. ^ 0
0,
если у2 - е.=1 х2р. <0;
Е.= 1 ХзР. - Уз , если Уз - Е.- х2Рг. < 0, 0,
Очевидно, что
(8) (9)
84
Щ =
если УЗ - Е3=1 < Р. > 0
7-1
Уз - X)хЗРг. = ^З - Щ
2 ■ з
г=1
и
з-1
(10) yjj - Е xjptJ = Zjj + Wj
i=1
Тогда задача будет состоять в том, чтобы найти минимум функции Si:
100 5
(и) s = EE(zr + W)'
r=1 j=2
с ограничениями: pij ^ 0 при любых i, j; ^4=i+1 pij = 1 при любом i; yr = j1 xjPij + Zr - Wj, j = 2,..., 5; Zr ^ 0, Wj ^ 0, r = 1,..., N. При этом неизвестными параметрами в задаче будут величины pij, Zj, WJ, где i < j ^ 5, i = 1,..., 4.
Минимума функция Si достигает при следующих значениях параметров: p12 = 0, 5; p13 = 0,49; p14 = 0, 01; p15 = 0; P23 = 0, 68; p24 = 0, 32; p25 = 0; P34 = 1; P35 = 0; P45 = 1- Для пассажиропотока из 11 человек средние значения пассажиропотоков распределены следующим образом (таблица 6):
Таблица 6. Информация о среднем значении пассажиропотоков, полученная при помощи метода минимизации абсолютных величин разностей
и
Среднее количе-
№ остановки 2 3 4 5 ство пассажиров, находящихся на остановке
1 2,5 2,45 0,05 0 5
2 0 2,04 0,96 0 3
3 0 0 1 0 1
4 0 0 0 2 2
Среднее ко- 2,5 4,49 2,01 2
личество
пассажиров,
вышедших на
остановке
Полученная информация о пассажиропотоках позволяет делать выводы относительно величин пассажиропотоков для различных направлений, а также внутри самого маршрута. Наибольшие пассажиропотоки наблюдаются между 1 и 2 остановочными пунктами, а также между пунктами 1 и 3, 2 и 3, 3 и 4. Аналогично можно определить участки маршрута, где наблюдается небольшой пассажиропоток.
3. Модель системы для К = 10
Рассмотрим модель, которая включает десять остановок и пассажиропотоки между остановками в прямом направлении. Определим характеристики пассажиропотоков параметрическими методами.
Было проведено сто экспериментов, по итогам которых составлены таблицы перемещений пассажиров (таблица 7, 8).
Таблица 7. Информация о количестве вошедших в автобус пассажиров
№ остановки № эксперимента
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
1 5 6 6 4 6 5 6 4 4 4 5 6 4 4 6
2 3 3 2 3 2 3 2 2 3 2 3 2 3 2 2
3 2 1 1 2 2 1 1 2 2 2 2 2 2 2 1
4 3 2 3 3 3 3 3 3 3 3 3 3 3 3 3
5 3 4 3 3 3 5 5 3 4 5 4 3 4 3 5
6 3 3 3 3 4 3 4 3 4 3 4 4 4 3 4
7 1 1 2 1 2 1 1 1 1 2 2 2 1 1 1
8 3 2 2 2 2 2 2 3 2 3 2 2 2 3 2
9 1 2 2 1 1 1 2 2 2 1 2 1 2 2 2
В таблицах 7, 8 представлена информация о полученных пассажиропотоках для первых 15 экспериментов.
На основе экспериментальных данных найдены значения неизвестных параметров методом наименьших квадратов и методом минимизации суммы абсолютных величин. В таблице 9 86
Сетевые модели в управлении Таблица 8. Информация о количестве вышедших из автобуса
пассажирах
№ остановки № эксперимента
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
2 3 3 3 2 3 3 3 3 2 2 2 3 2 3 3
3 5 6 5 5 4 4 5 3 5 4 6 4 5 3 5
4 2 1 0 1 2 1 1 2 1 1 1 2 1 2 1
5 1 2 3 1 4 2 3 2 2 2 3 4 2 2 3
6 5 3 3 6 2 7 4 3 6 6 4 2 6 3 4
7 1 2 2 3 3 1 3 2 3 2 4 3 3 2 3
8 3 3 4 1 4 3 3 3 2 4 3 4 2 3 3
9 3 1 2 2 2 1 2 3 2 3 1 2 2 3 2
10 1 3 2 1 1 2 2 2 2 1 3 1 2 2 2
представлена информация о найденных значениях параметров.
Значения параметров, полученные разными методами близки друг к другу, что позволяет делать вывод о возможности применения любого из этих методов при поиске значений пассажиропотоков.
Таким образом, на основе результатов, полученных в сочетании натурных экспериментов и статистических методов, возможно собрать информацию не только о наиболее и наименее популярных остановках среди пассажиров, но и информацию о передвижениях пассажиров, в том числе на каких остановках они садятся, куда хотят приехать и т.д. Полученная информация может быть полезна при планировании новых маршрутов или оптимизации действующих за счет увеличения или уменьшения интервалов движения автобусов.
Результаты исследования возможно применять не только для автобусных перевозок. Это могут быть пригородные, внутригородские, наземные, подземные, а также водные или монорельсовые перевозки. Однако наиболее актуальной на данный момент остается проблема наземного городского транспорта, о чем свидетельствуют многочисленные известия о ситуациях с пробка-
ми во многих городах, а также тендеры на официальном сайте zakupki.gov.ru на решение (исследование) проблем общественного транспорта, а также анализа и оценки пассажиропотоков [810]. Поэтому в работе такое большое внимание уделяется именно наземному городскому транспорту.
Таблица 9. Информация о распределении пассажиропотоков, полученная при помощи параметрических методов
МНК ММС МНК ММС МНК ММС
Р12 0,48 0,5 Р29 0 0,01 Р56 0,99 0,97
Р13 0,43 0,37 Р210 0 0 Р57 0 0,02
Р14 0,04 0,06 Р34 0,77 0,65 Р58 0,01 0,01
Р15 0 0,01 Р35 0,09 0,1 Р59 0 0,01
Р16 0,05 0,05 Р36 0,1 0,13 Р510 0 0
Р17 0 0,01 Р37 0,02 0,06 Р67 0,78 0,65
Р18 0 0,01 Р38 0,03 0,05 Р68 0,21 0,24
Р19 0 0,01 Р39 0 0,01 Р69 0,01 0,11
Р110 0 0 Р310 0 0,01 Р610 0 0
Р23 0,98 0,98 Р45 0,66 0,61 Р78 0,79 0,59
Р24 0,01 0,01 Р46 0,2 0,21 Р79 0,2 0,26
Р25 0,01 0 Р47 0,04 0,08 Р710 0,01 0,15
Р26 0 0,01 Р48 0,09 0,1 Р89 0,86 0,72
Р27 0 0 Р49 0,01 0 Р810 0,14 0,28
Р28 0 0 Р410 0 0 Р910 1 1
4. Выводы и перспективы
В работе проведено исследование статистической модели интенсивности пассажиропотоков между остановками. На основании натурных экспериментов построена математическая модель системы. На примерах для К = 5, К =10 продемонстрирована методика определения пассажиропотоков между остановками. Величины пассажиропотоков найдены при помощи методов наименьших квадратов и минимизации абсолютных величин разностей. Заметим, что предложенная процедура определения 88
пассажиропотоков легко реализуема на практике. Например, в автобусах в больших городах пассажиры входят на остановке в первую дверь и регистрируют свой билет в специальном автомате. Выход осуществляется через другие двери, где можно просто установить дополнительные фиксирующие устройства.
Литература
1. АРТЫНОВ А.П., ЕМБУЛАЕВ В.Н., ПУПЫШЕВ А.В. И
ДР. Автоматизация управления транспортными системами. - М.: Наука, 1984. - 272 с.
2. БАНДИ Б. Методы оптимизации. Вводный курс. - М.: Радио и связь, 1988. — 128 с.
3. БОРОВКОВ А.А. Теория вероятностей. - М.: Эдиториал УРСС. - 472 с.
4. БУРЕ В.М., ПАРИЛИНА Е.М. Теория вероятностей и математическая статистика: учебник. -- СПб.: «Лань», 2013.-416 с.
5. ВАРЕЛОПУЛО Г.А. Организация движения и перевозок. -М.: Транспорт, 1981. - 199 с.
6. ГИЛЛ Ф., МЮРРЕЙ У., РАЙТ М. Практическая оптимизация. - М.: Мир, 1985. — 510 с.
7. ДРЕЙПЕР Н., СМИТ Г. Прикладной регрессионный анализ. Т. 1 - 2-е изд., перераб. и доп. - М.: Финансы и статистика, 1986. - 366 с.
8. Официальный сайт Российской Федерации в сети Интернет для размещения информации о размещении заказов на поставки товаров, выполнение работ, оказание услуг http://zakupki.gov.ru/pgz/public/action/orders/info/common_-info/show?notificationId=6795209 (дата обращения: 20.12.2013).
9. Официальный сайт Российской Федерации в сети Интернет для размещения информации о размещении заказов на поставки товаров, выполнение работ, оказание услуг http://zakupki.gov.ru/pgz/public/action/orders/info/common_-
info/show?source=epz¬ificationId=8275955 (дата обращения: 25.12.2013).
10. Официальный сайт Российской Федерации в сети Интернет для размещения информации о размещении заказов на поставки товаров, выполнение работ, оказание услуг http://zakupki.gov.ru/pgz/public/action/orders/info/common_-info//show?source=epz¬ificationId=5629611 (дата обращения: 25.12.2013).
11. ПОСТНИКОВА Е. Квантильная регрессия. [Электронный ресурс]. - URL: http://allmath.ru/highermath/probability/probability48/proba-bility.htm (дата обращения: 12.09.2013).
ESTIMATING PASSENGERS TRAFFIC CHARACTERISTICS IN TRANSPORT SYSTEMS
Vladimir Bure, St.Petersburg State University, St. Petersburg, Doctor of Science, professor (vlb310154@gmail.com). Vladimir Mazalov, Institute of Applied Mathematical Research Karelian Research Center of RAS, Petrozavodsk, Doctor of Science, professor (vmazalov@krc.karelia.ru).
Nina Plaksina, Petrozavodsk State University, Petrozavodsk, post-graduate student (plaksina_nv@mail.ru).
Abstract: We consider a statistical model ofpassengers' flow between bus stops. Intensity of incoming and outcoming passenger flows is supposed to be known. We propose a model to estimate the function of distribution of passengers, who come in at a certain bus stop, among all successive bus stops. We also present the results of computer simulations based on empirical experiments.
Keywords: game theory, passenger traffic, transport system, natural experiment, equilibrium.
Статья представлена к публикации членом редакционной коллегии А. А. Печниковым
Поступила в редакцию 30.11.2013. Опубликована 31.01.2014.