河口海岸物質(zhì)輸運(yùn)時間尺度研究進(jìn)展

張 琴,楊 潔,陶建峰,張長寬,張 曉

(1.上海勘測設(shè)計研究院有限公司,上海 200434; 2.河海大學(xué)港口海岸與近海工程學(xué)院,江蘇 南京 210098)

河口海岸物質(zhì)輸運(yùn)時間尺度研究進(jìn)展

張 琴1,楊 潔2,陶建峰2,張長寬2,張 曉2

(1.上?？睖y設(shè)計研究院有限公司,上海 200434; 2.河海大學(xué)港口海岸與近海工程學(xué)院,江蘇 南京 210098)

從物質(zhì)輸運(yùn)時間概念、研究方法及應(yīng)用方面綜述了近年來河口海岸物質(zhì)輸運(yùn)時間尺度的研究進(jìn)展,論述了目前物質(zhì)輸運(yùn)時間尺度主要數(shù)學(xué)模型的優(yōu)點(diǎn)和局限性,著重討論了基于物質(zhì)組成的年齡和存留時間(CART)理論的研究進(jìn)展。針對當(dāng)前非保守性物質(zhì)輸運(yùn)時間尺度研究中存在的不足,指出易被泥沙吸附的污染物的輸運(yùn)時間尺度需要考慮其與泥沙的相互作用過程,加強(qiáng)非保守性物質(zhì)、泥沙和易被泥沙吸附的污染物輸運(yùn)時間尺度研究有助于全面了解復(fù)雜物質(zhì)的輸運(yùn)規(guī)律;CART理論應(yīng)拓展應(yīng)用于計算非保守性物質(zhì)、泥沙和考慮泥沙影響的污染物輸運(yùn)時間尺度的研究。

物質(zhì)輸運(yùn);污染物；泥沙；時間尺度;河口;海岸

在河口海岸生態(tài)系統(tǒng)中,流體介質(zhì)攜帶了大量的浮游生物、污染物、溶解相氣體、懸浮顆粒等物質(zhì),它們的輸運(yùn)過程不僅影響地形地貌,而且也顯著影響到水體系統(tǒng)的水質(zhì)和生態(tài)[1-2],因此研究物質(zhì)輸運(yùn)過程對河口海岸生態(tài)系統(tǒng)具有重要意義。水體中的這些物質(zhì),主要在水動力作用下進(jìn)行輸移擴(kuò)散,為了解其輸運(yùn)過程,針對不同的問題,在研究方法上從余流和物質(zhì)輸運(yùn)時間尺度方面開展研究,在數(shù)值手段上分別基于歐拉法和拉格朗日法開展研究。物質(zhì)輸運(yùn)時間尺度是評價水體環(huán)境質(zhì)量和環(huán)境容量的重要指標(biāo),海灣和近海海域的污染物通過對流擴(kuò)散等物理過程輸運(yùn)至外海海域,使污染物濃度降低,實(shí)現(xiàn)水體的自凈。不同區(qū)域的水體由于其水動力強(qiáng)弱不同,其自凈能力也就不同,因此物質(zhì)輸運(yùn)時間尺度可用來分析排入受納水體中污染物的輸運(yùn)過程,定量反映水體的自凈能力,并給出合理的生態(tài)保護(hù)建議,對水體系統(tǒng)的生態(tài)環(huán)境具有重要意義[1]。

目前對于物質(zhì)輸運(yùn)時間尺度的研究多集中于水體或溶解性物質(zhì),對非保守性物質(zhì)、泥沙或易被泥沙吸附的污染物研究較少。本文主要從物質(zhì)輸運(yùn)時間尺度概念、計算方法及應(yīng)用狀況等方面對河口海岸物質(zhì)輸運(yùn)時間尺度的研究情況進(jìn)行綜述,并對其未來的研究方向作了展望。

1 物質(zhì)輸運(yùn)時間尺度概念

從20世紀(jì)70年代開始,根據(jù)研究對象和問題的不同,不同研究者定義了不同的時間尺度概念來定量描述物質(zhì)輸運(yùn)過程,如Bolin等[3]在1973年根據(jù)海水交換率的定義引入沖淡時間(flushing time)的概念,即水體全部被置換的時間,主要用來表征河口或海灣的水體被更新的速率,反映的是研究水體整體的運(yùn)動特征。一些學(xué)者認(rèn)為研究區(qū)域內(nèi)物質(zhì)的衰減并非線性衰減,而是呈指數(shù)形式,于是提出轉(zhuǎn)運(yùn)時間(turn-over time)[4]的概念,即研究區(qū)域內(nèi)物質(zhì)總量減少至初始值的e-1(約0.37)所需的時間。此外,Bolin等[3]還引入了年齡(age)概念來描述質(zhì)點(diǎn)自進(jìn)入控制區(qū)域到某點(diǎn)所經(jīng)歷的時間。Luff等[5]提出了半交換時間(half-life time)的概念,其定義為研究海域內(nèi)濃度降低至初始值一半時所需的時間,是使用稀釋速率的快慢來反映海域的水體更新能力。1976年,Zimmerman[6]在年齡的概念基礎(chǔ)上定義了存留時間(residence time)的概念,將空間某點(diǎn)的存留時間定義為該點(diǎn)釋放的物質(zhì)第一次離開控制區(qū)域所需的時間。Zimmerman認(rèn)為在同一位置釋放后不同質(zhì)點(diǎn)移動路徑不同,存留時間依賴于質(zhì)點(diǎn)位置和釋放時間,需要對不同質(zhì)點(diǎn)存留時間進(jìn)行平均來衡量這一空間位置的水交換能力。在此基礎(chǔ)上,Takeoka[7]提出了剩余函數(shù)的概念(remnant function),并建立了剩余函數(shù)和平均存留時間之間的關(guān)系,剩余函數(shù)可以反映控制區(qū)域內(nèi)所研究物質(zhì)含量的下降,因此利用在某時刻剩余在控制區(qū)域內(nèi)物質(zhì)的比例可得到物質(zhì)的存留時間。

存留時間在一定程度上反映了污染物對一個特定水體生態(tài)系統(tǒng)的健康所造成的影響大小,是一個重要的物理變量,它表示物質(zhì)在水體中生命的剩余時間,與年齡互為補(bǔ)充[1,5]。物質(zhì)在離開控制區(qū)域后,由于受潮汐及擴(kuò)散的影響,可能再次甚至多次返回控制區(qū)域,因此在一些考慮生物和地球化學(xué)過程的研究中,物質(zhì)在控制區(qū)域內(nèi)存留的總時間顯得更為重要,如受潮汐影響較大的浮游植物,其生長與其在控制區(qū)域存留的總時間密切相關(guān),因此定義暴露時間為物質(zhì)在控制區(qū)域內(nèi)存留的總時間[1]。在衡量一個海域受污染物質(zhì)影響的時間時,如果這種污染物質(zhì)為保守物質(zhì),并且在控制區(qū)域的邊界處動力條件或生物化學(xué)條件沒有劇烈變化時,就需要考慮其再次進(jìn)入控制海域的時間；如果離開控制區(qū)域后釋放的物質(zhì)的特性會發(fā)生劇烈變化,那么存留時間的概念則更適用[8]。暴露時間在數(shù)值上比存留時間大,它們之間的差值與近海海域的環(huán)流有關(guān)[9],例如在強(qiáng)潮河口物質(zhì)更容易在潮汐作用下再次返回控制區(qū)域,使得強(qiáng)潮河口物質(zhì)的暴露時間與存留時間差異較大。

以上研究表明,物質(zhì)輸運(yùn)時間尺度可以從系統(tǒng)整體考慮,如沖淡時間、轉(zhuǎn)運(yùn)時間、半交換時間,也可以從局部出發(fā),如物質(zhì)在空間各點(diǎn)的存留時間、暴露時間和年齡等。對于復(fù)雜海域的物質(zhì)輸運(yùn),更需要反映物質(zhì)輸運(yùn)時間尺度的時空分布特征,因此,年齡、存留時間、暴露時間的概念更適用于實(shí)際海域中反映物質(zhì)輸運(yùn)過程的時空變化特征。迄今為止,沒有一種輸運(yùn)時間尺度可以較好地完整刻畫物質(zhì)的整個輸運(yùn)過程[2],在使用過程中應(yīng)謹(jǐn)慎對待,針對不同的研究問題選擇合適的輸運(yùn)時間尺度。

2 物質(zhì)輸運(yùn)時間尺度計算方法及應(yīng)用

隨著水動力數(shù)值模型及物質(zhì)輸運(yùn)模型的發(fā)展,利用各種數(shù)學(xué)模型來進(jìn)行物質(zhì)輸運(yùn)時間尺度的計算越來越多。目前常用的計算物質(zhì)輸運(yùn)時間尺度的方法有箱式模型(box model)法、拉格朗日粒子追蹤法(Lagrange particle tracking method)、對流擴(kuò)散模型法和基于物質(zhì)組成的年齡和存留時間(constituent-oriented age and residence time,CART)理論等[10-11]。

2.1 箱式模型法

箱式模型法以物質(zhì)守恒為基礎(chǔ),利用物質(zhì)通量與水體體積的關(guān)系得到海水交換速率,實(shí)際應(yīng)用中可分為單箱模型和多箱模型,將研究水體作為一個整體或者分區(qū)研究。典型的箱式模型法有水體體積和水體通量比法、納潮量法(tidal prism method)和淡水組分法(freshwater fraction method)。

沖淡時間最初利用控制區(qū)域內(nèi)水體的總體積與水體通量的比值或者是物質(zhì)總量與物質(zhì)通量的比值來計算,但是水體體積和通量通常很難計算,于是出現(xiàn)了不同的沖淡時間的計算方法。如Dyer等[12]提出納潮量法,其原理是假設(shè)漲潮期間進(jìn)入河口的水體與河口內(nèi)原有水體完全混合,河口內(nèi)的原有水體隨落潮流出后不再返回控制區(qū)域,但該方法不考慮徑流影響。Fischer等[13]提出淡水組分法,考慮淡水流量后,沖淡時間可以表示為水體的總體積與潮流、徑流通量之和的比值。根據(jù)Knudsen提出的鹽度平衡理論[9],潮流和徑流通量之比可用海水鹽度和河口內(nèi)水體的鹽度來反映,于是利用不同的鹽度平衡理論可以得到不同的淡水組分法的表達(dá)形式[9]。淡水組分法通過隨時間變化的海水鹽度反映潮流邊界的水交換情況,但是需要詳盡的水流和鹽度資料。

在箱式模型的理論基礎(chǔ)上,國內(nèi)外眾多學(xué)者對不同河口海灣的水交換時間進(jìn)行了研究。高抒等[14]根據(jù)狹長形海灣的特點(diǎn),將狹長形海灣分成若干個區(qū)域,建立多箱物理模型研究了象山港的水交換時間,得出象山港灣頂80%水體被灣外水體替換所需時間長達(dá)1 a,而半交換時間需近0.5 a。潘偉然[15]將湄洲灣作為一個整體,采用單箱模型分別計算了灣內(nèi)海水的平均交換率和平均半交換時間,同時指出單箱模型由于采用灣外海水流入灣內(nèi)瞬時即與整個灣內(nèi)海水充分混合的假設(shè),與實(shí)際情況有一定的差距。箱式模型在估算時把海灣假設(shè)成一個或多個箱體，與灣外部水體在一個潮周期內(nèi)充分交換,并假定控制區(qū)域外海水進(jìn)入控制區(qū)域內(nèi)立即與整個計算區(qū)域內(nèi)水體完全混合,而實(shí)際上進(jìn)入箱內(nèi)的水體不能立即與原有水體發(fā)生混合,因此計算出來的物質(zhì)輸運(yùn)時間往往偏小[16]。此外,海區(qū)不同位置存在動力差異,水體交換并不是整體簡單的進(jìn)退,而是有先后,甚至有些地方小海灣可能出現(xiàn)滯留現(xiàn)象,輸運(yùn)時間延長,交換能力也相對較差,這類方法不能區(qū)分海灣內(nèi)不同區(qū)域水體的輸運(yùn)時間,不能反映水體或者其他物質(zhì)的輸運(yùn)特征,尤其是對于地形復(fù)雜的海灣[16]。

2.2 拉格朗日粒子追蹤法

拉格朗日粒子追蹤法在利用數(shù)值模擬或?qū)崪y手段得到流場的基礎(chǔ)上,在研究區(qū)域內(nèi)布置用于代表水體的粒子,通過計算粒子的運(yùn)動來研究水體的流通和交換情況。通過在某一瞬時釋放大量粒子,計算其濃度降為一半或e-1所需的時間,得到半交換時間或轉(zhuǎn)運(yùn)時間;通過統(tǒng)計流出研究區(qū)域邊界的粒子數(shù)目,除以初始粒子數(shù)目計算出研究區(qū)域的水交換率;還可以統(tǒng)計每個粒子在水體中的流動時間，計算水體在研究海域的存留時間或年齡。Anabela等[17]曾通過追蹤大量粒子,統(tǒng)計分析得到葡萄牙Tagus河口的水體平均存留時間,發(fā)現(xiàn)該河口的水體存留時間存在明顯的時空變化,并且指出拉格朗日粒子追蹤法需要足夠高的精度才能避免產(chǎn)生嚴(yán)重的錯誤。Bilgili等[18]利用拉格朗日粒子追蹤法計算了新罕布什爾州大灣河口的存留時間,并對比了有無疊加隨機(jī)游走模式的計算結(jié)果,結(jié)果表明即使很小的擴(kuò)散率也會導(dǎo)致粒子的存留時間減少,指出計算存留時間需考慮擴(kuò)散過程。Liu等[19]采用拉格朗日粒子追蹤法研究了臺灣淡水河口三維水體中的存留時間和年齡,結(jié)果表明存留時間和徑流量存在指數(shù)關(guān)系,同時密度流對滯留時間和水體年齡影響顯著。沈永明等[20]將大伙房水庫劃分為7個水區(qū),利用拉格朗日粒子追蹤法分別計算了各個水區(qū)的存留時間,這比單純的計算整個水區(qū)的平均存留時間更符合實(shí)際。Sandery等[21]通過追蹤大量粒子模擬了Bass Strait海域真實(shí)驅(qū)動力作用下的沖淡時間、存留時間和水體年齡的分布,結(jié)果表明該海域的物質(zhì)輸運(yùn)時間尺度存在季節(jié)性變化,同時指出未來可研究水體生態(tài)系統(tǒng)中的各種物質(zhì)的輸運(yùn)時間尺度,如溶解性氣體、懸浮顆粒物等,可有助于分析其輸運(yùn)特征。拉格朗日粒子追蹤法雖然簡單直觀,且能反映出粒子的運(yùn)動軌跡,通過分區(qū)釋放粒子也可以估算研究區(qū)域的輸運(yùn)時間尺度的空間分布特征,但要得到輸運(yùn)時間尺度更細(xì)致的時空分布特征無疑要釋放大量的粒子,計算量很大;另外,盡管拉格朗日粒子追蹤法能通過添加隨機(jī)游走模塊模擬擴(kuò)散過程,但在渦黏系數(shù)存在很大的空間變化時,這種模擬方法不能夠直觀地表達(dá)擴(kuò)散過程[22]。同時,若考慮非定常流動時,需考慮不同初始釋放時間對計算結(jié)果的影響。

2.3 對流擴(kuò)散模型法

對流擴(kuò)散模型法在數(shù)值模擬得到流場的基礎(chǔ)上求解對流擴(kuò)散方程,從而得到示蹤物的濃度分布,再通過示蹤物濃度計算海水或者其他物質(zhì)的半交換時間、平均存留時間等。彭輝等[10]利用三維潮流數(shù)學(xué)模型,以溶解態(tài)的保守性物質(zhì)為示蹤劑,建立對流擴(kuò)散型的海灣水交換數(shù)值模型,通過計算示蹤劑濃度,計算了象山港水體半交換時間和平均存留時間,指出對流擴(kuò)散模型計算結(jié)果與箱式模型得到的結(jié)果差異較大,箱式模型可能會高估海水的交換能力。劉哲等[11]在Tarkeoka提出的剩余函數(shù)概念的基礎(chǔ)上,通過對流擴(kuò)散模型得到了膠州灣水體平均存留時間的空間分布,并通過追蹤粒子運(yùn)動驗(yàn)證了對流擴(kuò)散模型計算平均存留時間的正確性。Braunschwei等[23]利用高精度二維水動力數(shù)值模型耦合示蹤物輸運(yùn)模型研究了地形對物質(zhì)存留時間的影響;也有一些學(xué)者將研究區(qū)域分成多個子區(qū)域,在剩余函數(shù)和年齡概念定義的基礎(chǔ)上,通過對流擴(kuò)散示蹤物濃度數(shù)值模型計算多個粒子的濃度分布,分區(qū)計算得到水體存留時間和年齡[2,20，24]。

基于歐拉方法的對流擴(kuò)散模型計算得到的輸運(yùn)時間尺度仍是整體的概念,雖然可以通過分區(qū)分別計算各個區(qū)域的輸運(yùn)時間,但是得不到輸運(yùn)時間的空間分布特征。

2.4 CART理論

CART理論是Delhez和Deleersnijder建立的一種基于歐拉方法的計算水體年齡和存留時間的基本理論。根據(jù)Bolin等[3]對年齡的定義,不考慮源匯項(xiàng),CART理論中年齡的計算方法是首先求解物質(zhì)濃度方程和年齡濃度的對流擴(kuò)散方程,然后將空間各點(diǎn)的年齡濃度除以物質(zhì)濃度值,得到物質(zhì)在空間各點(diǎn)的平均年齡[25]。

CART理論通過求解兩個對流擴(kuò)散方程,可以得到任意時間任意網(wǎng)格節(jié)點(diǎn)上的年齡值,即可以得到物質(zhì)輸運(yùn)年齡的時空分布。由于同一位置釋放的不同粒子移動路徑不同,因此不同粒子從同一釋放點(diǎn)到達(dá)同一位置所需的時間也不盡相同,即任意位置的粒子團(tuán)包含了具有不同年齡的粒子,當(dāng)不需要給出每個粒子的年齡時,可采用平均年齡來確定整個粒子團(tuán)的年齡。CART理論最關(guān)鍵之處在于引入平均年齡的假定[25],即假定在空間某點(diǎn)的一組粒子的平均年齡是各個粒子年齡的質(zhì)量加權(quán)平均數(shù),表達(dá)成歐拉形式即為年齡濃度與物質(zhì)濃度的比值,因此,CART理論計算的年齡實(shí)際上為平均年齡。

CART理論中存留時間的計算方法是通過引入伴隨變量C*(T,x),假定伴隨變量滿足一定條件,利用污染物濃度方程,通過數(shù)學(xué)推導(dǎo)得到平均存留時間的控制方程,加上存留時間的邊界條件,通過一次數(shù)值計算就可以得到存留時間的空間分布[26],該方法也稱為伴隨方法。但伴隨方法在時間上要進(jìn)行逆向計算,即從后一時刻往前一時刻積分,計算中時間步長為負(fù)數(shù),向后積分的物理意義在于t時刻x位置處釋放的物質(zhì)的存留時間取決于[t,T]這一時間間隔內(nèi)的動力作用，其中T為向后積分時的初始時刻,對應(yīng)于數(shù)值模擬的最后時刻[26]。該方法雖通過一次計算就可得到存留時間的時空分布,但存在一個缺點(diǎn)就是它需要逆向計算,這使得該方法較難實(shí)現(xiàn)[27]。伴隨方法需要使用預(yù)先計算好的流場數(shù)據(jù),同時,進(jìn)行逆向計算前,需對研究區(qū)域水體的傳輸時間尺度進(jìn)行預(yù)估,以保證足夠的計算時長,因而流場刻畫的精細(xì)程度需要在流場存儲空間與計算精度之間作平衡,選擇較為合適的時間步長及計算時長輸出流場數(shù)據(jù)。

隨后,Delhez[28-29]在平均存留時間的理論基礎(chǔ)上提出了暴露時間的計算方法,在邊界處理上考慮物質(zhì)的再返回過程,且暴露時間的計算域要比研究區(qū)域大(為了反映物質(zhì)的再返回過程),從而得到物質(zhì)在控制區(qū)域內(nèi)存留的總時間,即暴露時間。

綜上所述,從計算方法上看,物質(zhì)輸運(yùn)時間尺度計算方法主要可以分為歐拉方法(如CART理論)和拉格朗日方法(如拉格朗日粒子追蹤法)。拉格朗日法簡單直觀,不僅能得到各個粒子的輸運(yùn)時間,還能得到粒子的輸運(yùn)軌跡,但是若要得到輸運(yùn)時間的時空分布則需要大量的計算;歐拉法能較簡便地得到輸運(yùn)時間的時空分布特征,但得到的是年齡和存留時間等的平均值。各方法均基于一定的假設(shè),各有優(yōu)缺點(diǎn),在具體的運(yùn)用過程中應(yīng)根據(jù)研究需要、水體特征、適用范圍等具體條件加以選擇。

3 CART理論研究進(jìn)展

CART理論因其簡單方便,且能反映輸運(yùn)時間的時空變化特征,被廣泛應(yīng)用于河口、海灣、湖泊以及海峽的物質(zhì)輸運(yùn)研究中,如在分析水體或保守性物質(zhì)輸運(yùn)時間的影響因素和分布特征等方面。Shen等[30]曾利用CART理論中的年齡計算方法研究了美國York River 河口的水流輸運(yùn)年齡,發(fā)現(xiàn)在水平和垂直方向上,河口的年齡分布均與鹽度的分布類似,且與徑流流量有密切關(guān)系。Shen等[31]通過CART理論中年齡的計算方法,研究了Chesapeake海灣的溶解性物質(zhì)的運(yùn)輸過程,得出徑流、密度環(huán)流、風(fēng)應(yīng)力是影響該海灣物質(zhì)輸運(yùn)的重要因素。de Brye等[32]曾利用CART 理論研究了Scheldt河口的年齡、存留時間和暴露時間,表明潮汐和徑流是影響水體年齡和平均存留時間的重要因素,存留時間和暴露時間受物質(zhì)釋放時刻影響,不同釋放時刻引起的差異可達(dá)75%,并指出河口的物質(zhì)再返回控制區(qū)域的可能性較大,需要考慮物質(zhì)的再返回過程。K?rn?等[33]模擬了Columbia River河口的水體年齡分布,指出水體年齡除受徑流影響顯著外,也受主通道兩側(cè)環(huán)流影響,兩側(cè)環(huán)流較弱,導(dǎo)致水流主通道兩側(cè)水體年齡較大。張琴等[34]基于CART理論模擬了浙江臺州灣海域不同流量條件下的淡水輸運(yùn)年齡,并探討了岸線變化對淡水年齡的影響。Andutta等[27]運(yùn)用對流擴(kuò)散模型得到中潮河口的存留時間和暴露時間的分布特征和影響因素,并利用CART理論計算得到回復(fù)系數(shù)(暴露時間減去存留時間與存留時間的比值),指出CART理論能很好地計算出水體再返回的頻率。蔡忠亞[8]將該理論應(yīng)用于研究渤海灣水體平均存留時間及其季節(jié)變化上,并分析了潮流、風(fēng)、密度流對渤海水體平均存留時間的影響。這些研究均以淡水水體或溶解性物質(zhì)為研究對象,假設(shè)物質(zhì)只受水體對流擴(kuò)散的影響,忽略其源匯項(xiàng)。對于其他物質(zhì),如非保守性物質(zhì)、泥沙和易被泥沙吸附的污染物等物質(zhì)的輸運(yùn)時間尺度缺乏研究。

從CART理論中年齡和存留時間的計算方程推導(dǎo)過程可以看出,如果在污染物濃度方程中加入源匯項(xiàng),那么年齡和存留時間的計算方程也將發(fā)生改變,當(dāng)源匯項(xiàng)滿足一定條件時,可推導(dǎo)得到考慮源匯項(xiàng)后物質(zhì)的輸運(yùn)時間計算方程。因此,CART理論中年齡的計算方法在2007年首次被拓展應(yīng)用于研究泥沙的輸運(yùn)年齡[35]。Mercier等[35]基于CART理論中年齡的計算方法對比利時淺海海域的三維懸沙輸運(yùn)時間進(jìn)行了研究,定義了兩種泥沙運(yùn)輸時間:一是“再懸浮年齡”,指泥沙自床面掀起后在水中所經(jīng)歷的時間(不包括泥沙在床面所滯留的時間);二是“輸運(yùn)年齡”,指泥沙自進(jìn)入控制區(qū)域到空間某點(diǎn)所需的時間(包括泥沙在床面所滯留的時間)。與以往水體年齡或保守性物質(zhì)的輸運(yùn)年齡不同的是,該研究根據(jù)懸沙濃度方程(考慮泥沙的沉降和再懸浮過程,即源匯項(xiàng)),推導(dǎo)出懸浮泥沙的年齡濃度方程,進(jìn)而分別計算了泥沙的“再懸浮年齡”和“輸運(yùn)年齡”,并得出在床面剪切力較弱的海域,懸沙因長期滯留在床面,導(dǎo)致“再懸浮年齡”較小,而“輸運(yùn)年齡”較大。隨后,龔文平等[36]將該研究中“輸運(yùn)年齡”的概念應(yīng)用于美國York River 河口,研究上游來沙的輸運(yùn)特征，指出保守性污染物的年齡主要受對流擴(kuò)散的影響,而泥沙的“輸運(yùn)年齡”還與徑流和泥沙的沖淤運(yùn)動密切相關(guān),隨著泥沙沉速以及臨界沖刷剪切應(yīng)力的減小,泥沙的“輸運(yùn)年齡”也相應(yīng)減小。Delhez等[37]指出,對于易被泥沙吸附的污染物,它們的輸運(yùn)過程需要考慮泥沙對其的吸附解吸作用,于是將年齡概念進(jìn)一步拓展到考慮泥沙影響的污染物輸運(yùn)時間研究中;以一維河口為例,假設(shè)泥沙對污染物的吸附作用大小可用一個分配系數(shù)(反映懸移質(zhì)泥沙的吸附能力)來反映,將泥沙對污染物的吸附解吸作用體現(xiàn)在污染物濃度方程的源匯項(xiàng)中,并根據(jù)該污染物濃度方程推導(dǎo)出考慮泥沙影響的污染物年齡濃度方程,進(jìn)而計算出考慮泥沙影響的污染物年齡,并討論了不同分配系數(shù)(懸移質(zhì)泥沙的吸附能力大小)對年齡的影響,指出泥沙對污染物的吸附能力越強(qiáng),年齡越大,但當(dāng)吸附能力達(dá)到一定程度時,這種影響將不再增大。Delhez等[37]首次將污染物年齡和泥沙運(yùn)動結(jié)合在一起,具有十分重要的意義,但研究局限于一維河道,且假定泥沙對污染物的吸附能力不隨空間變化,存在一定的局限性。

到目前為止,CART理論用于計算水體年齡和保守性物質(zhì)的年齡已相對較成熟,對用于非保守性物質(zhì)、泥沙和易被泥沙吸附的污染物年齡的計算上還有較大的發(fā)展空間,缺乏非保守性物質(zhì)、泥沙和易被泥沙吸附的污染物的存留時間、暴露時間的研究。能否將伴隨方法與非保守性物質(zhì)的降解過程、泥沙的沉降和再懸浮過程、易被泥沙吸附的污染物吸附解吸過程相結(jié)合[38-39],計算這些物質(zhì)的存留時間和暴露時間,是有待探討的一個問題。

4 結(jié) 語

物質(zhì)輸運(yùn)時間尺度是研究物質(zhì)輸運(yùn)過程的有效手段之一,對生態(tài)環(huán)境研究具有重要意義。自20世紀(jì)70年代開始不同的物質(zhì)輸運(yùn)時間尺度概念被提出,概念眾多,存在定義不嚴(yán)謹(jǐn)?shù)那闆r。從各種輸運(yùn)時間尺度的定義可以看出,沒有一種時間概念可以完整地刻畫整個物質(zhì)輸運(yùn)過程,也沒有一種輸運(yùn)時間計算方法適用于所有的研究領(lǐng)域,必須根據(jù)研究需要、水體特征、適用范圍等具體條件加以選擇。

在計算方法上,基于歐拉法的CART理論計算年齡和存留時間簡單實(shí)用,物理意義明確,在河口海岸水體及保守性物質(zhì)研究中得到了廣泛應(yīng)用。目前,該理論運(yùn)用于水體年齡或保守性物質(zhì)年齡的計算研究已經(jīng)較為成熟,但用于輸運(yùn)過程更為復(fù)雜的非保守性物質(zhì)、泥沙和易被泥沙吸附的污染物年齡的計算上還有較大的發(fā)展空間;利用伴隨方法雖然可通過一次計算得到存留時間、暴露時間的時空分布,但是也存在一定的缺陷,一方面它較難實(shí)現(xiàn),另一方面該方法基于伴隨方程的推導(dǎo),在應(yīng)用于計算考慮源匯項(xiàng)的物質(zhì)存留時間和暴露時間時不具有普適性,需根據(jù)具體情況重新推導(dǎo)存留時間的計算方程。將非保守性物質(zhì)、泥沙或易被泥沙吸附的污染物輸運(yùn)過程與CART理論相結(jié)合,研究其輸運(yùn)時間尺度,能更全面地掌握這些在水體中較為復(fù)雜的物質(zhì)的輸運(yùn)機(jī)理,從而為環(huán)境保護(hù)及治理提供幫助。

[1] MONSEN N E,CLOERN J E,LUCAS L V,et al.A comment on the use of flushing time,residence time,and age as transport time scales[J].Limnology and Oceanography,2002,47(5):1545-1553.

[2] 王鵬程.大遼河口物質(zhì)輸運(yùn)時間的數(shù)值研究[D].青島:中國海洋大學(xué),2013.

[3] BOLIN B,RODHE H.A note on the concepts of age distribution and transit time in natural reservoirs[J].Tellus,1973,25(1): 58-62.

[4] PRANDLE D.A modelling study of the mixing of137Cs in the seas of the european continental shelf[J].Philosophical Transactions of the Royal Society of London,1984,310(1513):407-436.

[5] LUFF R,POHLMANN T.Calculation of water exchange times in the ICES-boxes with a eulerian dispersion model using a half-life time approach[J].Ocean Dynamics,1995,47(4):287-299.

[6] ZIMMERMAN J T F.Mixing and flushing of tidal embayments in the western Dutch Wadden Sea:part I distribution of salinity and calculation of mixing time scales[J].Netherlands Journal of Sea Research,1976,10(2):149-191.

[7] TAKEOKA H.Exchange and transport time scales in the Seto Inland Sea[J].Continental Shelf Research,1984,3(4):327-341.

[8] 蔡忠亞.渤海水體平均存留時間及其季節(jié)變化數(shù)值研究[D].青島:中國海洋大學(xué),2013.

[9] ANDUTTA F P,RIDD P V,DELEERSNIJDER E,et al.Contaminant exchange rates in estuaries:new formulae accounting for advection and dispersion[J].Progress in Oceanography,2014,120(2):139-153.

[10] 彭輝,姚炎明,劉蓮.象山港水交換特性研究[J].海洋學(xué)研究,2012,30(4):1-12.(PENG Hui,YAO Yanming,LIU Lian.Study on the features of water exchange in Xiangshangang Bay[J].Journal of Marine Science,2012,30(4):1-12.(in Chinese))

[11] LIU Zhe,WEI Hao,LIU Guangshan,et al.Simulation of water exchange in Jiaozhou Bay by average residence time approach[J].Estuarine,Coastal and Shelf Science,2004,61(1):25-35.

[12] DYER K R,TAYLOR P A.A simple,segmented prism model of tidal mixing in well-mixed estuaries[J].Estuarine & Coastal Marine Science,1973,1(4):411-418.

[13] FISCHER H B,LIST J E,KOH R C,et al.Mixing in inland and coastal waters[M].New York:Academic Press,1979.

[14] 高抒,謝欽春.狹長形海灣與外海水體交換的一個物理模型[J].海洋通報,1991(3):1-9.(GAO Shu,XIE Qinchun.Pysical model for water exchange process in Xiangshan Bay[J].Marine Science Bulletin,1991(3):1-9.(in Chinese))

[15] 潘偉然.湄洲灣海水交換率和半更換期的計算[J].廈門大學(xué)學(xué)報(自然科學(xué)版),1992,31(1):65-68.(PAN Weiran.Computation of the turn-over rate and the half exchange period of seawater in Meizhou Bay[J].Journal of Xiamen University (Natural Science),1992,31(1):65-68.(in Chinese)).

[16] 陳益明.規(guī)劃環(huán)評中海灣水環(huán)境容量的計算方法[J].水資源保護(hù),2014,30(3):70-75.(CHEN Yiming.A method for calculation of environmental capacity of bay water during planning environmental impact assessment[J].Water Resource Protection,2014,30(3):70-75.(in Chinese))

[17] ANABELA O,BAPTISTA A M.Diagnostic modeling of residence times in estuaries[J].Water Resources Research,1997,33(8):1935-1946.

[18] 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(4):607-624.

[19] LIU Wencheng,CHEN Weibo,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 & Geosciences,2011,37(8):1148-1161.

[20] SHEN Yongming,WANG Jinhua,ZHENG Binghui,et al.Modeling study of residence time and water age in Dahuofang Reservoir in China[J].Science China Physics Mechanics & Astronomy,2011,54(1):127-142.

[21] SANDERY P A,KMPF J.Transport timescales for identifying seasonal variation in Bass Strait,South-Eastern Australia[J].Estuarine,Coastal and Shelf Science,2007,74(4):684-696.

[22] DELEERSNIJDER E,DELHEZ E J M.Timescale and tracer-based methods for understanding the results of complex marine models[J].Estuarine,Coastal and Shelf Science,2007,74(4):5-7.

[23] BRAUNSCHWEI G F,MARTINS F,CHAMBEL P,et al.A methodology to estimate renewal time scales in estuaries: the Tagus estuary case[J].Ocean Dynamics,2003,53(3):137-145.

[24] BRAUWERE A D,BRYE B D,BLAISE S,et al.Residence time,exposure time and connectivity in the Scheldt estuary[J].Journal of Marine Systems,2011,84(3/4):85-95.

[25] DELHEZ E J M,CAMPIN J M,HIRST A C,et al.Toward a general theory of the age in ocean modeling[J].Ocean Modeling,1999,99(1):17-27.

[26] DELHEZ R J M,HEEMINK A W,DELEERSNIJDER R.Residence time in a semi-enclosed domain from the solution of an adjoint problem[J].Estuarine,Coastal and Shelf Science,2004,61(4):691-702.

[27] ANDUTTA F P,HELFER F,MIRANDA L B D,et al.An assessment of transport timescales and return coefficient in adjacent tropical estuaries[J].Continental Shelf Research,2016,124:49-62.

[28] DELHEZ E J M.Transient residence and exposure times[J].Ocean Science,2005,2(1):1-9.

[29] DELHEZ E J M.On the concept of exposure time[J].Continental Shelf Research,2013,71(6):27-36.

[30] SHEN Jian,HAAS L.Calculating age and residence time in the tidal York River using three-dimensional model experiments[J].Estuarine,Coastal and Shelf Science,2004,61(3):449-461.

[31] SHEN Jian,WANG H V.Determining the age of water and long-term transport timescale of the Chesapeake Bay[J].Estuarine,Coastal and Shelf Science,2007,74(4):585-598.

[32] DE BRYE B,DE BRAUWERE A,GOURGUE O,et al.Reprint of water renewal timescales in the Scheldt Estuary[J].Journal of Marine Systems,2013,128(1):3-16.

[34] ZHANG Qin,TAO Jianfeng,YANG Jie.Numerical study on the transport timescale in a river-influenced macro-tidal estuary[J].Journal of Coastal Research,2016(Sup75):193-197.

[35] MERCIER C,DELHEZ E J M.Diagnosis of the sediment transport in the Belgian Coastal Zone[J].Estuarine,Coastal and Shelf Science,2007,74(4):670-683.

[36] GONG Wenping,SHEN Jian.A model diagnostic study of age of river-borne sediment transport in the tidal York River estuary[J].Environmental Fluid Mechanics,2010,10(1/2):177-196.

[37] DELHEZ E J M,WOLK F.Diagnosis of the transport of adsorbed material in the Scheldt estuary:a proof of concept[J].Journal of Marine Systems,2013,128:17-26.

[38] 馬俊超,李瓊芳,陸國賓，等. 三峽-葛洲壩梯級水庫不同蓄水階段對其下游泥沙特性的影響分析[J].水資源保護(hù),2016,32(1):75-78. (MA Junchao,LI Qiongfang,LU Guobin,et al. Analysis of impacts of different impound phases of Three Gorges-Gezhouba cascade reservoirs on sediment regime in downstream[J].Water Resources Protection, 2016,32(1):75-78.(in Chinese))

[39] 王鵬,李瓊瓊,蔡赟杰.太湖沉積物再懸浮對雙酚A吸附性能的影響[J].水資源保護(hù),2015,31(5):35-41.(WANG Peng,LI Qiongqiong,CAI Yunjie. Effect of sediment resuspension on BPA adsorption in Taihu Lake[J]. Water Resources Protection, 2015,31(5):35-41. (in Chinese))

Advancesinstudyofsubstancestransporttimescalesinestuarineandcoastalareas//

ZHANG Qin1,YANG Jie2,TAO Jianfeng2,ZHANG Changkuan2,ZHANG Xiao2

(1.ShanghaiInvestigation,Design&ResearchInstituteCo.,Ltd.,Shanghai200434,China; 2.CollegeofHarbour,CoastalandOffshoreEngineering,HohaiUniversity,Nanjing210098,China)

Research progresson the transport timescales of substances in estuaries and coastal areas in recent years was reviewed from three aspects, including the concepts, the research methods and the applications of transport timescales. The advantages and limitations of the current mathematical models for the computation of substance transport timescales, especially the progress of the constituent-oriented age and residence time (CART) theory are discussed emphatically in this paper. Aiming at the shortcomings of the current non-conservative substance transport timescales, it is pointed out that the interaction between adsorptive pollutants and suspended particles should be considered to describe the transport timescales completely. Strengthening the study on the transport timescales of non-conservative materials, sediments and adsorptive pollutants will help to reveal the transport characteristics of complex substances. The CART theory should be applied to calculate the transport timescales of non-conservative substances, sediments and the pollutant transport process that considers the interaction with sediments.

transport of substances;pollutants; sediments; transport timescales; estuaries; coastal areas

國家自然科學(xué)基金(51409093)

張琴(1992—),女,助理工程師,主要從事水環(huán)境數(shù)值模擬及環(huán)境影響評價研究。E-mail: 1104264732@qq.com

楊潔(1980—),女,講師,博士,主要從事海岸水環(huán)境與泥沙輸運(yùn)研究。E-mail: jie_yang@hhu.edu.cn

10.3880/j.issn.1006-7647.2017.06.015

P731.26

A

1006-7647(2017)06-0088-07

2016-12-14 編輯：熊水斌)