EFDC模型在水動(dòng)力環(huán)境影響評(píng)價(jià)中的應(yīng)用

福建省環(huán)境科學(xué)研究院 郭素銘

?

EFDC模型在水動(dòng)力環(huán)境影響評(píng)價(jià)中的應(yīng)用

福建省環(huán)境科學(xué)研究院 郭素銘

運(yùn)用EFDC模型在擬建的福建東山華浮碼頭進(jìn)行模型構(gòu)建，通過(guò)東山大澳漁港的3個(gè)站位實(shí)測(cè)潮位潮流數(shù)據(jù)驗(yàn)證模型模擬結(jié)果，表明模擬結(jié)果可信；進(jìn)一步模擬福建東山華浮碼頭建設(shè)前后，碼頭所在海域及周邊海域的潮流流向和流速的變化情況，模擬結(jié)果顯示，福建東山華浮碼頭的建設(shè)對(duì)所在海域水動(dòng)力環(huán)境的影響較小，也體現(xiàn)了該模型對(duì)東山灣海域的適用性。

福建東山華浮碼頭 EFDC 水動(dòng)力 模擬

福建東山灣海洋環(huán)境自然資源非常豐富，地理位置優(yōu)越，是天然的旅游勝地。但是，隨著經(jīng)濟(jì)不斷發(fā)展，越來(lái)越多的碼頭出現(xiàn)在東山灣海岸，改變了東山灣海岸的岸線分布情況，對(duì)航道的正常通行產(chǎn)生了影響，也改變了局部的沖淤環(huán)境，造成局部水域泥沙淤積嚴(yán)重，失去原使用功能。因此，通過(guò)模型對(duì)碼頭建設(shè)前后水動(dòng)力環(huán)境的模擬，可以較好地掌握工程建設(shè)前后水動(dòng)力環(huán)境可能產(chǎn)生的變化，從而采取相應(yīng)的措施，從海洋環(huán)境保護(hù)角度對(duì)工程可行性做出明確結(jié)論，為管理部門決策、建設(shè)單位海洋環(huán)境管理提供科學(xué)依據(jù)。

1 研究區(qū)域概況

東山縣是福建省第二大島，位于廈、漳、泉閩南三角經(jīng)濟(jì)區(qū)的南端，東瀕臺(tái)灣海峽，西臨詔安灣與詔安一水之隔；康美鎮(zhèn)地處東山縣東北部，東接銅陵鎮(zhèn)，西連樟塘鎮(zhèn)，地理位置優(yōu)越。

擬建東山華浮碼頭位于東山縣康美鎮(zhèn)城垵村東北側(cè)，地處東山港區(qū)城垵作業(yè)區(qū)，其地理坐標(biāo)為東經(jīng)117°30′、北緯23°44′。東北向與廈門經(jīng)濟(jì)特區(qū)毗鄰，南與廣東省汕頭市接壤，東瀕臺(tái)灣海峽，與臺(tái)灣省隔海相望；水路離廈門77海里、距汕頭73海里、距廣州332海里，陸路距東山縣城約10km，距漳州市約160km，規(guī)劃的廈深鐵路東山鐵路支線的終點(diǎn)站緊鄰港區(qū)，水陸交通十分方便。具體地理位置見圖1。

2 模型簡(jiǎn)介

環(huán)境流體動(dòng)力學(xué)模型，簡(jiǎn)稱EFDC模型(Environmental Fluid Dynamics Computer Code)是由美國(guó)Virginia海洋研究所的Hamrick等根據(jù)多個(gè)數(shù)學(xué)模型集成開發(fā)研制的綜合模型，現(xiàn)在是美國(guó)環(huán)保署(EPA)推薦使用的模型。該模型是一個(gè)多任務(wù)、高集成的環(huán)境流體動(dòng)力學(xué)模塊式計(jì)算程序包，用于模擬水系統(tǒng)一維、二維和三維流場(chǎng)、物質(zhì)輸送（包括溫、鹽、非粘性和粘性泥沙的輸送）、生態(tài)過(guò)程及淡水入流。其模擬范圍為：河口、河流、湖泊、水庫(kù)、濕地以及自近岸到陸架的海域。可以同時(shí)考慮風(fēng)、浪、潮、徑流的影響，并可同步布設(shè)水工建筑物。該模型到目前為止已經(jīng)用于幾十個(gè)海域的相關(guān)計(jì)算，得到廣泛的應(yīng)用[1-10]，被譽(yù)為21世紀(jì)最有發(fā)展前途的環(huán)境流體動(dòng)力學(xué)模型。采用該數(shù)學(xué)模型對(duì)本工程海域潮流場(chǎng)進(jìn)行模擬計(jì)算，計(jì)算中采用水平方向上的變笛卡爾正交坐標(biāo)與垂直方向上的Sigma坐標(biāo)相結(jié)合以及三維數(shù)學(xué)模型二維化的方法。

圖1 擬建東山華浮碼頭地理位置圖

動(dòng)力學(xué)方程是基于三維水動(dòng)力學(xué)方程組，在水平方向上采用曲線正交坐標(biāo)變化和在垂直方向上采用Sigma坐標(biāo)變換得到的，經(jīng)過(guò)兩種變換后的流體動(dòng)力學(xué)方程組分別為：

式(1)～(6)中u和v分別為坐標(biāo)x和y方向上的水平速度分量；mx和my為水平坐標(biāo)變換因子；經(jīng)坐標(biāo)變換后垂直方向z方向的速度w與坐標(biāo)變換前的垂直速度w*間的關(guān)系為：

H = h＋δ為總水深；p為壓力；動(dòng)量方程(1)和(2)中，f為Coriolis系數(shù)，Av為垂直紊動(dòng)粘性系數(shù)；Qu和Qv為動(dòng)量源匯項(xiàng)；QS和QT為溫鹽源匯項(xiàng)；ρ為海水密度；S為海水鹽度；T為海水溫度；b為浮力；連續(xù)方程(4)是在區(qū)間(0，1)對(duì)z積分并用垂直邊界條件當(dāng)z = (0，1)時(shí)，w = 0，運(yùn)動(dòng)邊界條件和方程(7)以得到深度積分連續(xù)方程(5)。

給出垂向紊動(dòng)和擴(kuò)散系數(shù)，方程(1)～(8)則給出了一個(gè)求解變量u，v，w，p，S，T和ζ的封閉的系統(tǒng)。紊動(dòng)粘性和擴(kuò)散系數(shù)采用的是Mellor和Yamada(1982)模型，模型相關(guān)的參數(shù)由下式確定：

以上各式中，q為紊動(dòng)強(qiáng)度，l為紊動(dòng)長(zhǎng)度，Rq為Richardson數(shù)，v和b是穩(wěn)定函數(shù)，以分別確定穩(wěn)定和非穩(wěn)定垂向密度分層環(huán)境的垂直混合或輸送的增減。

紊動(dòng)強(qiáng)度和混合長(zhǎng)度由下列方程確定：

式中B1、E1、E2和E3均為經(jīng)驗(yàn)常數(shù)；Qq和Ql為附加源匯項(xiàng)；例如子網(wǎng)格水平擴(kuò)散；垂直耗散系數(shù)Aq一般取與垂直紊動(dòng)粘性系數(shù)Av相等；上述式中m = mxmy。

式(1)～(8)與Mellor和Yamada(1982)紊動(dòng)模型(10)～(13)一起及適當(dāng)?shù)某踹呏禇l件給出了一個(gè)求解u，v，w，p，S，T，ρ和ζ的封閉的系統(tǒng)。

動(dòng)力學(xué)方程采用有限體積法和有限差分結(jié)合的方法來(lái)求解，水平方向采用交錯(cuò)網(wǎng)格離散。數(shù)值解分為沿水深積分長(zhǎng)波重力波的外模式和與垂直流結(jié)構(gòu)相聯(lián)系的內(nèi)模式求解。

3 模型構(gòu)建

3.1 模型模擬邊界

計(jì)算區(qū)域包括整個(gè)東山灣，見圖2，海岸線主要采用岸線修測(cè)成果，并結(jié)合歷史海圖和遙感圖確定；采用變迪卡爾正交網(wǎng)格，工程附近網(wǎng)格最密50m×50m，最大網(wǎng)格間距為200m×200m，水平網(wǎng)格數(shù)為224×189，總網(wǎng)格數(shù)24139；設(shè)東和南兩個(gè)開邊界，用實(shí)測(cè)潮位數(shù)據(jù)驅(qū)動(dòng)，由于實(shí)測(cè)潮流潮位數(shù)據(jù)的時(shí)間為1月，故不考慮漳江的徑流；垂向分為1層，計(jì)算時(shí)間步長(zhǎng)1s；由于該海區(qū)潮差較大，采用動(dòng)邊界。

圖2 東山灣網(wǎng)格圖

3.2 模擬運(yùn)算過(guò)程及檢驗(yàn)

模擬采用零初始條件，為了保證計(jì)算的穩(wěn)定性，強(qiáng)迫的邊界潮位從零開始逐步增加，經(jīng)過(guò)兩個(gè)潮周期后達(dá)到正常變化，第三個(gè)潮周期后形成穩(wěn)定的潮波，選擇大潮周期計(jì)算結(jié)果做分析。利用東山大澳中心漁港3個(gè)站位實(shí)測(cè)潮位潮流數(shù)據(jù)作為對(duì)比，結(jié)果見圖3～圖6，從驗(yàn)證結(jié)果看，潮位和流速驗(yàn)證較好，比較準(zhǔn)確地反映了潮汐特征，可以認(rèn)為模擬結(jié)果是可信的。

圖3 潮位驗(yàn)證曲線

圖4 A站位流速流向驗(yàn)證曲線

圖5 B 站位流速流向驗(yàn)證曲線

圖6 C 站位流速流向驗(yàn)證曲線

4 水動(dòng)力變化模擬結(jié)果

4.1 流場(chǎng)變化

工程前后數(shù)學(xué)模型陸地邊界：考慮到大澳中心漁港已經(jīng)批建，而且據(jù)了解，對(duì)面島附近的漁港防波堤正在建設(shè)中，故把中心漁港防波堤等填海區(qū)加入原始海岸線邊界，作為工 程前陸地邊界，然后進(jìn)行數(shù)學(xué)模型的模擬計(jì)算；工程后既包括本工程填海區(qū)，又包括東側(cè)臨近的已批未建的旗濱碼頭填海區(qū)，作為工程后的陸地邊界進(jìn)行數(shù)值模擬；對(duì)工程前后的數(shù)值模擬結(jié)果進(jìn)行比較分析。工程附近海域工程實(shí)施前后海域的大潮漲急、落急的流場(chǎng)見圖7～圖10。

圖7 工程前附近海域漲急流場(chǎng)圖

圖8 工程附近海域工程后漲急流場(chǎng)圖

圖9 工程前附近海域落急流場(chǎng)圖

圖10 工程后附近海域落急流場(chǎng)圖

從圖7～圖10可以看出，工程附近海域工程施工后的潮流場(chǎng)趨勢(shì)仍與工程前相同，工程附近流態(tài)略有變化。本工程的填?；疚挥跒┩恐?，水深較淺，本身流速較小。工程后，本工程對(duì)潮流有阻擋和挑流作用，漲潮時(shí)，潮流從旗濱碼頭和本工程?hào)|側(cè)連線往北流動(dòng)，到達(dá)本工程碼頭后，沿著本工程的東北側(cè)填海區(qū)延線往西北流動(dòng)，在本工程北端往西流動(dòng)；另有一小股漲潮流沿著旗濱碼頭前緣外西北流動(dòng)，然后往西流入旗濱碼頭北側(cè)與本工程形成的呈“凹入”地形的灘涂中，流速較小；落潮，潮流趨勢(shì)相反。

4.2 流速變化

為了進(jìn)一步分析工程實(shí)施前后流速的變化規(guī)律，選擇了工程附近36個(gè)試驗(yàn)點(diǎn)（位置見圖11），進(jìn)行工程前后流速比較，結(jié)果見表1。

圖11 工程前后模擬流速比較點(diǎn)位置圖

表1 工程前后試驗(yàn)點(diǎn)流速比較表（單位m/s）

根據(jù)表1和圖7～圖10中對(duì)試驗(yàn)點(diǎn)進(jìn)行工程前后流速的比較分析，可以看出：工程實(shí)施前，漲潮時(shí)，一部分漲潮流先往西流入工程西南側(cè)的小澳中，然后受地形影響往西北流出，到達(dá)小澳北端的突出角后轉(zhuǎn)往西流動(dòng)；落潮時(shí)相反，小澳內(nèi)潮水匯入落潮流，流速較小。

工程實(shí)施后，流速變化的點(diǎn)主要有14、15、19、20、21、22、25。14點(diǎn)位于本工程與旗濱碼頭形成的凹入地形中，流速變化率最大，工程后該處的流速明顯減小，漲潮減少為60%，落潮減少為67%；根據(jù)前面的分析，漲、落潮流流向變化也很大，漲落潮流從凹入地形口出入。15點(diǎn)位于本工程西側(cè)，由于受工程挑流的影響，漲、落潮最大減少變化量都為0.07m/s。19、25點(diǎn)位于旗濱玻璃廠東側(cè)，受填海區(qū)阻擋的作用，漲落潮的流速都有所減小。點(diǎn)20、21、22漲潮時(shí)，由于工程的挑流束窄作用，流速變大，落潮時(shí)20點(diǎn)由于工程的挑流作用及填海區(qū)匯入20點(diǎn)潮水量的減少，流速變小，21、22在落潮時(shí)由于束窄作用流速略微變大。本工程的建設(shè)對(duì)工程西側(cè)的潮流基本沒有影響，1~13試驗(yàn)點(diǎn)流速變化最大僅為0.01m/s。工程對(duì)北側(cè)較遠(yuǎn)的點(diǎn)17、18、23、24、29、30、35、36也基本沒有影響，流速變化最大也僅為0.01m/s。東側(cè)較遠(yuǎn)處的26、27、28、31、32、33、34變化很小，流速變化均小于等于0.01m/s。

總之，本工程填海會(huì)使碼頭前沿極小段航道區(qū)垂直潮流方向靠碼頭一側(cè)潮流流速增大約0.04m/s，遠(yuǎn)離碼頭一側(cè)潮流流速減小約0.05m/s，對(duì)航道區(qū)其他區(qū)域的潮流流速基本沒有影響；會(huì)使旗濱碼頭附近潮流流速有所減小，減小約0.01m/s左右；對(duì)工程附近其他碼頭及整個(gè)東山灣的潮流流速?zèng)]有影響。

5 結(jié)論

EFDC模型對(duì)于東山灣海域具有較強(qiáng)的適用性，模擬結(jié)果較為準(zhǔn)確，預(yù)測(cè)值與實(shí)測(cè)值的擬合程度較好。模擬結(jié)果顯示，東山華浮碼頭的建設(shè)對(duì)所在海域水動(dòng)力環(huán)境影響較小。該模型有較大的應(yīng)用和推廣價(jià)值，還可以對(duì)未來(lái)發(fā)展進(jìn)行定量的預(yù)測(cè)研究，為實(shí)際的決策過(guò)程提供科學(xué)依據(jù)。

[1] 王翠,孫英蘭,張學(xué)慶.基于EFDC模型的膠州灣三維潮流數(shù)值模擬[J].中國(guó)海洋大學(xué)學(xué)報(bào):自然科學(xué)版,2008,38(5):833-840.

[2] Jin K R.Hamrick J H.Todd T Application of three-dimensional hydrodynamic model for Lake Okeechobee 2000.

[3] Jin K R.Hamrick J H.Todd T Application of three-dimensional hydrodynamic model for Lake Okeechobee 2000.

[4] 陳異暉.基于EFDC模型的滇池水質(zhì)模擬[J].云南環(huán)境科學(xué),2005,24(4):28 – 30.

[5] Lewis C Linker.Gary W Shenk.Robin L Dennis.Jeffry S.Sweeney Cross-Media Models of the Chesapeake Baywatershed and AirshedJ 2000(1-4)

[6] 姜恒志,沈永明,汪守東.甌江口三維潮流和鹽度數(shù)值模擬研究[J].水動(dòng)力學(xué)研究與進(jìn)展A輯,2009,24(1):63 – 70.

[7] JI Z G.MORTON M R.HAMRICKJM Three-dimensional hydrodynamic and water quality modeling in a reservoir 2004

[8] JIN K R.JI Z G.JAMES R T Three-dimensional water quality and SAV modeling of a large shallow lake 2008,(1).

[9] JI Z G.HU G D.SHEN J Three-dimensional modeling of hydrodynamic processes in the St.Lucie Estuary 2007,(1-2).

[10] 李瑞杰,諸裕良,鄭金海.虎門太平水道航道整治工程潮流數(shù)模計(jì)算[J].水運(yùn)工程,2003,(3):38 – 42.