Моделирование и оптимизация управления основной деятельностью порта

 

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

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

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

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

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

Процесс выполнения плана грузопереработки может рассматриваться как не стационарный процесс, так как если, в первом приближении предположить что плановое задание по объему грузопереработки Qпл (t) на момент t равно математическому ожиданию выполненного объема грузопереработки к моменту времени t , то

mq (t) = Qпл (t ) ¹ const (1)

где mq ( t ) - математическое ожидание выполненного объема грузопере­работки к моменту времени t .

Ожидаемое состояние выполнения плана грузопереработки к момен­ту времени t можно характеризовать вероятностью выполнения объема работы Q к этому моменту - P(Q). Система оперативного слежения фиксирует в момент времени t фактически выполненный объем грузо­переработки Qt , а вероятность выполнения планового задания определяется вероятностью выполнения оставшегося объема грузопереработки за оставшееся время tпл – t .

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

1. Своевременно не будут завезены или вывезены грузы из порта.

2. Произойдет отказ перегрузочного оборудования.

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

2. Отставание от заданного плана-графика грузопереработки в момент с оказалось настолько большим, что не может быть скомпен­сировано ускорением работ (путем увеличения концентрации и произ­водительности механизированных линий) за остаток планового периода tпл – t (2)

Рис.1

Введем обозначения:

D Qt - отклонение от плана грузопереработки в момент времени t (см. рис. 1).

P (D Qt ) - вероятность компенсации отклонения D Qt за остаток планового периода tпл - t при условии наличия необходимого коли­чества трудовых и технических ресурсов, безотказной работы пере­грузочного оборудования и обеспеченности материалами для произ­водства вспомогательных работ;

Pt ( Nk) - вероятность отсутствия отказов перегрузочного оборудо­вания к моменту времени t , приводящих к невыполнению планово­го объема грузопереработки;

Po ( t, tпл ) - вероятность обеспечения материалами для производ­ства вспомогательных работ при переработке грузов к моменту вре­мени t1.

Все перечисленные причины в момент t являются взаимонезаменяемыми событиями, поэтому вероятность выполнения планового объема грузопереработки запишем выражением

Pt ( Qпл ) = P (D Qt ) * Pt ( Nk) * Po ( t, tпл ) (3)

где Po ( t, tпл ) - вероятность выполнения планового объема грузо­переработки, определяемая на момент t.

Остановимся на сомножителях, входящих в формулу № 3.

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

Uпл = Qпл / t. (4)

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

mq ( t ) = utпл * t. (5)

Если система работает с недогрузом или перегрузом, то интен­сивность грузовых работ величина переменная:

mq ( t ) = ut * t. (6)

Где ut - средняя интенсивность грузовых работ на интервале времени [0,t ].

В общем случае при задании интенсивности грузовых работ как функции времени u (t) величина mq (t ) определяется:

t

mq ( t ) = Qпл ( t ) = ò u ( t ) dt (7)

и в реальных условиях математическое ожидание выполненных работ к моменту t равно плановому объему грузопереработки на этот момент

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

D Qt == Qt - Qпл ( t ) (8)

где Qt - фактический объем грузопереработки к моменту времени t .

Для выполнения планового объема грузопереработки в срок в последующий период необходимо перейти к новой интенсивности гру­зовых работ:

u1 = uпл + Duп = uп -D Qt / (tпл – t ) (9)

Вероятность компенсации D Qt равна вероятности перехода к работе с интенсивностью u t .

Очевидно, что в каждый конкретный момент времени интенсивность грузовых работ может приобретать одну из множеств конкретных зна­чений .

Будем считать, что вид закона P ( u ) и его параметры известны. Тогда вероятность компенсации (выполнение плана) и перекомпенсации отклонения (перевыполнение плана) равна вероятности того, что фактическая интенсивность грузовых работ на оставшемся интервале времени t, tпл примет значение: umax

P ( Q ³ Qt ) = P (u ³ ut ) = ò P ( u ) du, (10)

ut

где umax - максимальная интенсивность грузовых работ.

Значительно сложнее обстоит дело с определением двух других сомножителей: Pt ( Nk ) и Po (t, tпл). Здесь несложно предложить формальный способ их определения. Так, для вероятности своевременного поступления материалов для производства вспомогательных работ, связанных с грузовыми Pt(Nk) если обозначить Pt (ni) - вероятность поступления n i материала до заданного момента времени t, Pi (qni) - вероятности выполнения плана поставок и Pt mp - вероятность безотказной работы транспорта, используемого для доставки материалов, можно с учетом независимости Pt ( qni ) и P imp найти

PI (ni)= Pt(Nk)* P imp .

В общем случае при независимости работы поставщиков вероят­ности своевременной поставки

k – го материала Pt(Nk) запишем::

k k

Pt(Nk)=Õ Pi (ni) = Õ Pt(Nk)* P imp (11)

i=1 i=1

где k

Nk = S ni. (12)

i=1

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

Исследуя сомножитель Po (t,tпл) выражения (2) очевидно, что если имел место отказ перегрузочного оборудования в начале вы­полнения месячного плана грузопереработки или такой же отказ прои­зошел в последний день, то по своим последствиям эти события явно неравноценны. В первом случае есть время произвести ремонт, и все таки выполнить план, а во втором случае при невозможности выпол­нения перегрузочных работ другими механизмами этой возможности уже, очевидно, нет.

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

Работая с максимальной интенсивностью грузовых работ umax план грузопереработки возможно выполнить за минимальное время:

Tmin = Qпл / umax (13)

Таким образом, если в начале месяца (работ) произошел отказ перегрузочного оборудования, продолжавшийся до момента времени t0, то при работе с umax план грузопереработки может быть выполнен при задержке на время

t0 = tпл – tmin = t – Qпл / umax. (14)

Если в момент времени t выполнен объем грузопереработки Qt мы можем задержать выполнение работ на время t1 = t, а затем перейти на работу с umax.

Другими словами, если в момент времени t произошел отказ пе­регрузочного оборудования, то при его восстановлении за время t1 = t имеется возможность выполнения планового объема грузопере­работки.

Очевидно, что чем ближе момент t к концу планового периода, тем меньше интервал времени t1 = t остается для ремонта оборудова­ния..

Допустимое время ремонта при сохранении возможности выполнения плана грузопереработки Dtв = t1 –t легко определить из графика, представленного на рис.

Dtв = tпл – t – ( Qпл – Q1 ) / umax , при t = 0,Qt = 0, (15)

Dtв = tпл - Qпл / umax. = t0 , при Qt = Qпл , Dtв = 0. (16) Таким образом, допустимое время ремонта (восстановления) оборудования представляет собой функцию времени:

tв= f ( t ); 0 £ t £ tв (17)

В приведенных рассуждениях рассмотрена лишь принципиальная сторона вопроса и установлено наличие связи между текущим момен­том времени t и допустимым временем tв которое возможно, при определенных параметрах механизированных линий (интенсивности) обеспечить выполнение плана грузопереработки. Отклонение в выполнении планового объема грузопереработки (рис. 1 пунктир) однозначно связано с допустимым временем вос­становления и определяется текущими отклонениями D Qв. Считывая тот факт, что как плановую, так и максимальную интенсивность пере­работки грузов следует понимать как усреднение множества их значе­ний на заданном интервале, допустимое время восстановления также следует понимать как среднее значение Dtв .

Не останавливаясь на определении значений umax. и Dtв как так эта задача травиальна, покажем один из возможных путей нахождения Po(t,tпл). Будем считать известным среднее время ремонта перегрузочного оборудования Dtр , тогда при том условии, что вероятность одновременного появления двух и более отказов оборудо­вания мала, вероятность отказа перегрузочного оборудования на ин­тервале ( t1 , tпл).

m

Q(t,tпл) = S qi(t,tпл),

I=1

где qi(t,tпл) - вероятность возникновения отдельного вида отказа перегрузочного оборудования, m - число видов отказов. План грузопереработки не будет выполнен при появлении одного отка­за со временем ремонта больше

_ _

допустимого времени восстановления Dtр>Dtв или нескольких отказов с суммарным временем ремонта

_ _ _ _

Dtр = Dtр1 + Dtр2 + …+ Dtрh>в , (18)

_

Здесь Dtв - приведенное время восстановления, учитывающее, что отказы могут возникнуть в разные моменты времени.

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

qo(t,tпл) =Skiqi(t,tпл ) (19)

где ki - коэффициент, принимающий значения 0 и 1 в за­висимости от соотношения Dtв и Dtp .

_ _

ó 0, при Dtв < D tр

ki = í _ _ (20)

î 1, при Dtв и Dtp

На рис. 1 видно, что при малых значениях t допустимое время восстановления оборудования велико, а число видов отказов, приводящих к невыполнению плана грузопереработки, мало.

С ростом t число видов таких поломок оборудования растет, но вероятность их возникновения уменьшается, так как уменьшается интервал времени tпл - t , в течение которого они могут возник­нут ть

Таким образом, вид функции q(t,tпл) может в конкретных ус­ловиях оказаться весьма сложной.

Вероятность отсутствия отказов перегрузочного оборудования, приводящих к невыполнению планового задания

Po(t,tпл) = 1 - q(t,tпл) (21)

_ _

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

Нельзя ли путем определенных упрощений определить P(Nk) и Po(t,tпл) обойдя недостаток информации о выполнении плана поставщиками поставки крепежных материалов и сложность опреде­ления.. Анализируя причины возникновения отказов перегрузочного оборудования, следует, заметить что для морских портов величина вероятности отказов этого оборудования весьма мала, кроме того, почти всегда имеется резервное перегрузочное оборудование, которое при отказе основного быстро вводится в действие, или, например, при выходе из строя перегрузочного оборудования на обрабатываемом судне, такое судно можно перешвартовать на другой причал порта. Поэтому можно допустить, что Po(t,tпл)» 1 , а при выходе из строя перегрузочного оборудования и отсутствии

резерва нужно вво­дить ожидаемое время ремонта D tр по данным службы механизации порта.

За время ремонта перегрузочного оборудования возникает отри­цательное отклонение: _

DQp = uпл * D tр , (22)

_

а в ожидаемый момент окончания ремонта t + D tр ожидаемое отклонение составит величину DQt -DQp , где DQt - отклонение, имевшее место в момент отказа перегрузочного оборудования. Тогда требуемая интенсивность грузовых работ с момента DQt аналогично определится по выражению:

DQt - DQp DQt + uпл * D tр

ut = uпл + ---------------- = uпл + -------------------. (23)

_ _

tпл- t - D tр tпл- t - D tр

Принятое допущение окажется Po(t,tпл)» 1 несправедливым лишь в отдельных случаях, когда перегрузочный комплекс представляет собой систему взаимосвязанных элементов со сравнительно высокой надежностью, а также не имеет резервов, что не характерно для морских портов. В этом случае необходимо определить q0(t,tпл) по формуле (19), а при возникновении отказа перегрузочного обору­дования рассчитать ut по (23).

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

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

Если учесть, что время простоя D tn из-за несвоевременного обеспечения грузовых работ крепежными материалами и соответ­ствующее ему дополнительное отклонение DQp, возможно опреде­лить по данным службы материального обеспечения порта. Тогда тре­буемая интенсивность работ определяется аналогично (23) , где вместо DQp и D tр вводятся соответственно DQn и D tn.

В общем случае требуемая интенсивность работ будет определяться выражением: _ _

DQt -DQp - DQt

ut = uпл + ------------------- . (24)

_ _

tпл - D tp - D tn

Таким образом, вероятность выполнения и перевыполнения планового задания, найденная в момент времени t , с принятым допущением Po(t,tпл)» 1 будет определяться выражением, аналогичным выражению (21):

 

umax

P ( Qпл) = P (u ³ ut ) = ò P ( u ) du, (25)

U1

где нижний предел интеграла описывается (24) , a umax пара­метр известного закона распределения вероятностей интенсивности работы P ( u ).

Зная требуемую интенсивность работ u (t) в момент t и закон P ( u ) , можно указать величену ожидаемого отклонения выпол­ненного объема грузопереработки от планового DQпл к концу пла­нового периода:

 

tпл

DQпл= ò mu ( t ) dt – ut ( ut ( tпл – t ). (26)

t

где mu ( t ) - математическое ожидание интенсивности грузопепеработки в момент времени t.

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

1. Фиксируется момент времени t .

2. Фиксируется выполненный объем грузовых работ Qt .

3. Определяется остаток грузопереработки Qn = Qt.

4. Определяется отклонение выполненного объема грузовых работ от планового в момент t - ± DQt.

5. Определяется остаток времени грузовых работ tпл - t.

_ _

6. Если введено Dtp, Dtn корректируется остаток времени работ.

7. Определяется требуемая интенсивность грузопереработки.

8. Определяется вероятность выполнения и перевыполнения планового задания

P ( Q ³ Qпл ) = P (u >> ut ).

9. Определяется прогнозируемое отклонение от плана к моменту окончания работ ± D Qпл .

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