劉 清 董詩(shī)瑀 王 磊
(武漢理工大學(xué)交通學(xué)院1) 武漢 430063) (國(guó)家水運(yùn)安全工程技術(shù)研究中心2) 武漢 430063)
內(nèi)河水上交通安全事故大多發(fā)生在航道條件復(fù)雜的水域,從海事保障角度辨識(shí)航道“黑點(diǎn)”,既是滿足水域通航安全保障的需要,也是航道改善輔助決策的必要基礎(chǔ),更是從本質(zhì)上保障船舶通航安全的必然途徑[1-2].
目前道路交通“黑點(diǎn)”的有關(guān)研究已較為成熟,已有關(guān)于航道“黑點(diǎn)”的研究主要采用了模糊綜合評(píng)價(jià)及動(dòng)態(tài)聚類法[3]、基于網(wǎng)格理論的模糊綜合評(píng)價(jià)[4]、層次聚類算法[5]等方法,這些方法大多基于問(wèn)卷的數(shù)據(jù)處理,主觀性強(qiáng),對(duì)內(nèi)河航道“黑點(diǎn)”的定量研究比較缺乏,這一缺乏主要源自兩方面的難題:①如何實(shí)現(xiàn)復(fù)雜航道條件下船舶通航風(fēng)險(xiǎn)特征變量的定量描述;②如何用恰當(dāng)?shù)哪P图八惴▉?lái)驗(yàn)證航道“黑點(diǎn)”辨識(shí)的效果,并證明這種方法的合理性.這些難題也就是本文力圖解決的研究命題[6-7].
文中明確界定了內(nèi)河航道黑點(diǎn)的內(nèi)涵,建立了內(nèi)河航道“黑點(diǎn)”辨識(shí)的MEA-BP神經(jīng)網(wǎng)絡(luò)模型,運(yùn)用思維進(jìn)化算法優(yōu)化模型并迅速找出“黑點(diǎn)”位置,通過(guò)對(duì)“黑點(diǎn)”影響因子的敏感度分析,提取影響內(nèi)河航道“黑點(diǎn)”的關(guān)鍵影響因素.
根據(jù)內(nèi)河水上交通事故特點(diǎn),從環(huán)境因素出發(fā),將內(nèi)河“黑點(diǎn)”航道界定為:一個(gè)比較長(zhǎng)的時(shí)間范圍里或者規(guī)定的統(tǒng)計(jì)周期內(nèi),某些位置(內(nèi)河地點(diǎn)、航段、區(qū)域)上發(fā)生的水上交通事故數(shù)量或者事故特征明顯比其他航段異常,或者該位置航道條件威脅航行安全,即內(nèi)河航道上影響船舶通航安全因素集中的航段[8].
采用基于思維進(jìn)化算法優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)建模方法.BP神經(jīng)網(wǎng)絡(luò)揭示數(shù)據(jù)樣本中蘊(yùn)含的非線性關(guān)系,靈活方便地對(duì)多成因的復(fù)雜未知系數(shù)進(jìn)行建模[9].BP神經(jīng)網(wǎng)絡(luò)已經(jīng)在道路黑點(diǎn)辨識(shí)研究中廣泛應(yīng)用,并取得良好的效果.運(yùn)用BP神經(jīng)網(wǎng)絡(luò)進(jìn)行內(nèi)河“黑點(diǎn)”航道辨識(shí),只要收集輸入層和輸出層的數(shù)據(jù),不必考慮輸入層多個(gè)影響因素之間的權(quán)重,模型會(huì)根據(jù)數(shù)據(jù)特征,自動(dòng)調(diào)整權(quán)值與閾值,當(dāng)目標(biāo)誤差和檢驗(yàn)誤差率達(dá)到預(yù)測(cè)標(biāo)準(zhǔn)時(shí)網(wǎng)絡(luò)停止運(yùn)算,模型建立[10-11].
針對(duì)BP神經(jīng)網(wǎng)絡(luò)容易陷入局部極小,收斂速度慢,選取思維進(jìn)化算法(mind evolutionary algorithm,MEA)優(yōu)化BP神經(jīng)網(wǎng)絡(luò)的初始權(quán)值和閾值[12].MEA把群體劃分為優(yōu)勝子群體和臨時(shí)子群體,趨同和異化操作分別進(jìn)行局部搜索和全局搜索,這兩種功能相互協(xié)調(diào)且保持一定的獨(dú)立性,便于分別提高效率,任一方面的改進(jìn)都對(duì)提高算法的整體搜索效率有利[13].
影響內(nèi)河船舶通航安全的航道條件因素很多,在選擇評(píng)價(jià)指標(biāo)時(shí),要考慮其量化的可能性及相互之間的邏輯關(guān)系,否則,會(huì)失去辨識(shí)的客觀性與真實(shí)性.充分考慮內(nèi)河航道影響因素構(gòu)建的評(píng)價(jià)指標(biāo)體系見(jiàn)圖1.
圖1 內(nèi)河航道“黑點(diǎn)”辨識(shí)評(píng)價(jià)指標(biāo)體系
采取輸入層八個(gè)元素,隱含層七個(gè)元素,輸出層一個(gè)元素的BP神經(jīng)網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu).建立的模型結(jié)構(gòu)見(jiàn)圖2.
圖2 內(nèi)河航道黑點(diǎn)辨識(shí)模型結(jié)構(gòu)示意圖
采用思維進(jìn)化優(yōu)化BP網(wǎng)絡(luò)的基本步驟如下.
1) 產(chǎn)生初始結(jié)構(gòu)群體 設(shè)神經(jīng)網(wǎng)絡(luò)的最大隱含層數(shù)為n,每一個(gè)隱含層最大節(jié)點(diǎn)數(shù)為m.隨機(jī)產(chǎn)生N組數(shù)作為初始結(jié)構(gòu)群體,每組數(shù)中包含n個(gè)元素,其中每個(gè)元素均從隱含層節(jié)點(diǎn)數(shù)0,1,…,m中選取.
2) 評(píng)價(jià)函數(shù)的確定 選擇訓(xùn)練集的均方誤差的倒數(shù)作為各個(gè)個(gè)體與種群的得分函數(shù).
(1)
式中:yi為第i個(gè)訓(xùn)練樣本的網(wǎng)絡(luò)輸出值;ti為目標(biāo)輸出;p為訓(xùn)練樣本數(shù).
3) 訓(xùn)練權(quán)值和閾值 對(duì)于每個(gè)網(wǎng)絡(luò)結(jié)構(gòu),在(0,1)之間以均勻分布產(chǎn)生R組隨機(jī)數(shù),作為初始的權(quán)值和閾值群體.根據(jù)網(wǎng)絡(luò)計(jì)算規(guī)則,按照得分函數(shù)計(jì)算每個(gè)個(gè)體得分,誤差越小,則認(rèn)為個(gè)體所含權(quán)值和閾值信息越好,得分就越高.得分最高的q個(gè)個(gè)體被稱為優(yōu)勝者.
4) 子群體趨同過(guò)程 趨同操作是在每個(gè)子群體內(nèi)部進(jìn)行的.分別以每一個(gè)優(yōu)勝者為中心,服從正態(tài)分布產(chǎn)生個(gè)體,形成M個(gè)優(yōu)勝子群體和T個(gè)臨時(shí)子群體,每個(gè)子群體包含w個(gè)個(gè)體.正態(tài)分布的中心就是優(yōu)勝者的坐標(biāo),即優(yōu)勝者的權(quán)值.
5) 子群體異化過(guò)程 異化操作是全局空間內(nèi)各子群體為成為優(yōu)勝者的競(jìng)爭(zhēng)過(guò)程.全局公告板記錄了各子群體的得分及成熟度,若一個(gè)臨時(shí)子群體的得分高于某個(gè)成熟的一般子群體的得分,則該成熟的一般子群體被獲勝的臨時(shí)子群體替代,原成熟的一般子群體中的個(gè)體被釋放,形成新的臨時(shí)子群體.
6) 解析最優(yōu)個(gè)體 當(dāng)滿足迭代停止條件時(shí),思維進(jìn)化算法結(jié)束優(yōu)化過(guò)程.此時(shí),根據(jù)編碼規(guī)則,對(duì)尋找到的最優(yōu)個(gè)體進(jìn)行解析,從而得到對(duì)應(yīng)的BP神經(jīng)網(wǎng)絡(luò)的權(quán)值和閾值.
7) 訓(xùn)練MEA優(yōu)化的BP神經(jīng)網(wǎng)絡(luò) 用得到的最優(yōu)權(quán)值與閾值建立BP神經(jīng)網(wǎng)絡(luò).
敏感度分析是用于研究輸出結(jié)果的不確定性以及不同變量之間關(guān)系的模型.MEA-BP神經(jīng)網(wǎng)絡(luò)模型的預(yù)測(cè)過(guò)程屬于“黑箱”過(guò)程,無(wú)法得出各影響因素對(duì)“黑點(diǎn)”辨識(shí)結(jié)果的影響程度,而了解影響程度最大的影響因素對(duì)提升整治航道“黑點(diǎn)”的效率與效果具有重要的作用.利用MEA-BP模型,改變“黑點(diǎn)”航段各影響因素的輸入值后分析輸出結(jié)果,得到對(duì)輸出結(jié)果的變化影響最大的輸入變量,即對(duì)“黑點(diǎn)”航道影響因素進(jìn)行敏感度分析,提取“黑點(diǎn)”航道敏感因子,作為航道條件改善的依據(jù).
三峽庫(kù)區(qū)通航環(huán)境較為復(fù)雜,與其他航段相比,由航道條件因素導(dǎo)致的事故數(shù)量較多.同時(shí)受水庫(kù)蓄水泄洪影響,庫(kù)區(qū)通航環(huán)境并不穩(wěn)定,尤其在庫(kù)尾回水變動(dòng)區(qū),航道隨水位變化在天然河流航道與庫(kù)區(qū)航道間轉(zhuǎn)換,對(duì)船舶通航安全造成極大影響.故選取三峽庫(kù)區(qū)作為典型航段進(jìn)行模型驗(yàn)證.
將三峽庫(kù)區(qū)長(zhǎng)江航道里程上游50 km至上游700 km,設(shè)置5 km長(zhǎng)度為一個(gè)單元,可分為130個(gè)單元.收集自2011—2018年的三峽庫(kù)區(qū)水上交通事故數(shù)據(jù)共80組作為樣本數(shù)據(jù).樣本數(shù)據(jù)見(jiàn)表1.
表1 樣本數(shù)據(jù)
在樣本的訓(xùn)練過(guò)程中,隨機(jī)選取70組樣本數(shù)據(jù)為訓(xùn)練樣本,其余10組為測(cè)試樣本.取種群規(guī)模M=20,優(yōu)勝子種群個(gè)數(shù)a=5,臨時(shí)子種群個(gè)數(shù)b=5,初始權(quán)值與閾值取值范圍為(-1,1),最大迭代次數(shù)為100次.BP神經(jīng)網(wǎng)絡(luò)的最大訓(xùn)練次數(shù)為100次,訓(xùn)練目標(biāo)誤差為0.01,學(xué)習(xí)率為0.1.
利用MEA算法得到的最優(yōu)權(quán)值與閾值建立BP模型,訓(xùn)練完畢后,將三峽庫(kù)區(qū)所有單元的數(shù)據(jù)輸入網(wǎng)絡(luò),通過(guò)正向計(jì)算輸出各單元輸出值,輸出值統(tǒng)計(jì)圖見(jiàn)圖3.
圖3 三峽庫(kù)區(qū)實(shí)例驗(yàn)證模型輸出值統(tǒng)計(jì)圖
將輸出值大于1的單元記為內(nèi)河“黑點(diǎn)”航道,共計(jì)10處,見(jiàn)表2.
表2 基于MEA-BP神經(jīng)網(wǎng)絡(luò)模型的三峽庫(kù)區(qū)水上交通“黑點(diǎn)”辨識(shí)結(jié)果
對(duì)應(yīng)的長(zhǎng)江航道里程為上游60—65,65—70,95—100,105—110,125—130,155—160,545—550,670—675,765—680,690—695 km.
經(jīng)過(guò)對(duì)航段實(shí)際情況的調(diào)研以及與海事部門(mén)統(tǒng)計(jì)的事故集中多發(fā)位置比較,運(yùn)用水上交通“黑點(diǎn)”辨識(shí)模型所辨識(shí)到的“黑點(diǎn)”位置與實(shí)際相符.以上游670—680 km航段為例,長(zhǎng)江航道里程670—680 km屬于三峽庫(kù)區(qū)庫(kù)尾回水變動(dòng)區(qū),該航段航道彎曲程度大,航道寬度突然變窄,存在多處礁石淺灘,水流湍急,架設(shè)李家沱長(zhǎng)江大橋,兩岸有多處碼頭,屬于重點(diǎn)監(jiān)控航段.由此表明MEA-BP神經(jīng)網(wǎng)絡(luò)模型辨識(shí)水上事故“黑點(diǎn)”是可行的且具有一定的準(zhǔn)確性.
通過(guò)改變辨識(shí)出的“黑點(diǎn)”航道單元相應(yīng)指標(biāo)數(shù)值,再次正向輸入訓(xùn)練完畢的“黑點(diǎn)”辨識(shí)模型,分析不同輸入變量的改變對(duì)輸出結(jié)果的影響程度,排序得到與輸出結(jié)果的變化聯(lián)系最為密切的輸入變量,提取庫(kù)區(qū)“黑點(diǎn)”航道風(fēng)險(xiǎn)敏感因子.
改變后的輸入值與對(duì)應(yīng)的輸出值見(jiàn)表3.
表3 用于敏感因子提取的模型輸入值與輸出值
通過(guò)分析影響因子數(shù)值修改前后輸出值,可以得到各“黑點(diǎn)”航段影響因素對(duì)輸出值影響程度排序,從而提取出敏感因子,結(jié)果見(jiàn)表4.
表4 各“黑點(diǎn)”航段風(fēng)險(xiǎn)敏感因子
文中界定了航道“黑點(diǎn)”的內(nèi)涵,針對(duì)影響內(nèi)河船舶通航安全的諸多因子,從航道復(fù)雜性和事故多發(fā)兩方面提出了航道“黑點(diǎn)”辨識(shí)指標(biāo)體系,建立了內(nèi)河航道黑點(diǎn)辨識(shí)的MEA-BP神經(jīng)網(wǎng)絡(luò)模型.以三峽庫(kù)區(qū)2011—2018年的事故為例,對(duì)模型進(jìn)行驗(yàn)證,輸出結(jié)果表明所構(gòu)建的模型與算法與實(shí)際相符.在此基礎(chǔ)上,通過(guò)對(duì)辨識(shí)模型影響因子的敏感性分析,提取了“黑點(diǎn)”航道重點(diǎn)治理因素.研究結(jié)果可用于內(nèi)河航道整治的輔助決策,對(duì)保障內(nèi)河船舶通航安全具有應(yīng)用價(jià)值.
模型主要針對(duì)現(xiàn)有航道條件進(jìn)行分析,模型動(dòng)態(tài)性不足,在監(jiān)控航道條件動(dòng)態(tài)變化方面仍有研究的空間,如何獲取更多更精確數(shù)據(jù)也對(duì)模型穩(wěn)定性產(chǎn)生很大影響,在后續(xù)研究中,可考慮在靜態(tài)辨識(shí)基礎(chǔ)上考慮加入動(dòng)態(tài)因素,建立動(dòng)態(tài)辨識(shí)模型.