胡思骙 張詩軍 任晨媛 王慶國* 梁 英

1.四川大學(xué)建筑與環(huán)境學(xué)院

2.成都市濕地保護(hù)中心

3.成都大熊貓繁育研究基地

數(shù)學(xué)模型的發(fā)展從上世紀(jì)60年代開始已經(jīng)有眾多學(xué)者開展了相關(guān)的研究，從開始的單層、單室、單成分和零維的簡單模型到70年代后期的多層、多室、多成分和三維的復(fù)雜生態(tài)模型用來揭示湖泊水體生態(tài)系統(tǒng)的動力學(xué)演變，并逐步開始應(yīng)用在河湖水庫的污染管理和水生態(tài)恢復(fù)治理等水環(huán)境管理方面。根據(jù)對不同水生態(tài)模型建立框架和思路的不同，可將模型分為簡單水質(zhì)回歸模型、單一營養(yǎng)鹽負(fù)荷平衡模型、生態(tài)動力學(xué)模型和復(fù)雜水質(zhì)-生態(tài)-水動力耦合模型四類[1]。

加拿大學(xué)者Vollenweider 在1968 年首次提出簡單的總磷輸入輸出平衡模型，開始了營養(yǎng)鹽平衡模型的發(fā)展，該數(shù)學(xué)模型假定湖庫中的TP 濃度隨時間的變化值是單位體積水體內(nèi)輸入的磷與從湖庫中輸出的磷及沉淀到沉積物中的磷之間的差值。之后Laesen、Mercier 和Dillon 等學(xué)者在后續(xù)的研究中對Vollenweider模型進(jìn)行了相關(guān)改進(jìn)和修正，考慮湖泊上下分層和輸入輸出及沉積物之間磷的交換，應(yīng)用不同形態(tài)磷的微分方程將營養(yǎng)鹽濃度與生物量描述為隨時間和水深變化的連續(xù)函數(shù)。近年來，營養(yǎng)鹽平衡的水質(zhì)模型已經(jīng)從一維發(fā)展為三維，從單個線性方程到多個線性系統(tǒng)方程和初步的相互作用非線性系統(tǒng)模型[2]。水質(zhì)模型在應(yīng)用到湖泊水體時需要的較少且簡單的基礎(chǔ)數(shù)據(jù)，因此在預(yù)測湖泊水體的水質(zhì)演變是具有較高的靈活度和適用性。Dahl 等利用水質(zhì)模型對瑞典Vnern 湖的磷和懸浮顆粒物的相互關(guān)系進(jìn)行模擬，通過對目標(biāo)水體進(jìn)行分區(qū)提高了模型的準(zhǔn)確率[3]。Arhonditsis 等針對華盛頓湖開發(fā)了富營養(yǎng)化模型，對于包括磷酸鹽、總磷、總氮和溶解氧的大部分水質(zhì)參數(shù)都能單一進(jìn)行平衡模擬，且模型預(yù)測值的相對誤差低于20%[4]。然而目前的營養(yǎng)鹽平衡水質(zhì)模型并未徹底認(rèn)清污染物在湖泊水體中的遷移轉(zhuǎn)化規(guī)律，使模擬的數(shù)值與實際水質(zhì)偏差較大。為提高水生態(tài)模型對不同類型湖泊水體的適用性，逐步開發(fā)出一批用于湖泊水庫水生態(tài)系統(tǒng)演變規(guī)律的承受商用水生態(tài)模型，比如MIKE系列模型、EFDC 模型、WASP 模型、AQUATOX 模型和QUAL2K模型。

近年來，國內(nèi)外對于富營養(yǎng)化水體進(jìn)行的數(shù)模研究也越來越多，特別是針對國內(nèi)一些富營養(yǎng)化問題嚴(yán)重大型湖泊，比如劉玉生在1991 年建立最早的滇池生態(tài)系統(tǒng)動力學(xué)模型[5]、太湖流域藻類生長綜合數(shù)學(xué)模擬[6]、基于Ecopath 模型建立的巢湖生態(tài)系統(tǒng)模型[7]、武漢東湖三維水質(zhì)和水動力數(shù)值模擬[8]。對目前國內(nèi)對湖泊生態(tài)系統(tǒng)數(shù)值模擬進(jìn)行綜合分析，具有代表性的研究有：

楊具瑞等根據(jù)對滇池的實測水質(zhì)和邊界數(shù)據(jù)資料對其富營養(yǎng)化情況建立了Chl-a、透明度、藻量生物量和TN、TP 之間的相互關(guān)系，搭建了垂向平均化的水動力-水質(zhì)模型，對滇池的上述五個指標(biāo)的分布進(jìn)行了數(shù)值模擬計算[9]；李一平等利用EFDC 模型對太湖采用拉丁超立方（LHS）取樣的方法研究了風(fēng)拖曳、粗糙度、渦粘系數(shù)、紊流擴(kuò)散系數(shù)及風(fēng)遮擋系數(shù)對湖泊流場的影響模擬，發(fā)現(xiàn)風(fēng)場和床面粗糙度是主要影響參數(shù)[10]；謝彪以武漢大東湖生態(tài)水網(wǎng)的構(gòu)建為研究對象，首先使用MIKE21FM 軟件構(gòu)建東湖的水動力模型，并針對東湖的實測水質(zhì)數(shù)據(jù)利用MIKE ECO Lab的水質(zhì)模塊編寫使用的水質(zhì)模型，搭建東湖水動力-水質(zhì)耦合模型，對大東湖生態(tài)系統(tǒng)構(gòu)建的幾種工程措施進(jìn)行了數(shù)值模擬，分析其對水環(huán)境的影響[11]；計勇以紊動擴(kuò)散理論為基礎(chǔ)，在對流擴(kuò)散方程中引入了底泥污染項參數(shù)，構(gòu)建考慮了底泥污染影響的垂向水平三維水質(zhì)模型，并結(jié)合玄武湖底泥營養(yǎng)物釋放數(shù)據(jù)建立底泥污染物釋放通量方程，耦合為水動力-水質(zhì)-底泥水生態(tài)模型，對玄武湖的流場和濃度場進(jìn)行了較為精準(zhǔn)的模擬[12]；劉曉臣等使用生態(tài)動力學(xué)模型DELFT-3D對興凱湖的入湖污染物負(fù)荷和富營養(yǎng)化狀況進(jìn)行了數(shù)值模擬，根據(jù)實測的基礎(chǔ)數(shù)據(jù)校正模型確立了適用于興凱湖的水環(huán)境特征參數(shù)，并較好的模擬了水體的水質(zhì)演變情況和藻類濃度的動態(tài)變化，為興凱湖可能發(fā)生的藻類水花提供了預(yù)警[13]；姜毅等針對人工景觀湖泊比較封閉的情況，基于WASP 富營養(yǎng)模型對西安芙蓉湖進(jìn)行了水生態(tài)模擬研究，考慮了氨氮、硝態(tài)氮、無機(jī)磷、浮游植物、BOD、溶解氧、有機(jī)氮和有機(jī)磷等8個水質(zhì)指標(biāo)，系統(tǒng)的分析了浮游植物生長系統(tǒng)、氮循環(huán)、磷循環(huán)和溶解氧平衡，與實測數(shù)據(jù)對比擬合度較高[14]；陳彥熹針對有非常規(guī)水體作為主要補(bǔ)水來源的城市景觀水體進(jìn)行了水環(huán)境數(shù)值模擬研究，利用CFD 模型FLUENT對湖泊的水動力進(jìn)行分析，并結(jié)合湖泊的氣象資料和外源污染負(fù)荷分析利用水生態(tài)模型AQUATOX 進(jìn)行水質(zhì)和水生態(tài)數(shù)值模擬，分析了水體的流動摻混、營養(yǎng)鹽循環(huán)、顆粒物及附浮游生物的食物網(wǎng)之間的相互關(guān)系進(jìn)行了分析[15]。陳無歧等利用AQUATOX 模型對洱海的富營養(yǎng)化程度加劇及藍(lán)藻爆發(fā)等問題進(jìn)行了模擬研究，模擬報道了營養(yǎng)鹽的輸入對洱海富營養(yǎng)化狀態(tài)的響應(yīng)關(guān)系，結(jié)果顯示藻類對N、P的輸入都比較敏感[16]。

總的來說，國內(nèi)對于水環(huán)境模型的研究多是針對大型的湖泊或水庫，針對城市景觀水體的研究仍不多，且多為借鑒和引用已經(jīng)成熟的國外商用模型為主。