李睿雪，馬 良，劉 勇

（上海理工大學(xué) 管理學(xué)院，上海 200093）

隨著工業(yè)的不斷發(fā)展，工廠選址問題作為工業(yè)生產(chǎn)開展的第一環(huán)節(jié)，始終是工業(yè)中研究的重點(diǎn)。在以往的容量限制選址模型中通常將需求量作為確定值考慮，但是，在如今工業(yè)發(fā)展和市場環(huán)境不斷變化的情況下，需求點(diǎn)的需求不再是確定值。因此，考慮需求的不確定性對選址模型的影響顯得尤為重要。

近年來，國內(nèi)外學(xué)者在選址問題中針對需求不確定以及容量限制作了大量的研究。在不確定需求方面，Bayram 等[1]在應(yīng)急設(shè)施選址模型中考慮疏散需求的不確定性，建立了基于需求場景的兩階段隨機(jī)規(guī)劃疏散模型。張鑫等[2]通過區(qū)間灰數(shù)理論預(yù)測不確定需求，將差分進(jìn)化和遺傳算法混合后對應(yīng)急物資選址模型進(jìn)行求解，但是，求解算法依賴大量的數(shù)據(jù)。成律等[3]在傳統(tǒng)配送中心選址中，利用隨機(jī)規(guī)劃理論和模糊理論構(gòu)建了隨機(jī)客戶需求情況下配送中心選址的數(shù)學(xué)模型，并應(yīng)用遺傳算法通過算例對模型進(jìn)行驗證，但遺傳算法存在早熟等缺點(diǎn)。隨后國內(nèi)外學(xué)者對容量限制因素進(jìn)行研究，Jena等[4]發(fā)現(xiàn)容量限制對動態(tài)設(shè)施選址有極大影響。李建光等[5]通過離散場景表征避難需求的不確定性，考慮情形權(quán)重以及避難位置的規(guī)模，構(gòu)建了一個魯棒隨機(jī)規(guī)劃模型來確定避難場所選址。商麗媛等[6]在多情景下需求不確定的應(yīng)急物流配送中心選址模型中，將配送中心容量作為約束，對不確定需求用區(qū)間灰數(shù)表示，設(shè)計了免疫量子粒子群算法進(jìn)行求解。徐小峰等[7]構(gòu)建了考慮同時取送貨需求和帶模糊容量約束的在線設(shè)施動態(tài)選址模型。

然而，目前研究大多是考慮需求不確定或容量限制單一因素對工廠選址的影響，而實際中這兩種因素共同影響選址決策?；诖耍疚慕⒉淮_定性需求下分級生產(chǎn)規(guī)模容量限制工廠選址模型。將備選工廠設(shè)置為按生產(chǎn)規(guī)模分級分段以滿足需求的不斷變化，并通過區(qū)間聯(lián)系數(shù)對不確定的需求進(jìn)行預(yù)測。在算法上，之前的研究算法通常采用遺傳算法和粒子群算法等求解，存在早熟和收斂性差等缺點(diǎn)。由于本文提出的工廠選址模型考慮需求的不確定性和容量限制規(guī)模的變化，求解算法編碼和計算量相對之前研究更為復(fù)雜。而具有較好的適應(yīng)性和并行性的蜂群算法更適用于本文模型。因此，采用人工蜂群算法求解以保證能夠達(dá)到全局最優(yōu)解。

1 生產(chǎn)規(guī)模分階段容量限制工廠模型

1.1 工廠選址模型及參數(shù)

容量限制選址作為選址問題的擴(kuò)展模型，其含義為在備選點(diǎn)滿足需求點(diǎn)需求并且不超過其本身的容量限制前提下，選擇總距離成本和建設(shè)成本最低的位置作為選址方案。針對需求不確定性這一特征，將生產(chǎn)規(guī)模設(shè)置為分階段規(guī)劃。根據(jù)備選位置不同需求量之和選擇合適的規(guī)模，更有利于節(jié)約成本以及后續(xù)發(fā)展。為此建立了改進(jìn)的容量限制選址問題模型。

模型的選址目標(biāo)是在備選位置中選擇能使距離成本最小和建設(shè)成本最小的最優(yōu)位置。設(shè)I為備選工廠的位置集；J為需求點(diǎn)的位置集；K為時間集合；N為備選工廠不同生產(chǎn)階段規(guī)模的集合；ωj表示需求點(diǎn)j需求量歷年變化率的平均值；Rjk表示需求點(diǎn)j第k年的需求量變化率；dj為需求點(diǎn)j的需求量；tij表示從備選工廠位置i到需求點(diǎn)位置j之間的運(yùn)輸成本；ein表示備選工廠位置i實際選擇的第n階段生產(chǎn)規(guī)模；gi表示備選工廠位置i的最大生產(chǎn)規(guī)模；hi表示備選工廠位置i在該地點(diǎn)的基礎(chǔ)建設(shè)成本；kin表示備選工廠位置i在第n階段生產(chǎn)能力下每單位的生產(chǎn)建設(shè)成本；f為選擇方案的總成本。

在工廠的生產(chǎn)規(guī)模大于需求點(diǎn)的需求之和的前提下，選擇最優(yōu)的備選位置，模型如下：

模型中目標(biāo)函數(shù)式（1）表示工廠選址的運(yùn)輸成本和固定生產(chǎn)成本之和最??；式（2）表示需求點(diǎn)的需求量歷年變化率的平均值；約束式（3）表示備選工廠所服務(wù)的需求點(diǎn)的需求之和不超出該備選工廠的最大生產(chǎn)能力；式（4）表示確保需求點(diǎn)和開放的備選位置保持對應(yīng)；式（5）和式（6）為決策變量約束，xij為工廠i對需求點(diǎn)j服務(wù)，yi為在備選工廠位置i選擇建設(shè)工廠，其中，開放取值為1，關(guān)閉取值為0。

1.2 基于聯(lián)系數(shù)的區(qū)間灰數(shù)的不確定需求預(yù)測

針對事物不確定性和確定性關(guān)系，本文提出聯(lián)系數(shù)這一概念，表示事物確定性測度和不確定性測度以及兩者之間關(guān)系。不確定需求的區(qū)間數(shù)通過聯(lián)系數(shù)表示，則可以更大程度地保留原來需求的不確定性[8]。

假設(shè)u=[A,B]是一個正區(qū)間數(shù)，可將其轉(zhuǎn)化為聯(lián)系數(shù)，表示為u=A+(B-A)i，根據(jù)聯(lián)系數(shù)定義，i∈[-1, 1]。轉(zhuǎn)化后區(qū)間數(shù)的范圍擴(kuò)大了一倍，變?yōu)閇A-(B+A),A+(B-A)]。但由于需求的實際情況，本文不需要負(fù)范圍，故取i∈[0, 1]。

灰色GM(1, 1)模型是常用的對含有不確定因素進(jìn)行預(yù)測的模型，其中，灰色系統(tǒng)的基本元素為灰數(shù)，通常用“Gr”表示，指在某一個區(qū)間或者某個數(shù)集內(nèi)取值的不確定數(shù)。區(qū)間灰數(shù)則是既有下界c又有上界d的灰數(shù)。記作Gr∈[c,d]，c≤d，上界與下界的差值為區(qū)間灰數(shù)的長度，記作l=d-c。在預(yù)測模型中需要將一些列區(qū)間灰數(shù)組成一個區(qū)間灰數(shù)序列，記作Gr{Gr1,Gr2, ···,Grn}[9]。

設(shè)區(qū)間灰數(shù)Gr∈[c,d]，則區(qū)間灰數(shù)序列為L{Gr}={Gr1,Gr2, ···,Grn}，可將區(qū)間灰數(shù)白化，其白化值為Gri=c+μd，μ=[0, 1]，μ為0～1的隨機(jī)數(shù)。同樣可以對區(qū)間灰數(shù)序列白化，其中，C{Gr}為白部，D{Gr}為區(qū)間序列白化后的灰部，如式（7）所示。

分別對區(qū)間灰色序列白部和灰部構(gòu)建預(yù)測模型，將累加式還原到原始序列，得到式（8）和式（9）。

式中，β為最小二乘法的參數(shù)列中的參數(shù)。

結(jié)合聯(lián)系數(shù)和區(qū)間數(shù)之間的關(guān)系，將不確定需求的區(qū)間轉(zhuǎn)化為聯(lián)系數(shù)，再利用灰色預(yù)測模型將聯(lián)系數(shù)作為預(yù)測灰數(shù)，對需求進(jìn)行預(yù)測。基于聯(lián)系數(shù)的區(qū)間灰數(shù)預(yù)測模型既適用于小樣本的數(shù)據(jù)，又能挖掘出樣本之間隱含的不確定性。

設(shè)區(qū)間灰數(shù)Gr∈[a,b]，U=A+Bi，A=a,B=b-a,i∈[0, 1]，則U=A+Bi為其對應(yīng)的聯(lián)系數(shù)。U={U1,U2, ···,Un}，Uk=Ak+Bki，Ak=ak，Bk=bk-ak，其中，A和B分別為聯(lián)系數(shù)U的同部和異部，A={A1,A2, ···,An}和B={B1,B2, ···,Bn}為聯(lián)系數(shù)序列U的同步序列和異部序列。根據(jù)區(qū)間灰數(shù)預(yù)測模型分別對A和B構(gòu)建預(yù)測模型，可得同部和異部的響應(yīng)函數(shù)。

通過A和B值還原成區(qū)間灰數(shù)預(yù)測，即可得到需求量的模擬值。

2 改進(jìn)的人工蜂群算法

2.1 人工蜂群算法

人工蜂群算法（artificial bee colony algorithm，ABC）是Karaboga等[10]在蜜蜂自組織采蜜的群體智能行為上提出的一種非數(shù)值優(yōu)化計算方法，其中，蜂群實現(xiàn)采蜜的智能行為主要包含3個部分：蜜源、采蜜蜂、待工蜂，蜂群的采蜜過程主要包括：搜索蜜源、招募蜜蜂和放棄蜜源。

蜜源：蜜源表示存在的可能的解。通常為了方便求解，將收益度作為衡量蜜源優(yōu)劣的標(biāo)準(zhǔn)；采蜜蜂：表示在蜜源中采蜜的蜜蜂，采蜜蜂通常攜帶蜜源信息，并且以一定比例成為引領(lǐng)蜂，引領(lǐng)蜂攜帶蜜源信息通過搖擺舞與其他蜜蜂分享這些信息；待工蜂：表示等待尋找蜜源采集的蜜蜂，他們可以分為兩種，偵察蜂和追隨蜂，偵查蜂搜索蜂巢附件的新蜜源，追隨蜂在蜂巢內(nèi)等待,通過分享采蜜蜂的信息, 尋找蜜源。

2.2 人工蜂群算法的步驟

在ABC算法中，將蜜源的位置抽象為解空間中的點(diǎn)，即待求解問題的可行解。蜜源i的質(zhì)量用收益度fit表示，np為蜜源的數(shù)量。ABC算法中引領(lǐng)蜂和待工蜂各占蜂群的一半，且同一時刻下同一蜜源中只有一只引領(lǐng)蜂采蜜[11]。

a. 設(shè)求解問題的維度為D，蜜源i的位置表示為X，其中，t為當(dāng)前的迭代次數(shù)。蜜源的初始位置根據(jù)式（12）在搜索空間中隨機(jī)產(chǎn)生。

式中：xij代表第i個蜜源xi的第j維度值，i取值{1, 2, ···,sn}，j取值{1, 2, ···,d}；xmax,ij和xmin,ij分別為維度j的最大值和最小值。

初始化蜜源就是對每個蜜源的所在維度通過式（12）選取一個隨機(jī)值。

b. 引領(lǐng)蜂階段，在引領(lǐng)蜂階段，引領(lǐng)蜂用式（13）更新新蜜源。

式中：xkj代表領(lǐng)域蜜源，k取值{1,2, ···,sn}，且k不等于i；φ為取值在[-1, 1]之間的隨機(jī)數(shù)。

通過式（13）更新得到新蜜源，再借助貪婪算法的原則，比較更新前、后的蜜源適應(yīng)度值，選擇最優(yōu)者。

c. 追隨蜂階段，引領(lǐng)蜂搜索結(jié)束后，追隨蜂開始搜索，在這一階段，在舞蹈區(qū)接受引領(lǐng)蜂分享的蜜源信息，追隨蜂采用輪盤賭策略來選擇其所追隨的引領(lǐng)蜂。為了保證適應(yīng)度更高的個體可以得到保留。同時對每個蜜源設(shè)置參數(shù)R，當(dāng)蜜源更新后被保留，R為0；反之，則R加1。從而記錄一個蜜源沒有被更新的次數(shù)。

d. 偵察蜂階段。如果一個蜜源經(jīng)過多次開采后，R值過高，并且超過了設(shè)定的閾值T，則需要放棄這個蜜源，同時進(jìn)入偵察蜂階段。這一階段，偵察蜂利用式（12）隨機(jī)尋找新的蜜源，以代替被拋棄的蜜源。

e. 評估所有蜜蜂搜索得到食物源的適應(yīng)度并進(jìn)行排序，進(jìn)行下一次迭代。

2.3 改進(jìn)人工蜂群算法

基本的人工蜂群算法是為解決連續(xù)空間的問題而設(shè)計，解是由蜜蜂在蜜源附近搜索生成的，通常都是在一個連續(xù)的范圍空間。但是，在處理調(diào)度安排、選址方案等離散型問題時，蜜蜂的鄰域搜索從連續(xù)變?yōu)殡x散，這就使得算法原來的公式不再有效。為了使這一算法可以應(yīng)用在離散問題上，本文基于人工蜂群算法的基本原理以及相關(guān)研究[12]，依照本文模型對公式進(jìn)行改進(jìn)，以適應(yīng)離散化的選址模型，提出了改進(jìn)的離散人工蜂群算法（new discrete artificial bee colony algorithm,NDABC），使得人工蜂群算法在保留原有并行性的優(yōu)點(diǎn)下，對離散問題也有良好的適應(yīng)性。同時也對其他算法的研究提供思路。

a. 個體編碼。建立蜜蜂個體的編碼形式是選址方案表示的關(guān)鍵，本文構(gòu)建的工廠選址問題主要包括兩方面：每個備選點(diǎn)是否被選中和每個需求點(diǎn)和備選點(diǎn)的對應(yīng)關(guān)系。本文將備選點(diǎn)按位置進(jìn)行排序，設(shè)定標(biāo)號（1, 2, ···,m），需求點(diǎn)同樣是按位置排序后形成一列n個元素的數(shù)列。

例如，有8個需求點(diǎn)，4個備選點(diǎn)，在4個備選點(diǎn)中選取2個位置作為工廠位置，解X為：

其中，X表示在4個備選位置中選取3和1號位置，3號位置為需求點(diǎn)1，3，4，6服務(wù)，1號位置為需求點(diǎn)2，5，7，8服務(wù)。

b. 初始化種群。作為整個算法計算的起點(diǎn)，優(yōu)秀的起點(diǎn)可以更好地將蜜蜂分散在可行解中。為了更好地讓蜜蜂分布在可行范圍內(nèi)，初始解采用隨機(jī)生成的辦法。備選點(diǎn)的選址數(shù)量和位置采用隨機(jī)生成方式，以擴(kuò)大選址方案。在確定備選位置后，隨機(jī)在備選位置上分配需求點(diǎn)，并且將選址位置的方案作為蜜源，保證在初始生成的種群中個體分布在不同蜜源上。

為了保證蜜蜂群體的隨機(jī)搜索能力，在蜜蜂群每次迭代更新中將群體的一小部分待工蜂作為偵察蜂。通過數(shù)據(jù)實驗發(fā)現(xiàn)，根據(jù)適應(yīng)度值進(jìn)行排序后的45%個體作為追隨蜂；選取總體種群中5%的個體作為偵察蜂。這種比例分配可以達(dá)到較好的優(yōu)化效果，保證進(jìn)行隨機(jī)搜索，防止陷入局部最優(yōu)。

c. 蜜蜂更新規(guī)則。

針對每一個個體都在自身的蜜源內(nèi)進(jìn)行搜索更新，本文參考了文獻(xiàn)[13]的更新策略，提出了一種多方式的領(lǐng)域更新策略。其中，包括3種領(lǐng)域更新策略：領(lǐng)域交換更新策略，記作U1，在待更新個體中隨機(jī)選取2個位置，并交換其值；領(lǐng)域插入更新策略，記作U2，在該更新個體中隨機(jī)選取2個位置，插入1個備選位置的序號，后面的元素依次后延；領(lǐng)域內(nèi)倒轉(zhuǎn)更新策略，記作U3，待更新個體中隨機(jī)選取2個位置，將2個位置間的位置進(jìn)行倒置更新。表1列出各策略的具體示例。

表1 蜜蜂領(lǐng)域更新策略Tab.1 Strategy of bee field renewal

d. 觀察蜂更新策略。

在人工蜂群算法中追隨蜂是在所追隨的采蜜蜂的蜜源附近進(jìn)行搜索采蜜，所以，使用式（13）確定觀察蜂的蜜源Y=[ya,yb, ···]，其中，a,b表示備選位置的編號，Sf(i)表示在整個選址方案中選取備選位置i所產(chǎn)生的成本，同時保留其追隨的采蜜蜂的選址策略。

e. 追隨更新。在傳統(tǒng)蜂群算法中，蜜蜂是在不同蜜源間進(jìn)行搜索，針對選址問題，首先將工廠位置按照地理位置進(jìn)行排序，將選址方案中的備選點(diǎn)位置作為蜜源。因此，追隨蜂則會在所追隨的選址方案中隨機(jī)地在選取位置的前后進(jìn)行變化，以模擬追隨蜂在蜜源附近搜索的行為。當(dāng)追隨蜂表示的備選點(diǎn)確定后，則開始進(jìn)行備選點(diǎn)與需求點(diǎn)間的對應(yīng)。其中，更新原則按照上文的更新策略進(jìn)行更新。

f. 偵察蜂更新策略?？傮w種群中選取總量5%的個體作為偵察蜂，進(jìn)行隨機(jī)搜索，以廣大搜索范圍，防止陷入局部最優(yōu)。

g. 局部搜索。在傳統(tǒng)人工蜂群算法中，為了防止單個蜜蜂多次在同一蜜源進(jìn)行搜索，設(shè)置了最大搜索次數(shù)，在本算法中為了進(jìn)一步擴(kuò)大算法的搜索范圍，防止陷入局部最優(yōu)。每個蜜蜂個體均在個體更新過程中更新R次，保留R次中適應(yīng)度值最大的個體。

h. 選擇策略。采取類錦標(biāo)賽，根據(jù)模型實際情況，與傳統(tǒng)錦標(biāo)賽算法中選擇策略不同，本策略中在采蜜蜂個體中隨機(jī)選擇n個個體，找出這n個個體中的最大值，作為觀察蜂的追隨個體。這樣可以擴(kuò)大隨機(jī)搜索方位，提高算法的搜索能力。

2.4 算法過程

a. 初始化蜜蜂種群，初始時刻所有蜜蜂都是偵察蜂，全局搜索蜜源并計算收益度；

b. 對步驟a中所有蜜蜂根據(jù)收益度的相對大小，蜜蜂轉(zhuǎn)化為采蜜蜂(收益度靠前)和觀察蜂，觀察蜂根據(jù)采蜜蜂的舞蹈來接收招募；

c. 每個采蜜蜂都會在原蜜源附近采蜜(局部搜索)；

d. 觀察蜂按照蜜源以一定概率選擇蜜源，并在附近局部搜索；

e. 采蜜蜂和觀察蜂都在其蜜源附近使用局部搜索策略，保留最優(yōu)的個體，進(jìn)入下一次更新；

f. 記錄當(dāng)前所有蜜蜂找到的最優(yōu)蜜源，執(zhí)行步驟b。直至算法迭代結(jié)束。

3 仿真實驗與分析

為了驗證改進(jìn)算法的有效性，以及模型的可行性，由于真實數(shù)據(jù)難以獲取，本文通過模擬文獻(xiàn)[14]的數(shù)據(jù)進(jìn)行仿真實驗。依照文獻(xiàn)[14]的數(shù)據(jù)特征，通過Matlab在其基礎(chǔ)上模擬生成更大規(guī)模的數(shù)據(jù)，保證算例可以更好地模擬實際環(huán)境并能更好地驗證算法和模型。本文中備選點(diǎn)為20個，需求點(diǎn)為150，每個備選點(diǎn)可分為4階段的生產(chǎn)規(guī)模，每階段的單位生產(chǎn)成本為[16，14，12，10]（元/萬件），數(shù)據(jù)如表2所示。對每個需求點(diǎn)記錄近5年的需求，如表3所示。

表2 備選工廠位置參數(shù)Tab.2 Alternative plant location parameters

表3 需求點(diǎn)位置參數(shù)Tab.3 Location parameters of demand points

為了方便對比不同算法的迭代效果，本文加大了文中算法的迭代次數(shù)，減少群體的規(guī)模數(shù)量。保證得到更明細(xì)的算法優(yōu)化結(jié)果曲線。因此，改進(jìn)蜂群算法參數(shù)設(shè)置：初始種群個體數(shù)為40，迭代500次。同時通過遺傳算法（genetic algorithm, GA）、粒子群算法（particle swarm optimization, PSO）、海鷗算法（seagull optimization algorithm, SOA）和文獻(xiàn)[15]提出的改進(jìn)螢火蟲算法（firefly algorithm, FA）作為對比組，改進(jìn)蜂群算法參數(shù)設(shè)置：初始種群個體數(shù)為40，迭代500次；遺傳算法的初始種群為40，迭代500次，采用單點(diǎn)交叉和單點(diǎn)變異，其中交叉概率為0.5，變異概率為0.05；粒子群算法的初始種群為40，最大迭代次數(shù)500次，慣性因子為0.6，C1=C2=2，C1，C2為學(xué)習(xí)因子；螢火蟲算法的初始種群為40，迭代500次，其中，移動概率0.7，變異概率0.1。海鷗算法的初始種群為40，最大迭代次數(shù)為500次，控制因子Fk的值從2線性降低到0。分別對每個算法進(jìn)行20次獨(dú)立計算后對比每種算法的平均值和標(biāo)準(zhǔn)差。每個需求點(diǎn)近5年的需求量使用式（10）和式（11）進(jìn)行預(yù)測，將其代入算法中進(jìn)行計算。選址工廠數(shù)為2。計算后算法迭代如圖1所示。

圖1 不同算法迭代對比圖Fig.1 Iterative comparison of different algorithms

根據(jù)圖1結(jié)果分析，本文算法相較于其他算法在迭代中不僅種群整體更新，而且加入了防止陷入局部最優(yōu)的自身更新，由此提高了搜索效率，能夠快速穩(wěn)定地得到最優(yōu)解，迭代次數(shù)遠(yuǎn)小于其他算法，收斂速度有了明顯的提高。另外，從表4中也可以看出，相對其他4種算法，本文提出的改進(jìn)的人工蜂群算法，在20次獨(dú)立計算的平均值以及標(biāo)準(zhǔn)差明顯優(yōu)于其他算法。在選址方案中，本文算法并結(jié)合其他算法得出的方案，1號位置及附近和12號位置及附近都是選址方案中較優(yōu)的位置。本文提出的算法得到的選址方案剛好選擇1號和12號位置，并且在幾種算法對比中得到最優(yōu)的總成本。根據(jù)備選工廠的數(shù)據(jù)可以看出，12號位置的生產(chǎn)規(guī)模比1號位置更大。本文的方案12號位置相對1號位置服務(wù)更多的需求點(diǎn)，可以更好地滿足需求點(diǎn)的需求。最后結(jié)合實際分析選址位置可以看出，本文算法所選的工廠位置更加分散，有利于實際情況下工廠的發(fā)展。

表4 不同算法計算結(jié)果對比Tab.4 Comparison of calculation results of different algorithms

在驗證了改進(jìn)算法有效后，對模型本身進(jìn)行深入分析。開始考慮改變?nèi)萘肯拗乒S選址模型的選址數(shù)量對整個選址模型的影響，對于本文的工廠選址問題模型選址位置不易過多，因此，僅分別求解選址數(shù)量N=2，3，4這3種情況，每種情況獨(dú)立運(yùn)行20次，計算20次的平均值、標(biāo)準(zhǔn)差以及迭代情況，運(yùn)行結(jié)果如表5和圖2所示。

圖2 不同工廠選址數(shù)量的迭代曲線圖Fig. 2 Iterative curve of different numbers of plant sites

表5 不同選址數(shù)量下選址模型結(jié)果對比Tab.5 Comparison of location model results under different location numbers

通過分析不同選址數(shù)量的選址方案結(jié)果得出，當(dāng)N=2時，即僅選擇2個備選位置作為工廠位置，得到的成本最小，所選擇的工廠的位置趨近于中心位置，其生產(chǎn)規(guī)模更大，服務(wù)需求點(diǎn)的數(shù)量更多，每個備選位置的后續(xù)擴(kuò)建能力較差；而當(dāng)N=3時，選擇3個備選位置產(chǎn)生的成本相對2個備選位置有所增加，工廠有選擇更小的生產(chǎn)規(guī)模的情況，可以兼顧后續(xù)需求增加后的擴(kuò)建需求；最后N=4時，整個模型的選址策略所需的總成本最大，工廠的位置相比較其他策略更加分散，備選工廠僅需要選擇最小的生產(chǎn)規(guī)模就可以滿足需求。因此，如果選址方案的目標(biāo)是追求低成本，那么可以選擇更少的選址數(shù)量，以達(dá)到更低的選址方案成本。但是，如果為了后續(xù)的發(fā)展以及更大程度地滿足需求的不斷增長，可以在選址方案中適當(dāng)增加選址數(shù)量，以滿足后續(xù)發(fā)展中需求的增加以及服務(wù)尚未發(fā)掘的潛在需求點(diǎn)。

最后在選址位置上分析，3種不同選址數(shù)量的方案中均選擇到了2號位置附近、9號位置和12號位置附近。由此分析可得，這3個位置是后續(xù)選址決策中可以重點(diǎn)考慮分析的位置，為后續(xù)的選址方案提供參考。

4 結(jié)束語

分析了現(xiàn)有容量限制選址問題中缺少對需求不確定性考慮和容量限制約束多為固定值這兩部分。同時將這兩部分因素加入容量限制選址模型中，用聯(lián)系數(shù)和區(qū)間灰數(shù)將不確定的需求轉(zhuǎn)化為固定值以及將固定的容量限制約束設(shè)置為分階段，建立了全新的容量限制選址模型。并且采用了針對模型的特征而改進(jìn)的人工蜂群算法對模型進(jìn)行求解。通過與其他算法對比，驗證了算法在此類問題上更加高效。但本文的研究僅為單一產(chǎn)品，下一步可以拓展為多種產(chǎn)品下的選址問題。