,,
(1.國家海洋局第一海洋研究所,山東 青島 266061;2.海洋環(huán)境科學(xué)和數(shù)值模擬國家海洋局重點實驗室,山東 青島 266061;3.青島海洋科學(xué)與技術(shù)國家實驗室區(qū)域海洋動力學(xué)和數(shù)值模擬功能實驗室,山東 青島 266061)
海上溢油是造成海洋環(huán)境污染損害的主要因素之一。在各種海洋污染中,石油污染無論在發(fā)生頻率、分布廣度還是在危害程度上均居首位。隨著世界海洋運輸業(yè)的發(fā)展和海上油田不斷投入生產(chǎn),溢油事故發(fā)生頻率不斷升高,給海洋安全及人民健康帶來重大損害[1-2]。
我國近海發(fā)生過多次嚴重的海上溢油事故,如2003年11月在渤海灣發(fā)生的“塔斯曼?!陛喴缬褪鹿?,生態(tài)索賠額高達1.2億元人民幣。又如2006年2月發(fā)生的“長島海域油污染事件”,范圍波及山東、天津、河北沿岸,是有史以來渤海發(fā)生的影響范圍最大的油污染事件。在所有這些溢油事故中,對殘存油量的估計都是非常重要的一個問題。
在多種作用因素中,波浪破碎導(dǎo)致的混合攪拌效應(yīng)通過影響油滴垂向混合,從而對溢油海面殘存油量估計起到了重要作用。觀測表明,高海況下被破碎波卷入水下的溢油量可以達到溢油總量的50%以上[3],所以研究波浪破碎對油滴的垂直混合是研究海上溢油海面殘油量隨時間、空間變化的重要基礎(chǔ),只有充分了解并掌握這一物理過程,才能夠全面地分析和預(yù)報海上溢油的變化趨勢,可以對溢油分散劑的使用量有一定的指導(dǎo)意義。同時,也可以為溢油事故發(fā)生后產(chǎn)生的經(jīng)濟賠償和生態(tài)損失評估提供強有力的理論基礎(chǔ)和數(shù)據(jù)證據(jù)。
對波浪破碎與油滴垂直混合之間關(guān)系的研究,之前主要采用了觀測數(shù)據(jù)分析和理論研究兩種研究手段。Delvigneand等[4]、Tkalich等[5]和孫寶楠等[6]根據(jù)實驗室數(shù)據(jù)并結(jié)合動力推導(dǎo)給出了波浪破碎與油滴垂向運動之間的一些經(jīng)驗和理論關(guān)系表達形式。溢油事故往往發(fā)生突然,并沒有非常實時的現(xiàn)場觀測波浪數(shù)據(jù)以供分析采用。在這種情況下,數(shù)值模式是一種非常有效的應(yīng)急補充方法,可以由計算得到的各類波浪破碎參數(shù),結(jié)合理論分析方法,給出具體時段和地點的海面溢油事故發(fā)生后油滴垂直混合情況,進而最終得到海面殘油量變化趨勢。
之前的分析研究因為現(xiàn)場觀測數(shù)據(jù)的缺乏,導(dǎo)致計算結(jié)果的可靠性存在問題。隨著衛(wèi)星觀測技術(shù)手段的豐富,雖然定量估計殘油量仍然難度較大,但是定性評估海面殘油量已經(jīng)成為可能[7-8]。這使得結(jié)合理論分析和數(shù)值模擬的方法,針對某一具體海上溢油時間開展波浪破碎過程對油滴垂直混合影響研究具有一定可行性。
針對波浪破碎過程對油滴垂直混合影響開展研究,必須首先建立能夠描述包括油滴入水和上浮兩個垂向運動過程在內(nèi)的統(tǒng)一動力方程。同時,為了最終給出海面殘油量的估計結(jié)果,該動力方程必須結(jié)合溢油蒸發(fā)模型,共同構(gòu)成描述海面殘油量變化的數(shù)學(xué)模型。
根據(jù)Tkalich等[5]的一階模型,把油滴垂向運動方程描述為
(1)
式中:Ms為海面浮油的質(zhì)量;Me為海面以下殘油的質(zhì)量(都是在單位海表面積內(nèi));Qo為單位時間單位水體破碎卷入的油滴體積(即卷入率)。
入水油的總質(zhì)量為:Ma=Ms+Me。在溢油之初,Ma等于溢漏出來的油的初始質(zhì)量M0。Qo在卷入過程中起著決定性的作用。下面用相似性定理分析Qo的表達形式,其中參與量綱分析的變量有4組。
③Z,油滴卷入深度,量綱:L;
基本量綱為長度量綱L、時間量綱T和質(zhì)量量綱M,共3個,根據(jù)關(guān)鍵變量形成其他變量的量綱一的量方程,求解得到量綱相容方程。
(2)
因為只有一個量綱相容方程,即只得到一個量綱一的量組,因此,有函數(shù)形式如下。
(3)
式中:f(x)為待定的函數(shù)形式。
選取待定函數(shù)為冪函數(shù)形式,即令f(x)=cxk,可以得到
(4)
選取最簡單的形式,待定系數(shù)c=1,k=1。式中Z和E分別為與波浪破碎過程有關(guān)的卷入深度和破碎能量損耗率。最終可以得到作用在單位海表面的Qo表達形式為
(5)
取kb=0.5,g=9.8 m/s2,ρw=1 025 kg/m3,Low=1 m。
油滴上浮的動力過程可用一階方程描述。
(6)
取Lwo=20 m,rc=50 μm,rmax=500 μm,
ρ=900 kg/m3,υ=10-5m2/s,p=2。
溢油蒸發(fā)方案采用Stiver等[9]提出的分析法模型。
(7)
綜合以上3個子模型,根據(jù)Tkalich等[5]的簡化方法,海面殘油量變化模型控制方程可寫為
(8)
通過已建立的運動學(xué)方程發(fā)現(xiàn),波浪破碎主要通過影響Qo對油滴垂直運動產(chǎn)生影響,而該要素主要受與波浪破碎有關(guān)的卷入深度(Z)和破碎能量損耗率(E)這兩個變量共同控制。針對這兩個要素,有多種計算方案可供選擇。
針對Z,Delvigne等[4]通過試驗室擬合給出的經(jīng)驗表達形式:ZD&S=αHb,其中α=1.5,Hb為破碎波高,取Hb=0.7Hs,Hs為有效波高。Yuan等[11]根據(jù)實際海況導(dǎo)出了理論表達形式。
(9)
(10)
為了分析采用不同計算方案對結(jié)果的最終影響,設(shè)計3組對比實驗,見表1。
表1 數(shù)值試驗方案設(shè)計
2006年3月24日渤海灣灤河河口以南曹妃甸附近發(fā)生了一起海洋溢油事故。根據(jù)當日獲取的歐洲空間局環(huán)境衛(wèi)星(Envisat)合成孔徑雷達圖像判斷,該海域形成了大面積的油膜污染,總面積達400 km2。由2006年4月1日獲取的Envisat圖像判斷,渤海灣發(fā)生的溢油現(xiàn)象已經(jīng)減輕,沒有發(fā)現(xiàn)形態(tài)完整、溢出時間較短的油膜。由于波浪破碎攜油入水,蒸發(fā)等的共同作用,油膜已經(jīng)分散,質(zhì)量明顯減少,溢油面積減少到100 km2左右,殘油量減少了初始質(zhì)量的75%以上[7]。
為了再現(xiàn)溢油時段的渤海海區(qū)波浪平面分布,應(yīng)用球坐標系下MASNUM海浪數(shù)值模式[13]模擬渤海的波浪場。水平分辨率是(1/12)°×(1/12)°,溢油時段為6 d。風場資料采用的是QuickSCAT融合風場數(shù)據(jù),計算區(qū)域為37°N~41°N,117°E~123°E,模式每小時輸出全場波浪要素和破碎參數(shù)。
模擬區(qū)域內(nèi)有效波高Hs和破碎能量損耗率(基于EY計算得到)的空間分布見圖1,帶色標的底圖描繪的是有效波高分布,帶數(shù)值的等值線描繪的是破碎能量損耗率分布。
圖1 模擬得到的不同時段有效波高和破碎能量損耗率平面分布
從整個渤海海區(qū)的波高分布上來看,溢油時間段內(nèi)的日平均有效波高都在1 m以內(nèi),2006年3月27日、28日、31日在渤海海峽附近出現(xiàn)了幾次大浪過程。渤海灣曹妃甸海域的平均波高在0.3 m左右,在31日最大為0.5 m。從整個區(qū)域的破碎能量損耗率上來看,在形態(tài)上與有效波高基本相同,波高較高的區(qū)域,能量損耗也越大,波高小于0.1 m的時候,破碎能量很小。從溢油區(qū)域的波高分布來看,溢油時間段內(nèi)的日平均浪高在0.3 m左右,在30日、31日兩天有相對較大的波浪傳入渤海灣,最大波高達到0.5 m。從溢油區(qū)域的破碎能量損耗率上來看,25日到29日日平均破碎能量損耗率在0.1 kg/s3以內(nèi),30日和31日隨著波高的增大最大可達到0.24 kg/s3。
比較溢油點附近兩種計算方案所得到的破碎能量損耗率發(fā)現(xiàn)(見圖2),在有效波高較低的情況下,兩種計算方案結(jié)果差異不大。但是當有效波高超過1 m時,EY的值陡然上升,對應(yīng)結(jié)果超過EH&K近3倍。
圖2 溢油發(fā)生點不同計算方案得到的破碎能量損耗率隨時間變化情況
計算結(jié)果見圖3。
圖3 不同實驗中海面殘油量和卷入量比例隨時間變化情況
3個試驗在最初的十幾個小時內(nèi),海表殘油量都迅速減少,主要是蒸發(fā)過程造成的。方案1和方案2對破浪破碎不敏感,在輕組分蒸發(fā)完以后,殘油量不再變化。3種試驗方案在3月31號的殘油量分數(shù)分別為45%,41%、20%,見表2。
表2 3組數(shù)值試驗結(jié)果與衛(wèi)星觀測的比較 %
如前所述,歐空局Envisat的雷達數(shù)據(jù)顯示,這一時段的海表殘油量最后減少為原來的75%以上[7]。由此可見,方案3的殘油量分數(shù)最接近這個觀測值。
導(dǎo)致這一試驗結(jié)果的原因可能有以下幾點:
首先,方案1和方案2在形態(tài)上大致相同,是因為都采用了Delvigne等[4]的卷入深度公式,雖然能量損耗率所選方案不同,但是這一時段渤海灣的海況并不惡劣,兩種破碎方案得到的結(jié)果比較接近,日平均差大約在0.03 kg/s3左右。所以最終計算得到的結(jié)果差別較小。
其次,方案2和方案3相比較,形態(tài)差異較大,是因為采用不同的卷入深度方案。實際計算結(jié)果中顯示,破碎動力水深ZY的值比經(jīng)驗關(guān)系式ZD&S的值要小。根據(jù)卷入率Qo定義,相同的破碎能量,卷入深度越小,單位時間單位水體卷入的數(shù)量越多。所以方案3比方案1和2的卷入分數(shù)都要大。最終海表殘油量分數(shù)要比其他方案小得多,更接近實測值。
分析結(jié)果顯示,采用Yuan等[11]根據(jù)實際海況導(dǎo)出的破碎卷入深度和破碎能量損耗率表達式,可以更好地再現(xiàn)油滴在不同波浪破碎背景情況下下沉及上浮的整個過程。通過與Envisat SAR衛(wèi)星的觀測結(jié)果相對比,理論與數(shù)值模擬相結(jié)合的方法,可以很好地重現(xiàn)溢油入海后海面殘油量的長期變化趨勢。當考慮波浪效應(yīng)后,估計海面殘油量可由40%以上縮小至20%左右,與衛(wèi)星觀測結(jié)果中的約25%相比非常接近。
波浪對海面殘油量的影響絕非是可有可無的,尤其是在高海況情況下,該效應(yīng)對海面殘油量的估計準確性至關(guān)重要。前人根據(jù)經(jīng)驗給出的參數(shù)化波浪破碎要素計算方案雖然可以大致刻畫波浪破碎與油滴垂向運動的關(guān)系,但是其準確性和可靠程度低于基于動力學(xué)理論推導(dǎo)得到的各要素計算方案。
[1] 秦銳鋒,史文強.渤海灣海洋石油開采過程溢油風險預(yù)警管理[J].船海工程,2015,44(5):133-137.
[2] 李偉峰,肖中賢,史國友.未確知測度在海上溢油應(yīng)急能力評估中的應(yīng)用[J].船海工程,2016(增刊):253-256.
[3] LI, M, GARRETT, C. The relationship between oil droplet size and upper ocean turbulence[J].Marine Pollution Bulletin,1998,36:961-970.
[4] DELVIGNE G A L, SWEENEY C E. Natural dispersion of oil[J]. Oil and Chemical Pollution,1988(4):281-310.
[5] TKALICH P, CHAN E S. Vertical mixing of oil droplets by breaking waves[J]. Marine Pollution Bulletin,2002,44(11):1219-1229.
[6] 孫寶楠,袁業(yè)立,韓磊,等.波浪破碎下油滴粒徑分布理論模型研究[J].海洋學(xué)報,2014,36(9):30-36.
[7] 于五一,李進,邵蕓,等.海上油氣勘探開發(fā)中的溢油遙感監(jiān)測技術(shù):以渤海灣海域為例[J].石油勘探與開發(fā),2007,34(3):378-383.
[8] 吳曉丹,宋金明,李學(xué)剛,等.海上溢油量獲取的技術(shù)方法[J].海洋技術(shù)學(xué)報,2011,30(2):50-54.
[9] STIVER W, MACKAY D. Evaporation rate of spills of hydrocarbons and petroleum mixtures[J]. Environ Sci Tech,1984,18(11):834-840.
[10] BUCHANAN I, HURFORD N. Methods for predicting the physical changes in oil spilt at sea[J].Chem. Pollut,1988(4):311-328.
[11] YUAN Y, HAN L, HUA F, et al. The Statistical Theory of Breaking Entrainment Depth and Surface Whitecap Coverage of Real Sea Waves[J]. Journal of Physical Oceanography,2009,39(1):143-161.
[12] KOMEN G J, HASSELMANN S, HASSELMANN K. On the existence of a fully developed wind-sea spectrum[J]. Journal of Physical Oceanography,1984(14):1271-1285.
[13] 楊永增,喬方利,趙偉,等.球坐標系下MASNUM海浪數(shù)值模式的建立及其應(yīng)用[J].海洋學(xué)報,2005,27(2):1-7.