УДК 519.681
В.Н. Никишов, Л.А. Сараев *
ПРОГНОЗИРОВАНИЕ ПОТЕРЬ И УПРАВЛЕНИЕ РИСКАМИ ГРУЗОПЕРЕВОЗОК
В публикуемой статье предложена математическая модель управления рисками грузоперевозок, которая представляет собой обобщение классической транспортной задачи, учитывающей риски потери грузов во время их доставки. Рассмотрены рыночные ситуации превышения предложения над потребностями. Получен численный анализ результатов расчетов в виде зависимости решения от параметров задачи.
Ключевые слова: поставщик, потребитель, целевая функция, ограничения, стохастическая транспортная задача, управление рисками, бета-распределение.
Пусть некоторый груз, сосредоточенный у т поставщиков Д в количестве щ единиц (¡=1..т), необходимо доставить п потребителям В, в количестве Ь, @=1..п). Стоимость перевозки единицы груза с, от поставщика Д к потребителю В, известна. Очевидно, что при транспортировке груза вследствие различных факторов риска (например ДТП) возможны его потери. Эти потери при перевозке груза от Д к В, в размере х, можно представить в виде
= 1 и У Ч .
Здесь I, — индикаторная случайная величина, фиксирующая реализацию риска и принимающая значения 0 или 1, у, — размер убытка от гибели или повреждения груза.
Р( = 1) = Р( > 0) = д1}, Р(1 У = 0) = 1 - Р( > 0) = 1 - д1} = Рц , (1)
где qij — вероятность реализации риска потери груза при перевозке от поставщика Д к потребителю В,.
Очевидно, что распределения величин ц и у, связаны соотношением
Р(у ц < г) = Р( < г\11} > 0). (2)
В актуарной математике размер убытка у, выражается через степень ущерба х, и стоимость перевозимого груза х, с помощью формулы у,=х, X,. При этом случайные величины I, и х, являются независимыми, поскольку попадание в аварию зависит от одних факторов, а размер ущерба совершенно от других факторов [1; 2]. Степень ущерба х, определяется отношением размера потерь к стоимости груза и удовлетворяет неравенству 0<х, <1.
* © Никишов В.Н., Сараев Л.А., 2011
Никишов Виктор Николаевич (TSh-sea05@yandex.ru), Сараев Леонид Александрович (saraev@ssu.samara.ru), кафедра математики и бизнес-информатики Самарского государственного университета, 443011, Российская Федерация, г. Самара, ул. Акад. Павлова, 1.
Будем считать, что степени ущерба не зависят от маршрута перевозки и имеют одинаковое бета-распределение [3].
Р(Х <= +11 П а > 0; Ь > 0. (3)
Г (а ) Г (Ь) 0
Для определения параметров а и Ь в законе распределения (3) нужно задать математическое ожидание
Хо = Е (Хц)
и дисперсию случайной величины степени ущерба ху
Ъ(ХЦ) = °2(Ху) = Я2(Х). Тогда параметры а и Ь определятся соотношениями
а = -Хо + Хо2(1 -Хо) 1°2(Х); Ь = (1 - Хо)а/Хо .
Математическое ожидание и дисперсию случайной величины степени ущерба Ху можно определить по статистическим данным
1 ' 1 '"1 Хо =-XХ1}(к); &2(Х) = -— XХ(к)"Хо)2.
' к =1 11 к =1
Здесь Ху(к) — фактически наблюдаемая реализация величины хр принимающая дискретный набор значений с вероятностями Р(Ху = Х(к)) = г (к). В этом
случае математическое ожидание принимает вид
£
Хо = Е(Хг}) = 2Х(к)г(к).
к =1
Доставка грузов от поставщика к потребителю осуществляется на основе контрактов, при этом условия перевозки зависят от того, кто несет риск случайной гибели, повреждения груза. В международной практике условия перехода риска от продавца к покупателю регламентируются условиями «Инкотермс». Здесь можно выделить несколько вариантов [4].
В первом варианте, наиболее распространенном на практике, покупатель оплачивает стоимость груза и самостоятельно вывозит его. В этом случае риск случайной гибели, утраты, повреждения груза возложен на потребителя, он также несет все последствия недопоставки груза или его излишнего завоза. Этот вариант перевозки приводит к задаче потребителя. Потребитель принимает весь риск на себя и минимизирует потери груза при транспортировке и потери, связанные с хранением излишне завезенной продукции.
Во втором варианте доставка поставщиком груза потребителю оплачивается только по факту доставленных грузов. Все убытки от потерь груза, его недопо-
ставки и расходы по хранению невывезенной продукции несет поставщик. Этот вариант перевозки приводит к задаче поставщика. Поставщик в силу конкуренции вынужден доставить груз, принимает на себя его потери, штрафы за ненадлежащее исполнение обязательств по доставке и несет дополнительные расходы, связанные с хранением остатка продукции.
В третьем варианте доставка осуществляется транспортно-экспедиционной организацией, а перевозчик принимает на себя риск потери груза.
Потери груза при перевозке от поставщика от Д к потребителю В) в размере Хц есть случайная величина
1аХаХа .
Среднее значение и дисперсия величины потерь даются выражениями
Е( ) = Е(I и)Е(х1})Ху = ^ХоХа ; £>( ) = Хо2ХI (Ра + ^ (Х)/Хо2). В результате потерь от поставщика Д в пункт В будет доставлен груз в размере
Ха - £ а = Х I;(1 -1а Ха).
Общие потери по всем перевозкам есть случайная величина
т п
5 = Ц^Ха1 аХа ■ (4)
1=1 ]=!
В случае если окажется, что имеет место неравенство
т
ЦХа (1- 1аХа ) < Ь, 1=1
то спрос потребителя не будет удовлетворен, а нанесенный потребителю ущерб будет пропорционален объему неудовлетворенного спроса
ф = g), = Ьа - ±Ха(1 - 1аХа). (5)
=1
Здесь g}'' — штраф за нехватку единицы продукта, Н(х) — функция Хевисайда. Если будет иметь место неравенство
т
Ц Ха (1 - 1аХа) > Ьа, =1
то возникнет необходимость в хранении избыточного продукта, а дополнительные затраты системы будут пропорциональны объему избыточного продукта
Прогнозирование потерь и управление рисками грузоперевозок 139
0&=-8), (6)
где 8 ) — затраты на хранение единицы.
Транспортные расходы в целом выражаются соотношением [5]
т п
с = XXе
I =1 Ц=1
Здесь с. — транспортные расходы по перевозке единицы груза из пункта Д в пункт В. Тогда общая величина расходов на транспорт с учетом потерь и штрафов составит величину
2 = X X ( С + ) + Ф + с(+)). (7)
1=1 Ц=1
Сведем задачу потребителя к задаче линейного программирования методом, аналогичным задаче случайного спроса, рассмотренного в работе [6]. Введем вспомогательную случайную величину дефицита
и ()={( ^))
и вспомогательную случайную величину избыточного продукта
V ( ^ ^^))
в пункте потребления с номером у. Величина потерь при перевозке грузов из
всех пунктов в пункт потребления с номером у выражается формулой
>) = "
< )=X1 хх.
=1
При этом имеет место очевидное соотношение
и Ц)- V ( )= X Хц - Ь - ^ =1
Не ограничивая общности рассуждений, рассмотрим случай, когда число поставщиков т=2, число потребителей п=3, а случайная величина /ц принимает три значения Хд=х(к) с соответствующими вероятностями г(к). В этом случае 5 пробегает значения от 0 до 9, и случайная величина
)= Х1 ц11 X Ц + Х2,12Ц = (Х1 ц11 ц + Х2и12Ц )Х(к)
будет иметь следующее распределение
s 0 1, 2, 3 4, 5, 6 7, 8, 9
0 Х1 аХ(5) Х2Х - 3) (Х а + Х2 а )Х(5 - 6)
Р( ) } Р1 аР2 а 41 аР2 / (5) Р1 ¡42 / (5 - 3) 41 ]42 / (5 - 6)
Целевая функция стохастической транспортной задачи записывается в виде
Е (2) = ЦЦ^ + Ц К Ц К
I=1 а=1 а =1 х=о а=1 *=о
3 9 3 9
+ц *г Ц а и а )+х*г Ц р( V аи шт.
а=1 х=о а =1 х=о
Ограничения для этой целевой функции имеют вид
(8)
Ц Ха = Щ; г = 1..2;
а =1 2
Ц Хга - Ьа + и() - V() = Ц(); а = 1..3; 5 = о..9;
г=1
Ц Ц (V() - и(5) + ^^ ЦЦ Ц Ха) - Ьа = Ц Ца, - ЦЬ,
9 3^2
} 1 1 а =15=о 5=о а=о V г=1
Ха > о; иа > о; viJ > о;
2 3
у ^=о V г=1 а =1 у
т
и() < ) < Цаг - Ьа; г = 1..2;; = 1..3; 5 = о..9.
г=1
(9)
Следует отметить, что если величины а превышают суммарные потребности потребителей, то первое ограничение из (9) следует удалить. Если имеет место равенство
23
Ц Щ = I Ьа ,
г=1 а=1
то правая часть третьего ограничения из (9) равна нулю. Если имеет место неравенство
23
Ц а > Ц Ьа ,
г=1 а=1
то третье ограничение их (9) следует удалить. Оно сохраняется только в том случае, если имеет место неравенство
23
Ц а < Ц Ьа .
г=1 а=1
п
Математическое ожидание стоимости доставленных грузов с учетом потерь будет даваться выражением
2 9
ь;=I ^ -I ^ Р] ■ (10)
i=1 5=0
Проведем численный сравнительный анализ классической и стохастической моделей. Два поставщика располагают продукцией в размере а1=300; а2=700. Потребности потребителей составляют Ь1=250; Ь1=350; Ь1=400. Задана стоимостная матрица перевозок сг Классическое решение транспортной задачи имеет вид
хи ] = 1 ] = 2 ] = 3 Лп х н си ] = 1 ] = 2 ] = 3
г = 1 250 0 50 300 1 = 1 0,05 0,1 0,15
г = 2 0 350 350 700 г = 2 0,15 0,08 0,2
Итого 250 350 400 1000
Целевая функция (минимальная стоимость расходов на транспортировку) составляет
т п
С = Цсихи = 118.
i=l ]=1
Рассмотрим теперь несколько вариантов решения этой задачи с помощью целевой функции (8) с ограничениями (9). Предположим сначала, что вероятности утраты груза равны нулю (#¿,=0). Тогда при нулевых штрафных санкциях
* Г = * Г = 0
решение будет иметь вид
хи 1 = 1 1 = 2 1 = 3 Цп х н
1 =1 300 0 0 300
1 =2 0 700 0 700
Итого 300 700 0 1000
Целевая функция (9) и стоимость транспортных расходов С принимают одно и то же значение
тп
Е (2) = С = IIcjXj = 71. ¿=1 }=1
Для классической транспортной задачи полученные значения х, не являются решениями, поскольку не выполняется одно из основных условий
т
Е хи =ь} ■
¡=1
Эквивалентом этого условия в задаче (8), (9) является наличие штрафа за неполную поставку > 0. Например, если £/+) = 0, то стохастическое решение практически совпадает с классическим решением
т п
Е (2) = с = ЕЕсх =111 ■
¡=1 ]=1
1 = 1 1 = 2 1 = 3 1Мп II
1 = 1 250,05 0 49,95 300
1 = 2 0 700 0 700
Итого 250,05 700 49,95 1000
При больших значениях g/ ) > 0,1; §(+) = 0 стохастическое решение полностью совпадает с классическим решением
тп
Е(2) = с = ЕЕсх = 118 ■
¡=1 }=1
Можно показать, что при наличии штрафов за излишне поставленную продукцию (§(—) = 0; g/+) > 0,1) или (§(—) > 0; §(+) > 0) стохастическое решение и классическое решение будут полностью совпадать.
Рассмотрим теперь задачу управления риском потери грузов. Утрата груза на маршрутах характеризуется значениями qij и /ц = /(к) с вероятностями г(к), (к=1..3).
1 =1 1 = 2 1 = 3 Х(к )
1=1 0 0 0,6 х(к) 0,25 0,5 0,75
1=2 0 0 0,6 г(к) 0,5 0,3 0,2
Решение системы (9) при g/ ) = §(+) = 0,2 имеет вид
хи 1 = 1 1 = 2 1 = 3 1Мп II
1 = 1 300 0 0 300
1 = 2 0 700 00 700
Итого 300 700 0 1000
Целевая функция Е( 0=231, все маршруты в третий пункт исключены, потери груза нулевые. Если же увеличить штраф за недоставку груза в третий пункт ^12(—) = 5^1,2,з(+) = 0,2; gз(—) = 0,7), то решение опять становится классическим (тре-
бование максимально обеспечить поставку), при этом расходы возрастут до величины Е(0=291,4 за счет потери груза в размере 102 единиц.
Наконец, рассмотрим задачу управления риском в условиях избыточного предложения. В реальных рыночных условиях, как правило, предложение превышает спрос.
Пусть каждый из поставщиков располагает продукцией, достаточной для удовлетворения спроса всех потребителей (а1 = 1200; а2 = 1200), а потребности потребителей по-прежнему составляют значения (Ь1 = 250; Ь1 = 350; Ь1 = 400). Параметры д^; g/—); g/+) примем те же, что и в предыдущей задаче. В этом случае первое и третье ограничения из (9) необходимо убрать. Решение будет иметь вид
1 = 1 1 = 2 1 = 3 н
г=1 250,00 0,00 400,00 650,00
г=2 0,00 350,00 0,00 350,00
Итого 250,00 350,00 400,00 1000,00
Целевая функция Е( 0=231, потери груза составляют 102 единицы. Все грузы в пункт В3 доставлены из пункта А1, поскольку эти транспортные расходы меньше, чем при доставке из пункта А2.
Предложенная методика может быть использована в расширенном варианте для произвольных условий перевозки при наличии риска потерь и без применения классических условий соответствия запасов поставщиков потребностям потребителей.
Библиографический список
1. Фалин Г.И. Математический анализ рисков в страховании. М.: Российский юридический издательский дом, 1994. 130 с.
2. Современная актуарная теория риска / Р. Каас [и др.]; пер. с англ. А.А. Новоселова; под ред. В.К. Малиновского. М.: Янус-К, 2007. 372 с.
3. Крамер Г. Математические методы статистики. М.: Мир, 1975. 648 с.
4. Инкотермс. Международные правила толкования торговых терминов. М.: Омега-Л, 2009. 80 с.
5. Гольштейн Е.Г., Юдин Д.Б. Задачи линейного программирования транспортного типа. М.: Наука, 1969. 384 с.
6. Юдин Д.Б. Математические методы управления в условиях неполной информации. М.: Сов. радио, 1974. 400 с.
V.N. Nikishov, L.A. Saraev * LOSS AND RISK MANAGEMENT OF CARGO FORECASTING
In the article a mathematical model of cargo transportation risk management is considered. This model is a generalization of the classical transportation problem, taking into account the risks of loss of goods during delivery. The market situation, the excess of supply over demand have been examined. The numerical analysis of the results of calculations is in the form of the solution depending on the parameters of the problem.
Key words: supplier, consumer, objective function, constraints, stochastic transportation problem, risk management, beta distribution.
* Nikishov Viktor Nikolaevich (TSh-sea05@yandex.ru), Saraev Leonid Alexandrovich (saraev@ssu.samara.ru), the Dept. of Mathematics and Business-Informatics, Samara State University, Samara, 443011, Russian Federation.