溢油擴展、漂移及擴散預測技術研究進展

李歡，邵偉增，李程，張松，李文善，徐珊珊

（1.國家海洋信息中心，天津300171；2.河海大學江蘇省海岸海洋資源開發與環境安全重點實驗室，江蘇南京210098；3.浙江海洋大學，浙江舟山316022；4.國家海洋技術中心，天津300112）

溢油擴展、漂移及擴散預測技術研究進展

李歡1,，2，邵偉增3，李程1，張松4，李文善1，徐珊珊1

（1.國家海洋信息中心，天津300171；2.河海大學江蘇省海岸海洋資源開發與環境安全重點實驗室，江蘇南京210098；3.浙江海洋大學，浙江舟山316022；4.國家海洋技術中心，天津300112）

海上溢油災害會嚴重破壞中國近海海洋環境，還會直接危害我國經濟發展。因此，開展溢油預測預警技術研究是非常有必要的，能夠為海上溢油應急響應提供技術支撐。本文綜述了溢油擴展、漂移及擴散數值預測技術的發展過程及相關研究成果，包括溢油擴展模型、溢油漂移與擴散模型、溢油預測預警系統，為未來開展溢油應急工作提供了理論依據和參考。

海上溢油；溢油擴展模型；溢油漂移與擴散模型；溢油預測預警系統

我國能源需求的不斷擴大促使海洋石油運輸業和開采業飛速發展，油輪不斷增多，海上石油勘探、開發、海底管線鋪設規模不斷加大，同時海上溢油事故，尤其是重特大事故的風險亦隨之大大增加。各種海上溢油災害，不僅會嚴重破壞我國近海海洋生態環境，還會直接危害我國社會經濟的正常發展。因此，開展海上溢油應急關鍵技術尤其是溢油漂移與擴散預測預警技術研究，為溢油事故應急響應、處置提供決策支持的技術平臺，提升溢油應急響應能力和技術水平是非常緊迫和必要的。不僅為保護我國海洋環境安全提供技術支撐，也將促進我國海上航運業的安全穩定發展，為我國經濟社會的持續穩定保駕護航。

本文從溢油擴展、漂移及擴散的變化特點入手，描述已有的數值預測方法、模型系統，綜述溢油主要物理行為變化預測研究，為廣大溢油應急相關部門工作人員提供有益的參考。

1 溢油擴展、漂移及擴散預測

其中，D0為初始時刻油膜直徑；dw、do分別為水和油的比重；Kr為Blokker常數；V為溢油總體積；為時間。該模型雖然可以初步反映油膜的物理擴展狀態，但由于其忽略了表面張力和粘性力的作用，對于油膜擴展的刻畫較為片面。針對其缺陷，Fay（1969）提出了自己的油膜擴展理論。該理論成立的前提條件是海面平靜的狀態下，大體積溢油在擴展過程中油的性質不變且在垂向受力平衡，油膜以圓形擴展，擴展范圍可以用圓的直徑來衡量。在此

溢油的擴展指油膜由自身物理特性引發的面積增大；溢油的漂移與擴散指油在海面風、海流以及海浪等海洋環境動力要素的共同作用下發生的遷移輸運。對于溢油本身，以上幾個過程緊密相關，實際過程中因考慮到海洋環境的不斷變化而較為復雜，也給針對這些過程的數值預測帶來了較大的難度和挑戰性，目前國際上主流的方法是對上述過程進行階段性模擬，雖然可以在一定程度上完成溢油主要物理行為的預測，但是結果與真實情況仍存在較大的誤差。下文將針對溢油的幾種主要物理變化過程及預測技術進行重點闡述。

1.1 溢油擴展預測

油溢至海面時，受到表面張力、慣性力、重力和粘性力的共同作用，會向四周進行擴展，重力和表面張力為驅動力，慣性力和粘性力為阻力。重力在初始階段起主要作用，油的擴展情況很大程度上由溢出方式決定。靜止點源瞬時大量溢油的擴展速率遠快于靜止或者移動點源連續溢油的。隨著時間推移，油層厚度的逐漸減小，表面張力逐漸成為影響溢油擴展的主要因素，油在水面形成了中間比邊緣厚的膜。擴展后的油膜在海洋環境動力要素的非線性作用下，形態隨時間變化不斷改變，最終以碎片、油滴等形式存在于海水中。

早期，Blokker（1964）建立了自由面上油膜擴展的數學模型：前提下，油膜擴展過程根據各種力在不同階段所起的作用劃分為3個階段，即重力和慣性力共同作用階段、重力和粘性力共同作用階段以及表面張力和粘性力共同作用階段。

Fay三階段模型表達式如下：

上式中，D為油膜直徑；K1、K2、K3分別為各擴展階段的經驗系數；g為重力加速度；V為溢油體積；t為時間；ρw為海水密度；vw為水的運動粘滯系數；β=1-ρ0/ρw，ρ0為油的密度；δ為凈表面張力系數，即水與空氣間、油與空氣間、油與水間的表面張力系數。上述各階段的時間區分點根據相鄰兩階段擴展直徑相等的狀態來確定。當δ為0時，表示油膜擴展直徑達到穩定態，擴展過程結束，油膜面積由如下公式得出：

對于連續溢油的油膜擴展，仍然可以用Fay三階段模型表示：

上式中，Q為單位時間的溢油量；u為流速；δ為凈表面張力系數；ρw為海水密度；vw為水的運動粘滯系數；Δ=1-ρ0/ρw，ρ0為油密度。Fay三階段模型提出之后，諸多研究涉及對其的優化與改進，衍生出許多新的溢油擴展模型。Mackay（1980a，b）在Fay三階段模型中考慮了風的影響，并基于觀測數據分別建立了針對厚油膜和薄油膜的擴展模型。Lehr等（1984）針對Fay三階段模型全面考慮了海面風場、流場對油膜的影響，發現油膜在海面實際是以橢圓而非圓形的形式擴展，且長軸的方向與風向保持一致。趙文謙等（1992）基于Fay理論提出了的油膜擴展模型中考慮了油膜擴展和各項異性擴散作用以及油膜邊緣的消失過程，通過建立油膜在主要擴展方向和與之垂直的次要擴展方向上的尺度公式，從而獲得油膜的擴展面積。

1.2 溢油漂移與擴散預測

海面溢油的漂移與擴散過程主要受海洋動力環境要素變化的影響。溢油漂移與擴散預測的首要工作是獲得精確的海洋動力環境預報參數，重點考慮以下兩點：（1）表層海流和海面風場的直接影響；（2）非線性波產生的Stokes流的影響。因此，針對海流、海浪以及海面風場的精細化預報尤為重要，通過三維海流數值模型、海浪數值模型以及大氣數值模型獲得精準的海洋環境動力要素是溢油漂移與擴散預測的根本保證。

1.2.1 海洋和大氣數值模型

目前，國內外已開發了多種三維海流數值模型，根據空間離散方法的不同，可將這些模型分為有限差分類、有限元類和有限體積類。有限差分類模型發展較早，較為成熟，這一類模型通過有限差分將各控制方程中的導數用網格節點上的函數值的差商代替進行離散，從而建立以網格節點上的值為未知數的代數方程組，直接將微分問題變成代數問題的近似數值解法。這類模型計算速度快，應用廣泛，其中具有代表性的有POM（Mellor，1998）、ECOM（張越美等，2001）、ROMS（Moore et al，2001）等。然而，這類模型也存在明顯缺點，即對不規則區域的適應性較差，岸界彌合的能力明顯不足。與有限差分類模型不同，有限元類模型將計算域劃分為有限多個互不重疊的單元，在每個單元內，選擇一些合適的節點作為求解函數的插值點，將微分方程中的變量改寫成由各變量或其導數的節點值與所選用的插值函數組成的線性表達式，借助于變分原理或加權余量法，生成有限元離散方程，最后結合邊界條件和初始條件，求解由各個單元的離散方程組成的總體方程組的近似解。這類模型最大的優點是網格剖分靈活，能夠精確擬合復雜岸線，進行局部網格加密，計算精度比有限差分類模型高，具有代表性的有SHYFEM（Umgiesser，2014），但這類模型計算所需存儲量大，計算效率相對較低（圖1所示）。與有限元類模型類似，有限體積類模型的設計思路是將計算區域劃分為一系列互不重疊的控制體，這些控制點均以網格點為中心，將待解的微分方程對每一個控制體積分得到離散方程，結合邊界條件和初始條件求得數值解。由于這類模型在水平方向上采用了非結構化三角形網格，吸收有限差分和有限元類模型的優點，不僅可以對復雜的岸線精確的擬合，進行局部網格加密，還因為便于對原始方程組進行離散差分，很大程度上確保了較高的計算效率，這類模型中具有代表性的有FVCOM（Chen et al，2003；Zheng et al，2003）、SELFE（Zhang et al，2008）等。有限體積類模型的諸多優點使其發展成為精細化海流預報的主要模型之一，目前已經廣泛應用于溢油漂移與擴散預測中。

圖1 SHYFEM模型的非結構化網格配置（Umgiesser，2014）

對于海浪數值模型，目前國際上主流的有WaveWatch III和SWAN等（徐麗麗等，2015），這類模型能夠對海浪進行精細化預報。尤其是非結構化網格的SWAN，同樣確保了海浪預報的精準度。

對于大氣數值模型，目前主流的仍然是WRF（章國材，2014）。該模型的WRF-ARW版本，擁有先進的資料三維變分同化技術，能充分有效地將各種資料信息同化到模型初始場中，為模型提供更高質量的初始值，達到明顯改進數值預報質量的目的。

1.2.2 溢油漂移動態預測模型

溢油漂移動態預測模型發展至今，只存在兩類模型：一是基于歐拉-拉格朗日漂移理論建立的軌跡追蹤模型。這類模型中，油膜在表面流場和風的共同作用下發生平移輸運，通過歐拉-拉格朗日方法可以預測水面溢油的全動力軌跡。這類模型發展較早，具有代表性的有Miller等（1977）研發的Coast Guard模型、Williams等（1975）研發的SEADOCK模型以及Webb等（1970）研發的美國海軍Navy模型等。二是基于蒙特卡羅方法的溢油漂移與擴散軌跡預測模型。這類模型中，油膜被劃分為多個質點，質點的擴散看作是在湍流作用下的隨機擴散，表現為質點群沿著復雜軌跡進行的布朗運動，其運動方向在某一時刻是隨機的，但這種隨機性又受到整個運動場的控制?？梢栽诮o予隨機數的同時，給定湍流強度、時間尺度和質點數，可以獲取標識質點的位移（婁安剛等，2000）。這類模型中最有代表性的為油粒子模型。

油粒子模型中，油膜被離散化成大量的粒子，每個粒子代表一定油量，在海洋環境動力要素的作用下漂移，剪切流和湍流對溢油運動產生的作用通過粒子的隨機運動實現。而油的蒸發、乳化等風化過程則由油粒子的質量損失來體現；油膜厚度分布通過一定水面面積內的油粒子的個數、體積、質量來體現。這些質點相互獨立的運動，在進行下一步前所有粒子團都必須跟蹤到同一時刻。油膜的擴散即為粒子群的輸運擴散，受到油粒子尺寸分布、剪切流和湍流過程的控制，也可看作粒子群在湍作用下進行的拉格朗日運動，通過粒子的隨機運動來實現。使用粒子概念研究對流擴散過程，系統總質量是守恒的，計算原理是跟蹤粒子的路徑，而不是求解方程組，因此油粒子模型可直接模擬擴散方程的實際物理現象（牟林等，2011）。

油粒子模型將油膜漂移和擴散聯系在一起，通過油粒子質點坐標的追蹤即可計算溢油漂移軌跡和影響區域。溢油污染影響面積通過油膜掃海面積確定，即所有油粒子在一定時間內經過的海域面積。目前，已有研究將油膜擴展和漂移擴散幾個階段進行了銜接，預測結果更為合理（劉偉峰等，2009）。

1.3 溢油預測預警系統

目前國際上主流的溢油預測預警系統有美國的OILMAP（http：//www.asascience.com/software/ oilmap/）、英國的OSIS（http：//osisinternational. com/the-osis-oil-spill-solution）、挪威的OSCAR（http：//www.oscar-insurance.com/）、意大利的BOOM（http：//www.medess4ms.eu/）以及西班牙的MOHID-TESEO（http：//www.mohid.com/）等（De Dominicis et al，2013）。國內則有國家海洋信息中心的中國海洋溢油一體化多節點協同預測預警系統、中海石油環保服務有限公司的中國近海海上溢油預測與應急決策支持系統、大連海事大學的海上溢油應急預報信息系統、國家海洋環境預報中心的渤海溢油應急預報系統（http：//www.nmefc.gov.cn/ cgofs/oilspill.aspx）等。以上溢油預測預警系統的總體架構大同小異，具體表現在如下幾點：

（1）系統均實現了模塊化集成。將系統根據功能的不同分為若干模塊，每個模塊間既相互獨立又相互關聯，例如：海流模塊負責海流的預報，海浪模塊負責海浪的預報，大氣模塊負責海面風場的預報，溢油預測模塊負責預測溢油的行為預測，且各模塊間相互傳遞著信息，表現在：大氣模塊預報的風場驅動海流模塊和海浪模塊運行，預報的流場和風場信息同時傳遞給溢油預測模塊進行溢油行為預測。模塊化集成的最大優點可以使系統高效運行，因各模塊分工運作，且彼此間進行有序的信息傳遞，大大提高了溢油應急響應的速度。

（2）油粒子模型在各系統中的廣泛采用。油粒子模型已成為溢油預測的主流模型，并且在各系統中，預測溢油風化的不同經驗方程已與油粒子模型完全耦合，使溢油行為與歸宿預測結果與實際更為接近。

（3）系統均實現了可視化。由于各系統均以GIS平臺作為支撐，GIS技術可以向應急人員提供各類溢油數據的可視化信息，主要方式是將基礎地理數據、海洋環境動力要素數據、溢油漂移與擴散預測結果等信息進行分層疊加展示（見圖2），使應急人員可以更直觀地獲得溢油的現狀和未來發展趨勢。

圖2 中國海洋溢油一體化多節點協同預測預警系統展示界面

以上系統之間的區別則主要表現在海域針對性和附加功能上。對于OILMAP這樣的商業系統，或者GNOME這樣的開源系統，只要能夠提供相應的基礎地理信息，在任何海域均可以進行溢油應急預測預警。但對于OSIS、BOOM、MOHID這樣非商業和開源的系統，幾乎只用來針對其所屬國家近海海域的溢油預測預警任務，包括我國的渤海溢油應急預報系統、中國海洋溢油一體化多節點協同預測預警系統以及中國近海海上溢油預測與應急決策支持系統，也具有特定的針對性。在附加功能上，有些系統，除了能夠進行溢油預測預警，還增加了有關溢油環境敏感區域預警功能，例如中國海洋溢油一體化多節點協同預測預警系統。還有一些系統可依據應急計劃和應急反應專家的知識和經驗，根據溢油預測結果，形成應急反應決策輔助方案，例如OSIS、中國近海海上溢油預測與應急決策支持系統等。

綜上所述，不同的溢油預測預警系統在各沿海國家溢油應急響應工作中的作用極其重要，也代表了海上溢油應急預測預警領域的最高水準。

2 結論

本文從本文從溢油擴展、漂移及擴散的變化特點入手，描述已有的數值預測方法、模型系統，綜述溢油主要物理行為變化預測研究，得出以下結論：

（1）針對溢油物理擴展的預測，Fay理論及其數值模型起到了重要作用，眾多新擴展模型都是在其基礎上的優化與改進。

（2）針對溢油漂移和擴散預測，首要工作是海洋環境動力參數的精細化預報，選取合適的數值預報模型十分重要；其次是油粒子模型的采用，使溢油漂移與擴散的預測更加合理。

（3）國內外溢油預測預警模型系統大同小異，均實現了模塊化，幾乎都采用了油粒子模型且均以GIS平臺作為支撐平臺來展示溢油的現狀和未來的發展動態，不同的溢油預測預警系統在各沿海國家溢油應急響應工作中的作用極其重要，也代表了海上溢油應急預測預警領域的最高水準。

參考文獻

Blokker P C,1964.Spreadingand evaporation of petroleum products on water.Proc,4th International Harbour Congress,Antwerp,the Netherlands,911-919.

Chen C,Huang H,Beardsley R C,et al,2003.An Unstructured Grid,Finite-Volume,Three-Dimensional,Primitive Equations Ocean Model：Application to Coastal Ocean and Estuaries.Journal of Atmospheric&Oceanic Technology,20(1):159-186.

De Dominicis M,Pinardi N,Zodiatis G,et al,2013.MEDSLIK-II,a Lagrangian marine surface oil spill model for short-term forecasting-Part 1:Theory.Geosci.Model Dev.,6:1851-1869.

Fay J A,1969.Oil on the Sea.New York:Plenum Press,53-63.

Lehr W J,Cekirge H,Fraga R,et al,1984.Empirical studies of the spreading of oil spills.Oil and Petrochemical Pollution,2(1):7-11.

Mackay D,Buist I,Mascarenhas R,et al,1980b.Oil spill processes and models.Environment Canada Report EE-8.

Mackay D,Paterson S,Trudel K,1980a.A mathematical model of oil spill behavior.Environment Canada Report EE-7.

Mellor G L,1998.User guide for a three-dimensional,primitive equation, numerical ocean model.Program in Atmospheric and Ocean Science Princeton University,Princeton.

Miller and Stolzenback K D,1977.A review and evaluation of basic techniques predicting thebehavior of surface oil slicks.MIT Report.

MooreA M,ArangoH G,Broquet G,et al,2001.The Regional Ocean Modeling System(ROMS)4-dimensional variational data assimilation systems:Part II-Performance and application to the California Current System.Progress in Oceanography,91(1):74-94.

Umgiesser G,2014.Finite Element Model for Coastal Seas User Manual. Oceanography,ISMAR-CNR.

Webb L,Taranto R,Hashimoto E,1970.Operational Oil Spill Drift Forecasting.In:proceedings Of 7th Nary Symposium of Military O-ceanography,Annapolis,Maryland,114-119.

Williams G N,Hann R W,1975.Simulation Models for Oil Spill Transport and diffusion.Summer ComputerSimulation Conference,748-752.

Zhang Y,Baptista A M,2008.SELFE:a Semi-implicit Eulerian-Lagrangian Finite-Element model for cross-scale ocean circulation[J]. Ocean Modelling,21(3-4):71-96.

Zheng L,Chen C,Liu H,2003.A modeling study of the Satilla River estuary,Georgia.I:Flooding-drying process and water exchange over the salt marsh-estuary-shelf complex.Estuaries,26(3):651-669.

劉偉峰，孫英蘭，2009.海上溢油運動數值模擬方法的討論與改進.華東師范大學學報（自然科學版），（3）：91-97.

婁安剛，王學昌，于宜法，等，2000.蒙特卡羅方法在海洋溢油擴展預測中的應用研究.海洋科學，24（5）：7-9.

牟林，趙前，2011.海洋溢油污染應急技術.北京，科學出版社.

武周虎，趙文謙，1992.海面溢油擴展、離散和遷移組合模型.海洋環境科學，11（3）：35-42.

徐麗麗，肖文軍，石少華，等，2015.WaveWatch III和SWAN模型嵌套技術在業務化海浪預報系統中的應用及檢驗.海洋通報，2015（3）：283-294.

張越美，孫蘭英，2001.ECOM模型在丁字灣中的應用.青島海洋大學學報，31（5）：659-664.

章國材，2014.美國WRF模式的進展和應用前景.氣象，30（12）：27-31.

（本文編輯：袁澤軼）

Research on numerical prediction technologyof oil spill spreading，drift and diffusion

LI Huan1,2,SHAO Wei-zeng3,LI Cheng1,ZHANG Song4,LI Wen-shan1,XU Shan-shan1
（1.National MarineData and Information Service,Tianjin300171,China;2.Jiangsu Key Laboratory of Coast Ocean Resources Development and Environment Security,Hohai University,Nanjing210098,China;3.ZhejiangOcean University,Zhoushan 316022,China;4.National Ocean Technology Center,Tianjin 300112,China）

Marine oil spill may bring irreversible damage to the marine environmentof China Seas,and will cause immeasurable economical loss.Therefore,there is strong needs to study on oil spill prediction and early warning technology, providing technical support for marine oil spill response.This paper briefly introduced the development and achievements of the numerical prediction technology of oil spill spreading,drift and diffusion,including the oil spill spreading model,the oil spill drift and diffusion model,and the oil spill prediction and early warning system,which provide theoretic criterions for the future work on oil spill response.

Marine oil spill;Oil spill spreading model;Oil spill drift and diffusion model;Weathering Model;Oil spill prediction model system

X55

A

1001-6932（2017）04-0379-06

10.11840/j.issn.1001-6392.2017.04.003

2015-12-14；

2016-07-26

天國家科技支撐計劃項目（2014BAB12B02）；浙江省公益性技術應用研究工業項目（2015C31021）；天津市科技支撐計劃項目（14ZCZDSF00012）；江蘇省海岸海洋資源開發與環境安全重點實驗室開放基金資助項目（JSCE201505）

李歡（1983-），碩士，助理研究員，主要從事海洋防災減災研究。電子郵箱：usher02@126.com。

張松，碩士，助理研究員。電子郵箱：zhangsong1006@163.com