Вероятностно-статистический метод расчета вместимости склада морского порта

ВЕСТНИК«)

ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА ^^

МОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА

DOI: 10.21821/2309-5180-2018-10-3-507-519

STOCHASTIC METHOD FOR SEA PORT STORAGE CAPACITY CALCULATION

S. S. Valkova

Admiral Makarov State University of Maritime and Inland Shipping, St. Petersburg, Russian Federation

Deals with the development and expanding of the analytical calculation method of the assessment of the characteristics of the sea port's cargo storage. The idea of the method is based on the transfer from determined values to stochastic ones. Traditional methods provide the resulting parameters in the form of constants, which explains why the correspondent approach is referred to as "static technique". The static values do not permit to form the perception of the possible range of dispersion of the characteristics around the mean values caused by random fluctuation of the input datafor calculation. Accordingly, it does not provide any ways to study the sensitivity of the end results to the inaccuracy and indeterminacy of the primary values. The growth of the demands for the quality of the results of the technologic design procedure translates into the needs to receive all the parameters' assessments not in the form of constant values, but functional dependencies, describing possible methodological variety and stochastic fluctuations of input data. In the first head, this method should be practicedfor calculation of the primal basic values which serve as the input for following flow calculation procedure, thus determining all main operational characteristics of the infrastructural object under design. Aiming at this target, in the beginning of the study main cargo operational functions of the sea port are analyzed, which enables to identify the cargo storage as the main element involved in their realizations. Based on this reasoning, the main operational lows govern the sea port's cargo storage operations. These lows take form of differential and integral equations and serve as the main tools in the study of the behavior of the cargo volume stored in the storage facilities of the sea trade port. This study is described in small details in the paper. The discussion of the results allows to conclude that as a result it is possible to obtain not only more accurate mean values, but also to estimate the spread around them Keywords: sea port, storage, simulation, Monte-Carlo technique.

For citation:

Valkova, Svetlana S. "Stochastic method for sea port storage capacity calculation." Vestnik Gosudarstven-nogo universiteta morskogo i rechnogo flota imeni admirala S. O. Makarova 10.3 (2018): 507-519. DOI: 10.21821/2309-5180-2018-10-3-507-519.

УДК 656:6

ВЕРОЯТНОСТНО-СТАТИСТИЧЕСКИЙ МЕТОД РАСЧЕТА ВМЕСТИМОСТИ

СКЛАДА МОРСКОГО ПОРТА

С. С. Валькова

ФГБОУ ВО «ГУМРФ имени адмирала С. О. Макарова», Санкт-Петербург, Российская федерация

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

_Oi

507

ЛВЕСТНИК

............ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

Х^ОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА

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

Ключевые слова: морские порты, склады, моделирование, методы Монте-Карло.

Для цитирования:

Валькова С. С. Вероятностно-статистический метод расчета вместимости склада морского порта / С. С. Валькова // Вестник Государственного университета морского и речного флота имени адмирала С. О. Макарова. — 2018. — Т. 10. — № 3. — С. 507-519. DOI: 10.21821/2309-5180-2018-10-3-507519.

Введение

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

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

Методы и материалы

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

- передача грузопотоков между видами транспорта (функция интерфейса);

- согласование размеров транспортных партий (технологическое хранение);

- сглаживание неравномерности движения транспортных средств (буферная функция);

- коммерческое хранение грузов (складская функция);

- преобразование грузопотоков (логистическая функция). г В то же время единым интегральным параметром, характеризующим работу склада, является объем хранимого в нем в каждый момент груза [5] - [6]. Изменение этого параметра во времени, т. е. то, что называется поведением в системном смысле, является совокупным результатом действия всех факторов, относящихся к перечисленным ранее функциям. В этом смысле очевидной становится задача идентификации и оценки степени влияния всех факторов, включенных в соответствующие выполняемым функциям кластеры.

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

E

CO

oo о

где £(г) — объем хранения груза на складе в момент времени г; Щ) — объем завезенного на склад груза к моменту времени г; 0(г) — объем вывезенного со склада груза к моменту времени и

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

Если обозначить как Де(t), Д/(t), Д0(t) приращения объема хранимого, поступающего и убывающего груза, то соотношение (1) очевидным образом сохранится, преобразуюсь в его конечно-разностный аналог:

Де

Действительно, изменение количества хранимого на терминале груза в сторону увеличения или уменьшения определяется разностью поступившего и убывшего на терминал груза не за весь период [0, г], а за произвольный интервал Дt е [0, г].

Переход от значений Д £ (г), Д 1( г), Д0( г) к относительным величинам в рассматриваемом случае близок к операции дифференцирования, т. е.

AE(t) AI(t) AO(t) At At At

и при A(t) ^ 0 действительно им является:

dE0) _ Ш

& dt dt

Интерпретация выражений (3) и (4) не вызывает сложности: они означают, что скорость изменения объема складирования е(г) есть разница между скоростями роста /(г), т. е. скоростью поступления груза на склад и скоростью убывания 0(г) — скоростью вывоза груза со склада:

г

Разница в уравнениях (3) и (5) состоит лишь в том, что в первом случае используются разностные, а во втором — дифференциальные уравнения.

Дифференцированием уравнения (1) получим уравнение (5), которому, как известно, соответствует целое семейство первообразных функций, т. е. уравнений вида

t t t Е

0 0 0

где С — произвольная константа.

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

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

2 О

8

со г

оо о

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

Пусть груз, составляющий некоторую произвольную партию п объема Уп, хранится на складе неравномерно: часть партии ах хранится в течение времени Тр часть партии а2 хранится в течение времени Т2 и т. д. Хранение на складе партии объемом V требует, таким образом, определенное число место-суток, которое назовем работой склада по хранению партии:

а, =±К а,Т, = К^ТаТ.

г=1 г=1

I

Величина представляет собой средневзвешенное значение длительности хранения от-

i=1

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

Т, = ^Т,.

1=1

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

а. = = Г„Т„.

I=1

Общий объем работы склада при хранении всех партий Vn , 1, ... , N есть

А = 2Х .

п =1 п =1

Среднее значение объема хранения груза Е есть величина, которая за время Т дает то же значение объема работы склада, т. е.

ЕТ = А.

Если рассматриваемый период Т составляет год, т. е. Т = 365 сут, все партии равны между собой: Vп= V, и средний срок хранения Т для всех партий одинаков: Т = Тхр, то имеем

Е • 365 = УУПТ„ = МУлТп С12)

или

ЫУТ

Е —

365

Величина NV = QT есть объем груза, проходящего через склад за период Т = 365 сут, откуда имеем

ЫУТ Q Т

Е _ _^ _ год хр

" 365 " 365

г т Т 365

С другой стороны, величина Тинт есть средний интервал между поступлениями

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

и-1

NVT VT

E =-XL = —XL.

365 T V '

инт

Обе выведенные формулы (14), (15) допускают различную интерпретацию. Формула (14) может быть переписана в виде

е = А

•^ГОД гр

ТХР

365

Величина- в логистике известна как «оборачиваемость склада», она показывает сколько

хр

раз склад полностью возобновляет свое содержимое за рассматриваемый период времени (в данном случае за год). Умножение среднего объема хранения склада Е на этот показатель дает грузопоток за этот период 2год.

Тх

Не менее очевидна интерпретация второй полученной формулы (15): величина показы-

инт

вает, сколько партий в среднем прибывает на склад за время хранения одной партии, т. е. сколько партий должно одновременно находиться на складе. Умножив это число на объем партии V получим среднюю величину объема единовременного хранения Е.

Итак, совпадение значений полученных соотношений (8) и (10) при внесении в них упрощающих предположений (12) с известными и используемыми в расчетно-аналитических методах формулами (15) и (16) подтверждает корректность сделанного вывода. В то же время зависимости соотношений (8) и (10) для целей, поставленных в настоящем исследовании, оказываются намного более продуктивными. В частности, введенные понятия (срок хранения груза на складе и интервал его поступления) позволяют объективно разграничить технологическое хранение, необходимое для согласования равномерной интенсивности обработки в порту смежного транспорта и пиковой обработки морских судов, и внешнее коммерческое хранение.

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

т тхр

инт инт

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

2 2

говоря, является дополнительным хранением.

Результаты

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

Если эта партия является единственной в году, то среднее значение загрузки склада в этом

NVTx Tx

случае составляет значение E = = (эквивалентное хранение 1/365 доли от работы VTxp

по хранению одной партии).

2 О

8

ЛВЕСТНИК

............ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

Х^ОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА

Рис. 1. Эпюра простейшей единичной партии

Если через склад проходят грузовые партии с интервалом поступления Тинт = 2Тхр, то дина

мика изменения объема хранения будет такова, как показано на рис. 2. Как следует из этого рисун-

УТХ V

ка и формулы (15), средний объем запаса на склад составит величину Е — -у-^ = — •

инт 2

Рис. 2. Изменение объема складирования при Т = 2Т

"г Г инт хр

Если уменьшить интервал между поступлением партий до значения граничного = Тх ,

УТх

то средний объем хранения на складе составит Е = — У, как это показано на рис. 3.

со г

Рис. 3. Изменение объема складирования при Тинт = Тхр: а — треугольная форма партии; б — прямоугольная форма партии

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

и равной — (рис. 4).

ВЕСТНИК«)

ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА ^^

МОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА

/

V 2

V 2

Рис. 4. Изменение объема складирования при Т < Т

Таким образом, для случаев простейшего треугольного закона формирования партий (1), регулярного интервала их поступления (2) и равенства их объема (3) максимальный объем хранения может быть получен добавлением к среднему значению величины указанной амплитуды, т. е.

УТ

Е =_ХР

тах Т

У

+ — . 2

Для рассматриваемого случая треугольного закона формирования партий время хранения является половиной времени формирования партии (времени от начала поступления груза данной партии до убытия ее со склада), т. е. Тх = 2Тформ. В принятых обозначениях формула (17) может

быть записана более компактно:

форм

VT

Е __ХР

тах т

V _ + 2 _

форм V V

2Т + 2 _ 2

Тф

форм

V

инт

+1

Именно такой вид формулы предлагается в нормах технологического проектирования морских портов [7] - [8]. В то же время на форму результирующей кривой оказывают более существенное влияние другие факторы из числа идентифицированных ранее.

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

а)

Рис. 5. Динамика изменения объема хранения при колебаниях интервала поступления: а — изменение регулярных и иррегулярных поступлений; б — эпюра для регулярного интервала поступления; в — эпюра для иррегулярного интервала поступления; г — суммарный объем хранения склада для регулярного интервала поступления; д — суммарный объем хранения склада для иррегулярного интервала поступления

ЛВЕСТНИК

............ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

Х^ОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА

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

Рис. 6. Динамика изменения объема хранения при колебаниях размера партий

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

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

со г

оо о

514J

Обсуждение

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

Q т vT

•^год xp -ч xp

ся лишь две величины, а именно -, в правой части уравнения (15) — три, а именно у—.

365 инт

Следовательно, определяемая этими уравнениями новая случайная величина E может иметь иной статистический характер, и в этом случае уравнения (14) и (15) перестают быть тождественными.

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

Иными словами, при инженерных расчетах уравнение (15) по умолчанию интерпретируется

[V]-Гтх ]

как M [E] =- . Математические ожидания исходных величин в правых частях уравнений,

[уинт]

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

обратной величины — , т. е. M

. вовсе не есть M

1

, поскольку величина — является иной X

случайной величиной. Как следствие, использование арифметической операции деления в уравне -VT

нии —ХР будет приводить к погрешности оценки искомого значения даже при известном значе-

инт

нии M[T ].

инт

Таким образом, использование традиционных методов расчета по рекомендуемым нормами формулам (14) и (15) имеет методическую погрешность в расчете предполагаемого центрального значения, не позволяя судить о характере случайных разбросов оцениваемой величины вокруг этого значения. Если, к примеру, в 10 случаях некоторых испытаний будет наблюдаться значение некоторой величины 50, и в 10 случаях — 950, то средним значением будет 500, т. е. величина, равноудаленная от реально наблюдаемых величин.

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

VTx

Txp = 10 сут и одинаковым интервалом поступления Гнт = 5 сут. Формула E = тХР для выбранно-

10 инт

го примера дает значение E = 1000 • — = 2000.

Пусть теперь интервал поступления партий Тинт является случайной величиной, среднее значение которой составляет те же 5 сут. Объем хранения на складе E также будет являться случайной величиной, конкретные значения которой будут, как представляется, колебаться вокруг значения E = 2000. Рассмотрим это на числовых примерах.

Пример 1. Интервалы между поступлением партий составляют значения 2, 3, 5, 10. Сред-

2 + 3 + 5 +10

ний интервал составляет величину-4-- 5. Расчет размера склада для каждого интервала

по приведенной выше формуле дает значения 5000, 3333, 2000, 1000. Среднее значение склада

в этом случае составит величину 5000 + 3333 + 2°00 +1000 _ 2833, что значительно отличается от ожидаемых 2000. 4

Пример 2. Интервалы поступления составляют значения 1, 2, 8, 9. В среднем интервал составляет ту же величину 1 + 2 + 8 + 9 _ 5. Расчет размера склада дает значения 10000, 5000, 1250,

^ 4 10000 + 5000 +1250 +1111 1111. Среднее значение склада в этом случае составит величину ---= 4340.

Полученное значение более чем в два раза выше ожидаемого.

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

В проведенных рассуждениях предполагалось, что M

1

T

инт

. В действительно-

' M TJ

8

1 й й 1 1 1 Е

сти, —— есть новая случайная величина, значения которой составляют —, —,..., — . Ее среднее 2

инт T T TN

1 1 1

— +-+ ... +-

T T T

значение (математическое ожидание) есть величина —---—, вовсе не равная значению

N ^

-1-=-—-. Действительно, значения математического ожидания т в примерах

T + T2+...+ TN T + T2+_+T— инт

— 1 + 1 + + 1 1111 1 + 1 + + 1 1111

—1---H...H----i___i___i____—1---ъ...н--_+---1---1--

ГТ1 ГТ1 ГТ1 I I I ГТ1 ГТ1 ГТ1 l l l

1 и 2 есть 71 72_^ = 2 3 5 10 = 0,28 и _hL = 1289 = 0,43,

N 4 N 4

что при умножении на VTxp = 10000 дает те же результаты, что и в приведенных примерах 1 и 2,

что снимает кажущееся противоречие.

Г515

ЛВЕСТНИК

............ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

Х^ОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА

СО

г

Если к точности расчетных значений не предъявляются особых требований и с методической погрешностью можно мириться, то подобный подход может быть использован [9] - [13]. В инженерной практике результаты получают расчетно-аналитическими методами этого класса в попытке избежать всевозможных ошибок и отклонений путем умножения на некоторый коэффициент запаса К > 1, который в различных дисциплинах, предметах, объектах и формулах называют по-разному, но который, по сути, является коэффициентом незнания. В тех случаях, когда требуется более адекватное суждение о случайном поведении оцениваемой величины, в науке используют иной подход, в основе которого лежит метод статистических испытаний, или метод Монте-Карло.

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

Основным исполнительным механизмом метода Монте-Карло является генерация значений случайной величины по заданному закону ее распределения, что обычно выполняется с помощью метода обратной функции [14]. Функция распределения F случайной величины X есть вероятность того, что X < х, т. е. функция вида F(x) = Р(X < х). Областью определения функции F являются значения случайной величины, областью значений — интервал [0,1], областью определения обратной F- — интервал [0,1], областью значений — вероятные значения случайной величины. Если генерировать равномерно распределенные на интервале [0,1] случайные числа, то значения F- будут группироваться плотнее там, где крутизна этой кривой выше, т. е. вероятность появления значений будет выше там, где выше плотность функции F, т. е. в статистическом смысле будет соответствовать заданному распределению случайной величины.

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

а)

б)

в)

Рис. 7. Функции распределения исходных расчетных величин: а — объема партии; б — срока хранения; в — интервала поступления

Ш

В таблице приведен пример массива статистических испытаний — сгенерированные значения исходных значений и вычисленные в каждом статистическом испытании значения объема хранения на складе Е:

N исп 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17

к 960 760 860 760 960 960 560 1060 1260 560 760 1060 760 1160 460 960 1260

Т хр 9,4 9,9 9,9 11,1 10,1 9,9 9,6 11,1 9,4 9,6 11,1 9,9 10,1 9,9 9,6 9,1 11,1

Т инт 4,6 3,6 4,6 4,9 4,9 4,9 4,9 4,4 4,1 5,1 4,9 4,6 5,6 5,4 5,6 4,9 3,9

Е 1945,9 2070,3 1836,2 1734,4 1993,8 1944,6 1105,6 2695,4 2863,6 1051,7 1734,4 2263,2 1368,0 2131,2 787,1 1796,9 3617,4

ВЕСТНИК«)

ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА ^^

МОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА

На рис. 8 для данного примера приведен результат статистической обработки массива значений — гистограмма плотности распределения вероятности значений объема хранения на складе Е.

Рис. 8. Гистограмма плотности распределения объема хранения груза на складе

Математические ожидания объема партии, срока хранения и интервала прибытия партий для примера на рис. 7 составляют, соответственно: M [ V] = 1000, M [ Гр] = 10 сут и M [ Тинт] = 5 сут, т. е. совпадают с данными рассмотренных выше примеров.

Как видно из рис. 8, значения E действительно располагаются в районе 2000 единиц, т. е. совпадают с ожидаемыми, полученными в расчетном примере. В то же время точное среднее значение, подсчитанное по массиву статистических испытаний, составляет детерминированную величину, равную 1900, а кривая на рис. 8 позволяет судить о возможных отклонениях и их вероятности, что дает возможность принимать более обоснованные инженерные и предпринимательские решения.

Выводы

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

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

3. Средние значения объема хранения при детерминированных исходных данных определяются несколькими инвариантными алгебраическими соотношениями.

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

5. При рассмотрении всех значений как случайных величин эти поправки становятся мало- оо

о

значимыми, а сами математические методы утрачивают свою применимость. р

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

Заключение

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

_«

ГйТ

ЛВЕСТНИК

............ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

Х^ОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА

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

СПИСОК ЛИТЕРАТУРЫ

со г

e>j

ТЛЕ

1. Thorensen C. A. Port designer's handbook / C.A. Thorensen. — London: Thomas Telford Limited,

2010. — 554 p.

2. UNCTAD Monographs on Port Management. Monograph № 9. Multipurpose port terminals. Recommendations for planning and management. — March 1991.

3. Кузнецов А. Л. Классификация и функциональное моделирование эшелонированных контейнерных терминалов / А. Л. Кузнецов, А. В. Кириченко, А. А. Давыденко // Вестник Государственного университета морского и речного флота имени адмирала С. О. Макарова. — 2015. — № 6 (34). — С. 7-16. DOI: 10.21821/2309-5180-2015-7-6-7-16.

4. Current issues in shipping, ports and logistics / T. Notteboom (ed.). — Asp/Vubpress/Upa,

2011. —608 p.

5. Елисеева А. С. Принципы организации грузового терминала / А. С. Елисеева. — Ростов-на-Дону: Феникс, 2003. — 329 с.

6. Маликов О. Б. Склады и грузовые терминалы / О. Б. Маликов. — М.: АСТ, 2005. — 560 с.

7. РД 31.3.05-97. Нормы технологического проектирования морских портов. — М., 1998.

8. Кузнецов А. Л. О несовершенстве нормативной базы технологического проектирования морских портов / А. Л. Кузнецов, В. А. Погодин // Морские порты. — 2017. — № 6. — С. 18-22.

9. Щербакова-Слюсаренко В. Н. Разработка функциональной модели контейнерного терминала типа «сухой порт» и принципов ее использования в технологическом проектировании / В. Н. Щербакова-Слюса-ренко, В. А. Погодин, А. С. Ткаченко // Вестник Государственного университета морского и речного флота имени адмирала С. О. Макарова. — 2017. — Т. 9. — № 1. — C. 48-60. DOI: 10.21821/2309-5180-2017-9-1-48-60.

10. Михаэль Д. Складская логистика. Новые пути системного планирования / Д. Михаэль; пер. с нем.; под ред. Г. П. Манжосова. — М: КИА центр, 2004. — 136 с.

11. International Handbook of Maritime Economics / K. Cullinane (ed.). — Edward Elgar Publishing,

2011. — 520 p.

12. European Comission. Directorate-General for Mobility and Transport. White Paper on Transport: Road-map to a Single European Transport Area: Towards a Competitive and Resource-efficient Transport System. — Publications Office of the European Union, 2011.

13. Памбухчиянц В. К. Организация, технология и проектирование торговых предприятий / В. К. Пам-бухчиянц. — М.: ИВЦ «Маркетинг», 1999. — 320 с.

14. Кобзарь А. И. Прикладная математическая статистика / А. И. Кобзарь. — М.: ФИЗМАТЛИТ,

2012. — 813 с.

REFERENCES

1. Thorensen, C.A. Port designer's handbook. London: Thomas Telford Limited, 2010.

2. UNCTAD Monographs on Port Management. Monograph № 9. Multipurpose port terminals. Recommendations for planning and management. March 1991.

3. Kuznetsov, Aleksandr Lvovich, Aleksandr Viktorovich Kirichenko, and Aleksandr Aleksandrovich Davydenko. "Classification and functional modeling of echeloned container terminals." Vestnik Gosudarstvennogo universiteta morskogo i rechnogoflota imeni admirala S.O. Makarova 6(34) (2015): 7-16. DOI: 10.21821/2309-51802015-7-6-7-16.

4. Notteboom, Theo, ed. Current issues in shipping, ports and logistics. Asp/Vubpress/Upa, 2011.

5. Eliseeva, A.S. Printsipy organizatsii gruzovogo terminala. Rostov-na-Donu: Feniks, 2003.

6. Malikov, O.B. Sklady i gruzovye terminaly. M.: AST, 2005.

ВЕСТН1

ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

МОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА

7. Russian Federation. Guidance Document RD 31.3.05-97. Engineering Design Standards for Seaports. M.,

8. Kuznetsov, A.L., and V.A. Pogodin. "O nesovershenstve normativnoi bazy tekhnologicheskogo proektirovaniya morskikh portov." Morskie porty 6 (2017): 18-22.

9. Shcherbakova-Slyusarenko, Victoria N., Vladimir A. Pogodin, and Andrei S. Tkachenko. "The development of the functional model for the "dry port" type container terminal and principles of its use in the technologic design." Vestnik Gosudarstvennogo universiteta morskogo i rechnogo flota imeni admirala S.O. Makarova 9.1 (2017): 48-60. DOI: 10.21821/2309-5180-2017-9-1-48-60.

10. Mikhael', D. Skladskaya logistika. Novye puti sistemnogo planirovaniya. Edited by G.P. Manzhosov. M: KIA tsentr, 2004.

11. Cullinane, Kevin, ed. International handbook of maritime economics. Edward Elgar Publishing, 2011.

12. European Comission. Directorate-General for Mobility and Transport. White Paper on Transport: Roadmap to a Single European Transport Area: Towards a Competitive and Resource-efficient Transport System. Publications Office of the European Union, 2011.

13. Pambukhchiyants, V. K. Organizatsiya, tekhnologiya i proektirovanie torgovykh predpriyatii. M.: IVTs «Marketing», 1999.

14. Kobzar', A.I. Prikladnaya matematicheskaya statistika. M.: FIZMATLIT, 2012.

1998.

ИНФОРМАЦИЯ ОБ АВТОРЕ

INFORMATION ABOUT THE AUTHOR

Валькова Светлана Сергеевна — соискатель ФГБОУ ВО «ГУМРФ имени адмирала С. О. Макарова»

198035, Российская Федерация, Санкт-Петербург, ул. Двинская, 5/7

e-mail: vlvalkov@yandex.ru, kaf_pgt@gumrf.ru

Valkova, Svetlana S. — Applicant Admiral Makarov State University of Maritime and Inland Shipping

5/7 Dvinskaya Str., St. Petersburg 198035, Russian Federation

e-mail: vlvalkov@yandex.ru, kaf _pgt@gumrf.ru

Статья поступила в редакцию 15 мая 2018 г.

Received: May 15, 2018.

[5191