滿苗苗, 柳 林, 2, 程 鵬, 魏國忠, 張 省

?

基于元胞自動機的海洋溢油模型

滿苗苗1, 柳 林1, 2, 程 鵬1, 魏國忠3, 張 省3

(1. 山東科技大學(xué) 測繪科學(xué)與工程學(xué)院, 山東 青島 266590; 2. 海島(礁)測繪技術(shù)國家測繪地理信息局重點實驗室, 山東 青島 266590; 3. 山東省國土測繪院, 山東 濟南 250102)

本文采用三維元胞自動機模型, 對海洋溢油過程進行了模擬, 對風(fēng)流、水流等水平方向的影響系數(shù)進行了改進?？紤]了蒸發(fā)、垂直擴散、岸邊附著、溶解、乳化等因子的影響, 并且引入乳化含水率將乳化因子進行實際量化, 構(gòu)建了海洋溢油模型。并以“DeepSpill”的溢油實驗為基礎(chǔ)進行了模擬, 實驗結(jié)果Kappa系數(shù)達到0.902, 與實際相比具有較好的一致性。本文對海洋溢油事故進行預(yù)測模擬, 為提出合理有效的應(yīng)急預(yù)案提供科學(xué)依據(jù)。

元胞自動機; 海洋溢油; 模擬

2001年西班牙“威望”號油輪斷裂、2010年墨西哥灣石油鉆井平臺溢油事故[1-2], 對海洋生物、生態(tài)環(huán)境造成了巨大損失。利用元胞自動機對海洋溢油擴散進行模擬, 能夠更直觀、真實地顯示溢油的動態(tài)變化。近年來, 元胞自動機在海洋溢油方面得到越來越多的關(guān)注。Karafyllidis[3]綜合了風(fēng)、流、蒸發(fā)等因素, 采用元胞自動機研究了海洋溢油的運動和遷移[3]。王璐等[4]考慮了流場和風(fēng)場等因素, 采用二維元胞自動機進行了水體污染帶的模擬擴散研究。沈敬偉等[5]將元胞自動機與并行計算相結(jié)合, 綜合考慮了多種因素, 模擬了水體污染物擴散。李維乾等[6]結(jié)合3S技術(shù)將智能體引入元胞自動機模型, 進行了水污染模擬仿真。盡管有眾多學(xué)者早已利用元胞自動機研究海上溢油, 但考慮的影響因素不夠全面、影響系數(shù)比較復(fù)雜。因此, 本文綜合考慮了風(fēng)流、水流、蒸發(fā)、溶解、乳化、岸邊附著等影響因子, 對風(fēng)流、水流等影響系數(shù)進行了改進, 構(gòu)建了海洋溢油模型。

1 元胞自動機影響因子設(shè)定

本文不僅考慮了風(fēng)、流因子的作用, 還考慮了重力、浮力作用。另外綜合了岸邊附著、溶解、乳化等因子的影響。

1.1 風(fēng)流因子

風(fēng)、流是影響污染物擴散的最重要因素, 用影響系數(shù)表示風(fēng)、流狀態(tài)下不同鄰域元胞對中心元胞的影響, 該系數(shù)的計算由兩部分組成:

假設(shè)風(fēng)速的方向為從西向東, 則只有鄰域元胞(–1,,)對中心元胞(,,)產(chǎn)生影響, 此時在時風(fēng)產(chǎn)生的影響系數(shù)F中引入(,,)和(–1,,)處風(fēng)速的平均值, 并以該平均值與最大風(fēng)速的比值作為風(fēng)的影響系數(shù), 如公式(2)所示:

S為在時在某方向上水流產(chǎn)生的影響系數(shù),S修正為該方向上的水流速度與流域內(nèi)最大流速的比值。如公式(3)所示:

1.2 擴散因子

1.3 蒸發(fā)

其中,Z為蒸發(fā)率, 其方程表述為:

1.4 岸邊附著

其中,為油濃度;為吸附常數(shù),m為最大吸附量,為吸附指數(shù)。

1.5 溶解

溢油過程中有一部分會溶解于海水, 溶解于海水的溢油質(zhì)量r的公式為:

1.6 溢油乳化

乳化含水率W(%)利用公式(10)計算:

2 溢油模型的構(gòu)建

元胞自動機(Cellular Automaton, CA)是定義在一個有限的、離散狀態(tài)的元胞空間上并按照一定的局部轉(zhuǎn)化規(guī)則, 在離散的時間維上演化的動力學(xué)系統(tǒng)[12]。它由4部分構(gòu)成: 元胞、元胞空間、鄰域和狀態(tài)演化規(guī)則。用規(guī)則的幾何圖形將研究區(qū)域分為網(wǎng)格, 每一個網(wǎng)格就是一個元胞。所有的元胞按照一定規(guī)則排列組成的空間就是元胞空間。每個元胞都有自己的狀態(tài), 元胞的狀態(tài)是它周圍鄰域中的其他元胞的狀態(tài)共同作用決定的。鄰域類型一般有馮.諾伊曼鄰域[13](Von Neumann型)、摩爾(Moore)型、擴展摩爾型[14]三種, 如圖1所示。狀態(tài)演化規(guī)則是元胞自動機的核心, 是當前狀態(tài)元胞進行下一時刻元胞狀態(tài)轉(zhuǎn)化的變換函數(shù)。

其中,是4個正方向的擴散系數(shù);是4個斜角方向上的擴散系數(shù), 且當=0.16、=0.084時可得到最佳模擬效果。

考慮風(fēng)流對溢油的影響, 參照文獻[6]可將公式(11)修正為:

因此, 通過對以上各因子進行分析, 綜合考慮構(gòu)建了海洋溢油模型, 其公式為:

圖1 鄰域類型圖

圖2 質(zhì)量傳遞圖

3 元胞自動機擴散模擬

3.1 實驗數(shù)據(jù)及參數(shù)選取

3.2 模擬結(jié)果

通過以上條件及演化規(guī)則得出模擬結(jié)果如圖3, 并與初始圖像、檢驗圖像、對比圖像進行對比, 其中對比圖像為彭曉鵑等[15]的結(jié)果。從圖4—圖6可以看出溢油的方向相同, 吻合度較高。溢油范圍向四周擴散, 整體朝東南方向擴散。與彭曉鵑等[15]實驗中所得到的模擬結(jié)果進行了對比, 結(jié)果具有一定的改進(圖7)。

圖3 模擬結(jié)果

圖4 初始圖像與模擬結(jié)果

3.3 模擬結(jié)果檢驗

模型構(gòu)建完成后需要進行檢驗, 本文通過混淆矩陣[17]計算了CA模型的模擬精度, 并通過Kappa系數(shù)a[18]定量地檢驗?zāi)Ｐ偷哪M精度。

圖5 模擬結(jié)果與及檢驗圖像

圖6 與檢驗、初始圖像對比圖

圖7 與彭曉鵑實驗結(jié)果對比圖

表1 精度模擬及Kappa系數(shù)

Tab.1 Accuracy simulation and kappa coefficient

4 結(jié)語

該模型模擬溢油的擴散漂移, 不僅考慮了風(fēng)、流等因素, 而且考慮了溢油消失過程, 綜合了蒸發(fā)、乳化、溶解、岸邊附著等溢油行為, 可以更加準確地模擬出溢油的時空變化過程。模擬了“Deepspill”海上溢油進行驗證, 實驗表明模擬結(jié)果的總精度、Kappa系數(shù)均較高, 因此認為該模型可以用來模擬海洋溢油的過程。

[1] Abascal A J, Castanedeo S, Mendez F J, et al. Calibration of a lagrangian transport model using drifting buoys deployed during the prestige oil spill[J]. Journal of Coastal Research, 2009, 25(1): 80-90.

[2] Mishra D R, Cho H J, Ghosh S, et al. Post-spill state of the marsh: Remote estimation of the ecological impact of the gulf of Mexico oil spill on Louisianan Salt Marshes[J]. Remote Sensing of Environment, 2012, 118(6): 176-185.

[3] Karafyllidis I. A model for the prediction of oil slick movement and spreading using cellular automata[J]. Environment International, 1997, 23(6): 839-850.

[4] 王璐, 謝能剛, 李銳, 等. 基于元胞自動機的水體污染帶擴散漂移仿真[J]. 水利學(xué)報, 2009, 4: 481-485. Wang Lu, Xie Nenggang, Li Rui, et al. Simulation of drift-diffusion of water pollution zone based on cellular automata[J]. Shuili Xuebao, 2009, 4: 481-485.

[5] 沈敬偉, 彭安琪, 周廷剛, 等. 基于并行元胞自動機的水體污染物擴散模擬[J]. 測繪科學(xué)技術(shù)學(xué)報, 2016, 1: 105-110. Shen Jingwei, Peng Anqi, Zhou Tingkang, et al.Water Pollutant Spreading Simulation Based on Parallel Cellular Automata[J]. Journal of Geomatics Science and Technology, 2016, 1: 105-110.

[6] 李維乾, 解建倉, 李建勛, 等. 基于元胞自動機與智能體的水污染可視化模擬仿真[J]. 西北農(nóng)林科技大學(xué)學(xué)報(自然科學(xué)版), 2013, 3: 213-220. Li Weiqian, Jiancang, LiJianxun, et al. Visualization simulation of water pollution based on Cellular Automata and intelligent agent[J]. Journal of Northwest A & F University(Nat.Sci.Ed), 2013, 3: 213-220.

[7] 劉彥呈, 殷佩海, 林建國, 等. 基于GIS的海上溢油擴散和漂移的預(yù)測研究[J]. 大連海事大學(xué)學(xué)報, 2002, 3: 41-44. Liu Yancheng.Yiin Peihai, Lin Jianguo, et al. Prediction of oil spill spreading and transport over the sea[J]. Journal of DaLian Maritime University, 2002, 3: 41-44.

[8] 張存智, 竇振興, 韓康, 等. 三維溢油動態(tài)預(yù)報模式[J]. 海洋環(huán)境科學(xué), 1997, 1: 26-33. Zhang Cunzhi, Dou Zhengxing, Han Kang, et al. A Three Dimensional Model to Predict the Behavior of Oil Spills[J]. Marine Environmental Science, 1997, 1: 26-33.

[9] 李崇明, 趙文謙, 羅麟. 河流泥沙對石油的吸附、解吸規(guī)律及影響因素的研究[J]. 中國環(huán)境科學(xué), 1997, 1: 25-28. Li Chongming, Zhao Wenqian, Luo Lin. Study on characteristics and effect factors of the absorption desorption of oil by sediment in rivers[J]. China Environmental Science, 1997, 1: 25-28.

[10] 莊學(xué)強, 陳堅, 孫倩. 海面溢油數(shù)值模擬及其可視化實現(xiàn)技術(shù)[J]. 中國航海, 2007, 1: 97-100. Zhuang Xueqiang, Chen Jian, Sun Qian. Numerical Simulation and Visualization Technology of Marine Spilled Oil[J]. Navigation of China, 2007, 1: 97-100.

[11] 過杰, 孟俊敏, 何宜軍. 基于二維激光觀測的溢油及其乳化過程散射模式研究進展[J]. 海洋科學(xué), 2016, 2: 159-164. Guo Jie, Meng Junmin, He Yijun. Scattering model research based on two-dimensional laser observation of spilled oil and emulsification processes[J]. Marine Sciences, 2016, 2: 159-164.

[12] 周成虎, 歐陽, 馬廷, 等. 地理系統(tǒng)模擬的CA模型理論探討[J]. 地理科學(xué)進展, 2009, 6: 833-838. Zhou Chenghu, Ou Yang, Ma Ting, et al. Theoretical Perspectives of CA-based Geographical System Modeling[J]. Progress in Geography, 2009, 6: 833-838.

[13] 歐敏, 張永興, 胡居義, 等. 基于Geo-CA和GIS的滑坡穩(wěn)定性分析[J]. 中國地質(zhì)災(zāi)害與防治學(xué)報, 2004, 4: 4-9. Ou Min, Zhang Yongxing, Hu Juyi, et al.The analysis of landslide stability based on Geo-CA and GIS [J].HONGGUO DIZHIZAIHAI YU FANGZHI XUEBAO, 2004, 4: 4-9.

[14] 羅平, 耿繼進, 李滿春, 等. 元胞自動機的地理過程模擬機制及擴展[J]. 地理科學(xué), 2005, 6: 6724-6730. Luo Ping, Geng Jijin, Li Manchun, et al. M echan ism of Smi u lating Geographic Process and Extension of Cellular Automata[J].Scientia Geographica Sinica, 2005, 6: 6724-6730.

[15] 彭曉鵑, 張亦漢. 基于元胞自動機的海上溢油擴散模擬[J]. 海洋通報, 2015, 4: 415-422. Peng Xiaojuan, Zhang Yihan. Simulation of marine oil spill diffusion based on cellular automata[J].Marine Science Bulletin, 2015, 4: 415-422.

[16] 楊慶霄, 徐俊英, 李文森. 海上溢油溶解過程的研究[J].海洋學(xué)報(中文版), 1994, 3: 50-56. Yang Qingxiao, Xu Junying, Li Wenshen. Study of dissolution of marine oil spill[J]. Acta Qceanologica Sinica, 1994, 3: 50-56.

[17] 李宜展, 潘耀忠, 朱秀芳, 等. 土地覆蓋類別面積混淆矩陣校正與回歸遙感估算方法對比[J]. 農(nóng)業(yè)工程學(xué)報, 2013, 11: 115-123. Li Yizhan, Pan Yaozhong, Zhu Xiufang, et al. Comparison analysis on land cover area estimators: confusion matrix calibration and regression[J]. Transactions of the Chinese Society of Agricultural Engineering, 2013, 11: 115-123.

[18] 全泉, 田光進, 沙默泉. 基于多智能體與元胞自動機的上海城市擴展動態(tài)模擬[J]. 生態(tài)學(xué)報, 2011, 10: 2875-2887. Quan Quan, Tian Guangjin, Sha Moquan.Dynamic simulation of Shanghai urban expansion based on multi-agent system and cellular automata model[J].Acta Ecologica Sinica, 2011, 10: 2875-2887.

(本文編輯: 李曉燕)

Marine oil spill model based on cellular automata

MAN Miao-miao1, LIU Lin1, 2, CHENG Peng1, WEI Guo-zhong3, ZHANG Sheng3

(1. Geomatics College, Shandong University of Science and Technology, Qingdao 266590, China; 2. Key Laboratory of Surveying and Mapping Technology on Island and Reef, National Administration of Surveying, Mapping and Geoinformation, Qingdao 266590, China; 3. Shangdong Provincial Institate of Land Surveying and Mapping, Jinan 250102, China)

Based on the three-dimensional cellular automata model, this paper studies oil spill. The influence coefficients of the wind and flow in the horizontal direction are improved, in addition to considering the influence of evaporation and vertical diffusion in the vertical direction; the shore attachment, dissolution factor, and the emulsified water content parameters are ased in order to achieve a practical quantitative emulsifying factor. Based on the experiment “DeepSpill”, the simulation precision is considerably improved.

cellular automata; marine oil spill; simulation

Apr. 13, 2017

滿苗苗(1992-), 女, 山東濟寧人, 在讀碩士研究生, 研究方向: 地理信息處理及應(yīng)用, 電話: 13206427649, E-mail: 1509443062@qq.com; 柳林,通信作者, 副教授, 主要從事LBS、海洋GIS、智慧城市、3S集成等方面的研究, 電話: 0532-80681183, E-mail: liulin2009@126.com

X55

A

1000-3096(2017)09-0021-06

10.11759/hykx20170413003

2017-04-13;

2017-08-08

山東省自然科學(xué)基金(ZR2012FM015); 海島(礁)測繪技術(shù)國家測繪地理信息局重點實驗室資助項目(2013B08); 衛(wèi)星測繪技術(shù)與應(yīng)用國家測繪地理信息局重點實驗室經(jīng)費資助項目(KLAMTA-201407)

[Shandong Provincial Natural Science Foundation, No.ZR2012FM015; the Key Laboratory of Surveying and Mapping Technology on Island and Reed, State Bureau of Surveying, Mapping and Geoinformation, No.2013B08; the Key Laboratory of Satellite Mapping Technology and Application, National Administration of Surveying, Mapping and Geoinformation , No. KLAMTA-201407]