鄭文艷

德州學(xué)院信息管理學(xué)院，山東德州 253023

基于顏色Petri網(wǎng)的不確定庫存模型性能仿真

鄭文艷

德州學(xué)院信息管理學(xué)院，山東德州 253023

1 引言

分散型供應(yīng)鏈是基于局部信息的，并且具有這樣的特點，在供應(yīng)鏈的參與者之間缺少協(xié)調(diào)，所做出的決策經(jīng)常和不同類型的訂貨策略即向上游節(jié)點訂貨的數(shù)量相關(guān)聯(lián)。而供應(yīng)鏈管理的主要目標是整合供應(yīng)鏈中上下游之間的產(chǎn)品和信息流以避免庫存數(shù)量的過度波動。因此，對供應(yīng)鏈中的參與者來說產(chǎn)品和信息是必須要考慮的方面。仿真工具為供應(yīng)鏈的分析和驗證提供了一種有效的方法。

Lionel等[1]提出使用多目標優(yōu)化算法與遺傳算法相結(jié)合去優(yōu)化供應(yīng)鏈庫存策略。實驗表明所提出的庫存策略方法比目前現(xiàn)有的方法在減少庫存成本和提高服務(wù)水平方面有著更好的效果，但是最小化庫存的總成本以及運輸成本仍是需要研究的問題。曾強等[2]基于VB開發(fā)了實現(xiàn)(t,s,S)庫存優(yōu)化的軟件，模擬結(jié)果符合實際，由于顧客的實際需求不都是整數(shù)隨機分布，而可能是連續(xù)隨機分布，因此，軟件有待進一步研究與改進。王鵬飛[3]運用系統(tǒng)動力學(xué)方法建立了一對一的從制造商到終端客戶的兩極供應(yīng)鏈庫存控制模型，通過對其進行仿真分析和優(yōu)化，提出了參數(shù)及結(jié)構(gòu)優(yōu)化策略，而本文并未涉及一對多及多對多的多級庫存控制系統(tǒng)。Yang H C等[4]開發(fā)了一個虛擬生產(chǎn)控制系統(tǒng)，該系統(tǒng)位于制造商端，制造商根據(jù)短期需求計劃去調(diào)整供應(yīng)商的產(chǎn)品速度以便于降低庫存水平，還可以遠程監(jiān)控供應(yīng)商的生產(chǎn)進程。該系統(tǒng)尤其適合用戶控制他們的外包生產(chǎn)。Makaji?等[5]使用時間，分層的顏色Petri網(wǎng)建立了一個包括零售商、批發(fā)商和制造商的三級供應(yīng)鏈模型，模擬并進行了諸如牛鞭效應(yīng)，庫存水平的實時變化等方面的性能分析。Papanagnou等[6]提出了分析供應(yīng)鏈系統(tǒng)的模擬框架，使用顏色Petri網(wǎng)分析了三種不同的訂貨策略，由于供應(yīng)鏈是復(fù)雜的動態(tài)系統(tǒng)，因此需要進一步考慮現(xiàn)實供應(yīng)鏈的所有特征。

基于此，本文使用顏色Petri網(wǎng)結(jié)合不同的需求預(yù)測方法及(s,Q)庫存策略，設(shè)計決策過程并用仿真結(jié)果展示供應(yīng)鏈中各環(huán)節(jié)庫存實時狀態(tài)，庫存補充，配送及缺貨等各種現(xiàn)實系統(tǒng)中產(chǎn)生問題的主要部分，以及不同參數(shù)的設(shè)定與目標費用函數(shù)的關(guān)系，為供應(yīng)鏈的模擬和未來的仿真工作提供了新思路。

2 基礎(chǔ)知識

2.1 著色Petri網(wǎng)的定義

2.2 CPN的ML語言

在CPN Tools中使用CPN ML語言進行變量的聲明和Petri網(wǎng)的描述及定義。此處僅介紹函數(shù)類型。

函數(shù)的聲明語法：

其中exp1和exp2類型要相同。

3 供應(yīng)鏈的CPN模型及其描述

3.1 供應(yīng)鏈的功能描述

供應(yīng)鏈有五個參與者，顧客（customer），零售商（retailer），批發(fā)商或供應(yīng)商（supplier），分銷商（distributor）和制造商（manufacturer）。每一個參與者根據(jù)自己的局部情況做出決策，基于預(yù)測技術(shù)并遵循某一訂貨策略向上游產(chǎn)生訂貨訂單，向下游運輸配送產(chǎn)品同樣遵循特定的規(guī)則和條件。參與者在減少本身庫存（對下游提交的訂單產(chǎn)品進行配送）的情況下同時也在增加庫存量（根據(jù)下游訂單向上游提交新訂單）以便維持本身庫存水平。

圖1 供應(yīng)鏈流程圖

對每個參與者的幾個主要活動描述如下，限于篇幅此處僅以圖2所示的模型中retailer子頁為例進行說明：

（1）訂單的產(chǎn)生

分兩種情況，一是終端顧客產(chǎn)生的初始訂單；二是Retailer，supplier和distributor利用需求預(yù)測技術(shù)產(chǎn)生的訂單：

①終端顧客的需求行為從兩個方面進行限定，一是產(chǎn)生訂單的數(shù)量，二是前后兩次訂單間隔的時間。形成的訂單由三元組組成（訂單號，訂單產(chǎn)生的時間，訂單數(shù)量）。這些參數(shù)的值是可以按照實際需求進行修改的，并且很容易檢測當一個單一事件發(fā)生后對供應(yīng)鏈中其他節(jié)點諸如庫存水平等因素的影響。另外模型還允許其他類型的需求，比如季節(jié)性的浮動等均可以實現(xiàn)。

②根據(jù)下游提交的訂單形成新訂單即需求預(yù)測：由于配送下游訂單的產(chǎn)品會減少庫存量，為了保持庫存水平，減少由于缺貨造成的增加成本，參與者需要及時向上游訂貨以補充庫存。Retailer，supplier和distributor新訂單的數(shù)量均按以下兩種方式之一進行預(yù)測：各個階段加權(quán)平均值(MA)以及平滑指數(shù)分布(ES)[9-10]。而customer按照正態(tài)分布隨機產(chǎn)生訂單數(shù)量，manufacturer按固定量產(chǎn)生對原材料的需求。

（2）產(chǎn)品的配送

接收到訂單時，需要檢測當前庫存量是否滿足訂單需求，如果庫存顯示有足夠量的產(chǎn)品，那么訂單產(chǎn)品被一次性配送；否則，如果庫存量不足就有可能出現(xiàn)兩種情況，一種情況等待庫存量達到訂單數(shù)量后再進行配送，另一種情況是進行不完全配送，庫存量只要不為0就配送，直到訂單產(chǎn)品全部完成配送。即一個訂單的配送次數(shù)是多次的。本文CPN模型采取多次配送的方式。

圖2 retailer子頁

（3）庫存控制策略

為保障供應(yīng)鏈的正常運行，除customer外每個參與者均需要按照庫存控制策略保持自身的庫存水平。Retailer，supplier和distributor采取的庫存控制策略是根據(jù)上游提交訂單產(chǎn)品的配送及時補充；而manufacturer采取了(s,Q)[11]的庫存策略，當庫存量低于訂貨點s時，及時補充固定量Q的原材料。

3.2 目標函數(shù)說明

在庫存系統(tǒng)，涉及到訂貨費用（用CD表示），存儲費用（用CP表示），缺貨損失費用（用CS表示），目標實施總費用（用CT表示），計算公式表示如下[12-13]：

目標是使總費用CT最少。

具體說明如下：

（1）訂貨費用CD，包括兩部分，一部分是每次訂貨必須支付的與訂貨量無關(guān)的手續(xù)費用，用K表示，發(fā)生的訂貨次數(shù)用s表示；一部分是每件貨物需要支付的費用，用m表示，假定訂貨量為Z，則CD=K×s+m×Z。

（2）存儲費用CP，是指維持庫存所產(chǎn)生的費用，如果庫存量為0，則不需要計算該部分存儲費用。假定每件貨物單位時間內(nèi)存儲費用為h，每件貨物在倉庫存儲的時間為t，貨物數(shù)量為n，則CP=h×t×n。

（3）缺貨損失費用CS，是指供貨方發(fā)生缺貨行為時所產(chǎn)生的費用，假定每件貨物的缺貨損失費用為p，缺貨貨物數(shù)量為l，則CS=p×l。綜上，得：CT=CD+CP+CS=K×s+m×Z+h×t×n+p×l。

3.3 功能函數(shù)說明只配送訂單的一部分，配送數(shù)量為當前庫存的數(shù)量，同時庫存量更新為0。

（2）函數(shù)fun disp(i,j)=if i＞=j then j else i

函數(shù)功能描述：配送產(chǎn)品的數(shù)量。根據(jù)庫存量j和等待配送訂單隊列中隊首元素i判斷配送數(shù)量，如果庫存量j比等待配送的訂單數(shù)量i多，則此次配送數(shù)量為訂單數(shù)量i；否則，配送數(shù)量為當前的庫存量j。

函數(shù)功能描述：該函數(shù)用于更新列表信息，模型中為保證顧客訂單按產(chǎn)生的先后順序進行配送，并且保證一個訂單全部配送完后才開始配送下一個訂單，在兩個地方用到了該函數(shù)及列表，一個是等待配送的訂單號，另一個是等待配送的訂單數(shù)量。這樣才能保證訂單的正確配送。如果當前等待配送的訂單數(shù)量k小于當前的庫存量n，那么把該訂單號從列表中移除，并把訂單數(shù)量從等待配送的列表中移除；否則，如果庫存量不足以完成一次配送，那么訂單號繼續(xù)留在列表的表頭，同時更新配送列表中表頭的數(shù)量為k-n。這樣從訂單順序和訂單數(shù)量兩方面保證了配送的正確進行。

3.4 隨機變量

功能函數(shù)均使用ML語言實現(xiàn)：

（1）函數(shù)fun rest(i,j)=if(i＞=j)then i-j else 0

函數(shù)功能描述：根據(jù)接收的下游訂單數(shù)量j和自身的庫存量i去更新庫存，如果庫存量i大于等于訂單數(shù)量j，那么將對訂單進行配送，庫存量更新為i-j，否則

假設(shè)模型沒有信息交換和信息處理的延遲，即信息交換和處理是及時的，是不消耗時間的。整個模型中只有產(chǎn)品的配送運輸是要消耗時間的。

本模型涉及以下幾個隨機變量：

（1）顧客需求量隨機變量：即顧客產(chǎn)生的訂單數(shù)量，服從均值μ=100.0標準偏差σ2=8.0的正態(tài)分布，其概率密度函數(shù)為：

（3）貨物配送延遲時間隨機變量，指從訂貨到貨物入庫所需要的時間，根據(jù)統(tǒng)計數(shù)據(jù)，各階段的訂單配送延遲時間服從[1.4,4.2]之間的均勻分布，其分布函數(shù)為：

4 模型仿真及性能分析

模型模擬仿真1 500步，耗時51個時間單位?，F(xiàn)從以下幾個方面展示實驗結(jié)果。

（1）采取MA或ES需求預(yù)測技術(shù)前后的對比圖

圖3 未采取預(yù)測策略的牛鞭效應(yīng)圖

圖4 采取MA或ES預(yù)測技術(shù)后參與者產(chǎn)生訂單數(shù)量的分布圖

同樣的顧客需求條件下，沒有采取任何策略而產(chǎn)生的牛鞭效應(yīng)[14-15]圖如圖3所示，而采取了MA或ES策略后訂單向上游傳遞的信息如圖4所示。從前后的對比分析圖可以看出，牛鞭效應(yīng)最高可達到8倍，并且整個波動幅度比較大且頻繁。而采取了需求預(yù)測技術(shù)之后，消除了牛鞭效應(yīng)，需求預(yù)測量在逐漸降低。以參與者之一Retailer為例，進一步展示修改前后對比圖如圖5所示。隨著訂單向上游的傳遞，修改前后的差距在逐漸減小。

圖5 retailer采用預(yù)測技術(shù)前后產(chǎn)生的訂單數(shù)量對比圖

（2）各參與者從訂單接收、配送以及庫存缺貨狀態(tài)三方面的分析對比圖

以retailer為例展示從接收的下游訂單，對下游訂單的配送以及庫存缺貨三個方面做詳細分析，如圖6和圖7所示。圖中x軸表示訂單號，從每個訂單號可以看出該訂單配送的情況以及庫存缺貨的情況。訂單號1至6均未發(fā)生庫存缺貨的情況，并且訂單的配送均是一次性完成，訂單號7雖然也沒有缺貨，但是訂單的配送是分2次完成的。訂單號8至17均不同程度上發(fā)生了缺貨，并且訂單配送是分多次完成的。庫存缺貨和訂單配送次數(shù)是密切相關(guān)的，如果庫存未缺貨但訂單卻多次配送，說明庫存補充是及時的。如果庫存補充及時，訂單配送雖然多次完成，但缺貨狀態(tài)會減少，因為庫存補充是降低缺貨發(fā)生的關(guān)鍵因素。

圖6 retailer訂單接收、配送及庫存缺貨情況分布圖

圖7 retailer訂單接收配送及庫存缺貨及庫存補充情況分布圖

（3）初始庫存量與存儲費用，訂貨費用及缺貨延遲費用的對比分析圖

根據(jù)模型初始假定得知下游顧客提交的訂單數(shù)量在[94，106]之間，向上游提交的訂單數(shù)量在[37，106]之間，庫存補充時間延遲在[0，3]之間。根據(jù)目標實施總費用CT的計算公式：CT=CD+CP+CS=K×s+m×Z+h×t×n+p×l得知總費用以及需求預(yù)測與初始庫存的關(guān)系。

根據(jù)顧客提交的訂單數(shù)量范圍，選取五種初始庫存數(shù)量，50，100，200，300和500。圖8中x軸幾個坐標代表的含義分別為：（1）采取ES的需求預(yù)測方法；（2）采取MA的需求預(yù)測方法；（3）存儲費用；（4）缺貨延遲費用。

圖8 初始庫存分別為50，100，200，300，500時產(chǎn)生的總費用情況

從圖9可以看出如果要使總費用最低：

①存貨費用與初始庫存的關(guān)系并不嚴格成遞增或遞減關(guān)系，而為下游提交訂單數(shù)量的2倍左右為最佳。

②ES和MA兩種需求預(yù)測方法相對來說MA這種方式產(chǎn)生的訂貨費用要少一些。

③存儲費用與庫存量的關(guān)系基本符合正比關(guān)系；在200處為最低。

④缺貨延遲費用最低的為200，如果庫存比200大，則兩者關(guān)系成反比，反之兩者關(guān)系成正比關(guān)系。

圖9 不同初始庫存量訂貨費用、存儲費用及缺貨延遲費用分析對比圖

綜上，該模型直觀顯示了初始庫存量與訂貨費用，存貨費用以及延遲費用之間的關(guān)系，并且通過修改變量很容易變更訂單數(shù)量，前后訂單間隔時間，貨物配送延遲時間等，從而在各種因素間找出使得總費用最低的最優(yōu)庫存管理策略。

（4）不同(s,Q)庫存策略與初始庫存，存儲費用，訂貨費用及缺貨延遲費用的對比分析圖

①初始庫存與總費用及各項費用的關(guān)系

表1 (s,Q)固定情況下初始庫存與費用關(guān)系表元

下游提交的訂單數(shù)量在[10，106]之間，(s,Q)=(80,100)的情況下，訂貨費用最低的是初始庫存為100，存儲費用和缺貨費用最低均出現(xiàn)在初始庫存為0，而總費用最低也出現(xiàn)在初始庫存為0時。

②初始庫存為0時，訂貨點s與費用的關(guān)系

表2 初始庫存為0時訂貨點s與費用的關(guān)系表元

訂貨費用不變，存儲費用隨著訂貨點的增加而增大，缺貨費用與訂貨點成反比關(guān)系，總費用與訂貨點成正比關(guān)系，因此如果要使目標函數(shù)最小，則訂貨點也應(yīng)該越小。

③初始庫存為0時，訂貨點s=20時固定訂貨量Q與費用的關(guān)系

表3 初始庫存為0，訂貨點s=20時固定訂貨量Q與費用的關(guān)系表元

訂貨點固定時，如果固定訂貨量比下游訂單的上界數(shù)量大，那么存儲費用及訂貨費用與Q成正比，而缺貨費用與Q成反比，當Q足夠大時缺貨費用為0，如果Q比較小，缺貨費用的增長幅度是比較大的。因此如果要使目標函數(shù)最小，則固定訂貨量Q不能小于訂貨點s，并且也應(yīng)該越小越好。

5 結(jié)論

本文以顏色Petri網(wǎng)為工具，構(gòu)造了不同需求預(yù)測技術(shù)及庫存控制策略的CPN模型，在初始庫存量，訂貨點及訂貨量等參數(shù)未知情況下，通過實驗仿真數(shù)據(jù)，清晰準確地實現(xiàn)了使目標函數(shù)最小的最優(yōu)庫存管理策略。

[1]Amodeo L，Chen H，El Hadji A.Supply chain inventory optimisation with multiple objectives：an industrial case study[M]//Advances in computational intelligence in transport，logistics，and supply chain management.Berlin Heidelberg：Springer，2008：211-230.

[2]曾強，潘啟東，沈玲.(t,s,S)庫存系統(tǒng)的計算機優(yōu)化仿真[J].工業(yè)工程，2007（3）：123-127.

[3]王鵬飛.基于系統(tǒng)動力學(xué)的產(chǎn)銷兩級供應(yīng)鏈庫存控制仿真研究[D].重慶：重慶大學(xué)，2010.

[4]Yang H C，Chen Y L，Hung M H，et al.Virtual production control system[C]//2010 IEEE Conference on Automation Science and Engineering（CASE）.IEEE，2010：984-989.

[5]Makaji?-Nikoli? D，Pani? B，Vujo?evi? M.Bullwhip effect and supply chain modelling and analysis using CPN tools[C]// Proceedings of the Fifth Workshop and Tutorial on Practical Use of Coloured Petri Nets and the CPN Tools，Aarhus， 2004：219-234.

[6]Papanagnou C I，Halikias G D.Analysing different ordering Policies in a series supply chain by using Coloured Petri Nets[C]//Proceedings of the 20th European Conference on Modelling and Simulation，2006：28-31.

[7]鄭文艷.基于CPN的物流配送系統(tǒng)模型正確性的證明[J].計算機工程與應(yīng)用，2013，49（11）：228-231.

[8]鄭文艷.分層有色Petri網(wǎng)在物流配送系統(tǒng)仿真建模中的應(yīng)用[J].計算機系統(tǒng)應(yīng)用，2013（4）：164-168.

[9]Chen F，Ryan J K，Simchi-Levi D.The impact of exponential smoothing forecasts on the bullwhip effect[J].Naval Research Logistics，2000，47（4）：269-286.

[10]Chen F，Drezner Z，Ryan J K，et al.Quantifying the bullwhip effect in a simple supply chain：The impact of forecasting，lead times，and information[J].Management Science，2000，46（3）：436-443.

[11]Caplin A S.The variability of aggregate demand with（S，s）inventory policies[J].Econometrica：Journal of the Econometric Society，1985：1395-1409.

[12]盧劍峰，聶鐵鑄.隨機庫存系統(tǒng)最優(yōu)訂貨策略計算機模擬研究[J].計算機技術(shù)與發(fā)展，2013（6）：50-53.

[13]Brunner K，Cukierman A，Meltzer A H.Money and economic activity，inventories and business cycles[J].Journal of Monetary Economics，1983，11（3）：281-319.

[14]Papanagnou C，Halikias G D.A state-space approach for analysing the bullwhip effect in supply chains[J].Proceedings of ICTA，2005，5：79-84.

[15]Lee H L，Padmanabhan V，Whang S.Information distortion in a supply chain：the bullwhip effect[J].Management Science，2004，50（12_supplement）：1875-1886.

ZHENG Wenyan

College of Information Management,Dezhou University,Dezhou,Shandong 253023,China

Because of distorted information and poor product management often lead to uncertainty and to instability phenomena.This paper considers the case of the supply chain model using Timed Colored Petri nets and analyses the impact of various continuous inventory policies and known forecasting methods followed by supply chain participants.CPN tools are used for the design of decision-making processes and simulation results are presented to highlight the main issues arising in real systems and to provide insights for future work on modeling and simulation of supply chains.

supply chain;forecasting;inventory control;colored Petri nets

供應(yīng)鏈中由于信息傳遞過程中出現(xiàn)的信息膨脹引起牛鞭效應(yīng)造成各種成本的急劇增長，為使目標函數(shù)費用最小，基于顏色Petri網(wǎng)建立了不同需求預(yù)測方法及庫存策略的CPN模型。在訂單數(shù)量等變量隨機產(chǎn)生以及訂貨點等因素不確定的情況下，通過實驗仿真數(shù)據(jù)確定了存儲、訂貨及缺貨費用與各種不確定變量的關(guān)系，從而確定了最優(yōu)的庫存策略。通過對比實驗證明了該方法的有效性及正確性。

供應(yīng)鏈；需求預(yù)測；庫存策略；顏色Petri網(wǎng)

A

TP391

10.3778/j.issn.1002-8331.1402-0200

ZHENG Wenyan.Performance simulation of uncertainty inventory based on colored petri nets.Computer Engineering and Applications,2014,50（22）：250-255.

德州市社會科學(xué)研究基地項目；2013年山東省自然科學(xué)基金計劃（No.ZR2013GL001）；德州市軟科學(xué)研究計劃項目（No.（2013）德科軟第43號）；2013年校級科技發(fā)展計劃項目（No.311678）。

鄭文艷（1980—），女，講師，主要研究領(lǐng)域為Petri網(wǎng)應(yīng)用。

2014-02-20

2014-04-22

1002-8331（2014）22-0250-06

CNKI網(wǎng)絡(luò)優(yōu)先出版：2014-07-11,http://www.cnki.net/kcms/doi/10.3778/j.issn.1002-8331.1402-0200.html