張宇銘， 張淑芳， 韓成偉

(國家海洋環(huán)境監(jiān)測中心，遼寧 大連 116023)

近年來，隨著中國經(jīng)濟和社會的高速發(fā)展，沿海地區(qū)人口不斷聚集，城鎮(zhèn)規(guī)模不斷擴大，其負面效應(yīng)也日益凸顯，企業(yè)和生活污水的大量入海和海洋開發(fā)力度的加大導(dǎo)致了近岸海洋環(huán)境壓力持續(xù)增加。根據(jù)《中華人民共和國海洋環(huán)境保護法》第三十條規(guī)定，“入海排污口位置的選擇，應(yīng)當(dāng)根據(jù)海洋功能區(qū)劃、海水動力條件及有關(guān)規(guī)定，經(jīng)科學(xué)論證后，報設(shè)區(qū)的市級以上人民政府環(huán)境保護行政主管部門審查批準(zhǔn)”。而現(xiàn)有諸多入海排污口布設(shè)時間較早，并未經(jīng)過科學(xué)論證。因此，對沿岸直接入海排污口進行布局合理性評估并優(yōu)化，對于充分利用海洋物理自凈能力，保障海洋資源開發(fā)利用的同時改善近岸海洋環(huán)境非常重要。

目前，對入海排污口選址適宜性或排污布局合理性評價主要采用兩種方法：一是基于水動力數(shù)值模擬方法評估納污海域的水交換能力或納污能力[1-5]，并以此為依據(jù)對排污口布局合理性進行評估。該方法的優(yōu)點是充分考慮了海洋水動力條件對污染物擴散的影響，更符合海洋實際情況。缺點是評價因子單一，忽略了海洋敏感功能區(qū)等其他要素在排污布局中的重要性；二是構(gòu)建基于多因子評估的排污口布局優(yōu)化方法[6-10]。與前一方法相比，該方法同時考慮了生態(tài)敏感區(qū)、工程穩(wěn)定性、工程風(fēng)險和工程經(jīng)濟等多種要素對排污布局的影響，評價因子選取更為全面。然而，該方法由于評價因子眾多導(dǎo)致其對基礎(chǔ)數(shù)據(jù)要求高，實用性和針對性較差，同時基于該方法的研究對于評價排污對敏感海洋功能區(qū)的影響主要以排污口到敏感功能區(qū)的直線距離作為評價指標(biāo)，忽略了海水可流動性的客觀事實，可能造成評價結(jié)果缺乏科學(xué)性。

綜合考慮前人研究的優(yōu)點和不足，本文為滿足防治海洋污染的需要，以充分利用海洋物理自凈能力為出發(fā)點，構(gòu)建綜合考慮入海排污口鄰近海域水交換能力和排污對敏感海洋功能區(qū)環(huán)境影響的入海排污口布局合理性評估方法。考慮到渤海灣灣頂海域承接了環(huán)渤海經(jīng)濟區(qū)和海河水系的大量污染，且海域水交換和自凈能力較差，海洋環(huán)境現(xiàn)狀不容樂觀。因此，本文以渤海灣灣頂海域為例，應(yīng)用構(gòu)建的入海排污口布局合理性評估方法，對沿岸所調(diào)查的直接入海排污口進行了布局合理性評估。

1 入海排污口布局合理性評估方法建立

1.1 評估指標(biāo)

根據(jù)《中華人民共和國海洋環(huán)境保護法》第三十條規(guī)定，評判入海排污口空間布局科學(xué)合理的指標(biāo)應(yīng)包含兩個方面：一方面入海排污后需有利于污染物的稀釋擴散；另一方面需注意減輕陸源污染物對海洋敏感功能區(qū)的影響。陸源污染物入海后的輸運擴散狀況主要取決于排污口鄰近海域的水交換能力，若排污口鄰近海域水動力條件差，水交換能力弱，則污染物易于排污口鄰近海域堆積，從而加重局地海域的污染程度；若排污口鄰近海域水動力條件好，水交換能力強，則污染物排海后可以快速向外輸運擴散，污染物濃度得以稀釋，從而緩解局地海域的污染狀況。但若在水動力作用下大量陸源污染物向海洋生態(tài)保護區(qū)等水質(zhì)要求較高的海洋功能區(qū)輸運，則易引發(fā)嚴重的海洋生態(tài)環(huán)境問題。因此，本文將用兩項指標(biāo)來評估入海排污口的布局合理性：一是入海排污口鄰近海域水交換能力；二是入海排污對敏感海洋功能區(qū)的環(huán)境影響程度。

1.2 入海排污口鄰近海域水交換能力評估方法

本文中入海排污口鄰近海域水交換能力以采用拉格朗日粒子追蹤法[11-12]計算半交換時間[13]作為評判指標(biāo)。具體方法如下：首先，根據(jù)入海排污口的地理分布及鄰近海域水體輸運狀況，將空間位置相鄰較近、對區(qū)域水質(zhì)影響具有累加效應(yīng)的入海排污口劃分為同一個排污分區(qū)，并假定分區(qū)內(nèi)污染物混合均勻，每個分區(qū)分別作為所含入海排污口鄰近海域水交換能力的評價區(qū)域；其次，基于三維水動力數(shù)值模型和拉格朗日粒子追蹤模型，在各排污分區(qū)內(nèi)均勻釋放保守粒子，考察保守粒子在水動力作用下的運移狀況，進一步計算排污分區(qū)的半交換時間，以表征入海排污口鄰近海域水交換能力，定義初始位于某分區(qū)的粒子在水動力作用下輸運50%至分區(qū)外部所用的時間為該分區(qū)的粒子半交換時間；最后，為了對各排污分區(qū)的水交換能力進行分級排序，將各分區(qū)的粒子半交換時間從短到長排列，按照每25%累積率來劃分，將各排污分區(qū)的水交換能力劃分為強、較強、較弱、弱4個等級。水交換能力等級評判方法見表1。

表1 水交換能力評估指標(biāo)Table 1 Evaluation criteria for water exchange ability

1.3 入海排污對敏感海洋功能區(qū)環(huán)境影響程度評估方法

由于海洋特別保護區(qū)、農(nóng)漁業(yè)區(qū)、濱海旅游區(qū)等海洋功能區(qū)對入海污染物較為敏感，水質(zhì)要求通常為一、二類，因此本文將需滿足一、二類水質(zhì)要求的海洋功能區(qū)統(tǒng)稱為敏感海洋功能區(qū)?；谌S水動力模型和拉格朗日粒子追蹤法，在各排污分區(qū)內(nèi)釋放保守粒子來表征陸源污染物，并計算逐時刻粒子進入各敏感海洋功能區(qū)內(nèi)的比例，以此來表征入海排污口在各時刻對各敏感海洋功能區(qū)的影響程度。當(dāng)計算時間足夠長，入海污染物在輸移擴散作用下離開近岸海洋功能區(qū)范圍時，入海排污口對敏感海洋功能區(qū)的影響減小以至于可以忽略不計。以模型計算中粒子釋放時刻至計算時間足夠長以至于大部分粒子離開近岸海洋功能區(qū)范圍時刻的平均粒子比例表征各排污分區(qū)對敏感海洋功能區(qū)的環(huán)境影響程度，該比例越大，說明排污分區(qū)內(nèi)釋放的污染物在敏感海洋功能區(qū)內(nèi)的滯留量越大，即對敏感海洋功能區(qū)可能造成的環(huán)境影響越大。按計算所得各排污分區(qū)在各敏感海洋功能區(qū)內(nèi)的平均粒子比例，將入海排污的環(huán)境影響劃分為4個等級，分別為小、較小、較大、大，評判標(biāo)準(zhǔn)見表2。

1.4 排污布局合理性評估方法

從海洋環(huán)境保護角度出發(fā)，入海排污口鄰近海域水交換能力和排污對敏感海洋功能區(qū)的環(huán)境影響程度兩項指標(biāo)的相對重要性相當(dāng)，因此在評估排污布局合理性時取兩者權(quán)重相同，即排污布局合理性評估分值為水交換能力等級分值與敏感海洋功能區(qū)環(huán)境影響程度等級分值的等權(quán)重加和。

當(dāng)水交換能力等級為“弱”且敏感海洋功能區(qū)環(huán)境影響程度等級為“大”或“較大”時，或當(dāng)水交換能力等級為“較弱”且敏感海洋功能區(qū)環(huán)境影響程度等級為“大”時，說明污染物入海后易在排污口鄰近海域堆積，且排污極易造成敏感海洋功能區(qū)的環(huán)境影響，說明排污布局不合理，亟需進行排污口空間布局的調(diào)整。當(dāng)水交換能力等級為“強”且敏感海洋功能區(qū)環(huán)境影響程度為“小”或“較小”時，或當(dāng)水交換能力等級為“較強”且敏感海洋功能區(qū)環(huán)境影響程度為“小”時，說明污染物入海后能夠較為迅速地擴散，且擴散過程中不易出現(xiàn)污染物在敏感海洋功能區(qū)堆積的情況，即排污布局合理。排污布局合理性評估方法見表3。

表2 排污環(huán)境影響等級評估指標(biāo)Table 2 Evaluation criteria for environmental effects of pollution

表3 排污布局合理性評估方法Table 3 Evaluation criteria for the suitability of drain outlets

2 評估方法應(yīng)用——以渤海灣灣頂海域為例

2.1 渤海灣三維水動力模型的建立與率定

對入海排污口鄰近海域水交換能力和排污對敏感海洋功能區(qū)環(huán)境影響評估均建立在三維水動力模型基礎(chǔ)上。本文以區(qū)域海洋模式Regional Ocean Modeling System(ROMS)[14-17]為基礎(chǔ)，構(gòu)建渤海灣三維水動力模型。模型岸線和水深地形數(shù)據(jù)來自于2014年5月止的最新版海圖，模擬區(qū)域的地形設(shè)置見圖1。模型網(wǎng)格采用曲線正交網(wǎng)格，分辨率為450～600 m，垂向采用S坐標(biāo)[18]，均勻分為10層。東邊界為開邊界，其他均為閉邊界。選取了6個主要分潮(M2、S2、K1、O1、N2、Q1)的合成水位作為開邊界條件，潮汐調(diào)和常數(shù)由國家海洋環(huán)境監(jiān)測中心渤海大區(qū)三維水動力模型提供[19]。海氣界面數(shù)據(jù)(包括風(fēng)應(yīng)力和海氣通量數(shù)據(jù)等)由NCEP再分析資料根據(jù)塊體公式[20]計算。模型初始場包括水位、流場、溫度場和鹽度場?？紤]到水位和流速在模型運行過程中會迅速調(diào)整，因此初始條件中水位和流速均設(shè)置為0；溫度和鹽度的初始場取自《渤海、黃海、東海海洋圖集(水文分冊)》[21]中的氣候態(tài)月平均的溫度和鹽度場。

數(shù)學(xué)學(xué)習(xí)是很有吸引力的事，我們應(yīng)發(fā)揮其自身的魅力，通過理解算理、探索算法感受數(shù)學(xué)帶來的樂趣，激發(fā)學(xué)生對數(shù)學(xué)學(xué)習(xí)的興趣。小學(xué)生以形象思維為主，基于此，我們設(shè)計一些生活情境讓學(xué)生真正理解算理，由算理向算法過渡。

為了驗證渤海灣三維水動力模型模擬結(jié)果的可靠性，將模擬結(jié)果分別與該海域內(nèi)可獲得的水動力實測數(shù)據(jù)進行了比對分析。圖2為用于比對驗證的水動力實測站位分布圖。比對數(shù)據(jù)來自2014年春季4個調(diào)查站位、夏季6個調(diào)查站位以及塘沽水文站，比對要素為水位、潮位、海流流速和流向。

圖1 渤海灣模型地形Fig.1 Bathymetry of the model of the Bohai Bay

本文采用相關(guān)系數(shù)、均方根誤差和技術(shù)評分[22-23]

三項指標(biāo)評估數(shù)值模擬效果。相關(guān)系數(shù)用以評估數(shù)值模擬在趨勢上對于實際情況的反映能力；均方根誤差用以評估模擬與實測數(shù)據(jù)之間的差異性；技術(shù)評分用以評估數(shù)值模擬效果等級，其分級評估標(biāo)準(zhǔn)見表4[22]。三項指標(biāo)的具體計算方法如下：

圖2 渤海灣灣頂海域?qū)崪y站位分布Fig.2 Location of monitoring sites on top of the Bohai Bay

相關(guān)系數(shù)CC(Correlation Coefficient)

CC=

(1)

均方根誤差RMSE(Root-Mean-Square Error)

RMSE=

(2)

技術(shù)評分SS(Skill Score)

(3)

渤海灣三維水動力模型的模擬數(shù)據(jù)與海水實測數(shù)據(jù)的比對結(jié)果見圖4和表5?？梢钥闯?，在數(shù)據(jù)量為23 239個的大樣本統(tǒng)計情況下，實測數(shù)據(jù)與模擬數(shù)據(jù)的相關(guān)系數(shù)總體較高，均方根誤差較小，模擬效果技術(shù)等級為極好，說明采用渤海灣三維水動力模型模擬計算得到的海域水動力場的可靠性較高，可以滿足渤海灣灣頂海域入海排污口布局合理性評估指標(biāo)的計算需求。

表4 數(shù)值模型模擬效果技術(shù)等級標(biāo)準(zhǔn)Table 4 Evaluation criteria for numerical modeling

圖3 模擬(黑線)與實測(紅線)水位對比圖Fig.3 The comparison between simulated (black line) water level and the observations (red line)

2.2 入海排污口分布現(xiàn)狀及排污分區(qū)劃分

根據(jù)2014年渤海灣灣頂海域入海排污口調(diào)查結(jié)果，沿岸有排污布局規(guī)劃需求的入海排污口共21個。為了便于評估現(xiàn)有入海排污口對環(huán)境的影響，根據(jù)排污分區(qū)劃分原則，將渤海灣灣頂海域共劃分為20個排污分區(qū)，分區(qū)外緣線取自岸線向海延伸5 km(見圖4)。其中分區(qū)Ⅳ、Ⅵ、Ⅷ、Ⅹ、ⅩⅢ、ⅩⅣ、ⅩⅦ、ⅩⅩ內(nèi)目前暫無排污口，可作為排污布局優(yōu)化的備選分區(qū)。每個分區(qū)內(nèi)所包含的排污口情況詳見表6。

2.3 渤海灣灣頂海域入海排污口布局合理性評估結(jié)果

渤海灣灣頂海域需滿足二類水質(zhì)要求的海洋功能共有7個(見圖5)，由北向南依次為A：漢沽淺海生態(tài)系統(tǒng)海洋特別保護區(qū)、B：北塘濱海旅游休閑娛樂區(qū)、C：漢沽農(nóng)漁業(yè)區(qū)、D：高沙嶺旅游休閑娛樂區(qū)、E：天津東南部農(nóng)漁業(yè)區(qū)、F：大港濱海濕地海洋特別保護區(qū)、G：馬棚口農(nóng)漁業(yè)區(qū)。在進行排污對敏感海洋功能區(qū)的環(huán)境影響評估時重點針對以上7個敏感海洋功能區(qū)。

采用渤海灣三維水動力模型和粒子追蹤方法，計算得到各排污分區(qū)水交換能力和敏感海洋功能區(qū)環(huán)境影響程度，進而評估現(xiàn)有入海排污口的布局合理性。渤海灣灣頂海域入海排污口均采用閘門控制排放量，夏季為主要排污時間，其他季節(jié)排污極少，因此僅對夏季進行數(shù)值模擬計算并進行排污布局合理性評估。由于排污Ⅲ區(qū)外邊緣建有防波堤，口門很小，幾乎封閉，與外部水交換能力極弱，因此在數(shù)值模擬計算中將此港池區(qū)域視為全封閉區(qū)域，不參與計算。

表5 2014年模擬效果評估Table 5 Evaluation results of simulations in 2014

表6 排污分區(qū)與所含的入海排污口信息表Table 6 Pollution zones and their drain outlets

渤海灣灣頂海域入海排污口空間布局合理性評估結(jié)果見表7和圖6。Ⅲ區(qū)較為封閉，難以與外海進行水交換；Ⅴ區(qū)和Ⅵ區(qū)水交換能力均為較弱，且極易造成敏感海洋功能區(qū)(北塘濱海旅游休閑娛樂區(qū)和漢沽農(nóng)漁業(yè)區(qū))的污染，因此Ⅲ區(qū)、Ⅴ區(qū)和Ⅵ區(qū)排污布局均不合理。渤海灣灣頂海域共有2個入海排污口位于上述布局不合理區(qū)，分別為18號閘涵、19號閘涵(中心漁港)，需優(yōu)化空間布局；海濱浴場沉淀池排水閘、塘鹽揚水站和新馬棚口村北排污口I等3個排污口位于布局合理區(qū)，在滿足達標(biāo)排放前提下無需進行空間布局優(yōu)化；其余16個排污口位于布局較合理區(qū)，在適當(dāng)限制排污總量基礎(chǔ)上暫不需優(yōu)化位置。

圖4 渤海灣灣頂海域排污分區(qū)劃分Fig. 4 Pollution zones of the drain outlets on top of the Bohai Bay

圖5 渤海灣灣頂海域敏感海洋功能區(qū)分布Fig.5 Distribution of sensitive marine functional areas on top of the Bohai Bay

表7 渤海灣灣頂海域各排污分區(qū)排污布局合理性評估Table 7 The evaluation results of suitability of the pollution zones on top of the Bohai bay

Note:①Pollution zone；②Water exchange ability；③Evaluation of water exchange ability；④Environmental influences；⑤Evaluation of environmental influences；⑥Total score；⑦Suitability

圖6 渤海灣灣頂海域各排污分區(qū)排污布局合理性評估結(jié)果Fig.6 The evaluation results of suitability of the pollution zones on top of the Bohai Bay

3 結(jié)語

本文從近岸海洋環(huán)境保護角度出發(fā)，建立了入海排污口布局合理性評估方法。該評估方法以水動力數(shù)值模擬技術(shù)中的拉格朗日粒子追蹤法為基礎(chǔ)，通過充分考慮水文動力條件對入海污染物的輸運擴散作用，對入海排污口鄰近海域水交換能力和排污對敏感海洋功能區(qū)的環(huán)境影響程度兩項指標(biāo)分別進行分級評估，進而綜合評估入海排污口的布局合理性。

文中采用建立的評估方法對渤海灣灣頂海域入海排污口空間布局合理性進行了評估。將該海域劃分為20個排污分區(qū)，對其中包含的21個入海排污口的布局合理性分別評估。結(jié)果表明，現(xiàn)有入海排污口中有2個排污口布設(shè)不合理，需優(yōu)化空間布局；3個排污口位于布局合理區(qū)，在滿足達標(biāo)排放前提下無需進行空間布局優(yōu)化；16個排污口位于布局較合理區(qū)，在適當(dāng)限制排污總量基礎(chǔ)上暫不需優(yōu)化位置。

本文建立的入海排污口布局合理性評估方法是針對有限目標(biāo)的評估方法，重點考慮入海排污對海洋環(huán)境的影響因素，而不考慮經(jīng)濟成本、工程風(fēng)險等其他非海洋因素的影響，主要適用于以改善和保護近岸海洋環(huán)境為出發(fā)點的入海排污口布局評估，為科學(xué)論證區(qū)域入海排污口位置提供一種方法。

[1] 孫長青, 趙可勝, 孫英蘭. 八所近岸海域排污口選址及水質(zhì)預(yù)測[J]. 青島海洋大學(xué)學(xué)報(自然科學(xué)版), 1994, 8: 209-216.

Sun C, Zhao K, Sun Y. Selection of pollutant discharge outlet and prediction of water quality in the offshore area of Basuo [J]. Periodical of Ocean University of Qingdao, 1994, 8: 209-216.

[2] 孫英蘭, 孫長青, 趙可勝. 青島東部開發(fā)區(qū)排污口優(yōu)選[J]. 青島海洋大學(xué)學(xué)報(自然科學(xué)版),1994(s1): 134-141.

Sun Y, Sun C, Zhao K. Optimum location of pollution discharge in the development district of Qingdao [J]. Periodical of Ocean University of Qingdao, 1994(s1): 134-141.

[3] 孫英蘭, 王麗霞, 史峰巖. 洋浦經(jīng)濟開發(fā)區(qū)污水排?？偭靠刂芠J]. 青島海洋大學(xué)學(xué)報(自然科學(xué)版), 2000, 30(3): 383-389.

Sun Y, Wang L, Shi F. The total discharging control at Yangpu Economic Development Zone [J]. Periodical of Ocean University of Qingdao, 2000, 30(3): 383-389.

[4] 王學(xué)昌, 婁安剛, 鄭丙輝. 不同方式污水排海對海水水質(zhì)的影響[J]. 海洋環(huán)境科學(xué), 2002, 21(8): 57-60.

Wang X C, Lou A G, Zheng B H, et al. Effect of sewage discharge on sea water quality at different cases [J]. Marine Environmental Science, 2002, 21(8): 57-60.

[5] 張宇銘, 宋朝陽, 吳克儉, 等. 環(huán)渤海排污口臨近海域水交換能力研究[J]. 中國海洋大學(xué)學(xué)報(自然科學(xué)版), 2014, 44(5): 1-7.

Zhang Y M, Song Z Y, Wu K J, et al. The water exchange of the adjacent seas of the drain outlets into the Bohai Sea [J]. Periodical of Ocean University of China, 2014, 44(5): 1-7.

[6] 汪亮, 張剛, 解建倉, 等. 基于多因素的排污口優(yōu)化布設(shè)模型及其動態(tài)實現(xiàn)[J]. 沈陽農(nóng)業(yè)大學(xué)學(xué)報, 2012, 43(3): 374-377.

Wang L, Zhang G, Xie J C, et al. Drain layout optimization model and dynamic complement based upon multi factors [J]. Journal of Shenyang Agricultural University, 2012, 43(3): 374-377.

[7] 石紅, 翟雅男, 李媛, 等. 入海排污口選址生態(tài)適宜性分析與比選研究——以深圳近岸海域為例[J]. 安全與環(huán)境工程, 2016, 23(3): 75-82.

Shi H, Zhai Y, Li Y, et al. Ecology suitability evaluation and comparison for the location of marine outfalls - A case study of coastal areas in Shenzhen city[J]. Safety and Environmental Engineering, 2016, 23(3): 75-82.

[8] Zhao Y W, Qin Y, Chen B, et al. GIS-based optimization for the location of sewage treatment plants and sewage outfalls—A case study of Nansha District in Guangzhou City, China [J]. China Communications in Nonlinear Science and Numerical Simulation, 2009, 14(4): 1746-1757.

[9] Alvarez-Vázquez L J, Martínez A, Rodríguez C, et al. Numerical optimization for the location of wastewater outfalls [J]. Computational Optimization and Applocations, 2002, 22(3): 399-417.

[10] Alvarez-Vázquez L J, Martínez A, Rodríguez C, et al. Mathematical analysis of the optimal location of wastewater outfalls [J]. Ima Journal of Applied Mathematics, 2002, 67(1): 23-39.

[11] Bilgili A, Proehl J A, Lynch D R, et al. Estuary ocean exchange and tidal mixing in a Gulf of Maine Estuary: A Lagrangian modeling study [J]. Estuarine, Coastal and Shelf Science, 2005, 65: 607-624.

[12] Liu W C, Chen W B, Hsu M H. Using a three-dimensional particle-tracking model to estimate the residence time and age of water in a tidal estuary [J]. Computers and Geosciences, 2011, 37: 1148-1161.

[13] 魏皓, 田恬, 周鋒, 等. 渤海水交換的數(shù)值研究——水質(zhì)模型對半交換時間的模擬 [J]. 青島海洋大學(xué)學(xué)報(自然科學(xué)版), 2002, 32(4): 519-525.

Wei H, Tian T, Zhou F, et al. Numerical study on the water exchange of the Bohai Sea: Simulation of the half-life time by dispersion model [J]. Periodical of Ocean University of Qingdao, 2002, 32(4): 519-525.

[14] Shchepetkin A F, McWilliams J C. Quasi-monotone advection schemes based on explicit locally adaptive dissipation [J]. Monthly Weather Review, 1998, 126: 1541-1580.

[15] ShchepetkinA F, McWilliams J C. A method for computing horizontal pressure-gradient force in an oceanic model with a nonaligned vertical coordinate [J]. Journal of Geophysical Research, 2003, 108(C3): 3090, doi: 10. 1029/2001JC001047.

[16] Moore A M, Arango H G, Miller A J, et al. A comprehensive ocean prediction and analysis system based on the tangent linear and adjoint components of a regional ocean model [J]. Ocean Modelling, 2004, 7: 227-258.

[17] Shchepetkin A F, McWilliams J C. The regional ocean modeling system: A split-explicit, free-surface, topography followingcoorddinates ocean model [J]. Ocean Modelling, 2005, 9: 347-404.

[18] Song Y, Haidvogel D B. A semi-implicit ocean circulation model using a generalized topography-following coordinate system [J]. Journal of Computational Physics, 1994, 115(1): 228-244.

[19] 張宇銘, 張淑芳, 宋朝陽, 等. 基于質(zhì)點追蹤方法的渤海水交換特性[J]. 海洋環(huán)境科學(xué), 2014, 33(3): 412-417.

Zhang Y M, Zhang S F, Song Z Y, et al. Numerical simulation of water exchange characteristics in Bohai Sea based on particle tracing method [J]. Marine Environmental Science, 2014, 33(3): 412-417.

[20] Fairall C W, Bradley E F, Rogers D P, et al. Bulk parameterization of air-sea fluxes for tropical ocean-global atmosphere coupled-ocean atmosphere response experiment [J]. Journal of Geophysical Research, 1996, 101: 3747-3764.

[21] 海洋圖集編委會. 渤海、黃海、東海海洋圖集(水文分冊) [M]. 北京: 海洋出版社, 1992.

Editorial Board for Marine Atlas. Marine Atlas of Bohai Sea, Yellow Sea, East China Sea: Hydrology [M]. Beijing: Ocean Press, 1992.

[22] Allen J I, Somerfield P J, Gilbert F J. Quantifying uncertainty in high-resolution coupled hydrodynamic-ecosystem models [J]. Journal of Marine Systems, 2007, 64： 3-14, doi: 10. 1016/j. jmarsys. 2006. 02. 010.

[23] Liu Y, Mac Creany P, Hickey B M, et al. Evaluation of a coastal ocean circulation model for the Columbia River plume in summer 2004 [J]. Journal of Geophysical Research, 2009, 114, C00B04, doi: 10. 1029/2008JC004929.