何 健， 張大鵬， 文正中

(1. 軍事交通學(xué)院研究生管理大隊(duì)， 天津 300161； 2. 軍事交通學(xué)院軍用車(chē)輛系， 天津 300161)

基于系統(tǒng)動(dòng)力學(xué)的車(chē)輛器材庫(kù)存需求量預(yù)測(cè)

何 健1， 張大鵬2， 文正中1

(1. 軍事交通學(xué)院研究生管理大隊(duì)， 天津 300161； 2. 軍事交通學(xué)院軍用車(chē)輛系， 天津 300161)

針對(duì)車(chē)輛器材庫(kù)存需求量預(yù)測(cè)存在的影響因素考慮不全和預(yù)測(cè)精度不高的問(wèn)題，分析了影響庫(kù)存需求量的主要因素及各變量之間的因果關(guān)系，建立了車(chē)輛器材庫(kù)存需求量預(yù)測(cè)系統(tǒng)動(dòng)力學(xué)模型，并進(jìn)行了實(shí)例驗(yàn)證。結(jié)果表明：該模型考慮因素齊全，預(yù)測(cè)結(jié)果更接近實(shí)際，具有一定的參考和應(yīng)用價(jià)值，為科學(xué)預(yù)測(cè)車(chē)輛器材庫(kù)存需求量提供了新方法。

車(chē)輛器材；系統(tǒng)動(dòng)力學(xué)模型；庫(kù)存需求量預(yù)測(cè)

車(chē)輛器材是車(chē)輛裝備維修保障的物質(zhì)基礎(chǔ)，準(zhǔn)確預(yù)測(cè)車(chē)輛器材的庫(kù)存需求量對(duì)提高車(chē)輛裝備完好率具有重要意義[1]。然而，目前對(duì)車(chē)輛器材庫(kù)存需求量的研究主要停留在理論計(jì)算方法上，對(duì)裝備維修保障中庫(kù)存需求規(guī)律和發(fā)展趨勢(shì)研究相對(duì)較少。徐飛等[2]在分析裝備使用情況和需求數(shù)據(jù)波動(dòng)的基礎(chǔ)上，運(yùn)用粗糙集理論和層次分析法修正了灰色預(yù)測(cè)模型，該模型所需樣本少，易挖掘出數(shù)據(jù)中隱藏的規(guī)律，但對(duì)需求波動(dòng)性大、離散性強(qiáng)的器材庫(kù)存需求量難以進(jìn)行有效預(yù)測(cè)；羅曉萌等[3]基于時(shí)間序列法對(duì)歷史數(shù)據(jù)進(jìn)行了分析，同時(shí)考慮價(jià)格對(duì)需求的影響，預(yù)測(cè)了在一定服務(wù)水平下的庫(kù)存需求量，該方法雖然得出了動(dòng)態(tài)庫(kù)存采購(gòu)策略，但考慮的庫(kù)存需求影響因素較少，預(yù)測(cè)偏差較大；樓增龍等[4]以實(shí)際統(tǒng)計(jì)數(shù)據(jù)為樣本，建立了多元線(xiàn)性回歸預(yù)測(cè)模型，該模型對(duì)于間斷性需求或當(dāng)需求數(shù)據(jù)量較小時(shí)，預(yù)測(cè)偏差較大。

系統(tǒng)動(dòng)力學(xué)預(yù)測(cè)方法能夠綜合考慮各種影響因素，對(duì)于歷史數(shù)據(jù)不足、復(fù)雜非線(xiàn)性問(wèn)題均具有良好的預(yù)測(cè)效果，且可預(yù)測(cè)系統(tǒng)的發(fā)展趨勢(shì)和波動(dòng)周期。因此，筆者在系統(tǒng)全面分析庫(kù)存需求影響因素的基礎(chǔ)上，建立了車(chē)輛器材庫(kù)存需求量預(yù)測(cè)的系統(tǒng)動(dòng)力學(xué)模型，并通過(guò)實(shí)例驗(yàn)證模型的有效性，為科學(xué)籌措器材、管理庫(kù)存提供支持。

1 車(chē)輛器材庫(kù)存需求量影響因素分析

我軍器材倉(cāng)庫(kù)正常周轉(zhuǎn)的關(guān)鍵在于器材入庫(kù)——保管——出庫(kù)過(guò)程能夠維持一定程度的動(dòng)態(tài)平衡，為維持該動(dòng)態(tài)平衡，必須準(zhǔn)確預(yù)測(cè)庫(kù)存需求量。庫(kù)存需求量影響因素較多，主要有器材可靠性、器材運(yùn)行環(huán)境、器材使用強(qiáng)度和維修保障水平等。筆者通過(guò)部隊(duì)實(shí)地調(diào)研和專(zhuān)家評(píng)定，選取車(chē)輛裝備規(guī)模、器材新增量和器材報(bào)廢量3個(gè)參數(shù)來(lái)衡量上述因素對(duì)庫(kù)存需求量的影響程度[5]。

1.1 車(chē)輛裝備規(guī)模

車(chē)輛裝備規(guī)模是指為完成某次任務(wù)所投入使用的車(chē)輛裝備數(shù)，是影響庫(kù)存需求量的基本因素。車(chē)輛裝備規(guī)模越大，器材需求量會(huì)越大，相應(yīng)的庫(kù)存需求量也越大。此外，車(chē)輛裝備規(guī)模還影響故障器材數(shù)量和器材新需求量，進(jìn)而影響器材報(bào)廢量；尤其是當(dāng)受到戰(zhàn)爭(zhēng)爆發(fā)和上級(jí)授予緊急任務(wù)等其他無(wú)法預(yù)料的隨機(jī)性因素影響時(shí)，車(chē)輛裝備規(guī)模會(huì)突然增大，庫(kù)存需求量會(huì)隨之迅速上升。

1.2 器材新增量

器材新增量是指隨著軍事訓(xùn)練任務(wù)或車(chē)輛裝備數(shù)量增多而需要增加的車(chē)輛器材數(shù)量，與車(chē)輛裝備規(guī)模增加量和庫(kù)存調(diào)整量有關(guān)，是影響庫(kù)存需求量的重要因素。隨著車(chē)輛裝備規(guī)模的增加，所需儲(chǔ)備的車(chē)輛器材就會(huì)增多，庫(kù)存需求量也隨之增大。近年來(lái)，我軍擔(dān)負(fù)的軍事訓(xùn)練、演習(xí)任務(wù)越來(lái)越頻繁，器材新增量不斷提升，但整體增幅逐漸減小并趨于穩(wěn)定。

1.3 器材報(bào)廢量

器材報(bào)廢量是指由于質(zhì)量不符合技術(shù)標(biāo)準(zhǔn)而發(fā)生故障且失去修復(fù)意義或價(jià)值較低直接進(jìn)行更換的車(chē)輛器材數(shù)量，其與故障器材數(shù)量和報(bào)廢速率有關(guān)，是影響庫(kù)存需求量的關(guān)鍵因素。隨著車(chē)輛裝備服役年限的增加，器材使用壽命將逐漸下降，器材報(bào)廢速率逐漸加快，同時(shí)器材報(bào)廢量會(huì)隨之增大，庫(kù)存需求量也隨之增大。

2 車(chē)輛器材庫(kù)存需求量預(yù)測(cè)系統(tǒng)動(dòng)力學(xué)模型構(gòu)建

2.1 預(yù)測(cè)流程

運(yùn)用系統(tǒng)動(dòng)力學(xué)解決實(shí)際問(wèn)題的流程為：分析系統(tǒng)信息——理清系統(tǒng)變量間的關(guān)系——解決系統(tǒng)問(wèn)題。這是一個(gè)由淺入深、波浪前進(jìn)的過(guò)程，它從多個(gè)角度剖析問(wèn)題，綜合考慮各種影響因素，進(jìn)而使預(yù)測(cè)結(jié)果更加科學(xué)、可靠。圖1為車(chē)輛器材庫(kù)存需求量系統(tǒng)動(dòng)力學(xué)模型預(yù)測(cè)流程。

圖1 車(chē)輛器材庫(kù)存需求量系統(tǒng)動(dòng)力學(xué)模型預(yù)測(cè)流程

2.2 變量之間因果關(guān)系分析

根據(jù)車(chē)輛裝備實(shí)際運(yùn)行狀況和庫(kù)存需求量影響因素分析結(jié)果，建立車(chē)輛器材庫(kù)存需求量預(yù)測(cè)變量之間的因果關(guān)系，如圖2所示。圖中：“+”、“-”分別表示兩變量間的關(guān)系為正相關(guān)、負(fù)相關(guān)。

圖2 車(chē)輛器材庫(kù)存需求量預(yù)測(cè)變量之間的因果關(guān)系

1) 器材報(bào)廢速率是指單位時(shí)間內(nèi)器材報(bào)廢的數(shù)量，其與維修水平和使用年限有關(guān)。隨著維修技術(shù)水平的不斷提高，可維修翻新的器材數(shù)量隨之增大，報(bào)廢速率隨之減小。

2) 器材新需求量是指在一定保障度的要求下，當(dāng)車(chē)輛裝備規(guī)模增大后所需增加的器材數(shù)量。新需求量主要受車(chē)輛裝備規(guī)模影響，其隨車(chē)輛裝備規(guī)模的增大而增大。

3) 故障器材數(shù)量與車(chē)輛裝備規(guī)模和器材故障率有關(guān)。車(chē)輛裝備規(guī)模越大，需要維修或淘汰報(bào)廢的器材基數(shù)就越大，故障器材數(shù)量也隨之增大。

4) 器材庫(kù)存需求量主要受新增量和器材報(bào)廢量的影響。隨著器材新增量和報(bào)廢量的增大，器材庫(kù)存需求量也隨之增大。

主要反饋回路為：新增量→(+器材報(bào)廢量→)+庫(kù)存需求量→-庫(kù)存調(diào)整→+新增量，即當(dāng)器材新增量增大時(shí)，器材報(bào)廢量增大，庫(kù)存需求量就增大，為了使器材庫(kù)存需求量達(dá)到理想狀態(tài)但又避免積壓浪費(fèi)，需要調(diào)整庫(kù)存，減少器材采購(gòu)數(shù)量，抑制新增量增大[6]。

2.3 系統(tǒng)流程圖

系統(tǒng)流程圖可全面地反映系統(tǒng)的物質(zhì)流、信息流和反饋?zhàn)饔谩８鶕?jù)車(chē)輛器材庫(kù)存需求量預(yù)測(cè)變量之間的因果關(guān)系，運(yùn)用Vensim軟件繪制車(chē)輛器材庫(kù)存需求量預(yù)測(cè)系統(tǒng)流程圖，如圖3所示。

圖3 車(chē)輛器材庫(kù)存需求量預(yù)測(cè)的系統(tǒng)流程圖

系統(tǒng)動(dòng)力學(xué)模型中各符號(hào)及其含義如表1所示。

表1 模型中各符號(hào)及其含義

系統(tǒng)動(dòng)力學(xué)模型中的主要方程如下：

N車(chē)輛裝備規(guī)模=INTEG(車(chē)輛裝備增加數(shù)量，X)。

T 車(chē)輛裝備增加數(shù)量=WITH LOOKUP(time,數(shù)組)。

T 維修水平影響因子=WITH LOOKUP(維修水平，數(shù)組)。

L 器材新增量=新需求量+器材調(diào)整量。

L 新需求量=車(chē)輛裝備增加數(shù)量×單車(chē)用量。

L 理想庫(kù)存=SMOOTH(需求量，移動(dòng)平滑時(shí)間)。

L 器材報(bào)廢量=故障器材數(shù)量×報(bào)廢速率。

L 報(bào)廢速率t=自然報(bào)廢速率t+人為報(bào)廢速率t。

其中：報(bào)廢速率包括自然報(bào)廢速率和人為報(bào)廢速率。報(bào)廢速率t為器材使用t年后的報(bào)廢速率。由于器材發(fā)生報(bào)廢后要求立即進(jìn)行更換和補(bǔ)充，則在第t年，整個(gè)車(chē)輛裝備系統(tǒng)中將包含使用1～t年的器材。一般運(yùn)用排隊(duì)論來(lái)求解第t年的報(bào)廢速率。

L 庫(kù)存需求量=器材新增量+器材報(bào)廢量。

A 庫(kù)存偏差=(理想庫(kù)存-器材需求量)×(1-保障率)。

A 庫(kù)存調(diào)整量=庫(kù)存偏差。

R 自然報(bào)廢速率t=β(t-1)-β(t)。

自然報(bào)廢指隨著使用年限的增加，器材壽命達(dá)到極限而產(chǎn)生的報(bào)廢，主要與使用時(shí)間有關(guān)。β(t)為器材使用t年后的自然完好率，一般用雙參數(shù)Weibull分布函數(shù)來(lái)計(jì)算[7]，即β(t)=exp((-t/T)k)。其中：T為某種器材的平均使用壽命；k為特征參數(shù)。

R 人為報(bào)廢速率t=σ×器材故障率。

人為報(bào)廢是指故障器材中因無(wú)法進(jìn)行人為修復(fù)而導(dǎo)致的報(bào)廢，主要受維修水平和器材故障率影響。當(dāng)某種器材的關(guān)鍵性、價(jià)值較低，對(duì)故障件直接進(jìn)行更換時(shí)，則σ=1,人為報(bào)廢速率=器材故障率；當(dāng)某種器材有修復(fù)價(jià)值時(shí)，則σ=1-維修水平影響因子，此時(shí)，人為報(bào)廢率與維修水平有關(guān)。

T器材故障率=故障率表(time)×使用強(qiáng)度影響因子×運(yùn)行環(huán)境影響因子×操作水平影響因子。

3 實(shí)例分析

3.1 參數(shù)及變量確定

以2005年列裝部隊(duì)的運(yùn)輸車(chē)解放CA1121的XX器材為研究對(duì)象，其單車(chē)用量為2件，且具有修復(fù)價(jià)值。以其10年的需求及故障數(shù)據(jù)為預(yù)測(cè)樣本，保障率為0.99。由于部隊(duì)保障任務(wù)的不間斷性，當(dāng)器材消耗后必須立即進(jìn)行補(bǔ)充，因此，不考慮時(shí)間延遲。部隊(duì)進(jìn)行器材采購(gòu)的周期通常為1年，在采購(gòu)周期內(nèi)，XX器材的理想庫(kù)存量即為年期望需求量。通過(guò)移動(dòng)平滑法，應(yīng)用前5年的庫(kù)存需求量來(lái)預(yù)測(cè)該器材下一年的理想庫(kù)存量[8]。解放CA1121的裝備規(guī)模及相關(guān)數(shù)據(jù)如表2所示，其中車(chē)輛裝備工作時(shí)間用來(lái)度量其使用強(qiáng)度影響因子的大小。車(chē)輛裝備數(shù)量增長(zhǎng)率的表函數(shù)如圖4所示，其中：橫坐標(biāo)為年份；縱坐標(biāo)為車(chē)輛裝備增加數(shù)量。

表2 2005-2015年某器材倉(cāng)庫(kù)解放CA1121的裝備規(guī)模及相關(guān)數(shù)據(jù)

圖4 車(chē)輛裝備數(shù)量增長(zhǎng)率表函數(shù)

器材的報(bào)廢速率不僅受車(chē)輛使用年限的影響，而且受維修水平和器材故障率的制約。根據(jù)張鵬等[9]的研究結(jié)果及該器材倉(cāng)庫(kù)的實(shí)際情況，該器材的維修水平影響因子的表函數(shù)如圖5所示，其中：橫坐標(biāo)為維修水平的模糊評(píng)價(jià)值；縱坐標(biāo)為維修水平影響因子。根據(jù)陶金宸[10]的研究結(jié)果，應(yīng)用模糊綜合評(píng)判法得到該器材倉(cāng)庫(kù)的維修水平為良好，評(píng)價(jià)值為0.75，由圖5可得到相應(yīng)的維修水平影響因子為0.7。器材自身故障率表函數(shù)如圖6所示，其中：橫坐標(biāo)為年份；縱坐標(biāo)為自身故障率。

圖5 XX器材維修水平影響因子表函數(shù)

圖6 XX器材自身故障率表函數(shù)

器材故障率不僅與器材自身的可靠性有關(guān)，還與使用強(qiáng)度影響因子、運(yùn)行環(huán)境影響因子和操作水平影響因子有關(guān)。因此，修正后的XX器材故障率表函數(shù)如圖7所示，其中：橫坐標(biāo)為年份；縱坐標(biāo)為器材故障率。

圖7 XX器材故障率表函數(shù)

根據(jù)《中國(guó)人民解放軍車(chē)輛保養(yǎng)規(guī)定》，車(chē)輛裝備每6年進(jìn)行一次大修，通常情況下，累積使用6年的器材也應(yīng)進(jìn)行更換[11]。因此，計(jì)算自然報(bào)廢速率時(shí)，設(shè)定該器材的平均使用壽命T=6年，特征參數(shù)k=3。運(yùn)用排隊(duì)論方法和MATLAB軟件計(jì)算第t年該器材的報(bào)廢速率，則XX器材的自然完好率和相關(guān)報(bào)廢速率如表3所示。其中：2016-2019年的相關(guān)數(shù)據(jù)通過(guò)雙參數(shù)Weibull分布函數(shù)和排隊(duì)論方法計(jì)算得出。

表3 XX器材的自然完好率及相關(guān)報(bào)廢速率

3.2 模擬仿真

通過(guò)該器材倉(cāng)庫(kù)的實(shí)際需求數(shù)據(jù)和構(gòu)建的數(shù)學(xué)模型方法，計(jì)算出各參數(shù)及變量的數(shù)值，應(yīng)用仿真軟件Vensim進(jìn)行系統(tǒng)模型的仿真預(yù)測(cè)，具體庫(kù)存需求量仿真結(jié)果如圖8所示。

圖8 XX器材庫(kù)存需求量的模擬仿真結(jié)果

由圖8可以看出：1)該器材庫(kù)存需求量每隔6年會(huì)發(fā)生突增，這是由于筆者設(shè)定器材的使用壽命T=6年，器材累積工作6年時(shí)必須進(jìn)行更換。2)該器材庫(kù)存需求量基本呈上升趨勢(shì)，這是由于軍事訓(xùn)練、演習(xí)等任務(wù)不斷增多，車(chē)輛裝備的使用強(qiáng)度和發(fā)生故障的概率隨之增加，庫(kù)存需求量也隨之增大。

車(chē)輛器材庫(kù)存需求量預(yù)測(cè)系統(tǒng)動(dòng)力學(xué)模型的仿真預(yù)測(cè)值與實(shí)際需求量對(duì)比如圖9所示。系統(tǒng)動(dòng)力學(xué)模型仿真預(yù)測(cè)結(jié)果及相對(duì)誤差如表4所示。

圖9 庫(kù)存需求量系統(tǒng)動(dòng)力學(xué)模型仿真預(yù)測(cè)值與實(shí)際需求量對(duì)比

表4 系統(tǒng)動(dòng)力學(xué)模型仿真預(yù)測(cè)結(jié)果及相對(duì)誤差

年份實(shí)際需求量/件仿真預(yù)測(cè)值/件相對(duì)誤差/%備注200524228．33200622209．09max200727263．70200832320．00200934365．88201081843．70201129276．89201230286．66201334332．94201437355．40201538372．63平均誤差/%5．522<10%

由表4可以看出：仿真預(yù)測(cè)值與實(shí)際需求量的最大誤差為9.09%(<10%)，表明所建立的庫(kù)存需求量預(yù)測(cè)系統(tǒng)動(dòng)力學(xué)模型的預(yù)測(cè)精度達(dá)到要求，可用于部隊(duì)車(chē)輛器材的庫(kù)存需求量預(yù)測(cè)。運(yùn)用該模型對(duì)2016-2019年的庫(kù)存需求量進(jìn)行仿真預(yù)測(cè)，結(jié)果分別為57、32、38、42件。

4 結(jié)論

筆者構(gòu)建的模型不僅考慮了外部使用環(huán)境、器材自身故障率等因素，而且通過(guò)表函數(shù)計(jì)算了車(chē)輛裝備數(shù)量增長(zhǎng)率等不易確定的變量，有效地彌補(bǔ)了傳統(tǒng)預(yù)測(cè)方法考慮影響因素較少、預(yù)測(cè)精度偏低的不足。

結(jié)合某器材倉(cāng)庫(kù)的實(shí)際需求數(shù)據(jù)進(jìn)行了仿真預(yù)測(cè)，并將得到的仿真預(yù)測(cè)值和實(shí)際需求量進(jìn)行了對(duì)比，結(jié)果表明：相對(duì)誤差滿(mǎn)足精度要求，驗(yàn)證了預(yù)測(cè)模型的有效性，為車(chē)輛器材庫(kù)存需求量預(yù)測(cè)提供了新方法。最后，對(duì)該倉(cāng)庫(kù)2016-2019年XX器材的庫(kù)存需求量進(jìn)行了預(yù)測(cè)，為部隊(duì)車(chē)輛器材籌措及庫(kù)存管理提供了支撐。下一步將進(jìn)一步研究各參數(shù)對(duì)所建模型的影響程度(即靈敏度)。

[1] 王亮． 軍用車(chē)輛器材保障概論[M]． 北京： 金盾出版社， 2014： 9-13．

[2] 徐飛， 楊士英， 魏祥生． 基于改進(jìn)灰色預(yù)測(cè)模型的雷達(dá)備件消耗預(yù)測(cè)[J]． 雷達(dá)科學(xué)與技術(shù)， 2015， 12(6)： 609-614．

[3] 羅曉萌， 李建斌， 胡鵬． 基于時(shí)間序列預(yù)測(cè)的電子商務(wù)庫(kù)存優(yōu)化策略[J]． 系統(tǒng)工程， 2014， 6(6)： 91-98．

[4] 樓增龍， 陳海松， 賀偉雄，等． 基于多元線(xiàn)性回歸的維修器材消耗預(yù)測(cè)[J]． 裝備制造技術(shù)， 2014， 28(6)： 94-96．

[5] 富盛． 基于系統(tǒng)動(dòng)力學(xué)的乘用車(chē)保有量和報(bào)廢量的預(yù)測(cè)模型研究[D]． 武漢： 武漢理工大學(xué)，2014．

[6] 王晚香．鐵路客運(yùn)量的系統(tǒng)動(dòng)力學(xué)預(yù)測(cè)[J]． 大連交通大學(xué)學(xué)報(bào)， 2014， 35(5)： 10-12，16．

[7] 唐旭， 張寶生， 鄧紅梅，等． 基于系統(tǒng)動(dòng)力學(xué)的中國(guó)石油產(chǎn)量預(yù)測(cè)分析[J]． 系統(tǒng)工程理論與實(shí)踐， 2010， 30(2)： 207-212．

[8] 許長(zhǎng)新， 嚴(yán)以新， 張萍． 基于系統(tǒng)動(dòng)力學(xué)的港口吞吐量預(yù)測(cè)模型[J]． 水運(yùn)工程， 2006， 5(5)： 26-28，40．

[9] 張鵬， 葛婷． 基于系統(tǒng)動(dòng)力學(xué)的航材庫(kù)存管理策略研究[J]． 測(cè)控技術(shù)， 2013， 32(11)： 135-138．

[10] 陶金宸． 部隊(duì)車(chē)輛周轉(zhuǎn)器材庫(kù)存結(jié)構(gòu)調(diào)整研究[D]． 天津： 軍事交通學(xué)院， 2014．

[11] 王云鵬， 楊志發(fā)， 李世武． 基于系統(tǒng)動(dòng)力學(xué)的道路運(yùn)輸量預(yù)測(cè)模型[J]． 吉林大學(xué)學(xué)報(bào)， 2005， 35 (4)： 426-430．

(責(zé)任編輯: 王生鳳)

Vehicle Material Inventory Demand Forecast Based on System Dynamics

HE Jian1, ZHANG Da-peng2, WEN Zheng-zhong1

(1.Brigade of Postgraduate Management, Academy of Military Transportation, Tianjin 300161, China; 2. Department of Military Vehicles, Academy of Military Transportation, Tianjin 300161, China)

In view of the problems of incomplete consinderation of influence factors and low accuracy of the vehicle material inventory demand forecast, the main factors affecting the inventory demand and the causal relationship between the variables are analyzed, and a system dynamics model of the vehicle material inventory demand forecast is set up. Examples are verified and the results show that this model considers factors completely and results are closer to the actual, which has reference and application value, and it proposes a new method of scientifically forecasting the inventory demand.

vehicle material; system dynamics model; inventory demand forecast

1672-1497(2016)05-0025-06

2016-08-24

何 健(1991-)，男，碩士研究生。

E92

A

10.3969/j.issn.1672-1497.2016.05.005