何斌斌 呂靖

【摘 要】 為實(shí)現(xiàn)對21世紀(jì)海上絲綢之路沿線區(qū)域應(yīng)急儲(chǔ)備庫的全面覆蓋和重點(diǎn)區(qū)域的多重覆蓋的目標(biāo)，運(yùn)用ArcGIS系統(tǒng)對所涉區(qū)域進(jìn)行柵格化及風(fēng)險(xiǎn)評價(jià)，以成本最小為目標(biāo)，構(gòu)建儲(chǔ)備庫選址模型。根據(jù)該模型設(shè)計(jì)遺傳算法，在217個(gè)備選港口中確定了37個(gè)港口作為應(yīng)急儲(chǔ)備庫。

【關(guān)鍵詞】 應(yīng)急儲(chǔ)備庫選址;風(fēng)險(xiǎn)評價(jià);遺傳算法

0 引 言

人類開發(fā)利用海洋時(shí)會(huì)對海洋環(huán)境造成威脅，為了有效應(yīng)對這些問題，國際社會(huì)加強(qiáng)了海事合作，如以海洋環(huán)保和海上安保為主要內(nèi)容的海事保障。與此同時(shí)，隨著21世紀(jì)海上絲綢之路倡議的提出和推進(jìn)，一系列的海上安全問題也隨之出現(xiàn)。為了應(yīng)對這些問題，我國與海上絲綢之路沿線國家不斷加強(qiáng)彼此間的國際海事合作，為沿線區(qū)域提供充分的海事保障。對于我國海事部門來說，在21世紀(jì)海上絲綢之路建設(shè)中可以從規(guī)范重塑與標(biāo)準(zhǔn)修訂、大國協(xié)調(diào)與機(jī)制完善、國際交流與隊(duì)伍建設(shè)等3個(gè)方面入手，尋求建立一個(gè)涵蓋沿線區(qū)域的應(yīng)急儲(chǔ)備庫網(wǎng)絡(luò)。這將有利于促進(jìn)我國周邊外交的有效開展，保障21世紀(jì)海上絲綢之路的順利建設(shè)，有助于推動(dòng)大國間的安全合作。應(yīng)急儲(chǔ)備庫選址問題通常可分為P-中值（P-Median Problem）模型、P-中位（P-Center Problem）模型、集合覆蓋模型和最大覆蓋模型。本文從217個(gè)沿線港口中選擇37個(gè)港口作為應(yīng)急儲(chǔ)備庫備選點(diǎn)進(jìn)行研究，以實(shí)現(xiàn)對所涉區(qū)域的全面覆蓋和重點(diǎn)區(qū)域的多重覆蓋。

1 風(fēng)險(xiǎn)分析

本文主要研究的問題是海上絲綢之路沿線海域應(yīng)急儲(chǔ)備庫備選點(diǎn)的全面覆蓋和重點(diǎn)海域的多重覆蓋，因而需要先對海上絲綢之路沿線海域的主干線進(jìn)行風(fēng)險(xiǎn)分析。

本文采用水域風(fēng)險(xiǎn)計(jì)算公式R=f（P，L），根據(jù)要求首先對影響風(fēng)險(xiǎn)大小的風(fēng)險(xiǎn)發(fā)生概率P、風(fēng)險(xiǎn)后果嚴(yán)重程度L兩個(gè)因素進(jìn)行分析[1]，最后綜合計(jì)算出水域風(fēng)險(xiǎn)值大小R。

運(yùn)用ArcGIS柵格計(jì)算器自帶的高級編程模式，在確定各影響因素的權(quán)重系數(shù)（見表1）后，得到風(fēng)險(xiǎn)發(fā)生概率P及其分布圖。本文參考相關(guān)文獻(xiàn)，采用專家打分法來確定各層指標(biāo)和權(quán)重系數(shù)。

由于越接近海上通道關(guān)鍵節(jié)點(diǎn)（如馬六甲海峽）的地方遭遇風(fēng)險(xiǎn)的后果越嚴(yán)重，因此，利用緩沖區(qū)柵格生成工具對關(guān)鍵節(jié)點(diǎn)相關(guān)數(shù)據(jù)進(jìn)行分析處理，在得到處理結(jié)果后，同樣通過ArcGIS高級編程模式對柵格文件的結(jié)果進(jìn)一步計(jì)算，并導(dǎo)出風(fēng)險(xiǎn)后果嚴(yán)重程度L的分布圖。

在計(jì)算影響風(fēng)險(xiǎn)發(fā)生的P和L兩個(gè)要素后，首先對所得到結(jié)果進(jìn)行歸一化處理，然后對圖層?xùn)鸥窕?，由此確定水域風(fēng)險(xiǎn)值的分布圖及重點(diǎn)水域。

在得到沿線水域風(fēng)險(xiǎn)分布圖后，為了能夠清楚地看到不同風(fēng)險(xiǎn)等級的分布情況，需要對分布圖進(jìn)行重新分類。具體而言，同一風(fēng)險(xiǎn)等級使用相同的顏色，圖片清晰度和精準(zhǔn)度可以通過調(diào)節(jié)單元像素大小來完成?？紤]到模型設(shè)計(jì)是點(diǎn)對點(diǎn)的覆蓋，本文首先將風(fēng)險(xiǎn)分布圖網(wǎng)格化，完成風(fēng)險(xiǎn)劃分;然后用重心法對每一個(gè)網(wǎng)格取點(diǎn)，該風(fēng)險(xiǎn)點(diǎn)的屬性定為面積大小和風(fēng)險(xiǎn)等級，并將其錄入到矢量點(diǎn)屬性表中;最后運(yùn)用計(jì)算面積與風(fēng)險(xiǎn)等級的相乘得到各個(gè)點(diǎn)的風(fēng)險(xiǎn)程度屬性值。

2 模型的建立

建立水上應(yīng)急儲(chǔ)備庫選址模型，以最小的成本實(shí)現(xiàn)對水域的全面覆蓋和重點(diǎn)水域的多重覆蓋。模型的目標(biāo)函數(shù)[1]為

式中： i為水域單元的編號， i=1，2，…，I;? j為應(yīng)急儲(chǔ)備庫候選點(diǎn)的編號， j=1，2，…，J; rj為應(yīng)急儲(chǔ)備庫j的覆蓋半徑; dij為水域單元i到應(yīng)急儲(chǔ)備庫j的距離，當(dāng)dij≤rj時(shí)，應(yīng)急儲(chǔ)備庫j可以覆蓋水域單元i; Ni={ j|dij≤rij|}，為能夠覆蓋水域單元i的集合; y為0或1變量，如果候選點(diǎn)j被選中，則yj=1，否則yj=0; wi1為水域單元i的重要度，重點(diǎn)水域單元的判斷通過重要度的大小來界定; W為重要度的一個(gè)閾值，若wi1≥w，那么就需要覆蓋水域單元i兩次以上; Cj為建設(shè)應(yīng)急儲(chǔ)備庫j的固定成本。

目標(biāo)函數(shù)（1）為應(yīng)急儲(chǔ)備庫建設(shè)的最小成本;約束條件（2）為所有水域單元被完全覆蓋;約束條件（3）為至少2次覆蓋重要度超過閾值的水域單元，即實(shí)現(xiàn)多重覆蓋重點(diǎn)水域的目標(biāo);約束條件（4）使主要參數(shù)變量為0或1。

3 求解算法的設(shè)計(jì)

3.1 染色體編碼

因?yàn)楸疚乃芯康倪x址問題是離散的，且候選點(diǎn)是已知的，因而使用二進(jìn)制編碼的方法。本文共有J個(gè)候選點(diǎn)，染色體的長度為J。具體而言，每一個(gè)基因位1，2，…，J分別對應(yīng)一個(gè)候選點(diǎn)，若選中一個(gè)候選點(diǎn)并建立應(yīng)急儲(chǔ)備庫，那么該基因位的取值為1，否則為0。

3.2 初始可行解

隨機(jī)產(chǎn)生遺傳算法的初始可行解，種群數(shù)量設(shè)為M，染色體隨機(jī)產(chǎn)生并判斷其是否可行，如果可行就將其保留，否則會(huì)出現(xiàn)一條新的染色體。當(dāng)染色體數(shù)量達(dá)到M時(shí)停止產(chǎn)生。

3.3 適應(yīng)度函數(shù)

基于本文的選址方案，計(jì)算應(yīng)急儲(chǔ)備庫建設(shè)成本，得到目標(biāo)函數(shù)值Z。在遺傳算法中，適應(yīng)度函數(shù)的選取關(guān)系到收斂的速度和能否確定最優(yōu)解，本文設(shè)計(jì)的適應(yīng)度函數(shù)為

式中： fn為染色體n的適應(yīng)度值; Zn為染色體n相應(yīng)的目標(biāo)值; I為一個(gè)相對大的數(shù)，同時(shí)適應(yīng)度函數(shù)值越大越好。

3.4 遺傳算子選擇

在遺傳算法中，下一代的個(gè)體的選擇采用輪盤賭模式，通過單點(diǎn)交叉、單點(diǎn)變異的方式迭代運(yùn)算。

3.5 精英保留策略

精英保留策略是指提取已經(jīng)得到適應(yīng)度最高的個(gè)體，不參與隨后的交叉和變異運(yùn)算，并且用適應(yīng)度最高的個(gè)體來代替本代群體中經(jīng)過遺傳操作后適應(yīng)度最低的個(gè)體。為了保護(hù)交叉、變異過程中的上一代最優(yōu)個(gè)體，在算法中將每一代染色體的適應(yīng)度值排序，并確定出精英個(gè)體，予以單獨(dú)保留。在完成一次遺傳操作后，精英個(gè)體代替適應(yīng)度值最低的進(jìn)入下一代種群。

3.6 懲罰函數(shù)

懲罰函數(shù)的設(shè)定是為了將有約束最優(yōu)化問題變?yōu)闊o約束最優(yōu)化問題。在種群每一次迭代過后會(huì)出現(xiàn)不可行的個(gè)體，那么就需要運(yùn)用懲罰函數(shù)來降低不可行個(gè)體進(jìn)入下一代的概率。通常，懲罰函數(shù)法有內(nèi)部懲罰函數(shù)和外部懲罰函數(shù)。本文運(yùn)用內(nèi)部懲罰函數(shù) fn'=fnHaP，即減去一個(gè)常數(shù)P，降低該個(gè)體的適應(yīng)度值，使其逐漸被淘汰。

3.7 終止條件

設(shè)置最大迭代次數(shù)，當(dāng)遺傳算法達(dá)到設(shè)定次數(shù)并且種群的適應(yīng)度值趨于穩(wěn)定后停止運(yùn)行。

4 實(shí)例分析

本文選取21世紀(jì)海上絲綢之路沿線港口共計(jì)217個(gè)，涉及國家21個(gè)。采用Matlab軟件進(jìn)行遺傳算法的相關(guān)編程，具體為：種群大小M=50，染色體長度J=217，交叉概率設(shè)定為0.9，變色概率設(shè)定為0.1，每次交叉保留2個(gè)個(gè)體。同時(shí)，在適應(yīng)度函數(shù)中I=300，懲罰函數(shù)P=200，迭代次數(shù)為100次。程序運(yùn)行共50次，41次得到相同結(jié)果，平均運(yùn)行時(shí)間為2.33 s。在眾多迭代中選取一次迭代計(jì)算（見圖1）。最終，通過建立模型和遺傳算法求解，在217個(gè)備選港口中確定了37個(gè)港口作為21世紀(jì)海上絲綢之路急儲(chǔ)備庫（見表2）。

5 結(jié) 語本文將選址模型運(yùn)用到21世紀(jì)海上絲綢之路中，將沿線國家的港口作為應(yīng)急儲(chǔ)備庫備選點(diǎn)，實(shí)現(xiàn)對沿線區(qū)域的全面覆蓋和重點(diǎn)區(qū)域的多重覆蓋，并通過模型算法得到了最優(yōu)解。在解決實(shí)際問題時(shí)，本文運(yùn)用ArcGIS對相關(guān)水域進(jìn)行風(fēng)險(xiǎn)劃分并柵格化，確定了重點(diǎn)區(qū)域，由此順利地與模型相結(jié)合。本文解決這一問題的思路，可為解決類似問題提供參考。

參考文獻(xiàn)：

[1] 艾云飛，呂靖，張麗麗，等. 水上應(yīng)急物資儲(chǔ)備庫選址――分配優(yōu)化模型[J]. 大連海事大學(xué)學(xué)報(bào)，2015（2）：62-66.