CC BY
52
Вестник СПбГУ. Сер. 8. 2007. Вып. 4
О. И. Заяц, В. М. Макаров, С. В. Такина
СОВЕРШЕНСТВОВАНИЕ МОДЕЛЕЙ УДЕЛЬНОГО ДЕФИЦИТА И УДЕЛЬНОГО ОСТАТКА ЗАПАСА НА СКЛАДЕ ДЛЯ РЕШЕНИЯ ЗАДАЧ ЛОГИСТИЧЕСКОГО МЕНЕДЖМЕНТА
В статье исследуется типовая стохастическая модель управления запасом, основанная на вычислении уровня дефицита. Выявлены ее принципиальные недостатки: независимость от коэффициента вариации потребления и уровня остаточного запаса. Предложена усовершенствованная модель, учитывающая указанные факторы. На ее основе разработан метод оптимизации резервного запаса на складе. Критерием оптимальности служит сумма потерь от дефицита и затрат на хранение. Решение доведено до числовых результатов, представленных в виде таблиц оптимальных параметров управления складом. В статье указаны границы применимости метода и приведен пример расчета.
Рациональное управление запасами является одним из ключевых вопросов производственного и операционного менеджмента. Запасы всегда возникают там и тогда, где и когда нарушается непрерывность движения материальных потоков. Исключить их появление в многономенклатурном дискретном материальном производстве различных отраслей промышленности, строительства и т. д., а также в каналах распределения произведенных товаров практически невозможно. Но и принижать их роль в качестве демпфера между смежными звеньями производства, имеющими различающиеся характеристики функционирования, нельзя. Вместе с тем хорошо известно, что запасы — слишком «дорогое удовольствие», так как их рост отрицательно сказывается на финансовых показателях деятельности предприятий и корпораций [Вумек, Джонс, 2004]. Налицо двойственность влияния запасов на эффективность производства, а следовательно, и актуальность отыскания их оптимального уровня.
Первая задача оптимизации применительно к управлению запасами была решена почти сто лет назад. Она позволяет на основе идеальной модели управления рассчитать наилучшие параметры пополнения запаса на
© О. И. Заяц, В. М. Макаров, С. В. Такина, 2007
складе. Формулу расчета оптимальной партии поставки для идеальной модели в зарубежной литературе обычно называют формулой Уилсона, или Economic Order Quantity (EOQ) [Джонсон и др., 2004; Чейз, Эквилайн, Якобс, 2001]. Одним из основных признаков идеальности (или условности) такой модели является то, что запас пополняется мгновенно при снижении остатка ресурса на складе до нуля. И хотя одна из разновидностей идеальной модели — модель планируемого дефицита [Долгов, Уваров, 2003; Модели и методы..., 2007] — позволяет рассчитать оптимальную партию даже в условиях разрешенного дефицита, тем не менее при таком допущении решение вопросов о времени заказа очередной партии и создании на складе резерва для недопущения дефицита не представляется возможным.
До 60-х гг. прошлого века считалось, что потери от дефицита в рыночной экономике слишком велики в сравнении с издержками на содержание запаса и он должен быть по возможности исключен. Но с развитием нового направления исследований в области менеджмента — логистического менеджмента (логистики) — стала очевидна несправедливость этого утверждения: оценка бизнесом затрат и потерь, обусловленных наличием запасов, поднялась до уровня, вполне сопоставимого с уровнем потерь от дефицита. Тогда возникла еще одна задача оптимизации управления запасами: определить допустимый уровень дефицита с точки зрения экономически эффективного функционирования склада. В настоящее время это одна из актуальных проблем теории и практики управления запасами.
Эта задача в принципе не разрешима методами детерминированной математики, в рамках которой выведена формула EOQ, если учесть случайный характер потребления ресурса и, соответственно, случайность возникновения дефицита. Только вероятностный подход позволяет перейти к моделям, на практике оптимизирующим соотношение «запас — дефицит». Вероятностная постановка задачи управления запасами отражает факт независимости спроса на складируемый ресурс от менеджмента предприятия или корпорации, которой принадлежит склад. Задача менеджера здесь — управлять поставками ресурса на склад, поддерживая на нем оптимальный уровень запаса.
ПОСТАНОВКА ЗАДАЧИ
В современных литературных источниках, где рассматриваются вопросы управления запасами (см., напр.: [Долгов, Уваров, 2003; Чейз, Эквилайн, Якобс, 2001; Букан, Кенигсберг, 1967] и др.) наибольшее распространение получила следующая реальная (в отличие от идеальной) постановка задачи управления запасами. Ресурс поступает на склад крупными партиями объемом ппост, а потребляется мелкими или поштучно. Время t меняется дис-
кретно, принимая целочисленные значения. Интенсивность потребления ресурса It для любого учетного интервала является случайной величиной. Значения It при различные t независимы и распределены по нормальному закону с одним и тем же математическим ожиданием MI и среднеквадрати-ческим отклонением öI. Новая партия ресурса заказывается, когда уровень наличного ресурса опускается до точки заказа НТЗ (модель с фиксированной партией поставки) или время заказа наступает с установленным ритмом RnocT (модель с фиксированным ритмом поставки). Время поставки партии TnocT фиксировано и не равно нулю.
Классические модели управления запасами в стохастической постановке базируются на задании такой характеристики, как вероятность бездефицитной работы склада Р0. На ее основе требуется рассчитать все параметры управления запасом: точку заказа, резервный запас, емкость склада, партию поставки и другие, характерные для различных моделей. Чем выше желаемая вероятность бездефицитной работы, тем большим должен быть резерв, а с ним и средний уровень хранящегося запаса, и затраты, обусловленные его хранением.
Вместе с тем с ростом Р0 снижаются объем возникающего дефицита и потери, обусловленные им. Наличие двух противоположных тенденций изменения затрат (потерь) при изменении величины Р0 приводит к необходимости решения оптимизационной задачи. Ее суть — поиск такой вероятности бездефицитной работы, при которой сумма затрат на хранение запаса и потерь ввиду его дефицита будет минимальна. Практически задача сводится к поиску оптимального уровня резервного запаса, а с ним и других параметров управления, зависящих от Р0. Как видим, задача оптимизации управления запасом получает экономическое содержание благодаря смыслу ее целевой функции.
Однако поставленная задача не может быть решена без дополнительного специального исследования целевой функции. Дело в том, что вероятность Р0 определяет только факт отсутствия (появления) дефицита, но не его величину, а следовательно, в целевой функции оказываются неопределенными потери, вызванные им. Таким образом, для решения задачи необходимо сначала установить зависимость между вероятностью бездефицитной работы склада, заложенной в основу управления запасом, и ожидаемым в результате такого управления дефицитом ресурса.
Вид этой зависимости в 60-е гг. XX в. предложил известный американский ученый Р.Браун [Brown, 1967]. Он выразил нормированный резервный запас Z (резервный запас, приходящийся на единицу среднеквадратичного отклонения), а также ожидаемый удельный дефицит E(Z) (средний дефицит, нормированный таким же образом) через вероятность бездефицитной работы склада Р0. В результате им была построена универсальная таблица,
связывающая Z и E(Z). Так как выразить Z через Р0 не составляет труда, это позволило на практике оценивать ожидаемый дефицит как функцию Р0. Указанный метод и сегодня рекомендуется без каких-либо изменений в наиболее авторитетных руководствах по логистике [Чейз, Эквилайн, Якобс, 2001; Ballou, 1991].
Чтобы решить, можно ли пользоваться этой таблицей, необходимо предварительно ответить на три вопроса:
1) при каких предположениях и каким образом была построена таблица?
2) все ли факторы, влияющие на решение, учтены при ее построении?
3) существуют ли какие-либо ограничения на область ее применения?
Отвечая на них, отметим следующее. Оригинальное решение [Brown,
1967] основано на гипотезе о нормальности It. Последняя допускает отрицательные значения It, причем вероятность того, что It < 0, зависит от коэффициента вариации у = ö ¡/M1 и быстро растет с ростом у. Это ставит под сомнение вывод выражения для E(Z), данный в [Brown, 1967]. Для его корректировки авторы предлагают, во-первых, с самого начала ограничить у сверху и, во-вторых, вести вычисления так, чтобы результат был более чувствителен к значениям у. С этой целью в основу всех расчетов закладывается распределение числа дней дефицита запаса на складе Тдеф, которое отсутствует в [Brown, 1967].
Сказанное вынудило авторов выполнить полный альтернативный вывод выражения для удельного дефицита. При этом было показано, что широко цитируемая таблица Брауна учитывает не все факторы, влияющие на решение. Далее в тексте приведена уточненная авторами таблица.
Что касается гипотезы нормальности потребления, выдвинутой Брауном, то, по мнению авторов, она действительно справедлива в некотором диапазоне изменения параметров спроса, границы которого также исследованы в настоящей статье. При этом показано, что таблица Брауна построена для недопустимого соотношения параметров спроса.
Вслед за уточнением расчета ожидаемого дефицита авторам для решения задачи оптимизации соотношения «запас — дефицит» пришлось также вывести выражение для оценки ожидаемого удельного остатка ресурса на складе к моменту поставки очередной партии (среднего остатка, приходящегося на единицу среднеквадратичного отклонения).
Конечная цель работы — доведение методики оптимизации резервного запаса до числовых результатов, позволяющих применять ее в практике конкретных логистических расчетов.
Отметим, что задачи оптимизации экономических показателей работы склада постоянно рассматриваются в научной литературе, но в иной поста-
новке. Так, в [Arrow, Harris, Marschak, 1951] строилось оптимальное управление штучным запасом, когда It является дискретной целочисленной случайной величиной. В 60-е гг. прошлого века было получено решение оптимизационной задачи, которую сегодня часто называют «проблемой разносчика газет» [Кофман, 1966; Черчмен, Акоф, Арноф, 1968]. Эта и подобные ей задачи в той или иной степени являются развитием идеальной модели управления запасами с разрешенным дефицитом. Сегодня это направление теории управления запасами активно развивается и в работах отечественных авторов [Долгов, Козлов, Уваров, 2004; Модели и методы..., 2007; Исследования операций в экономике, 2006].
Идеи классической постановки задачи и ее решения Р. Брауном используются в современных зарубежных разработках по логистике, таких как [Бауэрсокс, Клосс, 2001; Сток, Ламберт, 2005]. Однако развитие и уточнение решения классической задачи и по сей день остаются весьма актуальными.
МЕТОД РАСЧЕТА ОЖИДАЕМОГО ДЕФИЦИТА
Будем считать, что известны параметры случайного спроса на ресурс, хранящийся на складе, — MI и ар а также время поставки — Тпост и задана желаемая вероятность бездефицитной работы Р0. Тогда могут быть рассчитаны основные параметры управления запасом:
♦ точка заказа (Нтз) в модели с фиксированной партией поставки:
нтз=м 1тП0СТ.+)с /./тосг, (1)
где MI Тпост — ожидаемое потребление ресурса за время пополнения запаса; а /V ТП0СТ — среднеквадратичное отклонение этой случайной величины; P0 — вероятность того, что эта случайная величина примет любое значение, не превышающее Нтз, т. е. дефицита не будет; ^(P0) — квантиль порядка Р0 стандартного нормального закона распределения, который отыскивается по таблицам функции накопленной вероятности для нормального закона (см. формулу (7));
♦ условная емкость склада Н*кл в модели с фиксированным ритмом (обозначения даны по: [Макаров, 2006]):
Сл. = M I (пост. + «ПОСТ.) + ^(P0 )а ¡у1 тпост. + «пост.- (2)
где MI (ТП0СТ. + «ПОСТ. ) и а/V ТП0СТ. + «ПОСТ. — параметры распределения
случайной интенсивности потребления ресурса со склада за время (Т + R )
( пост. пост.).
Принципиальный подход к решению поставленной задачи одинаков для обеих моделей управления, поэтому рассмотрим его только для первой
из них — модели с фиксированной партией поставки. Более подробно с обеими моделями можно познакомиться в работах [Джонсон и др., 2004; Долгов, Уваров, 2003; Макаров, 2006; Чейз, Эквилайн, Якобс, 2001].
Обозначим через Тдеф ожидаемое число дней дефицита ресурса на складе. Это случайная величина, принимающая целочисленные значения в интервале [0, Тпост]. Значение Тдеф = 0 означает полное отсутствие дефицита, а Тдеф = Тпост говорит о том, что он возник сразу же в первый день (час, месяц, неделю — в зависимости от дискретности планирования; далее будем говорить — день) после пополнения запаса. Закон распределения этой случайной величины задается вероятностями:
С * = Р(ТДеф. = к },0 < к < Тпост, (3)
Однако их удобнее выразить через вероятности того, что дефицит продлится не более к дней:
Р* = Р{Тдеф.< к}, (4)
так как последние проще вычисляются и в то же время связаны с интересующими нас вероятностями очевидным соотношением:
Р = Е (5)
1=1
Случайное событие (Тдеф < к) Эквивалентно тому, что первые (Тпост - к) дней после заказа очередной партии для пополнения запаса потребление ресурса не достигло уровня Нтз, т. е. склад работал без дефицита. Поэтому вероятность отсутствия дефицита за к дней до поступления заказанной партии (4) выразится в следующем виде:
р = р(х) - ^ 1 , 0 < к < Тпост., (6)
где Р(х) есть функция накопленной вероятности:
* 12
Р(х) = ^= | е"т(7)
значения которой берутся из соответствующей таблицы. При к = 0 формула (6) задает вероятность бездефицитной работы склада Р0 , при к =1 — вероятность возникновения дефицита не ранее предпоследнего дня, т. е. Р1, и т. д.
Аргумент функции Р в выражении (6) — это квантиль ^ (Р*), фактически показывающий остаток ресурса на складе при его среднем потреблении в течение (Тпост - к) первых дней работы склада после заказа очередной партии, или, по определению, резервный запас, выраженный в долях среднеквадратичного отклонения за тот же срок. При к = 0 его выражение переходит в (1).
Искомые вероятности Ск легко выражаются через вероятности Рк с использованием соотношения (5):
С* =
Ро, к = 0,
Р* - Р*-1, 0 < к < (8)
1 - Рпост -1, к = ТПОСТ. •
Выражение (8) при к = Тпост следует из того, что Рк = 1 при к = Тпост, поскольку Тдеф. не превосходит Тпост..
По определению математическое ожидание времени дефицита ресурса на складе (в днях) может быть рассчитано следующим образом:
т
пост.
м = Е Юк. (9)
к=1
Это выражение нужно преобразовать, подставив в него правую часть (8). Результат после приведения подобных членов имеет вид:
Т -1 Т -1
пост . пост.
МТдаф = ТПОСТ.- Е Р* = К1 - Р*). (10)
*=о *=о
Ожидаемый размер дефицита Мн ^ при заданном Р0 тогда легко рассчитывается путем умножения ожидаемого времени дефицита на математическое ожидание интенсивности спроса:
МНдаф . = МI ■ МТ№ф = МI К1 - Р* ), (11)
*=о
а выражение для удельного дефицита Е^) может быть записано в виде:
Т -1
пост
М м м Т 1 Е (1 - Р*)
Е (I) = Мн = Тдеф . = 1 • *=0 , (12)
а IV тпост. аIV тпост. ^ \/тпост.
где у = а 1/М1 — коэффициент вариации интенсивности потребления товара со склада (или спроса на товар); характеризует степень близости фак-
тических значений к математическому ожиданию и определяется средне-квадратическим отклонением потребления от среднего, выраженным в долях этого среднего.
Численная зависимость удельного дефицита от Р0 и значения у, выступающего здесь в качестве параметра, приведена в табл. 1 [Макаров, Заяц, Семенова, 2005]. Дополнительно в таблицу введен столбец со значениями ^(Р0), (Z в обозначениях Р. Брауна).
Сравнение с оригинальным решением [Brown, 1967] показывает, что столбцы Z и E(Z) для у = 1 (табл. 1) в пределах точности вычислений повторяют таблицу Р. Брауна с той лишь разницей, что в [Brown, 1967] таблица составлена исходя из изменения с заданным шагом величины Z, а не Р0. Это подтверждает правильность предложенного авторами подхода.
Таблица 1
Удельный дефицит Епри различных значениях коэффициента
вариации у
Ро $(Ро) или Z Y
0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,00
0,10 -1,28 1,7887 1,5485 1,4670 1,4243 1,3970 1,3774 1,3621 1,3496 1,3389 1,3295
0,15 -1,04 1,5541 1,3229 1,2452 1,2046 1,1787 1,1602 1,1457 1,1339 1,1238 1,1149
0,20 -0,84 1,3725 1,1510 1,0772 1,0388 1,0143 0,9968 0,9832 0,9720 0,9625 0,9542
0,25 -0,67 1,2215 1,0102 0,9404 0,9041 0,8811 0,8646 0,8519 0,8414 0,8325 0,8247
0,30 -0,52 1,0905 0,8900 0,8242 0,7902 0,7686 0,7532 0,7413 0,7315 0,7232 0,7159
0,35 -0,39 0,9739 0,7846 0,7229 0,6911 0,6710 0,6567 0,6456 0,6365 0,6288 0,6221
0,40 -0,25 0,8680 0,6904 0,6330 0,6035 0,5848 0,5716 0,5613 0,5529 0,5458 0,5396
0,45 -0,13 0,7707 0,6052 0,5521 0,5249 0,5077 0,4955 0,4861 0,4784 0,4719 0,4662
0,50 0,00 0,6803 0,5273 0,4784 0,4537 0,4380 0,4270 0,4184 0,4113 0,4054 0,4002
0,55 0,13 0,5956 0,4556 0,4112 0,3885 0,3746 0,3646 0,3568 0,3505 0,3451 0,3405
0,60 0,25 0,5158 0,3891 0,3492 0,3292 0,3165 0,3077 0,3007 0,2951 0,2903 0,2863
0,65 0,39 0,4402 0,3273 0,2921 0,2744 0,2632 0,2555 0,2494 0,2445 0,2403 0,2367
0,70 0,52 0,3683 0,2697 0,2391 0,2238 0,2142 0,2075 0,2023 0,1981 0,1945 0,1915
0,75 0,67 0,2997 0,2158 0,1900 0,1772 0,1692 0,1636 0,1592 0,1557 0,1527 0,1502
0,80 0,84 0,2342 0,1655 0,1446 0,1342 0,1277 0,1232 0,1197 0,1169 0,1146 0,1125
0,85 1,04 0,1714 0,1186 0,1023 0,0947 0,0900 0,0864 0,0838 0,0817 0,0800 0,0784
0,90 1,28 0,1114 0,0750 0,0642 0,0588 0,0555 0,0532 0,0515 0,0501 0,0489 0,0479
0,95 1,64 0,0540 0,0350 0,0294 0,0267 0,0250 0,0239 0,0230 0,0223 0,0217 0,0212
0,99 2,33 0,0104 0,0063 0,0051 0,0046 0,0042 0,0040 0,0038 0,0037 0,0035 0,0034
ИССЛЕДОВАНИЕ ФАКТОРОВ, ВЛИЯЮЩИХ НА ВЕЛИЧИНУ ОЖИДАЕМОГО ДЕФИЦИТА
Расхождения в результатах выполненного расчета и первоисточника наблюдаются лишь в третьем-четвертом знаке после запятой. Можно предположить, что причина различий заключается в несовпадающей дискретности представления срока поставки. Дело в том, что Тпост — величина условная, которую можно представить разбитой на разное число учетных интервалов. Например, один месяц можно представить как четыре недели, двадцать рабочих дней, сорок рабочих смен или 320 рабочих часов.
Рост дискретности не меняет характера найденной зависимости, но в то же время позволяет повысить точность результатов, так как модель приближается к непрерывной [Джигарджян, Рубальский, 1988; Рубальский, 1977; Рыжиков, 2001]. Таким образом, отметим, что дискретность интервала управления — важный инструментальный, но не концептуальный фактор. Влияние количества учетных интервалов на результаты расчетов показано на рис. 1. Очевидно, что рост дискретности усложняет и удорожает выполнение учетно-плановых операций, но одновременно позволяет повысить точность управления запасами.
0,0245 ^ 0,0240 ^ 0,0235
е- 0,0230
V
£ 0,0225
g 0,0220
| 0,0215
¡5 0,0210>
0,0205 0,0200 0,0195
10 30 50 70 90 110 130 150 170 190
Количество учетных интервалов
Рис. 1. Зависимость удельного дефицита товара на складе от количества учетных
интервалов при Р0 = 0,95
Согласно формуле (12) результат находится в обратно пропорциональной зависимости от коэффициента вариации. Этот параметр в таблице Брауна отсутствует вовсе. Все расчеты выполнены им для у= 1.
Для практики производственно-логистических систем такое значение Y не характерно. Неопределенность потребления товара со склада в этом
случае оказывается очень большой, спрос — фактически непредсказуемым и не описываемым нормальным законом распределения. Тем не менее численно исследуем влияние на результат изменения параметра у в интервале [0,05; 1,0]. Оно показано в виде графика на рис. 2.
Влияние на решение задачи параметра у наиболее существенно на нижней границе интервала — в области, характерной для реальных структур. Это требует обязательного учета данного параметра при решении поставленной задачи. В то же время сказанное порождает новый вопрос: зависит ли от значения у величина МГдеф ? Если не зависит, то достаточно будет рассчитать и представить пользователю ряд значений функции МГ ^ = f(Pa) без учета изменения у по аналогии с тем, что было сделано Брауном для функции E(Z).
go,io
0,09
Я 0,08 К
f 0,07
st
« 0,06
а
g 0,05
а о,о4
>>
0,03 0,02 0,01 0,00
0,05 0,10 0,20 0,30 0,40 0,50 0,60 0,70 0,80 0,90 1,00
Коэффициент вариации
Рис. 2. Зависимость величины удельного дефицита от коэффициента вариации при
Р0 = 0,95
Для ответа на этот вопрос был проведен численный анализ зависимости МГ ^ = f(y) для ряда значений P0 как параметра [Заяц, Макаров, Такина, 2007]. Как видно из рис. 3, такая зависимость налицо, а, следовательно, без введения значений коэффициента вариации у в таблицу результатов не обойтись.
Таким образом, в ходе проведенного исследования подтверждена принципиальная правильность подхода к определению ожидаемого дефицита ресурса на складе в зависимости от вероятности его бездефицитной работы; предложены методы (либо один из методов) расчета величины такого
дефицита; уточнена таблица значений удельного дефицита путем ее существенного расширения за счет введения в рассмотрение еще одного параметра — коэффициента вариации спроса на ресурс.
ей
к
ЕГ а
■е-
ш ч
»Я ш № Ч О Ч и а ¡х1
4,5 4,0 3,5 3,0 2,5 2,0 1,5 1,0 0,5
Л, = 0,5
Л> = 0,6
Рот0'7
Рот 0,8
Ро^ 0,9
Ро у 0.99
0,0 0,2 0,4 0,6 0,8 1,0
Коэффициент вариации
Рис. 3. График зависимости МГдеф = /(у) для шести значений Р0
Остался открытым третий из поставленных вопросов — об области применения метода расчета удельного дефицита. Рассмотрим его подробнее.
О ГРАНИЦАХ ПРИМЕНИМОСТИ МЕТОДА РАСЧЕТА
ОЖИДАЕМОГО ДЕФИЦИТА
Р. Браун исходит из предположения о сохранении нормального закона потребления ресурса для любого временного интервала. Обоснованием этого служат центральные предельные теоремы (ЦПТ). Согласно ЦПТ, сумма случайных слагаемых при неограниченном увеличении их числа распределена по нормальному закону, каков бы ни был закон распределения каждого из них [Гнеденко, 2001]. Таким образом, применимость метода Брауна, а вместе с ним и метода, предлагаемого авторами, в первую очередь определяется выполнением условий ЦПТ в задачах управления запасами.
Поскольку ЦПТ требуют большого числа слагаемых, то на один учетный интервал (например, день) должно попадать достаточно много единичных заказов. Необходимое их число за день можно оценить, используя следующий факт: сумма двенадцати равномерно распределенных величин практически неотличима от нормальной величины [Ермаков, Михайлов, 1982]. Если распределение объема единичных партий отпускаемого ресурса отличается от равномерного, то указанное критическое число партий не-
сколько изменяется и колеблется от 10 до 15 штук в день. При меньшей активности клиентов для применения ЦПТ нужно укрупнять учетные интервалы до двух дней, недели и т. п.
Приведенное ограничение решает проблему нормальности, если ни одна из отпускаемых за день партий ресурса не преобладает в общей сумме выдачи за день. В ЦПТ для исключения этого ставятся специальные условия, например условие Линдеберга [Гнеденко, 2001]. Их физический смысл прост: дисперсия любого единичного заказа должна быть существенно меньше дисперсии всего ежедневного потребления.
В складских задачах контингент потребителей часто неоднороден (индивидуальные потребители, мелкие, средние, крупные фирмы и т. п.). Соответственно, не исключается доминирование небольшой части клиентов в общем обороте склада. Если численность доминирующей группы меньше указанного выше критического значения, то нормальность потребления следует проверить численно методами математической статистики. Используются выборка значений ежедневного потребления за предшествующий период и критерии согласия. Обычно применяют критерий % , однако в специальной литературе подробно описаны и более мощные критерии нормальности [Степнов, 1985].
Нормальность распределения потребления является необходимым, но отнюдь не достаточным условием применимости рассматриваемого метода. Ранее уже отмечалось, что прогнозирование спроса, а следовательно, и управление запасами может быть успешным только для сравнительно небольших значений коэффициента вариации потребления у. Докажем, что это условие принципиально важно при нормальной аппроксимации закона потребления. Для этого обратимся к рис. 4, на котором приведены графики плотности вероятности /(х) ежедневного потребления ресурса со склада.
Рисунок 4а отвечает большему коэффициенту вариации у= 1, для которого, как было показано выше, и строил свою таблицу Р. Браун. Видно, что в положительной области нормальная аппроксимация дает здесь приемлемый результат, однако она допускает отрицательные значения интенсивности потребления, причем отрицательный «хвост» распределения быстро увеличивается с ростом у.
Вероятность 5 отрицательных значений потребления определяется равенством:
5=р { < е
0 (х - М I
(х - МI )2
dx
(13)
и выражается через функцию накопленной вероятности в виде:
5 = Р(-1/у) = 1 - F(1/Y).
(14)
Рис. 4. Точный закон распределения потребления ресурса (сплошная линия) и его
аппроксимация (пунктир)
Численная зависимость S от у представлена в табл. 2.
Таблица 2
Зависимость S от у
Y 0,248 0,269 0,323 0,429 0,781 1
S 0,00001 0,0001 0,001 0,01 0,1 0,159
Составлено по: [Митропольский, 1972].
Согласно этой таблице, результаты, полученные Р. Брауном, относятся к случаю, когда отрицательные значения потребления ресурса допускаются с вероятностью почти 16%, что является слишком грубым приближением к реальности. Если ограничить погрешность модели десятыми долями процента, то коэффициент вариации у не должен превосходить 30-40%, что вполне соответствует действительности.
Тем не менее в практике логистических расчетов иногда встречаются задачи, когда потребление либо не является нормальным, либо хотя и может считаться таковым, но характеризуется слишком большим значением у. Как решать такие задачи?
Ответ можно найти в работах самого Р. Брауна, который рекомендовал в подобных случаях применять метод Монте-Карло [Brown, 1959]. В этом методе на основе выборки значений потребления за предшествующий период вначале проводится оценка закона распределения f(x). На основе выявленного закона методами статистического моделирования [Ер-
маков, Михайлов, 1982] разыгрывается выборка ежедневного потребления ресурса за достаточно длительный период. Полученные реализации характеристик работы склада затем обрабатываются статистическим путем так, словно они наблюдались на практике. Искомые значения ожидаемого дефицита Мн ф получаются численно с помощью соответствующих компьютерных расчетов. Они заменяют рассчитанные аналитически с использованием выражения (11).
МЕТОД РАСЧЕТА ОЖИДАЕМОГО ОСТАТКА РЕСУРСА НА СКЛАДЕ
При кажущейся простоте получения окончательного решения задачи мы столкнулись еще с одной методической проблемой: само понятие «резервный запас» и методы его расчета, широко используемые в теории управления запасами, оказались, на наш взгляд, не всегда адекватными действительности. Резервным по определению считается тот запас, который остается на складе к моменту поступления туда очередной партии поставки при средней интенсивности потребления ресурса и расходуется — при интенсивности, выше средней. Однако это определение применимо лишь для детерминированных моделей управления запасами. В стохастических моделях резервный запас представлен вторым слагаемым в формуле (1), откуда явно видно, что его уровень прямо пропорционален величине ^(Р0). Тогда, если задать вероятность бездефицитной работы Р0 < 0,5, то рассчитанный по этим формулам резервный запас окажется отрицательным. Теоретически объяснить его смысл в таком случае можно (см., напр.: [Чейз, Эквилайн, Якобс, 2001]), но практически, и особенно с точки зрения учета затрат на его хранение, — нельзя. Поэтому рассмотрим подробнее формирование и методы расчета среднего запаса Нср.
В идеальной модели управления запасом его расчет очевиден; он вытекает из рис. 5а:
Нср. = «пост/2 + Нрез.. (15)
Видно, что средний запас складывается из оборотной и резервной частей.
Это справедливо и для детерминированных моделей управления запасом (рис. 56), в которых интенсивность потребления ресурса переменна, но ограничена жесткими рамками, поэтому дефицит невозможен. Средний запас на складе, усредненный по нескольким циклам потребления, здесь также складывается из резервной и оборотной составляющих, а его величина рассчитывается по той же формуле (15). Следует особо отметить, что усредненный оборотный запас не зависит от значения Р0 и при решении поставленной задачи остается неизменным, следовательно, на оптимизацию управления запасом влияния не оказывает. Поэтому его можно опустить и перейти к учету затрат на хранение только усредненного резервного запаса.
Рис. 5. Графики, иллюстрирующие варианты расчета среднего запаса на складе: а) для идеальной модели; б) для детерминированной модели с фиксированной партией поставки; в) для стохастической модели с фиксированной партией поставки
При стохастической постановке задачи дефицит возможен, а если Р0 < 0,5, то он наблюдается даже чаще, чем бездефицитная работа склада (рис. 5в). В этой ситуации оборотная часть запаса остается той же, что и в предыдущих случаях, и так же она может быть исключена из рассмотрения, но резерва не существует (по расчету он отрицателен). Однако в бездефицитных циклах потребления существует некоторый остаточный запас Ност , который фактически поднимает Нср выше уровня усредненного оборотного запаса пПост /2, т. е. выполняет функцию среднего запаса, в то время как включение в расчет непосредственно Нрез < 0 снижает Нср . Зададимся целью исследовать переменную НосТ и ее зависимость от Р0.
Будем считать, что Ност — это остаточный запас ресурса на складе к моменту поступления туда очередной партии поставки. Переменная Ност может принимать любые неотрицательные значения и является случайной
величиной. Она определенным образом связана со случайной величиной Тдеф , а именно: если дефицит имеет место (Тдеф > 0), то Ност = 0; если же дефицита нет (Тдеф = 0), то переменная Ност может быть только положительной (Ност > 0). Для уяснения закона распределения случайной величины Ност введем вероятность:
чх = Р{Ност.> (16)
где х > 0 — произвольное неотрицательное число.
Случайное событие Ност > х Эквивалентно тому, что потребление за время поставки было меньше, чем (Нтз - х). Используя нормальность закона потребления, по аналогии с (6), получаем:
Чх = Р
'Нтз - тпостМ7 -
о
(17)
Сравнивая это выражение с (6), обнаруживаем, что q0 = Р0, т. е. вероятность фактического наличия ненулевого остатка совпадает с вероятностью отсутствия дефицита. Вероятность отсутствия остатка является дополнением по отношению к q0 и выражается в виде Р{Ност = 0} = 1 - Р0, т. е. совпадает с вероятностью наличия дефицита. Это полностью совпадает с рассуждениями, представленными выше при введении в рассмотрение переменной Ност .
Сказанное также свидетельствует о том, что случайная величина Ност относится к смешанному дискретно-непрерывному типу. Дискретное нулевое значение она принимает с ненулевой вероятностью 1 - Р0, а в области положительных значений ведет себя как непрерывная случайная величина, распределенная по нормальному закону с математическим ожиданием (Нтз - Тпост М1) и среднеквадратическим отклонением о^ ТПост. . Следовательно, остаток Ност распределен по усеченному нормальному закону. Ряд преобразований позволяет получить выражение для расчета его математического ожидания:
( Н тз-ТдостМ I )2
Мн„ =(ЯТЗ -ГпостМ1) е . (18)
ост .
Видно, что резервный запас, традиционно использующийся в расчетах, представлен здесь только в виде первого слагаемого и, более того, взят с вероятностью Р0. Это логично, так как с уменьшением вероятности бездефицитной работы значимость резервного запаса в формировании среднего остатка снижается. Второе слагаемое, формирующее величину остатка, имеет весьма оригинальный вид. Исследуем его поведение при изменении Р0, а также его влияние на величину остатка.
С учетом (1) выражение (18) можно записать в более удобном виде:
М
н „,
= а />/Тп(
>р
)
2 Л
(19)
Так как вынесенное за скобки выражение а / Тпост. не зависит от Р0 и не равно нулю, разделим на него правую и левую часть равенства (19). Получаемое в левой части отношение R(Z) = МНост/а / V Тпост. — это удельный остаток, выраженный числом среднеквадратичных отклонений, являющийся безразмерной функцией Z (или Р0). Результаты его расчета представлены в виде графиков на рис. 6.
3 3,5
- первый член второй член сумма
0,4 0,6 0,8 1,0 Вероятность бездефицитной работы
Рис. 6. Графики, иллюстрирующие формирование удельного остаточного запаса на складе согласно формуле (19) и последующим рассуждениям
Предложенный и рассчитанный нами удельный остаточный запас товара на складе R(Z) аналогичен в концептуальном смысле введенному и рассчитанному Р. Брауном удельному дефициту товара Е(2).
Графики на рис. 6 показывают: 1) в области высокой вероятности бездефицитной работы склада кривые резервного и остаточного запасов
практически совпадают; 2) с ростом вероятности дефицита остаточный запас все в большей степени превышает резервный; 3) остаточный запас, в отличие от резервного, ни при каких значениях Р0 не уходит в область отрицательных значений и только асимптотически приближается к нулю.
Таким образом, поведение смоделированного удельного остаточного запаса в поставленной задаче более адекватно действительности, чем поведение резервного запаса, и это еще раз подтверждает, что именно он должен фигурировать в методике оптимизации складского резерва.
МЕТОД ОПТИМИЗАЦИИ РЕЗЕРВНОГО ЗАПАСА
Теперь обратимся непосредственно к задаче оптимизации резервного запаса. Как отмечалось, ее решение требует минимизации суммарных затрат, обусловленных наличием запаса, и потерь от дефицита ресурса на складе.
Прямо зависящие от вероятности бездефицитной работы склада уровень резервного, а с ним и уровень остаточного запаса определяют ожидаемые затраты хранения:
Зхр= h Мн , (20)
ност.
где h — затраты, обусловленные хранением на складе единицы товара в
течение всего заданного срока.
Эти затраты, кроме обычных складских расходов, включают также: затраты на страхование запаса, потери ввиду естественной порчи и убыли, потери неликвидов из-за непредсказуемости изменения спроса, потери прибыли в результате «замораживания» в запасах оборотных средств предприятия и др.
Находящаяся в обратной зависимости от вероятности бездефицитной работы склада нехватка товара определяет ожидаемые потери от дефицита:
Здеф= ^МНдеф, (21)
где g — потери от нехватки на складе единицы товара.
Состав и величина этих потерь в значительной степени зависят от того, о каком товаре идет речь. Если это исходные ресурсы: сырье, материалы, комплектующие изделия, покупной инструмент и оснастка, ГСМ и т. п., то их дефицит приводит к остановке производства и в первую очередь к потерям ввиду простоев. Как следствие, это может привести и ко вторичным потерям прибыли из-за срыва поставок готовой продукции.
Если речь идет непосредственно о дефиците готовой продукции, то очевидны потери прибыли в краткосрочном плане, однако больший вред дефицит может нанести предприятию в долгосрочном плане в результате потери имиджа надежного поставщика.
В качестве аргумента задачи оптимизации запаса можно использовать как саму величину Р0 или ^(Р0), так и точку заказа Нтз. Они связаны моно-
тонной зависимостью (1), и между ними существует взаимно однозначное соответствие. С математической точки зрения эти подходы равнозначны.
Используем все полученные ранее результаты для записи целевой функции задачи оптимизации:
3х = 3 хр.+3,
= к ■ М н + % • М н
' ОСТ. Д1
= (к ■ R(^) + %■Е^)) ^ Шп.
(22)
Явное выражение удельного дефицита и удельного остатка следует из (12) и (19):
Е (г ) =
1
к=0
1 - F
Т
г I пост.
+ ■
Тпост. - к ' Ул/ТП
R (г ) = г ■ Ро +
1 -— 1 е 2 ,
(23)
Лп
где нормированный резервный запас г тождественен ^(Р0). Можно доказать, что при Р0 = 0 производная Зх по Р0 отрицательна, а при Р0 = 1 — положительна, т. е. в интервале 0 < Р0 <1 существует минимум функции (22), зависящий от у как от параметра. Он и определяет решение задачи оптимизации.
Графическое решение задачи имеет вид, показанный на рис. 7. Здесь h = 225 руб./ед. цикл, g = 450 руб./ед., у = 0,3, Тпост = 64 дн., а масштаб затрат и потерь определен значением М1 = 100 ед./дн.
Таким образом, задачу оптимизации соотношения «запас — дефицит» можно считать решенной в методическом плане. Далее был поставлен вопрос о создании единой таблицы оптимальных решений, которая обеспечила бы простое, быстрое и удобное практическое использование предлагаемой методики.
Для этого необходимо привести целевую функцию к виду:
зх = а 17ТП0Г ■ g ■ (ж ■ R (г) + е (г)) ^ шш,
(24)
где ж = Щg (при g Ф 0). Назовем ж коэффициентом безопасности хранения запаса. Он тем больше, чем меньше удельные потери от дефицита. Постоянный множитель, вынесенный за скобки, не зависит от Р0 и у и только сдвигает кривую суммарных затрат вверх или вниз по оси затрат, не меняя оптимального решения.
С учетом сделанных предположений была рассчитана таблица значений оптимальной вероятности бездефицитной работы склада Р0опт- для практически интересных значений коэффициентов ж и у (табл. 3).
Т -1
пост.
ю Рн и
Я &
ё И
3
У &
сп
— потери от дефицита
— затраты хранения
— суммарные затраты
0,4 0,6 0,8 1,0
Вероятность бездефицитной работы
Рис. 7. Пример решения оптимизационной задачи
Таблица 3
Оптимальные значения вероятности бездефицитной работы склада (Р0ОПТ) при различных сочетаниях коэффициентов ж и у
ж У
0,1 0,2 0,3 0,4 0,5
0,25 0,8902 0,8490 0,8323 0,8229 0,8167
0,50 0,7781 0,7230 0,7028 0,6919 0,6849
0,75 0,6808 0,6249 0,6053 0,5950 0,5883
1,00 0,5998 0,5480 0,5303 0,5210 0,5150
1,25 0,5331 0,4868 0,4711 0,4629 0,4576
1,50 0,4780 0,4373 0,4234 0,4162 0,4115
1,75 0,4322 0,3965 0,3843 0,3779 0,3738
2,00 0,3937 0,3624 0,3516 0,3460 0,3423
2,25 0,3611 0,3335 0,3240 0,3189 0,3157
2,50 0,3332 0,3088 0,3003 0,2958 0,2929
2,75 0,3090 0,2874 0,2798 0,2757 0,2731
3,00 0,2880 0,2687 0,2618 0,2582 0,2558
3,25 0,2696 0,2523 0,2460 0,2427 0,2406
3,50 0,2533 0,2377 0,2320 0,2290 0,2270
3,75 0,2388 0,2247 0,2195 0,2167 0,2149
4,00 0,2258 0,2130 0,2083 0,2057 0,2041
Использование табл 3 существенно упрощает процесс принятия управленческого решения о величине оптимального запаса на складе. На практике для этого необходимо собрать информацию о величине затрат на хранение товара и потерь ввиду его дефицита, соотношение которых формирует значение коэффициента ж. На основе статистики спроса можно сделать прогноз математического ожидания и среднеквадратичного отклонения будущей интенсивности потребления товара со склада. Соотношение двух последних величин формирует значение коэффициента у.
Значение вероятности бездефицитной работы Р0ОПГ, находящееся на пересечении строки и столбца, соответствующих полученным ж и у, и будет определять точку оптимума. Суммарные затраты, представленные в виде целевой функции задачи (22), при этом значении Р0 будут минимальны. На основе полученного значения Р0ОПГ- рассчитываются все параметры управления запасом.
Значения Zoпт в рабочем диапазоне изменения безразмерных (относительных) параметров работы склада ж и у приведены в табл. 4.
Таблица 4
Оптимальный нормированный резервный запас Zoпт при различных сочетаниях коэффициентов ж и у
ж У
0,1 0,2 0,3 0,4 0,5
0,25 1,2276 1,0322 0,9633 0,9265 0,9029
0,50 0,7658 0,5918 0,5325 0,5012 0,4814
0,75 0,4699 0,3184 0,2671 0,2404 0,2232
1,00 0,2528 0,1206 0,0760 0,0527 0,0376
1,25 0,0831 - 0,0331 - 0,0725 - 0,0931 - 0,1065
1,50 - 0,0552 - 0,1578 - 0,1932 - 0,2116 - 0,2237
1,75 - 0,1708 - 0,2624 - 0,2942 - 0,3110 - 0,3218
2,00 - 0,2697 - 0,3521 - 0,3810 - 0,3961 - 0,4062
2,25 - 0,3555 - 0,4303 - 0,4565 - 0,4708 - 0,4798
2,50 - 0,4311 - 0,4993 - 0,5235 - 0,5365 - 0,5449
2,75 - 0,4987 - 0,5610 - 0,5834 - 0,5957 - 0,6035
3,00 - 0,5592 - 0,6167 - 0,6378 - 0,6489 - 0,6563
3,25 - 0,6140 - 0,6673 - 0,6871 - 0,6976 - 0,7044
3,50 - 0,6641 - 0,7137 - 0,7323 - 0,7421 - 0,7488
3,75 - 0,7102 - 0,7564 - 0,7739 - 0,7834 - 0,7895
4,00 - 0,7528 - 0,7961 - 0,8123 - 0,8214 - 0,8271
ПРИМЕР ОПТИМАЛЬНОГО УПРАВЛЕНИЯ ЗАПАСОМ
Использование предложенной методики для управления запасами на практике позволит получить ощутимый экономический эффект. Покажем это на примере.
Пусть ежедневное потребление ресурса со склада характеризуется параметрами: М1 = 100 ед./дн., а1 = 30 ед./дн., т. е. у = 0,3; удельные затраты хранения составляют h = 225 руб./ед. цикл, потери от дефицита g = 450 руб./ед., т. е. ж = 0,5; Тпост = 64 дня. Зададим желаемое значение Р0 = 0,95 без учета возможности оптимизации этого параметра. Тогда по табл. 1 находим нормированный резервный запас г = 1,64 и удельный дефицит Е(г) = 0,0294, рассчитываем по формуле (1) точку заказа Нтз = 6794 ед. и по формуле (12) ожидаемый дефицит за один цикл потребления Мн~ 7 ед. По формуле (19) может быть рассчитан ожидаемый средний остаток ресурса на складе Мн ~ 400 ед. Тогда значение целевой функции составит 93,15 тыс. руб.
Оптимальная стратегия управления, представленная в табл. 4, дает гопт = 0,5325 (Р0опт- « 70%). По табл. 1 отыскиваем Е(г) = 0,2391 и рассчитываем Мн ~ 57 ед. Расчет среднего остатка дает Мн ~ 171,5 ед. Хорошо
-Гсдеф . ост .
видно, что со снижением вероятности бездефицитной работы склада ожидаемый дефицит растет, но средний остаток снижается. В результате значение целевой функции уменьшается до 64,24 тыс. руб. (рис. 7). Эффект оптимизации значителен.
Чтобы показать экономическое преимущество предлагаемой методики по сравнению с методом Брауна, представим, что менеджер-логист задался целью обеспечить тот же уровень покрытия спроса [Чейз, Эквилайн, Якобс, 2001], а значит, и тот же объем дефицита, что и в предыдущем расчете при Р0опт- ~ 70%. Но при этом он обратился к таблице Брауна, чтобы определить вероятность Р0, которую ему надо заложить в основу управления запасом. Формально эти действия могут быть записаны так: на основе МНд ф = 57 ед., рассчитано Е(г) = 0,2375. По этому значению в табл. 1 в столбце у = 1 (другого у Брауна не было) получено: г ~ 0,39, Р0 = 0,65.
Чтобы оценить реальный удельный дефицит, который будет получен при таком управлении, возьмем значение Е(г) = 0,2921 из столбца у = 0,3 той же строки табл. 1. Ему соответствует Мн ф ~ 70 ед. Вероятности бездефицитной работы 65% соответствует Мн ~ 149 ед. В результате значение целевой функции увеличилось до 65,02 тыс. руб. Увеличение незначительное, но и оно свидетельствует о неэффективности решения, принятого на основе таблицы Брауна.
Дальнейшее исследование показало, что с ростом значения ж преимущество предлагаемой методики по сравнению с методом Брауна заметно возрастает.
Литература
Бауэрсокс Д. Д., Клосс Д. Д. Логистика: интегрированная цепь поставок. М.: Олимп-Бизнес, 2001.
Букан Дж., Кенигсберг Э. Научное управление запасами. М.: Наука, 1967.
Вумек Д., Джонс Д. Бережливое производство: Как избавиться от потерь и добиться процветания Вашей компании. М.: Альпина Бизнес Букс, 2004.
Гнеденко Б. В. Курс теории вероятностей. М.: Эдиториал УРСС, 2001.
Джигарджян О. Л., Рубалъский Г. Б. Модели управления запасами с винеровским спросом. М.: ВЦ АН СССР, 1988.
Джонсон Д., Вуд Д., Вордлоу Д., Мерфи-мл. П. Современная логистика. 7-е изд. М.: Вильямс, 2004.
Долгое А. П., Козлов В. К., Уваров С. А. Логистический менеджмент фирмы: концепция, методы и модели. СПб.: СПбГУЭФ, 2004.
Долгов А. П., Уваров С. А. Логистический менеджмент. Управление запасами. СПб.: СПбГУЭФ, 2003.
Ермаков С. М., Михайлов Г. А. Курс статистического моделирования. М.: Наука, 1982.
Исследование операций в экономике / Под ред. проф. Н. Ш. Кремера. М.: ЮНИТИ, 2006.
Кофман А. Методы и модели исследования операций. М.: Мир, 1966.
Макаров В. М. Логистика. Управление запасами в логистических системах. СПб.: СПбГПУ, 2006.
Макаров В. М., Заяц О. И., Семенова С. В. Метод оптимизации резервного запаса // Логистика и управление цепями поставок. 2005. № 5.
Заяц О. И., Макаров В. М., Такина С. В. Методика оптимизации резервного запаса на складе // Сборник статей германо-российской конференции по логистике «Ло-гистичекое взаимодействие» / Под ред. Д. А. Иванова и др. СПб.: СПбГПУ, 2007.
Митрополъский А. К. Интеграл вероятностей. Л.: ЛГУ, 1972.
Модели и методы теории логистики / Под ред. В. С. Лукинского. СПб.: Питер, 2007.
Рубалъский Г. Б. Управление запасами при случайном спросе (модели с непрерывным временем). М.: Советское радио, 1977.
Рыжиков Ю. И. Теория очередей и управление запасами. СПб.: Питер, 2001.
Степнов М. Н. Статистические методы обработки результатов механических испытаний. М.: Машиностроение, 1985.
Сток Д. Р., Ламберт Д. М. Стратегическое управление логистикой. М.: ИНФРА-М, 2005.
Чейз Р. Б., Эквилайн Н. Д., Якобс Р. Ф. Производственный и операционный менеджмент. М.: Вильямс, 2001.
Черчмен У. , Акоф Р., Арноф Л. Введение в исследование операций. М.: Наука, 1968.
Arrow K. J., Harris T., Marschak J. Optimal Inventory Policy // Econometrica. 1951. Vol. 19. N 3. P. 250-272.
Ballou R. H. Business Logistics Management. N. Y.: Prentice Hall, 1991.
Brown R. G. Statistical Forecasting for Inventory Control. N. Y. McGrow-Hill, 1959.
Brown R. G. Decision Rules for Inventory Management. N. Y.: Holt, Rinehart&Winston, 1967.
Статья поступила в редакцию 8 июня 2007 г.
