海洋水質(zhì)模型研究進(jìn)展及發(fā)展趨勢(shì)

于 寒，楊 靜，劉桂梅

海洋水質(zhì)模型研究進(jìn)展及發(fā)展趨勢(shì)

于 寒，楊 靜，劉桂梅

（國(guó)家海洋環(huán)境預(yù)報(bào)中心國(guó)家海洋局海洋災(zāi)害預(yù)報(bào)技術(shù)研究重點(diǎn)實(shí)驗(yàn)室，北京100081）

介紹了水質(zhì)模型的發(fā)展階段，從一維、二維、三維的角度分析了海洋水質(zhì)模型應(yīng)用情況，展望了海洋水質(zhì)模型未來(lái)的發(fā)展趨勢(shì)，可為控制污染物排放、海洋水質(zhì)監(jiān)測(cè)及水資源管理提供系統(tǒng)參考。

海洋水質(zhì)模型；研究進(jìn)展；發(fā)展趨勢(shì)

1 引言

沿岸多為人口密集區(qū)，人為活動(dòng)產(chǎn)生的污染物加劇了沿岸海域水污染程度，生態(tài)環(huán)境也受到一定程度的影響，加劇了人類可利用水資源的流失，因此水污染的治理與水資源科學(xué)管理規(guī)劃，對(duì)人類水資源的可持續(xù)利用顯得尤為關(guān)鍵。水質(zhì)模型著重研究的是水環(huán)境中水質(zhì)狀態(tài)變量（溶解氧、營(yíng)養(yǎng)鹽等）的遷移、轉(zhuǎn)化規(guī)律，以及影響因素之間的相互關(guān)系，可用于預(yù)測(cè)該海域水質(zhì)的發(fā)展趨勢(shì)，所以水質(zhì)模型不僅是水環(huán)境科學(xué)研究的主要內(nèi)容，同時(shí)也是水環(huán)境研究中必不可少的工具[1]。水質(zhì)模擬可分為定性模擬與定量模擬，在實(shí)際水質(zhì)問題的研究中，定量模擬更能反應(yīng)水質(zhì)在空間和時(shí)間上的發(fā)展變化，還可以利用定量模擬評(píng)價(jià)、預(yù)測(cè)、選擇合適的污染控制方案，為水質(zhì)管理規(guī)定的制定以及污染物排放標(biāo)準(zhǔn)的確定提供依據(jù)[2]。

建立可靠的水質(zhì)模型，對(duì)水質(zhì)狀態(tài)變量進(jìn)行準(zhǔn)確模擬及預(yù)測(cè)，是科研工作者多年來(lái)致力研究的重點(diǎn)，伴隨著實(shí)際應(yīng)用的多樣化需求，水質(zhì)模型的應(yīng)用從最初河水自凈化研究到水質(zhì)評(píng)價(jià)、污染物行為預(yù)測(cè)及過(guò)程模擬、水環(huán)境綜合管理等方面的研究，理論從最初的氧平衡模型到現(xiàn)在的模糊理論、灰色理論、隨機(jī)理論，研究范圍從河流擴(kuò)展到湖泊、河口、水庫(kù)、海岸，這些都極大的推動(dòng)了水質(zhì)模型的發(fā)展[3]。本文綜述了水質(zhì)模型研究進(jìn)展，分析了不同維數(shù)海洋水質(zhì)模型及其應(yīng)用，并展望了海洋水質(zhì)模型的發(fā)展趨勢(shì)。對(duì)海洋水質(zhì)模型的研究應(yīng)用，可以為海洋水質(zhì)模擬、污染物排放計(jì)算、水質(zhì)監(jiān)測(cè)、水資源管理提供系統(tǒng)的參考，使用戶對(duì)模型有初步大概的認(rèn)識(shí)，能夠根據(jù)需要選擇適合的海洋水質(zhì)模型研究區(qū)域海洋的具體問題。

2 水質(zhì)模型研究進(jìn)展

氧平衡模型是第一個(gè)水質(zhì)模型，早在1925年，Streeter等[4]在對(duì)Ohio河污染及其自凈化的綜合研究中，提出了最初的氧平衡模型形式，該模型假定河流自凈化過(guò)程中，存在兩個(gè)維持水中溶解態(tài)氧平衡的相反過(guò)程，即消耗水中溶解氧的氧化反應(yīng)，以及補(bǔ)充水中溶解氧的復(fù)氧過(guò)程，模式方程如下：

式中：L為河水中有機(jī)物需氧量，D為河水中的氧虧值，K1河水脫氧速率系數(shù)，K2河水復(fù)氧速度系數(shù)，t為反應(yīng)時(shí)間。

基于最初單變量的S-P（Streeter-Phelps）水質(zhì)模型，水質(zhì)模型在后續(xù)的發(fā)展中考慮了更多的相關(guān)過(guò)程，引入了新的變量，研究對(duì)象多樣化，數(shù)學(xué)工具的選擇更為多樣，針對(duì)水質(zhì)模型的發(fā)展階段，不同學(xué)者做出了不同的歸納[5-13]，徐祖信等[6]按時(shí)間將水質(zhì)模型的發(fā)展大致分為3個(gè)階段（見表1）：1925—1980年為發(fā)展的第一階段，水質(zhì)本身為研究主體，研究中注重水體內(nèi)部水質(zhì)各組分之間的發(fā)展規(guī)律，污染負(fù)荷輸入停留在點(diǎn)源層面上，經(jīng)歷了水質(zhì)模型空間維數(shù)（一維、二維至三維）、研究范圍以及應(yīng)用領(lǐng)域上的發(fā)展，考慮影響生化需氧量-溶解氧（Biological Oxidation-Dissolved Oxygen，BOD-DO）的因素對(duì)SP模型進(jìn)行改進(jìn)；1980—1995年為發(fā)展的第二階段，狀態(tài)變量增多的同時(shí)，水動(dòng)力模型、底泥等作用被納入到多維模型系統(tǒng)內(nèi)部，并且連接流域模型，使面污染源可以作為初始輸入，由于增加了水質(zhì)組分的約束，可以大大減少模擬的主觀性，如美國(guó)國(guó)家環(huán)保局開發(fā)的WASP（The Water Quality Analysis Simulation Program）系列模型，多介質(zhì)模型[14-15]；1995年至今為水質(zhì)模型發(fā)展的第三階段，考慮到大氣中沉降污染物（如有機(jī)化合物、金屬以及氮化合物）的輸入對(duì)河流水質(zhì)的影響愈加顯著，所以建立了大氣污染模型將大氣沉降（動(dòng)態(tài)或靜態(tài)）連接到給定水域。

海洋水質(zhì)模型是水質(zhì)模型在河口、海岸、海灣等區(qū)域數(shù)值計(jì)算中的應(yīng)用，針對(duì)海灣、海岸等區(qū)域的重金屬污染以及富營(yíng)養(yǎng)化、水質(zhì)監(jiān)測(cè)、環(huán)境容量計(jì)算的需求，以基本水質(zhì)模型為基礎(chǔ)，用以解決海水有關(guān)水質(zhì)的問題。海洋水質(zhì)模型嚴(yán)格意義上講，是水質(zhì)模型的一個(gè)應(yīng)用領(lǐng)域，在水質(zhì)模型發(fā)展的第二階段逐漸得到應(yīng)用，伴隨實(shí)際需求的不斷增加，從最初對(duì)河流問題的研究發(fā)展到對(duì)海灣、近岸等水體的研究[5]。

從海洋水質(zhì)模型的研究來(lái)看，美國(guó)、丹麥對(duì)海洋水質(zhì)模型的研究走在世界前列，美國(guó)環(huán)保局開發(fā)了WASP系列、丹麥水動(dòng)力研究所研究開發(fā)了MIKE系列模型，并且不斷對(duì)模型修訂、增強(qiáng)，推動(dòng)了海洋水質(zhì)模型的發(fā)展。以WASP模型為例，自最原始的版本在1983年發(fā)布以來(lái)，經(jīng)過(guò)不斷的修訂完善，用戶界面更加友好，操作系統(tǒng)適用性逐步增強(qiáng)，運(yùn)算速度顯著提高，廣泛應(yīng)用于海灣、近岸的重金屬污染、富營(yíng)養(yǎng)化過(guò)程[16]。國(guó)內(nèi)對(duì)海洋水質(zhì)模型的研究起步較晚，也取得一些成果：江毓武等[17]基于三維變邊界σ坐標(biāo)下的污染擴(kuò)散模型建立了廈門水質(zhì)模型；倪晉仁等[18]建立基于沿深度平均的二維平流-擴(kuò)散模型的水質(zhì)模型；王昆等[19]建立非保守型水質(zhì)預(yù)測(cè)預(yù)報(bào)模型，考慮化學(xué)分解和生物的制約作用；沈永明等[20]將污染物擴(kuò)散輸移的湍流模型與多組分污染物生化轉(zhuǎn)化模型結(jié)合，建立了綜合考慮物理、生物、化學(xué)過(guò)程的近岸海域多組分三維水質(zhì)動(dòng)力學(xué)模型。在實(shí)際需求的推動(dòng)下，國(guó)內(nèi)學(xué)者不僅基于基本方程建立適用于海域的水質(zhì)模型，同時(shí)也引進(jìn)國(guó)外廣泛應(yīng)用的海洋水質(zhì)模型，為區(qū)域海洋污染控制方案的制定提供了一定的參考。林衛(wèi)青等[21]利用丹麥水動(dòng)力研究所開發(fā)的MIKE21水質(zhì)模型在流場(chǎng)、溫鹽場(chǎng)、泥沙場(chǎng)的基礎(chǔ)上，以高錳酸鹽、氨氮等變量為模擬指標(biāo)，建立該區(qū)域的二維水質(zhì)模型，為環(huán)境容量的計(jì)算與分配提供了理論依據(jù)。

水質(zhì)狀態(tài)變量的計(jì)算需要依托水動(dòng)力模型提供基本的物理場(chǎng)，近年來(lái)應(yīng)用較多的海洋水質(zhì)模型大致可以分為兩類：一類僅模擬水質(zhì)的狀態(tài)變量，需要耦合其它的水動(dòng)力模型進(jìn)行實(shí)際的數(shù)值模擬，這類模型比較典型的為WASP、CE-QUAL-ICM（Corps of Engineers Integrated Compartment Water Quality Model）等；另一類可同時(shí)模擬水動(dòng)力和水質(zhì)，應(yīng)用于海洋區(qū)域比較典型的模型為地表水模型系統(tǒng)（Surface-Water Modeling System，SMS）。國(guó)內(nèi)外學(xué)者利用海洋水質(zhì)模型針對(duì)不同區(qū)域的海洋水質(zhì)模擬做了大量研究，對(duì)該海域水質(zhì)狀態(tài)變量變化趨勢(shì)的了解、預(yù)測(cè)，污染物的控制、環(huán)境容量的計(jì)算等方面都提供了一定的參考。例如針對(duì)膠州灣海域：Li等[22]將三維水動(dòng)力模型與從膠州灣圍隔實(shí)驗(yàn)的氮磷動(dòng)力學(xué)模型改編得到的水質(zhì)模型耦合，利用三維水質(zhì)模型對(duì)膠州灣的氮、磷容量進(jìn)行評(píng)估；Zhang等[23]利用三維物理和水質(zhì)模型耦合，其中物理模型采用的ECOM-si（Estuarine，Coastal and Ocean Model（semi-implicit））的修改版本，包括潮間帶的干濕過(guò)程，水質(zhì)模型采用WASP5水質(zhì)模型，研究了營(yíng)養(yǎng)鹽、生化需氧量、溶解氧、浮游植物等狀態(tài)變量的時(shí)空變化情況。Khangaonkar等[24]為了解普及特海灣營(yíng)養(yǎng)負(fù)荷對(duì)水質(zhì)以及生態(tài)健康的影響，采用非結(jié)構(gòu)有限體積海岸模型（An Unstructured Grid,Finite-Volume Coastal Ocean Model，F(xiàn)VCOM）提供水動(dòng)力條件，用CE-QUAL-ICM水質(zhì)模型模擬生化過(guò)程，通過(guò)模擬19個(gè)狀態(tài)變量來(lái)計(jì)算藻類生產(chǎn)、消亡及對(duì)溶解氧的影響。

表1 水質(zhì)模型發(fā)展的3個(gè)階段

3 水質(zhì)模型分類及常用的海洋水質(zhì)模型

3.1 水質(zhì)模型的分類

根據(jù)不同的分類標(biāo)準(zhǔn)，可將水質(zhì)模型分為不同的種類：水質(zhì)模型在空間維數(shù)上可以分為零維、一維、二維、三維，在變量數(shù)量上可以分為單變量、多變量模型，按照水質(zhì)組分時(shí)間變化特征可以分為動(dòng)態(tài)模型、穩(wěn)態(tài)模型，按照數(shù)學(xué)方法可以分為隨機(jī)性模型、確定模型、線性模型、非線性模型，根據(jù)描述水體的不同可以分為河流、湖泊、河口、海灣、地下水模型，根據(jù)對(duì)象的不同可以分為溶解氧、溫度、重金屬、有毒有機(jī)物、放射性模型等不同的分類[12]。

3.2 海洋水質(zhì)模型

3.2.1 一維海洋水質(zhì)模型

一維水質(zhì)模型多適用于狹長(zhǎng)的河流、渠道，考慮狀態(tài)變量在一個(gè)方向上的變化，忽略另外兩個(gè)方向上的變化，一般常用的類型為水平一維和垂向一維水質(zhì)模型，適用于沿河道變化顯著和在垂直方向上變化顯著的情況。

基于質(zhì)量守恒原理推導(dǎo)得到的一維水質(zhì)模型的數(shù)學(xué)表達(dá)式如下[25]：

式中：C1為水質(zhì)濃度，A為過(guò)水面積，Q為流量，Ex為縱向分散系數(shù)，Sc為動(dòng)力學(xué)反應(yīng)源匯項(xiàng)，Wc是外源輸入項(xiàng)。

一維水質(zhì)模型在狹長(zhǎng)河道得到較好應(yīng)用[26-27]，在近岸、海灣等海域用一維水質(zhì)模型可以模擬簡(jiǎn)化的物理和生化過(guò)程。楊紅等[28]在WASP模型基礎(chǔ)上，運(yùn)用系統(tǒng)動(dòng)力學(xué)軟件Stella9.0.2建立了適用于海洋圍隔生態(tài)系統(tǒng)的一維箱式水質(zhì)動(dòng)力學(xué)模型，利用象山港圍隔生態(tài)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行了模型驗(yàn)證和參數(shù)率定工作，成功模擬了不同水溫條件下圍隔水質(zhì)情況。建立海洋圍欄研究水質(zhì)的變化情況，為一維海洋水質(zhì)模型在近岸、海灣的應(yīng)用提供了參考。

3.2.2 二維海洋水質(zhì)模型

二維海洋水質(zhì)模型多為水平方向的二維模型，在海岸、海灣等區(qū)域水深較淺，表層和底層的物理、生物化學(xué)性質(zhì)差別較小，在外部強(qiáng)迫下，較易達(dá)到均一狀態(tài)，在相對(duì)于變量的水平變化，垂直變化可以忽略的情況下，水平二維水質(zhì)變量是較為合理的選擇。

MIKE21模型由丹麥水動(dòng)力研究所（Danish Hydraulic Institude，DHI）研究開發(fā)，在海域、河口地區(qū)的污染物輸運(yùn)、擴(kuò)散研究等方面得到廣泛應(yīng)用。林衛(wèi)青等選用了MIKE21建立了長(zhǎng)江口及毗鄰海域的二維水質(zhì)模型，其基本的對(duì)流擴(kuò)散方程如下：

式中：C為物質(zhì)濃度，u、v為x、y方向流速，Dx、Dy為x、y方向紊動(dòng)擴(kuò)散系數(shù)，s為源匯項(xiàng)，K為綜合降解系數(shù)。針對(duì)該海域的氮磷超標(biāo)問題，利用二維水質(zhì)模型計(jì)算得到該海區(qū)排污口排放總量的控制方案，并利用三維生態(tài)模型進(jìn)行校驗(yàn)，結(jié)果顯示全海域?qū)嵤┛刂品桨钢螅喑北┌l(fā)的概率和面積都大大減少[21]。Babu等[29]利用MIKE21計(jì)算了Kochi的離岸排水口在3種不同排放速率的情況下的吸收能力，基于距排水口8 km一橢圓區(qū)域?qū)Ω連OD（生化需氧量）值的限制，來(lái)確定BOD的每天合適的排放量。以上國(guó)內(nèi)外學(xué)者利用MIKE21水質(zhì)模型計(jì)算污染物的濃度，考慮了點(diǎn)源負(fù)荷對(duì)研究區(qū)域的影響，根據(jù)具體的排放限制，確定了合理的排放方案，體現(xiàn)了MIKE21海洋水質(zhì)模型在水質(zhì)管理規(guī)劃的應(yīng)用，具有較好的效果。實(shí)際研究中，由于實(shí)時(shí)的水質(zhì)數(shù)據(jù)難以獲得，利用海洋水質(zhì)模型來(lái)掌握水質(zhì)變量的變化情況是必要的，可以根據(jù)預(yù)定基本功能所要求的水質(zhì)及水域自凈能力，確定總排放量、日排放量、排放分布情況，以較小的人為代價(jià)獲得較好的水質(zhì)效果[5,30-31]。

RM4模型是美國(guó)陸軍工程兵團(tuán)開發(fā)的基于有限元法的二維污染物輸移水質(zhì)模型，主要應(yīng)用于河口、河流、湖泊、海岸等地區(qū)的污染物輸移擴(kuò)散、排污口排放量等方面的計(jì)算[32-34]。馮媛媛利用水動(dòng)力數(shù)學(xué)模型RM2與水質(zhì)模型RM4，對(duì)甌江河口工程前后的水動(dòng)力及水質(zhì)變化情況作了研究，RM4的基本方程如下：

式中：h為水深，C為污染物濃度，t為時(shí)間，u、v為x、y方向流速，Dx、Dy為x、y方向紊動(dòng)擴(kuò)散系數(shù)，k為衰減系數(shù)，σ為污染物源匯項(xiàng)，R(C)為降雨蒸發(fā)。以化學(xué)需氧量（Chemical Oxygen Demand, COD）作為污染物進(jìn)行擴(kuò)散研究，結(jié)果顯示圍墾工程對(duì)甌江河口的水質(zhì)不會(huì)產(chǎn)生影響[34]。利用海洋水質(zhì)模型考慮污染物在水體中的行為，應(yīng)考慮該海域的主要生物化學(xué)過(guò)程以及點(diǎn)源、非點(diǎn)源污染物的輸入，對(duì)甌江河口水質(zhì)的模擬未考慮到COD的降解，所以可能引起與實(shí)際情況的差異。海洋水質(zhì)模型最主要的功能即是對(duì)污染物在水體中行為進(jìn)行模擬和預(yù)測(cè)，由于污染物進(jìn)入水體之后，隨水體遷移過(guò)程中會(huì)受到水文、理化、水力等因素的影響，會(huì)發(fā)生降解、混合等復(fù)雜的行為，在借助模型來(lái)了解污染物的遷移運(yùn)動(dòng)規(guī)律時(shí)，應(yīng)注意考慮影響該區(qū)域水質(zhì)狀態(tài)變量的主要過(guò)程[5]。

3.2.3 三維海洋水質(zhì)模型

三維海洋水質(zhì)模型適用于流場(chǎng)環(huán)境復(fù)雜，水質(zhì)狀態(tài)變量在水平方向和垂直方向均有顯著變化的情況。水質(zhì)狀態(tài)變量的模擬需要基本的水動(dòng)力場(chǎng)信息（水位、流等），根據(jù)關(guān)注區(qū)域的基本水文情況選擇海洋水質(zhì)模型的水動(dòng)力模塊，對(duì)于只模擬水質(zhì)狀態(tài)變量的海洋水質(zhì)模式而言，需要耦合其他水動(dòng)力模塊。

Cerco和Cole開發(fā)的CE-QUAL-ICM，最初是應(yīng)用于美國(guó)切薩皮克灣（Chesapeake Bay）的富營(yíng)養(yǎng)化模型，能夠模擬一維、二維、三維水體結(jié)構(gòu)，可以模擬多種水質(zhì)狀態(tài)變量（多種藻類、浮游動(dòng)物、溶解氧等），由于CE-QUAL-ICM模型只模擬水質(zhì)，所以需要物理模型提供相應(yīng)的流、水位、擴(kuò)散系數(shù)等水動(dòng)力信息[35]。

該模型在國(guó)內(nèi)外河口、海灣都得到應(yīng)用，Taeyun等將FVCOM與CE-QUAL-ICM耦合，方程如下：

式中：Vj是第j個(gè)控制體的體積，Cj是第j個(gè)控制體的濃度，Qk是流過(guò)第j個(gè)控制體的k流面的體積流量，Ck是流經(jīng)k流面的流體濃度，Ak是k流面的面積，Dk是k流面的擴(kuò)散系數(shù)，n是與第j個(gè)控制體相連的流面的數(shù)量，S是第j個(gè)控制體中的外部壓力和動(dòng)力源匯，t是時(shí)間坐標(biāo)，x是空間坐標(biāo)。對(duì)比4個(gè)試驗(yàn)的模式結(jié)果與解析解，結(jié)果顯示二者耦合可進(jìn)一步考慮用于評(píng)估人為引起的該海域營(yíng)養(yǎng)富集，也推動(dòng)了水動(dòng)力模型與水質(zhì)模型耦合的研究應(yīng)用[36]。Khangaonkar等[24]采用FVCOM與CE-QUAL-ICM耦合，由FVCOM提供水動(dòng)力輸出，包括：溫度、鹽度、自由表面高度、水位、速度分量及邊界通量等，通過(guò)對(duì)19個(gè)狀態(tài)變量的模擬，研究了普及特海灣營(yíng)養(yǎng)負(fù)荷對(duì)水質(zhì)以及生態(tài)健康的影響。Kim等[37]將ROMS與CE-QUAL-ICM模型耦合，考慮河流以及水閘的營(yíng)養(yǎng)鹽點(diǎn)源負(fù)荷，研究了新萬(wàn)金（Saemangeum）潮汐湖及毗鄰的韓國(guó)海域水質(zhì)狀態(tài)變量的季節(jié)變化。CE-QUAL-ICM模型可以耦合不同的水動(dòng)力模型，可以輸入點(diǎn)源、非點(diǎn)源負(fù)荷，靈活性較強(qiáng)，也可以與基于個(gè)體的魚類生物能學(xué)模型耦合[38]，應(yīng)用前景廣泛。

WASP由美國(guó)環(huán)保局開發(fā)，可模擬一維、二維、三維水質(zhì)問題的水質(zhì)模型，由水動(dòng)力模型程序（Hydrodynamic Model，DYNHYD）與水質(zhì)模型程序WASP組成，水動(dòng)力模型程序DYNHYD適用于一維水動(dòng)力情況，水質(zhì)模型程序WASP可以與其連接計(jì)算，也可以耦合其他水動(dòng)力模型[18]。其中水質(zhì)模型程序WASP由兩個(gè)子程序組成：有毒化學(xué)物模型（The Toxic Chemical Model，TOXI）和富營(yíng)養(yǎng)化模型（Eutrophication Model，EUTRO），EUTRO可以用于模擬傳統(tǒng)污染物的遷移轉(zhuǎn)化規(guī)律（如DO、BOD和富營(yíng)養(yǎng)化），TOXI用于有毒物質(zhì)遷移轉(zhuǎn)化規(guī)律（有機(jī)化學(xué)物、金屬、沉積物等）[39-44]。

Zheng等利用修訂的ECOM-si與WASP5耦合，研究了薩蒂拉（Satilla）河口的低溶解氧濃度的原因[45]，水質(zhì)成分濃度的質(zhì)量平衡式為：

式中：C為水質(zhì)成分的濃度，u、v、w笛卡爾坐標(biāo)系下x、y、z方向的流速分量，Ah、Kh分別為水平粘性擴(kuò)散系數(shù)和垂直渦動(dòng)擴(kuò)散系數(shù)，S為水質(zhì)變量的內(nèi)部源匯項(xiàng)，W0是點(diǎn)源和非點(diǎn)源的外部輸入，t表示時(shí)間。通過(guò)模式結(jié)果分析得到結(jié)論為：河口的溶解氧值低，原因是潮間帶的沉積需氧量高，是主要的DO（溶解氧）消耗者，而與人類活動(dòng)無(wú)關(guān)[45]。

WASP模型靈活性較強(qiáng)，可以輸入點(diǎn)源、非點(diǎn)源營(yíng)養(yǎng)鹽負(fù)荷，與其他模型的耦合能力較強(qiáng)，可以與環(huán)境流體動(dòng)態(tài)模型地表水模擬系統(tǒng)（Environmental Fluid Dynamics Code，EFDC）耦合，兼具了EFDC極強(qiáng)的問題適應(yīng)能力、主流的數(shù)值方法和系統(tǒng)開發(fā)方法、多種水動(dòng)力過(guò)程、多種模擬計(jì)算方案等優(yōu)點(diǎn)[33]。李林子等[46]耦合WASP和EFDC模型研究了突發(fā)水污染事故影響的預(yù)測(cè)預(yù)警，取得了較好的結(jié)果。Xiong[47]選取EFDC和WASP作為水動(dòng)力和水質(zhì)模型，將水動(dòng)力、水質(zhì)、沉積模塊耦合，在Mobile海灣的水動(dòng)力、水質(zhì)以及沉積過(guò)程的研究中驗(yàn)證了三者耦合的可行性。

MIKE3是由丹麥水動(dòng)力研究所（DHI）研究開發(fā)的一款用于模擬水動(dòng)力、水質(zhì)、泥沙的專業(yè)工程軟件，主要應(yīng)用于水庫(kù)、港口、河流、湖泊、河口、海岸，具有先進(jìn)的前后處理功能和友好的用戶界面，在大流域、長(zhǎng)時(shí)間的數(shù)值模擬方面有著突出優(yōu)點(diǎn)[48]。

MIKE3模型應(yīng)用廣泛，王勇智等[49]利用MIKE3水動(dòng)力模型耦合粒子追蹤模塊，計(jì)算了渤海灣環(huán)流系統(tǒng)及其驅(qū)動(dòng)下的水質(zhì)點(diǎn)運(yùn)移規(guī)律和變化特征，為污染物輸運(yùn)的研究提供了基礎(chǔ)。黃少彬等[50]利用MIKE3模型建立珠江口河口附近的三維水動(dòng)力模型，計(jì)算出珠江口水體交換時(shí)間，并分析河流流量、外海潮汐與珠江口水體交換的相關(guān)性，對(duì)水體交換的定量計(jì)算，有利于水質(zhì)環(huán)境治理和污染物擴(kuò)散研究。陳瑤泓伶等[51]利用MIKE3的三維水動(dòng)力模型和對(duì)流擴(kuò)散模塊（Advection-Dispersion Model，AD）對(duì)天津南港污水排海進(jìn)行污染物輸移擴(kuò)散數(shù)值模擬，排污口放置水下9 m左右，結(jié)果表明深海排放利于控制表層水體污染物含量的增量，為短期污染物排放控制提供了一定的參考。

地表水模型系統(tǒng)——SMS，由美國(guó)Brigham Young大學(xué)環(huán)境模型研究實(shí)驗(yàn)室開發(fā)，可用于一維、二維、三維水文建模、設(shè)計(jì)，可以模擬和分析地表水運(yùn)動(dòng)規(guī)律，包括預(yù)處理及后處理平臺(tái)[6,11]。SMS軟件界面高度可視化，其中程序以及代碼都是美國(guó)政府公開的，包含來(lái)自美國(guó)陸軍工程兵水道實(shí)驗(yàn)站（U. S.Army Corps of Engineers Waterways Experiment Station）與美國(guó)聯(lián)邦公路管理局的計(jì)算模塊[52-53]。

SMS模型系統(tǒng)主要包含水動(dòng)力和泥沙模型，水質(zhì)變化過(guò)程有限，但可通過(guò)耦合其他水質(zhì)模型模擬水質(zhì)狀態(tài)變量。SMS模型系統(tǒng)有較為廣泛的應(yīng)用，吳志易等[54]利用SMS建立杭州灣海域模型，實(shí)現(xiàn)海域流場(chǎng)的可視化，并結(jié)合水動(dòng)力學(xué)模型進(jìn)行潮流場(chǎng)模擬，結(jié)合波浪動(dòng)力學(xué)模型進(jìn)行波浪場(chǎng)模擬，得到不同條件下的水質(zhì)擴(kuò)散模型，對(duì)于分析不同假設(shè)條件下的海域水流場(chǎng)變化有一定的參考價(jià)值。陳志娟等[55]通過(guò)SMS模擬了黃河口周邊大范圍海域的流場(chǎng)分布情況，分析了該海域的泥沙擴(kuò)散輸運(yùn)情況。魯海燕[56]利用SMS在杭州灣建立大范圍平面二維水流模型,通過(guò)實(shí)測(cè)資料和水下地形的驗(yàn)證該模型滿足預(yù)測(cè)要求。SMS模型具有較強(qiáng)的計(jì)算功能，對(duì)復(fù)雜流場(chǎng)、岸線的模擬都較好，且精度也較高。

4 海洋水質(zhì)模型發(fā)展趨勢(shì)

4.1 水質(zhì)模型內(nèi)部發(fā)展趨勢(shì)

水質(zhì)模型在建立之初注重的是流域水體本身的發(fā)展，隨著現(xiàn)實(shí)應(yīng)用的需求，經(jīng)歷由點(diǎn)源向面源的過(guò)渡，逐漸加入了底泥模塊、水動(dòng)力模塊、大氣污染模塊[6]。最初的水質(zhì)模型只包含DO、BOD在內(nèi)的少數(shù)水質(zhì)指標(biāo)，隨著不同機(jī)構(gòu)對(duì)水質(zhì)模型的不斷開發(fā)、修正、增強(qiáng)，狀態(tài)變量明顯增加。例如MIKE模型系統(tǒng)的變量包括DO、BOD、水溫、細(xì)菌、氮、磷、水生動(dòng)物、底泥、巖屑、藻類、金屬以及用戶自定義物質(zhì)[9]，CE-QUAL-ICM中的狀態(tài)變量也增加到36個(gè)。水質(zhì)模型中水質(zhì)組分的增加，在模擬不同區(qū)域的不同水質(zhì)狀況時(shí)，用戶可以通過(guò)選擇合適的水質(zhì)組分達(dá)到較好的模擬效果。

海洋水質(zhì)模型在內(nèi)部組分、數(shù)學(xué)方法、前后期處理等方面不斷發(fā)展，并且更多的生化過(guò)程考慮進(jìn)來(lái)，能夠更客觀的反映實(shí)際情況，也增加了模型的復(fù)雜性，其復(fù)雜性體現(xiàn)在以下幾個(gè)方面：水質(zhì)模型污染源負(fù)荷從點(diǎn)源、面源，到大氣污染源的加入，體現(xiàn)了不同發(fā)展階段污染源輸入的主要特征，增加了模型結(jié)構(gòu)的復(fù)雜性；基本的水動(dòng)力模型與水質(zhì)模型耦合的基礎(chǔ)上，更多可能影響水質(zhì)狀態(tài)變量的過(guò)程被考慮進(jìn)來(lái)，比如泥沙模塊、沉積模塊[47]、底泥模塊，增強(qiáng)了模型對(duì)生態(tài)系統(tǒng)的客觀描述，使模型本身也更加復(fù)雜；河口地區(qū)一維與多維模型的耦合[25]，建立可以適用于不同特征的相鄰區(qū)域的耦合模型，有利于計(jì)算效率和計(jì)算精度的提高。

海洋水質(zhì)模型的應(yīng)用日益廣泛，最主要應(yīng)用于水體中污染物行為的預(yù)測(cè)，隨著實(shí)際需求的增加，模型功能更加細(xì)化。海洋水質(zhì)模型應(yīng)用于水質(zhì)管理規(guī)劃與預(yù)測(cè)，可以合理規(guī)劃分配削減各污染源的排放量，力求在規(guī)定期限內(nèi)使水質(zhì)達(dá)標(biāo)，林衛(wèi)青等[21]在長(zhǎng)江口以及鄰近海域建立了二維水質(zhì)模型，計(jì)算得出總量分配方案，同時(shí)建立三維生態(tài)動(dòng)力學(xué)模型校驗(yàn)水質(zhì)模型計(jì)算出的總量控制方案；海洋水質(zhì)模型也應(yīng)用于計(jì)算水環(huán)境容量，是在確定的水質(zhì)目標(biāo)條件下，計(jì)算可以允許的最大排放量，喻良等[57]采用SMS模擬水質(zhì)，采用影響系數(shù)法計(jì)算環(huán)境容量，得出目標(biāo)污染負(fù)荷，計(jì)算方便易于理解；海洋水質(zhì)模型在水質(zhì)預(yù)警方面也廣泛應(yīng)用，曹立杰等[58]研發(fā)基于傳感器網(wǎng)絡(luò)的海洋環(huán)境監(jiān)測(cè)與預(yù)報(bào)系統(tǒng)，通過(guò)建立傳感器網(wǎng)絡(luò)，獲取實(shí)時(shí)監(jiān)測(cè)的水質(zhì)數(shù)據(jù)，利用水質(zhì)模型預(yù)測(cè)、預(yù)報(bào)富營(yíng)養(yǎng)化趨勢(shì)，為國(guó)內(nèi)近海富營(yíng)養(yǎng)化的監(jiān)測(cè)預(yù)報(bào)系統(tǒng)提供了參考。

4.2 水質(zhì)模型與新技術(shù)（新方法）結(jié)合

4.2.1 不確定性分析

水環(huán)境充滿很多不確定性因素,這種不確定性包括對(duì)水環(huán)境理化性質(zhì)的未知，也包括數(shù)據(jù)的不確定性，Takyi等[59]分析了優(yōu)化水質(zhì)模擬模型中不確定性的來(lái)源，可以概括為3個(gè)方面:（1）來(lái)自污染物排放量和河流背景值的隨機(jī)性；（2）河流與水質(zhì)資料的不充分對(duì)參數(shù)估計(jì)的影響；（3）對(duì)污染物輸移認(rèn)識(shí)不足，對(duì)簡(jiǎn)化水質(zhì)管理系統(tǒng)仍然缺乏理論支持。

不確定性水質(zhì)模型主要有3種[60]：基于概率論與數(shù)理統(tǒng)計(jì)的隨機(jī)方法（概率表述）、模糊數(shù)學(xué)（模糊度）、灰色系統(tǒng)理論（灰色度）。張榮祥等[61]從數(shù)據(jù)輸入輸出、參數(shù)識(shí)別、構(gòu)建模型等方面說(shuō)明了整個(gè)水質(zhì)模型構(gòu)建中不確定性的來(lái)源，其中分析了水質(zhì)模型構(gòu)建過(guò)程中采用的黑箱水質(zhì)模型和隨機(jī)理論水質(zhì)模型。扈華等[62]采用T-S模糊神經(jīng)網(wǎng)絡(luò)模型對(duì)烏梁素海水質(zhì)進(jìn)行合理的量化評(píng)價(jià)，該模型兼并模糊系統(tǒng)與人工神經(jīng)網(wǎng)絡(luò)的優(yōu)點(diǎn)，一定程度上解決了隨機(jī)性、模糊性的問題，結(jié)果合理且與實(shí)際環(huán)境狀況相符，可為監(jiān)管烏梁素海水域的提供理論依據(jù)。田一梅等[63]建立一維不確定性水質(zhì)模型（內(nèi)嵌神經(jīng)網(wǎng)絡(luò)），采用遺傳算法（改進(jìn)適應(yīng)度函數(shù)）優(yōu)化水質(zhì)模型的參數(shù)解，結(jié)果驗(yàn)證不確定性水質(zhì)模型具有很強(qiáng)的優(yōu)勢(shì)，體現(xiàn)在擬合精度更高、對(duì)污染物波動(dòng)更敏感，預(yù)測(cè)河流水質(zhì)的平均準(zhǔn)確度普遍高于確定性水質(zhì)模型。不確定性水質(zhì)模型比確定性水質(zhì)模型更能體現(xiàn)水質(zhì)的實(shí)際變化的情況，減少模式的主觀性，是將來(lái)發(fā)展一個(gè)大趨勢(shì)。

4.2.2 基于GIS的水質(zhì)模型

地理信息系統(tǒng)（Geographic Information System，GIS）應(yīng)用于水質(zhì)模型，可以在一定程度上彌補(bǔ)數(shù)學(xué)模型在管理維護(hù)數(shù)據(jù)、表現(xiàn)模擬結(jié)果以及空間分析上能力有限的缺點(diǎn)，不僅可以存儲(chǔ)管理大量的數(shù)據(jù)，還能對(duì)水質(zhì)模擬的結(jié)果進(jìn)行空間分析，為水質(zhì)模擬提供了預(yù)處理、后處理，使輸出輸出更加便捷[8]。江毓武等將GIS的地圖處理、數(shù)據(jù)分析功能應(yīng)用于海洋水質(zhì)模型的前期建立以及后期數(shù)據(jù)輸出中，縮短了建模周期，同時(shí)也提高了模型精度，并且使在地圖上查詢模擬結(jié)果成為現(xiàn)實(shí)[17]。Hartnett等[64]將DIVAST（DepthIntegratedVelocityandSoluteTransport）與GIS框架耦合，力求實(shí)現(xiàn)在特定領(lǐng)域中非專家人員對(duì)水質(zhì)模型的建立、執(zhí)行和水質(zhì)數(shù)據(jù)分析。基于GIS的水質(zhì)模型不僅在數(shù)據(jù)的管理上更加便利省時(shí)，還可以強(qiáng)化模型概念。

應(yīng)用GIS對(duì)模式數(shù)據(jù)的規(guī)范化管理，對(duì)結(jié)果進(jìn)行空間分析，是水質(zhì)模型發(fā)展的一個(gè)大方向，現(xiàn)階段將GIS與河流、水庫(kù)集成較多，結(jié)果相對(duì)于單獨(dú)的水質(zhì)模型結(jié)果要更優(yōu)，證實(shí)了水質(zhì)模型與GIS集成具有很好的應(yīng)用前景，為GIS在河口、海岸、海灣等海洋區(qū)域的應(yīng)用奠定了良好基礎(chǔ)。利用實(shí)時(shí)數(shù)據(jù)建立海洋監(jiān)測(cè)預(yù)警體系也是海洋水質(zhì)模型發(fā)展的一個(gè)重要方面，在沿海地區(qū)建立水質(zhì)監(jiān)測(cè)自動(dòng)化平臺(tái)，對(duì)海洋漁業(yè)、污染物治理、環(huán)境監(jiān)測(cè)有重要的指導(dǎo)意義。

5 結(jié)語(yǔ)

海洋水質(zhì)模型為水質(zhì)模擬、污染物治理、水資源規(guī)劃等方面提供了有力的工具，也為水質(zhì)預(yù)警工作提供了有力的技術(shù)支撐，隨著新技術(shù)的蓬勃發(fā)展，水質(zhì)模型在原有的基礎(chǔ)上推陳出新，適用性和模型的模擬能力也都有新的提高。海洋水質(zhì)模型在應(yīng)用中得到發(fā)展，使水質(zhì)模型有了更廣泛的應(yīng)用領(lǐng)域，建立綜合性更強(qiáng)、包含多種生物化學(xué)模塊的模式，有利于更加精確的描述海洋區(qū)域的生化過(guò)程，使模擬結(jié)果有更高的參考價(jià)值。

在應(yīng)用水質(zhì)模型的同時(shí)，也應(yīng)該注意水質(zhì)模型的發(fā)展的中面臨的問題：

（1）實(shí)時(shí)觀測(cè)數(shù)據(jù)缺乏。海洋水質(zhì)模型需要完整、連續(xù)、系統(tǒng)的水文、水質(zhì)資料，這是建立模型的基礎(chǔ)，觀測(cè)數(shù)據(jù)可以為模式提供初始場(chǎng)、邊界條件、強(qiáng)迫場(chǎng)，也可以用于驗(yàn)證模式的計(jì)算結(jié)果的準(zhǔn)確性，而現(xiàn)在走航式觀測(cè)和浮標(biāo)觀測(cè)，在時(shí)間上的連續(xù)性和數(shù)據(jù)的質(zhì)量方面還有待提高。建議不針對(duì)某個(gè)項(xiàng)目建立觀測(cè)點(diǎn)，在沿海選取代表性的位置布放浮標(biāo)，觀測(cè)變量包括基本物理變量和水質(zhì)狀態(tài)變量，數(shù)據(jù)自動(dòng)化傳輸處理，同時(shí)實(shí)行數(shù)據(jù)共享；

（2）水質(zhì)污染機(jī)理還不完全清楚，難以用數(shù)學(xué)公式詳盡的表達(dá)，且模式中存在諸多假設(shè)，難免與實(shí)際有差異。建議注重基本理論的研究，必要時(shí)進(jìn)行實(shí)驗(yàn)驗(yàn)證假設(shè)是否合理，與觀測(cè)數(shù)據(jù)結(jié)合，確認(rèn)假設(shè)成立的條件；

（3）近岸水體的水質(zhì)變化受眾多因素的影響，涵蓋了生物、物理、化學(xué)等不同學(xué)科，如何提高學(xué)科間的交叉、協(xié)作，以及基于水質(zhì)模型建立一個(gè)考慮多學(xué)科多介質(zhì)的海洋生態(tài)系統(tǒng)模型，更有利于海洋水質(zhì)的客觀模擬。

海洋水質(zhì)模型自產(chǎn)生以來(lái)得到廣泛應(yīng)用，通用性不斷增強(qiáng)，逐漸成為決策者的得力工具，基于海洋水質(zhì)模型，建立完善、系統(tǒng)的沿海水質(zhì)預(yù)警體系，預(yù)測(cè)、預(yù)報(bào)富營(yíng)養(yǎng)化的趨勢(shì)，同時(shí)提高海洋水質(zhì)模型模擬能力，是海洋水質(zhì)模型應(yīng)用的熱點(diǎn)。

[1]傅國(guó)偉,程聲通.水污染控制系統(tǒng)規(guī)劃[M].北京:清華大學(xué)出版社,1985:7-159

[2]王玲杰,孫世群,田豐.河流水質(zhì)模擬問題的探討[J].合肥工業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版）,2005,28(3):260-265.

[3]郭勁松,李勝海,龍騰銳.水質(zhì)模型及其應(yīng)用研究進(jìn)展[J].重慶建筑大學(xué)學(xué)報(bào),2002,24(2):109-115.

[4]Streeter H W,Phelps E B.A Study of the Pollution and Natural Purification of the Ohio River,Ⅲ.Factors Concerned in the Phenomena of Oxidation and Reaeration[M].Washington D C:US Government Printing Office,1925:1-73.

[5]李繼選,王軍.水環(huán)境數(shù)學(xué)模型研究進(jìn)展[J].水資源保護(hù),2006, 22(1):9-14.

[6]徐祖信,廖振良.水質(zhì)數(shù)學(xué)模型研究的發(fā)展階段與空間層次[J].上海環(huán)境科學(xué),2003,22(2):79-85.

[7]陳英.基于河流總氮水質(zhì)模型的研究進(jìn)展及發(fā)展趨勢(shì)[J].宜春學(xué)院學(xué)報(bào),2010,32(4):108-110.

[8]姜云超,南忠仁.水質(zhì)數(shù)學(xué)模型的研究進(jìn)展及存在的問題[J].蘭州大學(xué)學(xué)報(bào)（自然科學(xué)版）,2008,44(5):7-11.

[9]王斌,萬(wàn)天杰.河流水質(zhì)模型研究進(jìn)展[J].河南化工,2010,27(3): 2-3.

[10]樊石磊,呂松力,周志慧.水質(zhì)模型研究進(jìn)展[J].內(nèi)蒙古石油化工,2014,(21):32-33.

[11]周雪麗,孫森林,張寬義,等.水質(zhì)數(shù)學(xué)模型的研究進(jìn)展及其應(yīng)用[J].科學(xué)觀察,2011,(2):87-88.

[12]李程.水質(zhì)模型的研究進(jìn)展及發(fā)展趨勢(shì)[J].現(xiàn)代農(nóng)業(yè)科技, 2013,(6):208-209.

[13]李君,徐俊紅,曹永梅.河流水質(zhì)數(shù)學(xué)模型的研究進(jìn)展及發(fā)展趨勢(shì)[J].科技·經(jīng)濟(jì)·市場(chǎng),2007,(9):35-36.

[14]葉常明,顏文紅.多介質(zhì)環(huán)境循環(huán)模型的研究進(jìn)展[J].環(huán)境科學(xué)進(jìn)展,1994,2(1):9-25.

[15]Mackay D,MacLeod M.Multimedia Environmental Models:the Fugacity Approach[M].The Practice Periodical of Hazardous, Toxic,and Radioactive Waste Management,2002,6(2):63-69.

[16]陳美丹,姚琪,徐愛蘭.WASP水質(zhì)模型及其研究進(jìn)展[J].水利科技與經(jīng)濟(jì),2006,12(7):420-422,426.

[17]江毓武,洪華生,張珞平.地理信息系統(tǒng)(GIS)在廈門海域水質(zhì)模型中的應(yīng)用[J].廈門大學(xué)學(xué)報(bào)(自然版),1999,(1):90-95.

[18]倪晉仁,秦華鵬,趙智杰.基于水質(zhì)模型的海灣填海岸線選擇[J].環(huán)境科學(xué)學(xué)報(bào),2001,21(6):684-688.

[19]王昆,王年斌,孫家文,等.基于遼東灣海域的非保守型水質(zhì)預(yù)測(cè)預(yù)報(bào)模型研究[J].水動(dòng)力學(xué)研究與進(jìn)展,2013,28(6):691-700.

[20]沈永明,鄭永紅,吳修廣.近岸海域污染物遷移轉(zhuǎn)化的三維水質(zhì)動(dòng)力學(xué)模型[J].自然科學(xué)進(jìn)展,2004,14(6):694-699.

[21]林衛(wèi)青,盧士強(qiáng),矯吉珍.長(zhǎng)江口及毗鄰海域水質(zhì)和生態(tài)動(dòng)力學(xué)模型與應(yīng)用研究[J].水動(dòng)力學(xué)研究與進(jìn)展,2008,23(5):522-531.

[22]Li K Q,Zhang L,Li Y,et al.A Three-Dimensional Water Quality Model to Evaluate the Environmental Capacity of Nitrogen and Phosphorus in Jiaozhou Bay,China[J].Marine Pollution Bulletin, 2015,91(1):306-316.

[23]Zhang Y,Sun Y L,Yu J,et al.A Three-Dimensional Water Quality Model and Its Application to Jiaozhou Bay,China[J].China Ocean Engineering,2012,26(4):669-684.

[24]Khangaonkar T,Sackmann B S,Long W,et al.Puget Sound DissolvedOxygenModelingStudy:Developmentofan Intermediate Scale Water Quality Model[R].Technical Report AC05-76RL01830,2012.

[25]胡嘉鏜,李適宇,裴木鳳,等.珠江三角洲一維河網(wǎng)與三維河口耦合水質(zhì)模型模擬與驗(yàn)證[J].海洋與湖沼,2012,43(1):1-9.

[26]譙程駿,張東輝,張敏.長(zhǎng)江水質(zhì)的評(píng)價(jià)預(yù)測(cè)模型[J].工程數(shù)學(xué)學(xué)報(bào),2005,22(7):41-46.

[27]周新明.一維河流水質(zhì)模型在應(yīng)急監(jiān)測(cè)中的應(yīng)用分析[J].綠色科技,2011,(6):10-11.

[28]楊紅,丁駿,王春峰,等.象山港圍隔生態(tài)系水質(zhì)模型研究[J].海洋科學(xué),2012,36(7):14-22.

[29]Babu M T,Das V K,Vethamony P.BOD-DO Modeling and Water Quality Analysis of a Waste Water Outfall off Kochi,West Coast of India[J].Environment International,2006,32(2):165-173.

[30]朱仲平.水質(zhì)規(guī)劃的建立[J].環(huán)境科學(xué),1981,2(3):56-63.

[31]Kazmi A A,Hansen I S.Numerical Models in Water Quality Management:A Case Study for the Yamuna River(India)[J]. Water Science and Technology,1997,36(5):193-200.

[32]郭安.長(zhǎng)江口水動(dòng)力及污水?dāng)U散影響的數(shù)值模擬[D].南京:河海大學(xué),2007.

[33]陳躍燕.SMS在鎮(zhèn)江內(nèi)江水質(zhì)模擬與預(yù)測(cè)中的應(yīng)用[D].鎮(zhèn)江:江蘇大學(xué),2007.

[34]馮媛媛.甌江河口灘涂圍墾對(duì)水動(dòng)力及水質(zhì)的影響研究[D].南京:河海大學(xué),2008.

[35]Cerco C F,Cole T.Three-Dimensional Eutrophication Model of Chesapeake Bay[J].Journal of Environmental Engineering,1993, 119(6):1006-1025.

[36]Kim T,Labiosa R G,Khangaonkar T,et al.Development and Evaluation of A Coupled Hydrodynamic(FVCOM)and Water Quality Model(CE-QUAL-ICM)[J].American Society of Civil Engineers,2009,388:373-388.

[37]Kim C S,Lim H S,Cerco C F.Three-Dimensional Water Quality Modelling for Tidal Lake and Coastal Waters with ROMS-ICM [J].Journal of Coastal Research,2011,SI 64:1068-1072.[38] Dalyander P S,Cerco C F.Integration of An Individual-Based Fish Bioenergetics Model into A Spatially Explicit Water Quality Model(CE-QUAL-ICM)[R].Technical NoteADA518934,2010.

[39]Di Toro D M,Fitzpatrick J J,Thomann R V,et al.Documentation for Water Quality Analysis Simulation Program(WASP)and Model VerificationProgram(MVP)[R].Duluth,MN:US Environmental ProtectionAgency,1983.

[40]Ambrose R B,Wool T A,Connolly J P，et al.WASP4,A Hydrodynamic and Water Quality Model-Model Theory,User's Manual,andProgrammer'sGuide[R].Athens,GA:US Environmental ProtectionAgency,1988.

[41]Ambrose R B,Wool T A,Martin J L,et al.WASP5.X,A Hydrodynamic and Water Quality Model-Model Theory,User's Manual,and Programmer's Guide[M].Draft,Athens,GA:US Environmental Protection Agency,1991.[42]Wool T A,Ambrose R B,Martin J L,et al.Water Quality Analysis Simulation Program(WASP),Version6.0-DRAFT:User'sManual[R]. Atlanta,GA:US Environmental ProtectionAgency,2001.

[43]孫學(xué)成,鄧曉龍,張彩香,等.WASP6系統(tǒng)在三峽庫(kù)區(qū)水質(zhì)仿真中的應(yīng)用[J].三峽大學(xué)學(xué)報(bào)(自然科學(xué)版),2003,25(2):185-188.

[44]Hao L F,Cheng L B,Bi X W.Apply With WASP Water Quality Model[C]//Proceedings of the 2nd International Conference on Computer Application and System Modeling.Paris,France: Atlantis Press,2012.

[45]Zheng L Y,Chen C S,Zhang F Y.Development of Water Quality Model in the Satilla River Estuary,Georgia[J].Ecological Modelling,2004,178(3-4):457-482.

[46]李林子,錢瑜,張玉超.基于EFDC和WASP模型的突發(fā)水污染事故影響的預(yù)測(cè)預(yù)警[J].長(zhǎng)江流域資源與環(huán)境,2011,20(8): 1010-1016.

[47]Xiong Y.Coupling Sediment Transport and Water Quality Models [D].Mississippi:Mississippi State University,2010.

[48]馬騰,劉文洪,宋策,等.基于MIKE3的水庫(kù)水溫結(jié)構(gòu)模擬研究[J].電網(wǎng)與清潔能源,2009,25(2):68-71.

[49]王勇智,吳頔,石洪華,等.近十年來(lái)渤海灣圍填海工程對(duì)渤海灣水交換的影響[J].海洋與湖沼,2015,46(3):471-480.

[50]黃少彬,李開明,姜國(guó)強(qiáng),等.基于MIKE3模型的珠江口水體交換研究[J].環(huán)境科學(xué)與管理,2013,38(8):134-140.

[51]陳瑤泓伶,袁洪濤.天津南港污水排海工程對(duì)海域環(huán)境影響的數(shù)值模擬[J].水道港口,2014,35(5):539-544.

[52]薛立梅.SMS系統(tǒng)軟件在水利工程中的應(yīng)用[J].東北水利水電, 2010,28(11):63,68.

[53]張明進(jìn),張華慶.SMS水動(dòng)力學(xué)軟件[J].水道港口,2006,27(1):57-59.

[54]吳志易,徐項(xiàng)煜.基于SMS的長(zhǎng)江口—杭州灣海域建模[J].中國(guó)水運(yùn),2013,13(12):125-126.

[55]陳志娟,拾兵,韓艷.SMS在黃河口水流數(shù)值模擬中的應(yīng)用[J].人民黃河,2008,30(8):47-49.

[56]魯海燕.SMS模型在杭州灣潮流模擬中的應(yīng)用[J].浙江水利科技,2003,(3):4-6.

[57]喻良,劉遂慶,王牧陽(yáng).基于水環(huán)境模型的水環(huán)境容量計(jì)算的研究[J].河南科學(xué),2006,24(6):874-876.

[58]曹立杰,郭戈,靳玉峰,等.基于傳感器網(wǎng)絡(luò)的海洋水質(zhì)監(jiān)測(cè)及赤潮預(yù)報(bào)系統(tǒng)的設(shè)計(jì)[J].大連海洋大學(xué)學(xué)報(bào),2014,29(6):664-668.

[59]Takyi A K,Lence B J.Surface Water Quality Management Using A Multiple-Realization Chance Constraint Method[J].Water Resources Research,1999,35(5):1657-1670.

[60]金明.一維穩(wěn)態(tài)河流水質(zhì)的隨機(jī)微分方程模型[J].水利學(xué)報(bào), 1991,(2):19-25.

[61]張榮祥,劉鍵.不確定性方法應(yīng)用于水質(zhì)模型中的研究[J].環(huán)境與生活,2014,(10):57,59.

[62]扈華,王冬青.T-S模糊神經(jīng)網(wǎng)絡(luò)在水質(zhì)評(píng)價(jià)中的應(yīng)用研究[J].內(nèi)蒙古農(nóng)業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,36(4):128-132.

[63]田一梅,劉揚(yáng),王彬蔚.不確定水質(zhì)模型在城市河流水質(zhì)模擬中的應(yīng)用[J].土木建筑與環(huán)境工程,2011,33(3):119-123.

[64]Hartnett M,Nash S,Mills P.The Development of a GIS-Based Water Quality Model[J].River Basin Management II,2003,60: 401-409.

Review on the marine water quality model

YU Han，LIU Gui-mei，YANG Jing

（National Marine Environmental Forecasting Center,Beijing 100081 China）

In this paper,the development stages and modeling steps of the water quality model,are reviewed. The applications of the marine water quality model are introduced from one-dimensional,two-dimensional and three-dimensional perspectives and the future development trend of marine water quality model is prospected.It is concluded that the marine water quality model can provide the systematic reference for the emission of pollutant,monitoring the water quality and managing the water resources.

marine water quality model;research progress;development

X834

A

1003-0239（2017）02-0088-09

10.11737/j.issn.1003-0239.2017.02.012

2016-08-19；

2016-11-28。

國(guó)家自然科學(xué)基金（41222038，41206023）；中國(guó)科學(xué)院戰(zhàn)略性先導(dǎo)科技專項(xiàng)(XDA11020104)；國(guó)家海洋局海洋公益性行業(yè)專項(xiàng)（01205018）。

于寒（1991-），女，碩士在讀，從事海洋生態(tài)動(dòng)力學(xué)模型研究。E-mail:yuhan3297@163.com