黃軼康,李一平,邱 利,薛偲琦,章雙雙

(1.河海大學(xué)環(huán)境學(xué)院,江蘇 南京 210098; 2.南京信息工程大學(xué)濱江學(xué)院,江蘇 南京 210044)

基于EFDC模型的長江下游碼頭溢油風(fēng)險預(yù)測

黃軼康1,李一平1,邱 利1,薛偲琦1,章雙雙2

(1.河海大學(xué)環(huán)境學(xué)院,江蘇 南京 210098; 2.南京信息工程大學(xué)濱江學(xué)院,江蘇 南京 210044)

為了更有效地預(yù)測溢油事故,建立了基于EFDC(environmental fluid dynamlics code)模型的長江下游南京段碼頭溢油事故影響的預(yù)測模型。研究結(jié)果表明:構(gòu)建的二維水動力模型能夠較為準(zhǔn)確地反映研究區(qū)域風(fēng)場、流場、復(fù)雜地形條件綜合影響下長江溢油擴(kuò)展和遷移運動的整體規(guī)律,油膜在研究區(qū)域降解較慢,溢油發(fā)生在落急時刻和西風(fēng)風(fēng)向時,油膜更早向下游擴(kuò)散,對下游保護(hù)區(qū)的污染更嚴(yán)重。

溢油風(fēng)險;EFDC模型;碼頭溢油;風(fēng)險預(yù)測;長江下游;碼頭;油膜

進(jìn)入21世紀(jì)以來,經(jīng)濟(jì)的迅猛發(fā)展導(dǎo)致了我國的能源需求越來越大,在加快了石油運輸行業(yè)蓬勃發(fā)展的同時,也帶來了海洋和江河溢油事故的風(fēng)險。長江是我國最重要的水源地,其水質(zhì)安全關(guān)系到千家萬戶,建立長江溢油應(yīng)急反應(yīng)系統(tǒng)非常必要。溢油應(yīng)急反應(yīng)系統(tǒng)的關(guān)鍵是對溢油軌跡的準(zhǔn)確預(yù)測,既縮小環(huán)境污染的范圍,又能保護(hù)海洋和江河的生態(tài)環(huán)境。多目標(biāo)評估模型是目前溢油事故風(fēng)險評價的主要模型,其基本思路是:首先確定溢油風(fēng)險評價指標(biāo),然后通過各種評估模型分析溢油風(fēng)險。但是這種多目標(biāo)評估模型比較適合溢油事故發(fā)生后的油污影響評價,對溢油事故發(fā)生前的預(yù)測則起不到借鑒作用。筆者首次嘗試以EFDC(environmental fluid dynamics code)模型預(yù)測長江下游碼頭溢油事故,建立長江下游南京段溢油事故影響的預(yù)測模型?；贓FDC模型,能高效整合多種邊界數(shù)據(jù),模擬油膜的遷移軌跡,并最終在地理信息系統(tǒng)的圖層里以動畫形式實時展示溢油油膜到達(dá)下游環(huán)境保護(hù)區(qū)的時間、油膜的遷移時間以及油膜的濃度情況。

1 研究區(qū)域概況

長江下游南京段自然條件優(yōu)越,是長江三角洲地區(qū)的重要組成部分。本河段水流路徑綿延,洲灘眾多,除局部河寬較窄外,河段寬度一般在1 km以上(圖1)[1]。本河段潮汐為非正規(guī)半日淺海潮,每日兩次漲潮、兩次落潮,且有日潮不等現(xiàn)象。在徑流與河床邊界條件阻滯下,潮波變形明顯,漲落潮歷時不對稱,自下而上漲潮歷時短,落潮歷時長,潮差沿程遞減,落潮歷時沿程遞增,漲潮歷時沿程遞減[1]。

圖1 研究區(qū)域

2 研究方法

2.1 預(yù)測模型簡介

長江下游南京段碼頭溢油事故影響的預(yù)測模型由EFDC模型模擬而成。目前EFDC是美國環(huán)保署TMDL模擬工具箱中的重要模擬軟件。根據(jù)對南京段碼頭溢油事故水環(huán)境預(yù)測特點的要求,考慮長江該段的地理氣候條件,對比了現(xiàn)行眾多的地表水模擬軟件[2],筆者選擇EFDC模型。EFDC模型優(yōu)勢在于:①EFDC模型是美國環(huán)保署首推的模擬軟件,在水動力模擬[3-4]和水質(zhì)模擬[5]方面的準(zhǔn)確性已得到學(xué)術(shù)界的認(rèn)可;②與 POM、ECOM 模型相比,EFDC模型能夠更加有效地處理變邊界條件,文件輸入方式更通用,更好地耦合水動力、水質(zhì)和泥沙模塊,不同模型接口間的轉(zhuǎn)換過程可以省略;③EFDC-Explorer軟件擁有友好的可視化操作界面,能夠完整地完成模擬的前后處理,快速生成網(wǎng)格數(shù)據(jù),最終形成實時圖像；④EFDC是開源軟件,有利于預(yù)測模型的集成和二次開發(fā)。

2.2 區(qū)域網(wǎng)格的構(gòu)建

根據(jù)研究區(qū)域長江地形實際變化情況,對干流江段選擇采用矩形網(wǎng)格,對南京長江大橋斷面到三江口斷面江段則采用曲線正交貼體網(wǎng)格,以便更好地擬合岸界情況。笛卡爾網(wǎng)格單元的大小為100 m×100 m,曲線正交網(wǎng)格單元沿河寬方向最小1.4 km,最大2.9 km,沿河流方向全程37 km,在研究區(qū)域共概化出9 174個網(wǎng)格單元。由于長江在研究江段的水深較淺(水深遠(yuǎn)小于河寬),可忽略垂向變化,概化為二維平面網(wǎng)格。通過EFDC自身的網(wǎng)格生成程序,進(jìn)行設(shè)置邊界控制點,生成研究區(qū)域的計算網(wǎng)格。如圖2所示,長江下游南京段碼頭溢油事故預(yù)測模型都將建立在該空間概化的結(jié)果之上。

圖2 研究區(qū)域計算網(wǎng)格

2.3 源強(qiáng)分析

有關(guān)統(tǒng)計資料表明,物料碼頭可能發(fā)生的泄漏有3種類型:①船舶碰撞發(fā)生的事故性溢油;②裝卸作業(yè)時管道泄漏事故性溢油;③裝卸過程集裝箱掉落入江溢油。

參照一般工程實際情況,事故溢油主要為船舶自身的燃料油,一般1 000 t級以上船載儲油量可達(dá)到10t以上,一旦發(fā)生船舶相撞導(dǎo)致漏油現(xiàn)象,會造成比較大的溢油事故。按照一次最大溢油量,1 000 t級輪出現(xiàn)撞船等事故導(dǎo)致的溢油量為10 t。

按照國家標(biāo)準(zhǔn)集裝箱規(guī)格, 20尺(5.69 m×2.13 m×2.18 m)集裝箱和40尺(11.8 m×2.13 m×2.18 m)集裝箱的配貨毛質(zhì)量分別為27 t和26 t,因此,在考慮裝卸過程集裝箱掉落事故時,液體化工污染品100%入江,污染物最大入江量為27 t。但實際情況下,集裝箱內(nèi)存放化工品桶,單個化工品桶為200 kg,需考慮掉落化工品桶個數(shù),從而確定污染源強(qiáng)。

管道裝卸的污染源量需按照各工程案例實際情況進(jìn)行確定,不做統(tǒng)一要求。

2.4 模型選取

本研究使用的EFDC模型不同于常見的溢油模型,是采用Fay公式擴(kuò)展后的油粒子模型。該模型同時考慮了溢油的自身擴(kuò)展和漂移階段,且使自身擴(kuò)展和漂移階段同時進(jìn)行[6]。該模型的優(yōu)點在于吸取了油粒子模型和Fay模型各自的長處,成功突破Fay模型局限于靜止水面的限制,而特別適用于水域狹窄且河道復(fù)雜、水流狀況較為多變的河流溢油模擬。

油粒子模型的中心思想是將溢油離散化為大量油粒子,每個粒子表征一定的油量[7]。在流場中追蹤油粒子的運動軌跡,得到每一時刻各個油粒子所處的空間位置,統(tǒng)計各時刻油粒子的位置,可得到各時刻溢油的空間分布[8]。目前,針對事故性溢油的預(yù)測,國內(nèi)外都取得了一系列研究成果。沈永明等[9]建立了統(tǒng)籌物理、化學(xué)和生物反應(yīng)綜合作用的近岸海域多組分三維水質(zhì)動力學(xué)模型;熊德琪等[10]將風(fēng)場和流場數(shù)據(jù)庫作為強(qiáng)制條件,建立了大連港溢油應(yīng)急預(yù)報系統(tǒng),并借助 HAMSOM 三維海洋模式進(jìn)行潮流分析,完成珠江口水域的溢油模擬工作;汪守東等[10-11]考慮了風(fēng)、浪、流等多要素的誘導(dǎo)作用,模擬溢油油膜的遷移過程,選用POM模式作為水動力模型,實現(xiàn)溢油遷移和歸宿模型的三維化;Mannul等[12]選擇拉格朗日油粒子模型(LPTM),將采用ROMS(regional ocean modeling system)模式模擬出的流場以及美國國家環(huán)境預(yù)報中心NCEP(National Center for Environmental Prediction)每6 h一次的風(fēng)場作為邊界數(shù)據(jù), 模擬了“威望號”船舶溢油事故?？傮w而言,現(xiàn)階段有關(guān)溢油的研究主要集中于開闊的海域,而對長江等淺窄江道領(lǐng)域,鮮見油粒子模型在此相關(guān)領(lǐng)域的研究。對溢油的漂移擴(kuò)散過程,研究主要在油粒子模型方面展開;而對油膜擴(kuò)展過程,研究更多地集中在 Fay提出的三階段擴(kuò)展理論[13]。Fay 模型的假設(shè)是,在自身擴(kuò)展階段油膜在靜止水面擴(kuò)展,不會碰到水陸邊界,所以開敞海域里Fay假設(shè)是基本適用的,但河流則相對狹窄,實際情況下溢油在自身擴(kuò)展階段多數(shù)是會接觸水陸邊界的,因此將Fay假設(shè)用于江、河溢油預(yù)測會出現(xiàn)比較多的問題[14],故采用Fay公式擴(kuò)展后的油粒子模型是非常有必要的。

3 模型的建立

3.1 模型的基本原理

EFDC模型垂向上采用σ坐標(biāo)變換,能較好地擬合近岸復(fù)雜的岸線和地形;采用Gelperin等修正的Mellor-Yamada2.5階湍封閉模式能較客觀地提供垂向混合系數(shù),避免人為選取造成誤差。動量方程、連續(xù)方程及狀態(tài)方程為

?t(mxmyHu)+?x(myHuu)+?y(mxHvu)+

?z(mxmywu)-femxmyHv=-myH?x(p+

(1)

(2)

mxmyfe=mxmyf-u?ymx+v?xmy

(3)

(τxz,τyz)=AvH-1?z(u,v)

(4)

?t(mxmyH)+?x(myHu)+?y(mxHv)+

?z(mxmyw)=QH+δ(0)(Qss+Qsw)

(5)

?zp=-gH(ρ-ρ0)ρ0-1=-gHb

(6)

ρ=ρ(p,S,T)

(7)

垂向紊動黏滯系數(shù)Av通過Mellor and Yamada(1982)提出的二階矩紊動閉合模型求得:

Av=φAA0ql

(8)

(9)

(10)

(11)

(12)

(13)

式中:φA為垂向穩(wěn)定分層環(huán)境下的傳輸衰減系數(shù);q為比流量,L/(s·m);l為長度,m;A0為初始垂向紊動黏滯系數(shù);A1、B1、A2、B2、C1為與主控的紊動特征長度和紊動結(jié)構(gòu)特征有關(guān)的實驗常數(shù),分別取0.92、16.6、0.74、10.1和0.08。

動量方程的垂向邊界層考慮了水表面的風(fēng)拖曳力和湖底的摩擦力,即河床剪切應(yīng)力τxz和τyz取決于速度分量,可以根據(jù)二次阻力公式求得,具體方程為

(14)

式中:Uw、Vw分別為x、y方向在水表面10 m高處的風(fēng)速；τsx和τsy分為x、y方向水表面的風(fēng)拖曳力分量；cs為風(fēng)拖曳系數(shù)。cs求法如下:

(15)

式中,ρa(bǔ)和ρw分別為空氣和水的密度。

按：歐公《集古跋尾》云：“《盤谷序》石本貞元中所刻，以集本校之，或小不同，疑刻石誤。然以其當(dāng)時之物，姑存之以為佳玩，其小失不足校也?！痹敼搜裕顬橥ㄕ?。近世論者專以石本為正，如《水門記》《溪堂詩》，予已論之，《南海廟》《劉統(tǒng)軍碑》之類亦然，其謬可考而知也。[注] 《原本韓集考異》卷四《汴州東西水門記》：“‘文’，方從石、閣、蜀本作‘醇’。今按：此記方氏多從石本，石本固當(dāng)據(jù)信，但上條用字大誤，而此‘醇’字亦未安耳?！贝颂幩^“石本固當(dāng)據(jù)信”，乃就通例而言，而就此文此句而言，則石本不可據(jù)信，只不過朱子表達(dá)比較委婉而已。

湖底摩擦力的計算方法為

(16)

式中:u,v的下標(biāo)1指河流底部的對應(yīng)流速；τbx,τby分別為x、y方向湖底部切應(yīng)力分量；cb為底摩擦系數(shù)；cb的求取方程為

(17)

方程組在離散求解過程中,空間上采用二階精度的中心差分格式,時間上采用三次時間步蛙跳式差分格式。使用質(zhì)量守恒格式的干濕網(wǎng)格處理方法,更好地模擬湖流的漫灘過程,保證淺水區(qū)域計算結(jié)果的準(zhǔn)確性。具體離散求解方法可參見文獻(xiàn)[15]。

3.2 模型的設(shè)置與率定

數(shù)值模型的參數(shù)取值一直是影響模型結(jié)果的關(guān)鍵,因此必須對模型進(jìn)行靈敏性分析。由于對水流運動物理特性的研究已經(jīng)十分成熟,因此EFDC 模型的大部分物理參數(shù)未作改動,如,Mellor-Yamada 紊流模型[16]的相關(guān)參數(shù)與其他水動力模型(如Princeton海洋模型[17]與河口海岸海洋模型[18])中設(shè)置的參數(shù)是一致的。此外,考慮燃油的理化性質(zhì),模型中燃油的密度取0.8 kg/m3,衰減系數(shù)取0.2。類似地,Smagorinsky 公式[19]中的量綱黏度系數(shù)是不變的,取 0.2[20]。因此,本文確定使用 3 s 為模型的時間步長,以保證模型運行的穩(wěn)定性。

目前使用1997年7月19日 0:00～7月26日11:00(共180h)、1996年7月1日1:00～7月9日2:00(共194h)的實測水文資料。模型驗證結(jié)果顯示,沿程水面線與實測水面線基本一致,南京站檢驗水位誤差不超過20 cm的時段占總驗證時段(1996年為194 h,1997年為180 h)的70%～90%之間(圖3)?？傮w而言,利用一維水流模型對區(qū)域的概化處理基本合理,選用的參數(shù)基本反映了河道的水力特性。

圖3 南京站計算斷面的實測水位與計算結(jié)果的比較

由率定驗證可知,EFDC具有較好的匹配性,能夠較真實地模擬出長江的水動力變化情況。

本研究以南京港西壩港區(qū)西壩作業(yè)區(qū)八期液體化工碼頭工程為實例。南京港西壩港區(qū)為5萬 t級液體化工品碼頭,3萬t級泊位1個,吞吐量140萬t/a。預(yù)測該碼頭發(fā)生的事故是船舶碰撞溢油10 t。用笛卡爾直角坐標(biāo)網(wǎng)格建立長江下游南京段模型,水平面總計 9 174 個網(wǎng)格,每個網(wǎng)格單元邊長相當(dāng)于實際長度 100 m。EFDC 模型以大氣、表面風(fēng)力和出入江流量為動力邊界條件。初始條件設(shè)置了潮位、流量、風(fēng)向風(fēng)速。通過EFDC模型在項目溢油事故發(fā)生地點連續(xù)釋放示蹤劑(LPT)100個,追蹤油膜漂移的軌跡,模擬油粒子過程。時間步長取為3 s,共計算8 d。

為了方案的合理性,考慮影響最大的情況,本研究使用枯水期小潮位數(shù)據(jù)輸入邊界條件作為入江口和出江口的潮位數(shù)據(jù),設(shè)置項目所在地位置、溢油的流量、持續(xù)溢油的時間、風(fēng)向和風(fēng)速。從偏安全角度考慮,根據(jù)氣象氣候資料,選取年主導(dǎo)風(fēng)向(東北風(fēng))、年平均風(fēng)速3.4 m/s;不利風(fēng)向(西風(fēng))、強(qiáng)風(fēng)速10 m/s兩種情況下事故分別發(fā)生在剛好漲潮和剛好落潮的潮位過程作為設(shè)計工況條件(表1)。同時考慮到兩種類型事故下的溢油量相近,影響效果也會相似,因此只對10t溢油量進(jìn)行預(yù)測,分析事故發(fā)生后對上下游水環(huán)境保護(hù)目標(biāo)的影響。

表1 模型預(yù)測方案

注：各方案溢油量均為10 t。

4 結(jié)果與討論

假設(shè)分別在剛好漲潮時(漲急)和剛好落潮時(落急)發(fā)生事故泄漏,模型模擬了此江段沿程的各個保護(hù)區(qū)受到溢油事故影響的時間、溢油到達(dá)和離開各保護(hù)區(qū)的時間,見表2。

研究區(qū)域位于長江下游感潮河段,水流既受上游下泄徑流的影響,又受下游潮汐的影響,水流極其復(fù)雜。從圖4的結(jié)果可以看出,落急時長江下游南京段的計算區(qū)域內(nèi)流場分布不均勻,東北風(fēng)時,流速在0.01～0.65 m/s之間;西風(fēng)強(qiáng)風(fēng)速10 m/s時,計算區(qū)域內(nèi)流速在0.01～0.72 m/s之間。漲急時計算區(qū)域內(nèi)流場分布也不均勻,東北風(fēng)時,流速在0.01～0.15 m/s之間;西風(fēng)強(qiáng)風(fēng)速10 m/s時,計算區(qū)域內(nèi)流速在0.01～0.25 m/s之間。從落潮向漲潮變化的過程不會出現(xiàn)往復(fù)流,流速逐漸變小。

表2 溢油到達(dá)和離開保護(hù)區(qū)時刻

注：深色為長江石油類污染物本底質(zhì)量濃度(0.05 mg/L),符合國際Ⅱ類標(biāo)準(zhǔn)；淺色為高濃度狀態(tài)。圖4 研究區(qū)域流場分布

4.1 溢油泄漏時刻對油膜擴(kuò)散的影響

如圖5所示,方案1和方案3的結(jié)果表明,落急時刻油膜到達(dá)和離開各保護(hù)區(qū)的時間都比漲急的提前0.04 d。主導(dǎo)風(fēng)向溢油油膜到達(dá)龍?zhí)讹嬘盟吹氐臅r間一樣,而離開該保護(hù)區(qū)的時間則落急時刻比漲急時刻提前了0.04 d;溢油油膜到達(dá)興隆洲濕地的時間落急時刻比漲急時刻提前了0.08 d,離開該保護(hù)區(qū)的時間落急時刻提前了0.17 d。落急時刻,整個研究區(qū)域的流速較快,油膜往下游的擴(kuò)散更快,下游保護(hù)區(qū)更早受到污染。

圖5 相同風(fēng)場下不同潮位油膜漂移時間對比

4.2 風(fēng)場對油膜擴(kuò)散的影響

以往研究表明,溢油長期擴(kuò)展的主要影響因素是油膜和水表面的風(fēng)力,油膜運動速度大約為水表面10 m以上風(fēng)速的3%[21],所以風(fēng)場對溢油油膜的影響非常重要。本研究選用東北風(fēng)3.4 m/s、西風(fēng)10 m/s以及西風(fēng)3.4 m/s 3種方案,分析風(fēng)速和風(fēng)向?qū)τ湍U(kuò)散的影響(圖6～8)。方案3和方案5的預(yù)測結(jié)果表明:相同風(fēng)速、不同風(fēng)向下,油膜到達(dá)和離開保護(hù)區(qū)的時間明顯不同。對于龍?zhí)讹嬘盟吹?西風(fēng)下油膜0.19d到達(dá)、0.48 d離開,分別比東北風(fēng)下提前0 d和0.04 d;對于興隆洲濕地,西風(fēng)下油膜0.44 d到達(dá)、1.48 d離開,分別比東北風(fēng)下提前0.08 d和0.29 d。

圖6 方案3(3.4 m/s東北風(fēng))結(jié)果

圖7 方案4(10 m/s西風(fēng))結(jié)果

而通過對比方案4和方案5可以看出,相同風(fēng)向、不同風(fēng)速下,對于龍?zhí)讹嬘盟吹?10 m/s和3.4 m/s風(fēng)速方案油膜到達(dá)和離開的時間一致;對于興隆洲濕地,10 m/s風(fēng)速下油膜0.40 d到達(dá)、1.15 d離開,分別比3.4 m/s風(fēng)速下提前0.04 d和0.33 d。龍?zhí)讹嬘盟吹氐挠湍U(kuò)散受到風(fēng)向的影響更加明顯,西風(fēng)風(fēng)向是對下游保護(hù)區(qū)不利的條件。而風(fēng)速的變化影響油膜到達(dá)和離開興隆洲濕地的時間變化更明顯。

圖8 方案5(3.4 m/s西風(fēng))結(jié)果

4.3 油膜停留時間的討論

從方案1～5的結(jié)果可以看出,溢油油膜停留在龍?zhí)端吹氐臅r間最大為0.33 d(東北風(fēng)向漲急時刻),停留在興隆洲濕地的時間最大為1.25 d(東北風(fēng)向漲急時刻)。不難看出在最不利的兩種情況綜合影響下,方案3導(dǎo)致油膜對保護(hù)區(qū)的污染最大。長江下游南京段該區(qū)域的流向是自西向東北方向流動的,因此油膜更傾向于聚集在河道的北岸,即興隆洲濕地一側(cè)。同時,從底高程數(shù)據(jù)來看,興隆洲濕地區(qū)域高程較高,出現(xiàn)了淺灘,流速較慢,停留時間也會較長。

5 結(jié) 論

a. 以枯水期水文觀測資料對預(yù)測模型中水動力模型的驗證結(jié)果表明,基于EFDC構(gòu)建的二維水動力模型能夠較為準(zhǔn)確地反映研究區(qū)域的水動力情況。

b. 通過對南京港西壩港區(qū)西壩作業(yè)區(qū)八期液體化工碼頭工程情景的模擬可知,油膜在研究江段的擴(kuò)散過程主要受到潮位、風(fēng)速、風(fēng)向的影響。燃油泄漏發(fā)生在落急時刻時,油膜向下游擴(kuò)散得較快,較早到達(dá)保護(hù)區(qū)。但風(fēng)場對油膜的影響更顯著,風(fēng)速變化更直接地表現(xiàn)為溢油油膜對長江北岸保護(hù)區(qū)的影響,西風(fēng)風(fēng)向?qū)е掠湍じ绲竭_(dá)長江南岸保護(hù)區(qū)。

c. EFDC預(yù)測模型能對長江下游南京段碼頭溢油事故后污染物的水環(huán)境行為進(jìn)行快速的模擬和預(yù)測,并且能在地理信息系統(tǒng)的圖層上以動畫的形式實時展示事故影響的范圍和程度,可以為長江下游南京段碼頭溢油事故的預(yù)警應(yīng)急提供決策依據(jù)。

[ 1 ] 聞云呈,夏云峰,張世釗,等.長江下游南京至瀏河口河段沿程設(shè)計最低通航水位分析[J].水運工程,2012(3):100-103.(WEN Yuncheng,XIA Yunfeng,ZHANG shizhao,et al.Design lowest navigation water level of river reach from Nanjing to Liu River Estuary on lower Yangtze River[J].Port & Waterway Engineering,2012(3):100-103.(in Chinese))

[ 2 ] 錢瑜.環(huán)境影響評價[M].南京:南京大學(xué)出版社,2009:115-119.

[ 3 ] 齊珺,楊志峰,熊明,等.長江水系武漢段水動力過程三維數(shù)值模擬[J].水動力學(xué)研究與進(jìn)展,2008,23(2):212-219.(QI jun,YANG Zhifeng,XIONG Ming,et al.Three-dimensional modeling of hydrodynamic processes in the Wuhan catchments of the Changjiang River[J].Journal of hydrodynamics,2008,23(2):212-219.(in Chinese))

[ 4 ] 姜恒志,沈永明,汪守東.甌江口三維潮流和鹽度數(shù)值模擬研究[J].水動力學(xué)研究與進(jìn)展,2009,24(1):63-70.(JIANG Hengzhi,SHEN Yongming,WANG Shoudong.Numerical simulation study on three-dimensional tidal flow and salinity in the Oujiang Estuary[J].Journal of Hydrodynamics,2009,24(1):63-70.(in Chinese))

[ 5 ] 苗曉雨.基于EFDC的尹府水庫水質(zhì)數(shù)值模擬及預(yù)測[D].青島:中國海洋大學(xué),2012

[ 6 ] 趙琰鑫,王永貴,張萬順,等.河道溢油污染事故二維數(shù)值模型研究 [J].人民長江,2012,43(15):81-84.(ZHAO Yanxin,WANG Yonggui,ZHANG Wanshun,et al.River oil spill pollution in accident two-dimensional numerical model for research[J].Yangtze River,2012,43(15):81-84.(in Chinese))

[ 7 ] 蔡忠清.海上船舶溢油粒子化及適應(yīng)性數(shù)值模擬研究[D].武漢:武漢理工大學(xué),2010.

[ 8 ] 陳麗娜.可降解地表水污染物隨機(jī)游走改進(jìn)模式研究[D].南京:河海大學(xué),2008.

[ 9 ] 沈永明,鄭永紅,吳修廣.近岸海域污染物遷移轉(zhuǎn)化的三維水質(zhì)水動力學(xué)模型[J].自然科學(xué)進(jìn)展,2000,14(6):664-670.(SHENG Yongming,ZHENG Yonghong,WU Xiuguang.Migration of pollutants in near shore three-dimensional dynamic model of water quality of water[J].Progress in Natural Science,2000,14(6):664-670.(in Chinese))

[10] 熊德琪,楊建立,嚴(yán)世強(qiáng).珠江口區(qū)域海上溢油應(yīng)急預(yù)報信息系統(tǒng)的開發(fā)研究[J].海洋環(huán)境科學(xué),2005,24(5):63-66.(XIONG Deqi,YANG Jianli,YAN Shiqiang.R & D on oil spill contingency forecast information system for Pearl River Estuary[J].Marine Environmenal Science,2005,24(5):63-66.(in Chinese))

[11] 汪守東,沈永明,鄭永紅.海上溢油遷移轉(zhuǎn)化的雙層數(shù)學(xué)模型[J].力學(xué)學(xué)報,2006,38(4):452-461.(WANG Shoudong,SHEN Yongming,ZHENG Yonghong.The migration of oil spill into the double mathematical model[J].Chinese Journal of Theoretical and Applied Mechanics,2006,38(4):452-461.(in Chinese))

[12] 汪守東.基于 Langrange追蹤的海上溢油預(yù)報模型研究[D].大連:大連理工大學(xué),2008.

[13] MANUEL G,LUIS F,ADOLFO U,et al.Operational ocean graphy system applied to the prestige oil-spillage event[J].Journal of Marine System,2008,72(1/2/3/4):178-188.

[14] FAY J A.The spread of oil slicks on a calm sea[M].New York:Oil on the Sea Plenum,1969.

[15] 張艷軍,雒文生,雷阿林,等.河道溢油模型在三峽水質(zhì)預(yù)警系統(tǒng)中的研究與運用[J].環(huán)境科學(xué)與技術(shù),2010,33(1):128-131.(ZHANG Yanjun,LUO Wensheng,LEI Alin,et al.Research and application of river oil spill model in advanced warning system of water quality of Three Gorges [J].Environmental Science and Technology,2010,33(1):128-131.(in Chinese))

[16]TETRA TEC,INC.The environmental fluid dynamics code theory and computation [M].Washington DC:US EPA,2007:233-235.

[17] MELLOR G L,YAMADA T.Development of a turbulence closure-model for geophysical fluid problems[J].Reviews of Geophysics,1982,20(4):851-875.

[18] YIN Xunqiang,QIAO Fangli,TANG Yongzeng,et al.An ensemble adjustment Kalman filter study for Argo data[J].Chinese Journal of Oceanology and Limnoligy,2010,28(3):626-635.

[19] 堵盤軍,胡克林,孔亞珍,等.ECOMSED 模式在杭州灣海域流場模擬中的應(yīng)用[J].海洋學(xué)報,2007,29(1):7-16.(DU Panjun,HU Kelin,KONG Yazhen,et al.Application of ECOMSED model to the simulation of Hangzhou Bay tide current [J].Acta Oceanologica Sinica,2007,29(1):7-16.(in Chinese))

[20] BERNTSEN J.Internal pressure errors in sigma-coordinate ocean models [J].Journal of Atmospheric and Oceanic Technology,2002,19(9):1403-1414.

[21] 吳曉丹,宋金明,李學(xué)剛,等.海洋溢油油膜厚度影響因素理論模型的構(gòu)建[J].海洋科學(xué),2010,34(2):68-74.(WU Xiaodan,SONG Jinming,LI Xuegang,et al.Construction of theoretical model for thickness of marine oil spilled [J].Marine Sciences,2010,34(2):68-74.(in Chinese))

Risk prediction on wharf oil spill in the lower reaches of Yangtze River based on EFDC

HUANG Yikang1, LI Yiping1,QIU Li1, XUE Siqi1,ZHANG Shuangshuang2

(1.CollegeofEnvironment,HohaiUniversity,Nanjing210098,China2.BinjiangCollege,NanjingUniversityofInformationScience&Technology,Nanjing210044,China)

In order to predict more effectively the oil spill accident, a model to predict the influence of wharf oil spill accident of Nanjing section in the lower reaches of Yangtze River was set up based on EFDC (environmental fluid dynamics code). The study results show that the two dimensional hydrodynamic model set up can more accurately reflect the overall rules of the spilled oil spreading and migration in Yangtze River under the comprehensive effects of wind field, flow field, complex terrain conditions of the study area, and the degradation of oil slicks in the study area is slow, and when the oil spill occurred in the ebb time and in the westerly wind, oil slicks diffuse to the downstream more early, resulting in a more serious pollution on the downstream protection zone.

oil spill risk; EFDC model; wharf oil spill; risk prediction; lower reaches of Yangtze River;wharf; oil slick

10.3880/j.issn.1004-6933.2015.01.015

黃軼康(1990—),男,碩士研究生,研究方向為水環(huán)境數(shù)學(xué)模型研究和環(huán)境影響評價。E-mail:kolowin@163.com

X507

A

1004-6933(2015)01-0091-08

2014-04-11 編輯：彭桃英)