亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        波浪破碎過(guò)程對(duì)油滴垂直混合影響的數(shù)值模擬及實(shí)際驗(yàn)證

        2018-04-26 09:31:43,,
        船海工程 2018年2期
        關(guān)鍵詞:殘油損耗率波高

        ,,

        (1.國(guó)家海洋局第一海洋研究所,山東 青島 266061;2.海洋環(huán)境科學(xué)和數(shù)值模擬國(guó)家海洋局重點(diǎn)實(shí)驗(yàn)室,山東 青島 266061;3.青島海洋科學(xué)與技術(shù)國(guó)家實(shí)驗(yàn)室區(qū)域海洋動(dòng)力學(xué)和數(shù)值模擬功能實(shí)驗(yàn)室,山東 青島 266061)

        海上溢油是造成海洋環(huán)境污染損害的主要因素之一。在各種海洋污染中,石油污染無(wú)論在發(fā)生頻率、分布廣度還是在危害程度上均居首位。隨著世界海洋運(yùn)輸業(yè)的發(fā)展和海上油田不斷投入生產(chǎn),溢油事故發(fā)生頻率不斷升高,給海洋安全及人民健康帶來(lái)重大損害[1-2]。

        我國(guó)近海發(fā)生過(guò)多次嚴(yán)重的海上溢油事故,如2003年11月在渤海灣發(fā)生的“塔斯曼?!陛喴缬褪鹿剩鷳B(tài)索賠額高達(dá)1.2億元人民幣。又如2006年2月發(fā)生的“長(zhǎng)島海域油污染事件”,范圍波及山東、天津、河北沿岸,是有史以來(lái)渤海發(fā)生的影響范圍最大的油污染事件。在所有這些溢油事故中,對(duì)殘存油量的估計(jì)都是非常重要的一個(gè)問題。

        在多種作用因素中,波浪破碎導(dǎo)致的混合攪拌效應(yīng)通過(guò)影響油滴垂向混合,從而對(duì)溢油海面殘存油量估計(jì)起到了重要作用。觀測(cè)表明,高海況下被破碎波卷入水下的溢油量可以達(dá)到溢油總量的50%以上[3],所以研究波浪破碎對(duì)油滴的垂直混合是研究海上溢油海面殘油量隨時(shí)間、空間變化的重要基礎(chǔ),只有充分了解并掌握這一物理過(guò)程,才能夠全面地分析和預(yù)報(bào)海上溢油的變化趨勢(shì),可以對(duì)溢油分散劑的使用量有一定的指導(dǎo)意義。同時(shí),也可以為溢油事故發(fā)生后產(chǎn)生的經(jīng)濟(jì)賠償和生態(tài)損失評(píng)估提供強(qiáng)有力的理論基礎(chǔ)和數(shù)據(jù)證據(jù)。

        對(duì)波浪破碎與油滴垂直混合之間關(guān)系的研究,之前主要采用了觀測(cè)數(shù)據(jù)分析和理論研究?jī)煞N研究手段。Delvigneand等[4]、Tkalich等[5]和孫寶楠等[6]根據(jù)實(shí)驗(yàn)室數(shù)據(jù)并結(jié)合動(dòng)力推導(dǎo)給出了波浪破碎與油滴垂向運(yùn)動(dòng)之間的一些經(jīng)驗(yàn)和理論關(guān)系表達(dá)形式。溢油事故往往發(fā)生突然,并沒有非常實(shí)時(shí)的現(xiàn)場(chǎng)觀測(cè)波浪數(shù)據(jù)以供分析采用。在這種情況下,數(shù)值模式是一種非常有效的應(yīng)急補(bǔ)充方法,可以由計(jì)算得到的各類波浪破碎參數(shù),結(jié)合理論分析方法,給出具體時(shí)段和地點(diǎn)的海面溢油事故發(fā)生后油滴垂直混合情況,進(jìn)而最終得到海面殘油量變化趨勢(shì)。

        之前的分析研究因?yàn)楝F(xiàn)場(chǎng)觀測(cè)數(shù)據(jù)的缺乏,導(dǎo)致計(jì)算結(jié)果的可靠性存在問題。隨著衛(wèi)星觀測(cè)技術(shù)手段的豐富,雖然定量估計(jì)殘油量仍然難度較大,但是定性評(píng)估海面殘油量已經(jīng)成為可能[7-8]。這使得結(jié)合理論分析和數(shù)值模擬的方法,針對(duì)某一具體海上溢油時(shí)間開展波浪破碎過(guò)程對(duì)油滴垂直混合影響研究具有一定可行性。

        1 研究方法

        針對(duì)波浪破碎過(guò)程對(duì)油滴垂直混合影響開展研究,必須首先建立能夠描述包括油滴入水和上浮兩個(gè)垂向運(yùn)動(dòng)過(guò)程在內(nèi)的統(tǒng)一動(dòng)力方程。同時(shí),為了最終給出海面殘油量的估計(jì)結(jié)果,該動(dòng)力方程必須結(jié)合溢油蒸發(fā)模型,共同構(gòu)成描述海面殘油量變化的數(shù)學(xué)模型。

        1.1 油滴入水模型

        根據(jù)Tkalich等[5]的一階模型,把油滴垂向運(yùn)動(dòng)方程描述為

        (1)

        式中:Ms為海面浮油的質(zhì)量;Me為海面以下殘油的質(zhì)量(都是在單位海表面積內(nèi));Qo為單位時(shí)間單位水體破碎卷入的油滴體積(即卷入率)。

        入水油的總質(zhì)量為:Ma=Ms+Me。在溢油之初,Ma等于溢漏出來(lái)的油的初始質(zhì)量M0。Qo在卷入過(guò)程中起著決定性的作用。下面用相似性定理分析Qo的表達(dá)形式,其中參與量綱分析的變量有4組。

        ③Z,油滴卷入深度,量綱:L;

        基本量綱為長(zhǎng)度量綱L、時(shí)間量綱T和質(zhì)量量綱M,共3個(gè),根據(jù)關(guān)鍵變量形成其他變量的量綱一的量方程,求解得到量綱相容方程。

        (2)

        因?yàn)橹挥幸粋€(gè)量綱相容方程,即只得到一個(gè)量綱一的量組,因此,有函數(shù)形式如下。

        (3)

        式中:f(x)為待定的函數(shù)形式。

        選取待定函數(shù)為冪函數(shù)形式,即令f(x)=cxk,可以得到

        (4)

        選取最簡(jiǎn)單的形式,待定系數(shù)c=1,k=1。式中Z和E分別為與波浪破碎過(guò)程有關(guān)的卷入深度和破碎能量損耗率。最終可以得到作用在單位海表面的Qo表達(dá)形式為

        (5)

        取kb=0.5,g=9.8 m/s2,ρw=1 025 kg/m3,Low=1 m。

        1.2 油滴上浮模型

        油滴上浮的動(dòng)力過(guò)程可用一階方程描述。

        (6)

        取Lwo=20 m,rc=50 μm,rmax=500 μm,

        ρ=900 kg/m3,υ=10-5m2/s,p=2。

        1.3 溢油蒸發(fā)模型

        溢油蒸發(fā)方案采用Stiver等[9]提出的分析法模型。

        (7)

        1.4 海面殘油量變化模型

        綜合以上3個(gè)子模型,根據(jù)Tkalich等[5]的簡(jiǎn)化方法,海面殘油量變化模型控制方程可寫為

        (8)

        1.5 波浪破碎過(guò)程對(duì)油滴垂直運(yùn)動(dòng)過(guò)程的影響

        通過(guò)已建立的運(yùn)動(dòng)學(xué)方程發(fā)現(xiàn),波浪破碎主要通過(guò)影響Qo對(duì)油滴垂直運(yùn)動(dòng)產(chǎn)生影響,而該要素主要受與波浪破碎有關(guān)的卷入深度(Z)和破碎能量損耗率(E)這兩個(gè)變量共同控制。針對(duì)這兩個(gè)要素,有多種計(jì)算方案可供選擇。

        針對(duì)Z,Delvigne等[4]通過(guò)試驗(yàn)室擬合給出的經(jīng)驗(yàn)表達(dá)形式:ZD&S=αHb,其中α=1.5,Hb為破碎波高,取Hb=0.7Hs,Hs為有效波高。Yuan等[11]根據(jù)實(shí)際海況導(dǎo)出了理論表達(dá)形式。

        (9)

        (10)

        為了分析采用不同計(jì)算方案對(duì)結(jié)果的最終影響,設(shè)計(jì)3組對(duì)比實(shí)驗(yàn),見表1。

        表1 數(shù)值試驗(yàn)方案設(shè)計(jì)

        1.6 溢油事件

        2006年3月24日渤海灣灤河河口以南曹妃甸附近發(fā)生了一起海洋溢油事故。根據(jù)當(dāng)日獲取的歐洲空間局環(huán)境衛(wèi)星(Envisat)合成孔徑雷達(dá)圖像判斷,該海域形成了大面積的油膜污染,總面積達(dá)400 km2。由2006年4月1日獲取的Envisat圖像判斷,渤海灣發(fā)生的溢油現(xiàn)象已經(jīng)減輕,沒有發(fā)現(xiàn)形態(tài)完整、溢出時(shí)間較短的油膜。由于波浪破碎攜油入水,蒸發(fā)等的共同作用,油膜已經(jīng)分散,質(zhì)量明顯減少,溢油面積減少到100 km2左右,殘油量減少了初始質(zhì)量的75%以上[7]。

        1.7 波浪數(shù)值模式

        為了再現(xiàn)溢油時(shí)段的渤海海區(qū)波浪平面分布,應(yīng)用球坐標(biāo)系下MASNUM海浪數(shù)值模式[13]模擬渤海的波浪場(chǎng)。水平分辨率是(1/12)°×(1/12)°,溢油時(shí)段為6 d。風(fēng)場(chǎng)資料采用的是QuickSCAT融合風(fēng)場(chǎng)數(shù)據(jù),計(jì)算區(qū)域?yàn)?7°N~41°N,117°E~123°E,模式每小時(shí)輸出全場(chǎng)波浪要素和破碎參數(shù)。

        2 結(jié)果分析

        模擬區(qū)域內(nèi)有效波高Hs和破碎能量損耗率(基于EY計(jì)算得到)的空間分布見圖1,帶色標(biāo)的底圖描繪的是有效波高分布,帶數(shù)值的等值線描繪的是破碎能量損耗率分布。

        圖1 模擬得到的不同時(shí)段有效波高和破碎能量損耗率平面分布

        從整個(gè)渤海海區(qū)的波高分布上來(lái)看,溢油時(shí)間段內(nèi)的日平均有效波高都在1 m以內(nèi),2006年3月27日、28日、31日在渤海海峽附近出現(xiàn)了幾次大浪過(guò)程。渤海灣曹妃甸海域的平均波高在0.3 m左右,在31日最大為0.5 m。從整個(gè)區(qū)域的破碎能量損耗率上來(lái)看,在形態(tài)上與有效波高基本相同,波高較高的區(qū)域,能量損耗也越大,波高小于0.1 m的時(shí)候,破碎能量很小。從溢油區(qū)域的波高分布來(lái)看,溢油時(shí)間段內(nèi)的日平均浪高在0.3 m左右,在30日、31日兩天有相對(duì)較大的波浪傳入渤海灣,最大波高達(dá)到0.5 m。從溢油區(qū)域的破碎能量損耗率上來(lái)看,25日到29日日平均破碎能量損耗率在0.1 kg/s3以內(nèi),30日和31日隨著波高的增大最大可達(dá)到0.24 kg/s3。

        比較溢油點(diǎn)附近兩種計(jì)算方案所得到的破碎能量損耗率發(fā)現(xiàn)(見圖2),在有效波高較低的情況下,兩種計(jì)算方案結(jié)果差異不大。但是當(dāng)有效波高超過(guò)1 m時(shí),EY的值陡然上升,對(duì)應(yīng)結(jié)果超過(guò)EH&K近3倍。

        圖2 溢油發(fā)生點(diǎn)不同計(jì)算方案得到的破碎能量損耗率隨時(shí)間變化情況

        計(jì)算結(jié)果見圖3。

        圖3 不同實(shí)驗(yàn)中海面殘油量和卷入量比例隨時(shí)間變化情況

        3個(gè)試驗(yàn)在最初的十幾個(gè)小時(shí)內(nèi),海表殘油量都迅速減少,主要是蒸發(fā)過(guò)程造成的。方案1和方案2對(duì)破浪破碎不敏感,在輕組分蒸發(fā)完以后,殘油量不再變化。3種試驗(yàn)方案在3月31號(hào)的殘油量分?jǐn)?shù)分別為45%,41%、20%,見表2。

        表2 3組數(shù)值試驗(yàn)結(jié)果與衛(wèi)星觀測(cè)的比較 %

        如前所述,歐空局Envisat的雷達(dá)數(shù)據(jù)顯示,這一時(shí)段的海表殘油量最后減少為原來(lái)的75%以上[7]。由此可見,方案3的殘油量分?jǐn)?shù)最接近這個(gè)觀測(cè)值。

        導(dǎo)致這一試驗(yàn)結(jié)果的原因可能有以下幾點(diǎn):

        首先,方案1和方案2在形態(tài)上大致相同,是因?yàn)槎疾捎昧薉elvigne等[4]的卷入深度公式,雖然能量損耗率所選方案不同,但是這一時(shí)段渤海灣的海況并不惡劣,兩種破碎方案得到的結(jié)果比較接近,日平均差大約在0.03 kg/s3左右。所以最終計(jì)算得到的結(jié)果差別較小。

        其次,方案2和方案3相比較,形態(tài)差異較大,是因?yàn)椴捎貌煌木砣肷疃确桨?。?shí)際計(jì)算結(jié)果中顯示,破碎動(dòng)力水深ZY的值比經(jīng)驗(yàn)關(guān)系式ZD&S的值要小。根據(jù)卷入率Qo定義,相同的破碎能量,卷入深度越小,單位時(shí)間單位水體卷入的數(shù)量越多。所以方案3比方案1和2的卷入分?jǐn)?shù)都要大。最終海表殘油量分?jǐn)?shù)要比其他方案小得多,更接近實(shí)測(cè)值。

        3 結(jié)論

        分析結(jié)果顯示,采用Yuan等[11]根據(jù)實(shí)際海況導(dǎo)出的破碎卷入深度和破碎能量損耗率表達(dá)式,可以更好地再現(xiàn)油滴在不同波浪破碎背景情況下下沉及上浮的整個(gè)過(guò)程。通過(guò)與Envisat SAR衛(wèi)星的觀測(cè)結(jié)果相對(duì)比,理論與數(shù)值模擬相結(jié)合的方法,可以很好地重現(xiàn)溢油入海后海面殘油量的長(zhǎng)期變化趨勢(shì)。當(dāng)考慮波浪效應(yīng)后,估計(jì)海面殘油量可由40%以上縮小至20%左右,與衛(wèi)星觀測(cè)結(jié)果中的約25%相比非常接近。

        波浪對(duì)海面殘油量的影響絕非是可有可無(wú)的,尤其是在高海況情況下,該效應(yīng)對(duì)海面殘油量的估計(jì)準(zhǔn)確性至關(guān)重要。前人根據(jù)經(jīng)驗(yàn)給出的參數(shù)化波浪破碎要素計(jì)算方案雖然可以大致刻畫波浪破碎與油滴垂向運(yùn)動(dòng)的關(guān)系,但是其準(zhǔn)確性和可靠程度低于基于動(dòng)力學(xué)理論推導(dǎo)得到的各要素計(jì)算方案。

        [1] 秦銳鋒,史文強(qiáng).渤海灣海洋石油開采過(guò)程溢油風(fēng)險(xiǎn)預(yù)警管理[J].船海工程,2015,44(5):133-137.

        [2] 李偉峰,肖中賢,史國(guó)友.未確知測(cè)度在海上溢油應(yīng)急能力評(pí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é)報(bào),2014,36(9):30-36.

        [7] 于五一,李進(jìn),邵蕓,等.海上油氣勘探開發(fā)中的溢油遙感監(jiān)測(cè)技術(shù):以渤海灣海域?yàn)槔齕J].石油勘探與開發(fā),2007,34(3):378-383.

        [8] 吳曉丹,宋金明,李學(xué)剛,等.海上溢油量獲取的技術(shù)方法[J].海洋技術(shù)學(xué)報(bào),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] 楊永增,喬方利,趙偉,等.球坐標(biāo)系下MASNUM海浪數(shù)值模式的建立及其應(yīng)用[J].海洋學(xué)報(bào),2005,27(2):1-7.

        猜你喜歡
        殘油損耗率波高
        基于FHDI-GNWM 數(shù)據(jù)的全球超越概率波高宏觀分布特征分析
        超聲波輔助提取油茶籽油的工藝優(yōu)化研究
        基于漂流浮標(biāo)的南大洋衛(wèi)星高度計(jì)有效波高研究
        非平整港池的多向不規(guī)則波試驗(yàn)研究
        殘油艙水分蒸發(fā)處理方式在IOPP證書附錄3.2.3項(xiàng)簽注問題和記錄注意事項(xiàng)
        損耗率高達(dá)30%,保命就是保收益!這條70萬(wàn)噸的魚要如何破存活率困局?
        淺談鋼煤斗制作原材料損耗率的研究
        脫水撬脫水效果分析與評(píng)價(jià)
        白天氣態(tài)亞硝酸(HONO)未知源對(duì)我國(guó)東部沿海地區(qū)HONO和OH自由基收支的影響
        飽和秋色
        国产系列丝袜熟女精品视频| 18禁黄污吃奶免费看网站| 精品国产一区av天美传媒| 国产欧美日韩专区| 青青草手机在线观看视频在线观看| 日本熟妇美熟bbw| 欧美人妻精品一区二区三区| 欧美精品久久久久久三级| 中文字幕熟女激情50路| 日韩精品无码一区二区三区四区 | 日本少妇按摩高潮玩弄| 一区二区三区四区国产亚洲| 亚洲一区精品无码| 亚洲精品无码乱码成人| 国产在线拍偷自拍偷精品| 大香蕉视频在线青青草| 午夜爽爽爽男女免费观看影院| 亚洲va中文字幕无码| 国产在线视频h| 蜜桃国产精品视频网站| 国产亚av手机在线观看| 草草网站影院白丝内射| 日本一区二区亚洲三区| 国产一区亚洲二区三区极品| 人人色在线视频播放| 婷婷丁香91| 熟女人妻一区二区中文字幕 | 97日日碰曰曰摸日日澡| 欧美日韩亚洲精品瑜伽裤| 色偷偷亚洲av男人的天堂| 人妻少妇满足中文字幕| 无遮挡h肉动漫在线观看| 久久狠狠高潮亚洲精品暴力打| 一区二区三区成人av| 女人的精水喷出来视频| 少妇寂寞难耐被黑人中出| 久久久久久无码AV成人影院| 亚洲不卡一区二区视频| 久激情内射婷内射蜜桃| 91久久精品无码人妻系列 | 国产精品黑色丝袜在线播放|