含水層中重金屬溶質(zhì)運(yùn)移動(dòng)力學(xué)模型研究進(jìn)展

陶 安,劉曙光,婁 廈,代朝猛

(同濟(jì)大學(xué)土木工程學(xué)院,上海 200092)

?

含水層中重金屬溶質(zhì)運(yùn)移動(dòng)力學(xué)模型研究進(jìn)展

陶 安,劉曙光,婁 廈,代朝猛

(同濟(jì)大學(xué)土木工程學(xué)院,上海 200092)

在總結(jié)國(guó)內(nèi)外重金屬污染研究的基礎(chǔ)上,綜述含水層中重金屬溶質(zhì)運(yùn)移動(dòng)力學(xué)模型的發(fā)展,歸納總結(jié)模型的計(jì)算方法。討論了模型參數(shù)的確定方法,以及溫度,pH值,孔隙水流速,有機(jī)質(zhì)與膠體,土壤非均勻性對(duì)參數(shù)確定的影響。最后,指出含水層中重金屬溶質(zhì)運(yùn)移過(guò)程中目前研究存在的主要問(wèn)題,并對(duì)將來(lái)在該領(lǐng)域需重點(diǎn)研究的內(nèi)容進(jìn)行展望。

地下水;含水層;重金屬運(yùn)移;數(shù)值模擬;動(dòng)力學(xué)模型;土壤非均勻性

由于地表水容易受到自然條件和人為因素的影響,地下水已經(jīng)逐漸成為城市供水的主要水源。重金屬溶質(zhì)以其不可降解性和生物積累性等特點(diǎn)成為威脅地下水資源的主要污染物之一。研究表明,我國(guó)地下水中80%的淺層地下水和30%的深層地下水已不同程度地受到了重金屬污染,其中最主要的重金屬污染物為汞、鎘、鉻和鉛[1]。根據(jù)《全國(guó)地下水污染防治規(guī)劃(2011—2020年)》,近十幾年中我國(guó)地下水重金屬污染呈現(xiàn)由點(diǎn)到面、由淺到深、由城市向農(nóng)村蔓延的趨勢(shì),污染程度日益嚴(yán)重。重金屬污染已成為我國(guó)地下水污染控制的重中之重,研究重金屬溶質(zhì)在地下水含水層中的運(yùn)移規(guī)律也受到了眾多研究人員的高度關(guān)注。

重金屬溶質(zhì)在含水層中的運(yùn)移受到對(duì)流、水動(dòng)力彌散、吸附解吸和氧化還原等多種復(fù)雜作用的影響[2-3]。應(yīng)用數(shù)學(xué)模型模擬重金屬溶質(zhì)在含水層中的運(yùn)移可以預(yù)測(cè)地下水中重金屬污染的瞬時(shí)情況和污染物的擴(kuò)展范圍。數(shù)學(xué)模型可分為確定型模型和隨機(jī)型模型。確定型模型描述含水層中污染物的運(yùn)移是以1961年Nielsen等[4-5]提出的對(duì)流彌散方程為基礎(chǔ)的,根據(jù)確定的參數(shù)及邊界條件可求得確定型模型的解。隨機(jī)型模型是從研究含水層非均質(zhì)性及其尺度效應(yīng)開(kāi)始的,依賴于多孔介質(zhì)性質(zhì)變異的隨機(jī)函數(shù)變量[6]。本文綜述含水層中重金屬溶質(zhì)運(yùn)移動(dòng)力學(xué)模型的發(fā)展,并根據(jù)國(guó)內(nèi)外的相關(guān)研究總結(jié)模型的計(jì)算方法。在此基礎(chǔ)上,進(jìn)一步討論模型參數(shù)的確定方法及影響參數(shù)確定的因素。最后,作者針對(duì)目前含水層中重金屬溶質(zhì)運(yùn)移研究過(guò)程中存在的主要問(wèn)題進(jìn)行探討。

1 重金屬溶質(zhì)運(yùn)移動(dòng)力學(xué)模型的發(fā)展

1.1 確定型模型

達(dá)西定律自提出以來(lái),由于其形式簡(jiǎn)單,一直被認(rèn)為是地下水動(dòng)力學(xué)中最為基本的定律，地下水流動(dòng)基本微分方程、裘布依穩(wěn)定井流模型等都是基于其發(fā)展的。但達(dá)西定律提出后的很長(zhǎng)時(shí)間里,國(guó)內(nèi)外學(xué)者對(duì)地下水研究都未考慮溶質(zhì)運(yùn)移問(wèn)題。1952年Lpaidus和Amundson首次將一個(gè)類似于“對(duì)流-擴(kuò)散”方程的數(shù)學(xué)模型應(yīng)用于溶質(zhì)運(yùn)移問(wèn)題[7]。1960年Nielsen等[4-5]基于質(zhì)量守恒推導(dǎo)了“對(duì)流-彌散”方程(CDE),其一般形式見(jiàn)式(1)所示。此后,該方程一直是研究多孔介質(zhì)中溶質(zhì)運(yùn)移的經(jīng)典方程和基本方程[8]。

(1)

式中:θ為土壤體積含水量,%;c為溶質(zhì)質(zhì)量濃度,mg/L;ρ為土壤干容重,g/cm3;S為溶質(zhì)在土壤上的吸附量,%;D(θ,v)ij為水動(dòng)力彌散系數(shù),cm2/min;v為對(duì)流速度,cm/min;qi為土壤水流通量,cm/min;t為時(shí)間,min;xi為空間坐標(biāo),x1,x2,x3分別代表x,y,z;φk為源項(xiàng),包括生物吸收、化學(xué)反應(yīng)、降解、沉淀等過(guò)程,對(duì)于重金屬溶質(zhì)來(lái)說(shuō),可以忽略其降解項(xiàng)。

式(1)中左邊第2項(xiàng)求解是很難的,為了便于研究,一般進(jìn)行轉(zhuǎn)化,見(jiàn)式(2)所示。

(2)

(3)

式中:Kd為分配系數(shù),L/g;λ為源項(xiàng)系數(shù)。

本文引入了描述污染物吸附解吸作用的1個(gè)重要參數(shù),即阻滯系數(shù)Rd,其表達(dá)式見(jiàn)式(4)。

(4)

此后,國(guó)外學(xué)者Beyer等[10]基于實(shí)測(cè)土壤屬性和鎘離子的土柱淋溶實(shí)驗(yàn),建立了德國(guó)諾爾登哈姆地區(qū)地下水中鎘離子運(yùn)移的對(duì)流-彌散方程。錢(qián)天偉等[11]認(rèn)為,使用吸附等溫線的概念并不能解釋“固-液”界面所發(fā)生復(fù)雜物理化學(xué)反應(yīng)的內(nèi)在規(guī)律。他利用表面絡(luò)合吸附理論研究了地下水中“固-液”表面的溶質(zhì)分布特征,并將其與傳統(tǒng)的“對(duì)流-彌散”模型相耦合,建立了一個(gè)能夠考慮表面絡(luò)合吸附影響的溶質(zhì)運(yùn)移模型,模型模擬結(jié)果與實(shí)際值吻合得較好。

對(duì)于非均質(zhì)含水層中的溶質(zhì)運(yùn)移,很多情況下會(huì)出現(xiàn)反常擴(kuò)散,因此無(wú)法應(yīng)用傳統(tǒng)對(duì)流-彌散理論對(duì)其進(jìn)行描述。許多學(xué)者開(kāi)始質(zhì)疑傳統(tǒng)“對(duì)流-彌散”模型所基于的理想假設(shè)給現(xiàn)實(shí)含水層溶質(zhì)運(yùn)移問(wèn)題的描述和模擬帶來(lái)了不可避免的誤差和缺陷[12]。近年來(lái),對(duì)傳統(tǒng)二階“對(duì)流-彌散”方程進(jìn)行改進(jìn)的主流方法之一就是分別對(duì)時(shí)間或空間進(jìn)行分?jǐn)?shù)微分,進(jìn)而依據(jù)非線性函數(shù)的關(guān)聯(lián)方程描述相應(yīng)的彌散系數(shù)。常福宣等[13]對(duì)彌散過(guò)程考慮時(shí)間和空間相關(guān)性,用非局域處理法,推導(dǎo)出對(duì)時(shí)間和空間分?jǐn)?shù)階微分的“對(duì)流-彌散”方程,見(jiàn)式(5)所示。

(5)

Xia等[14]通過(guò)對(duì)非均勻含水層中的反常擴(kuò)散數(shù)值試驗(yàn)進(jìn)一步說(shuō)明了分?jǐn)?shù)階“對(duì)流-彌散”方程在研究含水層污染物運(yùn)移中的重要意義。

傳統(tǒng)“對(duì)流-彌散”模型能很好地體現(xiàn)出含水層中溶質(zhì)運(yùn)移的過(guò)程機(jī)理,其參數(shù)也具有明確的物理意義,故目前的應(yīng)用以確定型模型為主。但由于地下水運(yùn)動(dòng)中的不確定性因素較多,包括含水層介質(zhì)的物理化學(xué)特性以及溶質(zhì)運(yùn)移過(guò)程和運(yùn)動(dòng)參數(shù)等,模型中參數(shù)隨時(shí)空的變異性會(huì)產(chǎn)生尺度效應(yīng)問(wèn)題[15]。即使是改進(jìn)的分?jǐn)?shù)階“對(duì)流-彌散”方程,在模擬較長(zhǎng)距離的污染物運(yùn)移時(shí),仍會(huì)出現(xiàn)尺度效應(yīng)。因此,對(duì)于確定型模型在不同的模擬尺度下如何消除或減少尺度效應(yīng)是一個(gè)需要研究的問(wèn)題。

1.2 隨機(jī)型模型

近幾十年來(lái),隨機(jī)理論迅速地應(yīng)用于地下水污染物的運(yùn)移研究中。隨機(jī)型模型依賴于含水層性質(zhì)變異的隨機(jī)函數(shù)變量,引入含水層空間隨機(jī)場(chǎng)的概念,在理論上將污染物的宏觀彌散與含水層介質(zhì)的統(tǒng)計(jì)特征相聯(lián)系,從而解決了宏觀彌散系數(shù)的求解問(wèn)題[16]。Acar等[17]提出了估計(jì)一個(gè)地下水和污染物排放量的隨機(jī)模型,并用蒙特卡羅方法估計(jì)了參數(shù)不確定性。Coppola等[18]應(yīng)用隨機(jī)方法研究了介質(zhì)垂向非均勻性對(duì)于溶質(zhì)運(yùn)移的影響。

對(duì)于隨機(jī)型模型,非均勻含水層統(tǒng)計(jì)參數(shù)的方法、模型的理論完善及野外試驗(yàn)驗(yàn)證方面仍需要進(jìn)一步的研究。目前許多研究者都很重視隨機(jī)理論和隨機(jī)型模型的應(yīng)用研究,雖然其仍處于發(fā)展的初期,距離實(shí)際應(yīng)用還有一定距離,但它仍然是一種比較有應(yīng)用前景的理論模型。

2 模型的計(jì)算方法

2.1 確定型模型計(jì)算方法

除了在條件比較簡(jiǎn)單、理想的情況下可求得解析解外,數(shù)值方法是解決地下水污染實(shí)際問(wèn)題最常用的方法。數(shù)值方法可歸納為3類方法,即歐拉方法、拉格朗日方法以及二者的混合方法。目前,歐拉方法應(yīng)用較多,而拉格朗日方法和歐拉-拉格朗日方法還只停留在研究階段,在實(shí)際中應(yīng)用較少[16]。歐拉方法主要包括有限差分法、有限元法、有限體積法和邊界元法等。國(guó)內(nèi)外學(xué)者的研究情況如表1所示。

表1 確定型模型歐拉求解方法

2.2 隨機(jī)型模型計(jì)算方法

地下水隨機(jī)型模型計(jì)算方法主要有矩方程法、攝動(dòng)待定系數(shù)法和蒙特·卡羅法[26]。矩方程方法是通過(guò)求解有關(guān)均值和協(xié)方差的隨機(jī)偏微分方程獲得隨機(jī)問(wèn)題的解;攝動(dòng)待定系數(shù)法是把隨機(jī)變量表示為確定性部分和由攝動(dòng)引起的隨機(jī)部分;蒙特·卡羅法則是通過(guò)平均一系列反映含水層實(shí)際性質(zhì)的確定性問(wèn)題來(lái)模擬隨機(jī)過(guò)程的一種計(jì)算機(jī)模擬方法[16]。

經(jīng)過(guò)十幾年的發(fā)展,近年來(lái)國(guó)內(nèi)外不少學(xué)者提出了自己關(guān)于模型計(jì)算和建模理念方面的新見(jiàn)解。Mazzia等[27]提出了一種新的求解高鹽分在地下水中運(yùn)移動(dòng)力學(xué)二維非線性方程的數(shù)值方法。該方法對(duì)時(shí)間和空間采用了不同的差分格式,保證了流速場(chǎng)和濃度場(chǎng)預(yù)測(cè)的高精度性和高穩(wěn)定性。王浩然等[28]利用區(qū)域分解法建立了王旺莊水源地地下水的邊界單元法(BEM)與有限單元法(FEM)耦合模型,在需要重點(diǎn)研究的地段采用有限單元法處理,而非重點(diǎn)研究的外圍地區(qū)則采用邊界單元法解決。O’Malley等[29]綜合了考慮地球化學(xué)反應(yīng)、對(duì)流、彌散及其他微尺度作用,利用概率論方法求解了反常擴(kuò)散大尺度隨機(jī)型模型的解析解,并對(duì)參數(shù)的不確定性進(jìn)行了分析。

3 模型參數(shù)的確定方法及其影響因素

3.1 模型參數(shù)的確定方法

在動(dòng)力學(xué)模型的建立過(guò)程中,參數(shù)確定的準(zhǔn)確性研究是很重要的,是污染物在含水層中運(yùn)移轉(zhuǎn)化模擬計(jì)算的關(guān)鍵問(wèn)題。

水動(dòng)力彌散系數(shù)是最主要的運(yùn)移參數(shù),主要通過(guò)控制彌散過(guò)程來(lái)影響運(yùn)移,但其測(cè)定具有較大難度,且一直存在較多爭(zhēng)議[30]。國(guó)內(nèi)外測(cè)定含水層中水動(dòng)力彌散系數(shù)的計(jì)算方法主要有公式法、穿透曲線法和水平土柱吸滲法、垂直土柱吸滲法、極大似然法、瞬時(shí)剖面法、水分特征曲線法和反函數(shù)法等[31-32]。水動(dòng)力彌散系數(shù)的測(cè)定方法及國(guó)內(nèi)外學(xué)者的研究情況如表2所示。

阻滯系數(shù)是用來(lái)表示含水層介質(zhì)對(duì)溶質(zhì)運(yùn)移阻滯能力的物理化學(xué)參數(shù),可以采用室內(nèi)試驗(yàn)或野外試驗(yàn)方法進(jìn)行確定。野外試驗(yàn)成本大,耗時(shí)長(zhǎng),且重金屬溶質(zhì)會(huì)造成污染,因此國(guó)內(nèi)外都傾向于通過(guò)室內(nèi)試驗(yàn)方式獲得阻滯系數(shù)。室內(nèi)測(cè)定阻滯系數(shù)主要方法有動(dòng)態(tài)柱法和批式法。通常情況下,批式法得出的阻滯系數(shù)值要比動(dòng)態(tài)柱法的結(jié)果大得多。相比之下,動(dòng)態(tài)柱法更接近野外實(shí)際情況,但動(dòng)態(tài)柱法實(shí)驗(yàn)周期太長(zhǎng)[41]。

3.2 模型參數(shù)的影響因素

含水層中重金屬溶質(zhì)運(yùn)移轉(zhuǎn)化動(dòng)力學(xué)模型的參數(shù)受多種條件的影響,因此,在確定模型參數(shù)時(shí)必須綜合考慮溫度、pH值、孔隙水流速、有機(jī)質(zhì)與膠體和土壤非均勻性等因素的影響。

表2 水動(dòng)力彌散系數(shù)測(cè)定方法

3.2.1 溫度

溫度的改變會(huì)影響離子運(yùn)動(dòng)的劇烈程度。Medina[42]等研究認(rèn)為在恒定壓強(qiáng)下,彌散系數(shù)隨著溫度的增加而增大,但隨著壓強(qiáng)的增大,彌散系數(shù)受溫度的影響減小。Do[43]研究了不同溫度對(duì)于Zn、Cd離子在黏土含水層中運(yùn)移的影響。研究表明,溫度從15℃上升到55℃,兩種重金屬離子的彌散系數(shù)都會(huì)增大,多達(dá)10倍。溫度對(duì)于阻滯系數(shù)的影響是基于重金屬離子的吸附機(jī)理,Cd主要以離子態(tài)存在,阻滯系數(shù)受溫度影響變化小;Zn主要以碳酸鹽相存在,阻滯系數(shù)隨溫度增加而增加,多達(dá)2倍。因此,在溫度從15℃上升到55℃時(shí),含水層中Cd離子的彌散鋒從7 mm處擴(kuò)大到了15 mm處,而Zn離子的彌散鋒幾乎不變。

3.2.2 pH值

pH值反映的是地下水中的H+濃度,H+濃度會(huì)對(duì)重金屬離子與含水層介質(zhì)發(fā)生的絡(luò)合作用產(chǎn)生影響。Moradkhani等[44]通過(guò)試驗(yàn)研究發(fā)現(xiàn),Zn、Cd和Mn 3種重金屬離子的最大分解百分比隨著pH值的增大而減小。Chotpantarat等[45]對(duì)4種重金屬離子溶液在pH值分別為4和5的條件下進(jìn)行淋溶試驗(yàn),研究表明,在溶液pH值為5時(shí),由于H+濃度的減小,重金屬離子的阻滯系數(shù)較溶液pH值為4時(shí)大。pH值所產(chǎn)生的不同運(yùn)移效果可能有以下原因:pH值較低時(shí),以靜電吸附為主,H+濃度達(dá)到破壞含水層介質(zhì)表面官能團(tuán)與重金屬的絡(luò)合,使吸附態(tài)重金屬離子全部交換解吸下來(lái);隨pH值升高,含水層介質(zhì)表面含氧基團(tuán)得到活化,更多的重金屬離子被絡(luò)合,不易解吸下來(lái);當(dāng)pH值超過(guò)一定值后,會(huì)產(chǎn)生沉淀反應(yīng),也導(dǎo)致解吸量急劇減少[46]。

3.2.3 孔隙水流速

含水層中孔隙水流速會(huì)隨時(shí)間和空間變化,但在求解溶質(zhì)運(yùn)移方程時(shí)通常忽略其對(duì)水動(dòng)力彌散系數(shù)的影響。Bond[47]求解了考慮孔隙水流速對(duì)水動(dòng)力彌散系數(shù)影響的非飽和非穩(wěn)定流的水動(dòng)力彌散方程,研究結(jié)果表明,考慮孔隙水流速的影響會(huì)使彌散寬度在短時(shí)間內(nèi)增大很多。Pang等[48]通過(guò)均一壓實(shí)消除物理非平衡影響的土柱進(jìn)行孔隙水流速對(duì)于Cd、Zn、Pb離子運(yùn)移影響的實(shí)驗(yàn),結(jié)果表明,分配系數(shù)與孔隙水流速呈正相關(guān),阻滯系數(shù)與孔隙水流速呈負(fù)相關(guān),這說(shuō)明低流速使重金屬離子與介質(zhì)接觸得充分,利于吸附達(dá)到動(dòng)態(tài)平衡,因而吸附量較大。

3.2.4 有機(jī)質(zhì)與膠體

含水層中的有機(jī)質(zhì)是一個(gè)不均勻的混合體,由于有機(jī)質(zhì)具有較多的含氧官能團(tuán)(羧基、酚基、羥基等),重金屬與有機(jī)物質(zhì)之間有很強(qiáng)的親和性。章明奎等[49]用室內(nèi)分析方法研究了土壤顆粒狀有機(jī)質(zhì)(POM)對(duì)Cu、Cd、Co、Pb、Ni和Zn等6種常見(jiàn)重金屬離子的吸附作用。研究表明,土壤中顆粒狀有機(jī)質(zhì)對(duì)重金屬離子有較高的吸附能力,吸附過(guò)程是一個(gè)快速反應(yīng),在100 min內(nèi)可接近平衡。余貴芬等[50]對(duì)含水層中鎘的有效性研究發(fā)現(xiàn),有機(jī)質(zhì)對(duì)含水層中重金屬的離子活性有很大的影響,有機(jī)質(zhì)的分子量、組分和環(huán)境條件均會(huì)對(duì)鎘的移動(dòng)性和生物活性造成影響。膠體是指粒徑范圍在1 nm～10 μm之間的能起到污染物載體作用的顆粒物,包括層狀硅酸鹽、鐵鋁三氧化物、有機(jī)大分子等[51-52],在土壤中普遍存在。膠體與含水層中污染物的相互作用是一個(gè)極為復(fù)雜的微觀過(guò)程[53-55]。Saiers等[56]進(jìn)行了銫離子和無(wú)機(jī)膠體的協(xié)同運(yùn)移室內(nèi)實(shí)驗(yàn),結(jié)果表明減小孔隙水離子強(qiáng)度可以使無(wú)機(jī)膠體促進(jìn)銫離子的運(yùn)移。Baumann等[57]通過(guò)對(duì)垃圾填埋場(chǎng)的現(xiàn)場(chǎng)試驗(yàn)研究發(fā)現(xiàn),膠體會(huì)從垃圾填埋場(chǎng)與含水層的交界面滲出,并形成一層致密層。此致密層減小了交界面的水力滲透系數(shù),也增加了含水層介質(zhì)對(duì)重金屬離子的吸附能力。

3.2.5 土壤非均勻性

室內(nèi)測(cè)定水動(dòng)力彌散系數(shù)一般是在均勻土中進(jìn)行的,然而現(xiàn)實(shí)的土壤卻是非均勻的。土壤的非均勻性會(huì)導(dǎo)致水流在含水層中流向流速等發(fā)生變化,進(jìn)而影響水動(dòng)力彌散系數(shù)。尚熳廷等[58]以Br-作示蹤劑,通過(guò)室內(nèi)實(shí)驗(yàn)測(cè)定含大孔隙的原狀土和均質(zhì)擾動(dòng)土的水動(dòng)力彌散系數(shù)和彌散度。結(jié)果表明,大孔隙的存在使原狀土水動(dòng)力彌散系數(shù)比擾動(dòng)土大2個(gè)數(shù)量級(jí)。土壤的非均勻性和空間變異性導(dǎo)致了滲透系數(shù)存在尺度效應(yīng),而滲透性能的變異性直接影響到含水層中溶質(zhì)的運(yùn)移過(guò)程[59]。李國(guó)敏等[60]認(rèn)為產(chǎn)生尺度效應(yīng)主要原因就是含水層介質(zhì)具有非均質(zhì)性,他通過(guò)解析方法和數(shù)值模型求出了尺度效應(yīng)的分維數(shù),結(jié)果表明數(shù)值模型所得出的縱向彌散度隨研究尺度增加而增大的速度小于用解析模型所求出的相應(yīng)值。王永森等[61]通過(guò)現(xiàn)場(chǎng)觀測(cè)數(shù)據(jù)應(yīng)用濃度方程進(jìn)行了觀測(cè)井彌散系數(shù)的反求,應(yīng)用分維理論確定了彌散系數(shù)尺度效應(yīng)的分維數(shù)為1.82567,即彌散系數(shù)與尺度之間呈冪相關(guān)關(guān)系。

4 結(jié)論與展望

基于以上對(duì)含水層中重金屬溶質(zhì)運(yùn)移動(dòng)力學(xué)模型研究進(jìn)展的評(píng)述,今后該領(lǐng)域的研究應(yīng)重點(diǎn)關(guān)注以下幾個(gè)方面:

a.重金屬污染物在地下水系統(tǒng)運(yùn)移的整體性研究。重金屬污染物一般是先通過(guò)非飽和帶再進(jìn)入地下水,如垃圾填埋場(chǎng)滲漏的重金屬污染物等。由于非飽和帶的特殊性,水流在非飽和帶的運(yùn)動(dòng)與飽和帶不同,溶質(zhì)運(yùn)移過(guò)程則更為復(fù)雜。然而以往的研究主要是對(duì)非飽和帶與飽和帶進(jìn)行獨(dú)立研究,缺少符合實(shí)際情況的整體性研究,未來(lái)在該方向應(yīng)作深入的探索。

b.含水層的空隙。對(duì)于一般的含水層,空隙很小,地下水流速緩慢,均符合達(dá)西滲流定律。而當(dāng)含水層介質(zhì)中出現(xiàn)較大裂隙甚至是空洞時(shí),地下水流速會(huì)變得很大,此時(shí)達(dá)西定律不再適用,如喀斯特溶洞、寬大的構(gòu)造裂隙等。此時(shí)溶解于地下水中的重金屬溶質(zhì)隨著地下水的流動(dòng)而快速運(yùn)移,造成的危害也更大。目前,國(guó)內(nèi)外對(duì)于此問(wèn)題尚處在探索階段,機(jī)理研究、數(shù)值模擬方法等方面有待進(jìn)一步發(fā)展。

c.模型參數(shù)的確定及尺度效應(yīng)。重金屬溶質(zhì)的運(yùn)移參數(shù)受到多方面的影響。如重金屬離子的吸附與解吸,不僅與含水層介質(zhì)的粒徑有關(guān),還與原先吸附在介質(zhì)表面的重金屬離子種類、濃度和價(jià)位有關(guān)。這是一個(gè)復(fù)雜的物理化學(xué)微觀過(guò)程,而目前的數(shù)學(xué)模型對(duì)于該過(guò)程的研究極為薄弱。因此今后的研究應(yīng)加強(qiáng)學(xué)科間交叉,宏觀結(jié)合微觀來(lái)探索重金屬溶質(zhì)在含水層中的運(yùn)移。此外,模型在不同的模擬尺度下如何消除或較少尺度效應(yīng)有待于進(jìn)一步探索。

d.邊界條件的研究。隨著全球溫室效應(yīng)的加劇,極端天氣事件頻出。暴雨、干旱等通過(guò)影響含水層中流場(chǎng)的變化來(lái)影響重金屬的運(yùn)移規(guī)律。另外,在沿海地區(qū),復(fù)雜水動(dòng)力條件如波浪、潮汐和風(fēng)暴潮等,均會(huì)影響近岸含水層中重金屬的運(yùn)移,進(jìn)而威脅到海岸帶資源與生態(tài)環(huán)境。然而,國(guó)內(nèi)外對(duì)于這些方面的研究甚少。

含水層中重金屬溶質(zhì)運(yùn)移動(dòng)力學(xué)模型是在考慮主要影響因素的前提下對(duì)實(shí)際過(guò)程進(jìn)行一定簡(jiǎn)化的理論研究。實(shí)際過(guò)程非常復(fù)雜,某些因素的不當(dāng)簡(jiǎn)化往往會(huì)造成模擬結(jié)果與實(shí)際過(guò)程的較大偏差,因此對(duì)該過(guò)程的研究還有待進(jìn)一步的探索和細(xì)化。

[1] 王焰新.地下水污染與防治[M].北京:高等教育出版社,2007:7-10.

[2] CETIN B,AYDILEK A H,LI L.Experimental and numerical analysis of metal leaching from fly ash-amended highway bases[J].Waste Management,2012,32(5):965-978.

[3] KESKIN T E,TOPTAS S.Heavy metal pollution in the surrounding ore deposits and mining activity:a case study from Koyulhisar (Sivas-Turkey)[J].Environmental Earth Sciences,2012,67(3):859-866.

[4] NIELSEN D,BIGGAR J.Miscible displacement in soils:I.experimental information[J].Soil Science Society of America Journal,1961,25(1):1-5.

[5] NIELSEN D,BIGGAR J.Miscible displacement:III.theoretical considerations[J].Soil Science Society of America Journal,1962,26(3):216-221.

[6] FETTER C W,FETTER J C.Contaminant hydrogeology[M].New Jersey:Prentice Hall New Jersey,1999:73-80.

[7] OBERDORFER J A.Hydrogeologic modeling of submarine groundwater discharge:Comparison to other quantita-tive methods[J].Biogeochemistry,2003,66:159-169.

[8] 孫訥正.地下水污染:數(shù)學(xué)模型和數(shù)值方法[M].北京:地質(zhì)出版社,1989:38-50.

[9] ROWE R K,BOOKER J R.A finite layer technique for modeling complex landfill history[J].Canadian Geotechnical Journal,1995,32(4):660-676.

[10] BEYER C,ALTFELDER S,DUIJNISVELD W H M,et al.Modelling spatial variability and uncertainty of cadmium leaching to groundwater in an urban region[J].Journal of Hydrology,2009,369(3/4):274-283.

[11] 錢(qián)天偉,陳繁榮,武貴斌,等.一種耦合表面絡(luò)合吸附作用的地下水溶質(zhì)遷移模型初探[J].輻射防護(hù)通訊.2003 (3):14-18.(QIAN Tianwei,CHEN Fanrong,WU Guibin,et al.An initial research on solute migration model coupled with adsorption of surface complexation in groundwater[J].Radiation Protection Bulletin,2003(3):14-18.(in Chinese))

[12] 覃榮高,曹廣祝,仵彥卿.非均質(zhì)含水層中滲流與溶質(zhì)運(yùn)移研究進(jìn)展[J].地球科學(xué)進(jìn)展.2014(1):30-41.(QIN Ronggao,CAO Guangzhu,WU Yanqing.Review of the study of groundwater flow and solute transport in heterogeneous aquifer[J].Advances in Earth Science,2014 (1):30-41.(in Chinese))

[13] 常福宣,吳吉春,薛禹群,等.考慮時(shí)空相關(guān)的分?jǐn)?shù)階對(duì)流-彌散方程及其解[J].水動(dòng)力學(xué)研究與進(jìn)展:A輯，2005(2):233-240.(CHANG Fuxuan,WU Jichun,XUE Yuqun,et al.The fractional order advection-dispersion equation withtemporal and spatial correlation and its solution[J].Journal of Hydrodynamics,2005(2):233-240.(in Chinese))

[14] XIA Y,WU Jichun.Numerical simulation of fractional advection-diffusion equation:a method to anomalous diffusion[C]//Calibration and Reliability in Groundwater Modeling:Managing Groundwater and the Environment.Wuhan:China Univ Geosciences Press,2009:433-436.

[15] 陳崇希,李國(guó)敏.地下水溶質(zhì)運(yùn)移理論及模型[M].武漢:中國(guó)地質(zhì)大學(xué)出版社,1996:26-32.

[16] 王俊,張津濤,王莉靜.地下水污染數(shù)學(xué)模型綜述[J].天津城市建設(shè)學(xué)院學(xué)報(bào),2006(4):273-277.(WANG Jun,ZHANG Jintao,WANG Lijing.Development of mathematic model of groundwater pollution[J].Journal of Tianjin Institute of Urban Construction,2006(4):273-277.(in Chinese))

[17] ACAR O,KLAMMLER H,HATFIELD K,et al.A stochastic model for estimating groundwater and contaminant discharges from fractured rock passive flux meter measurements[J].Water Resources Research,2013,49(3):1277-1291.

[18] COPPOLA A,COMEGNA A,DRAGONETTI G,et al.A stochastic texture-based approach for evaluating solute travel times to groundwater at regional scale by coupling GIS and transfer function[C]//Four decades of progress in monitoring and modeling of processes in the soil-plant-atmosphere system:Applications and challenges.Naples:Procedia Enuironmental Sciences,2013,19:711-722.

[19] 張洋.重金屬污染物在多孔介質(zhì)中的遷移模型與仿真[D].重慶:重慶大學(xué),2012.

[20] TANTEMSAPYA N,NAKSAKUL Y,WIROJANAGUD W.Mathematical modeling of heavy metals contamination from MSW landfill site in Khon Kaen,Thailand[J].Water Science and Technology,2011,64(9):1835-1842.

[21] GALITSKAYA I,POZDNIAKOVA I,TOMS L.Simulation of contaminant transport for contamination risk assessment[C]//Calibration and reliability in groundwater modeling:managing groundwater and the environment.Wuhan:IAHS Publication,2011:172-178.

[22] ZHANG Q,VOLKER R,LOCKINGTON D.Influence of seaward boundary condition on contaminant transport in unconfined coastal aquifers[J].Journal of contaminant hydrology,2001,49(3):201-215.

[23] 孫志浩,王開(kāi)章.濟(jì)寧市湖東平原淺層地下水溶質(zhì)運(yùn)移數(shù)值模擬[J].人民黃河.2010,32(4):83-84.(SUN Zhihao,WANG Kaizhang.Numerical simulation of shallow groundwater solute transportin Jining koto plain[J].Yellow River,2010,32(4):83-84.(in Chinese))

[24] 呂歲菊.有限體積法在土壤水-鹽運(yùn)動(dòng)規(guī)律數(shù)值模擬中的應(yīng)用[D].銀川:寧夏大學(xué),2005.

[25] 田玲玲.直接邊界元法及其在地下水滲流問(wèn)題中的應(yīng)用[J].西華大學(xué)學(xué)報(bào):自然科學(xué)版.2006,25(1):44-45.(TIAN Lingling.Boundary element method and its applicationin seepage of groundwater problem[J].Journal of Xihua University:Natural Science,2006,25(1):44-45.(in Chinese))

[26] 金忠.多孔介質(zhì)中水分及溶質(zhì)運(yùn)移的隨機(jī)理論[M].北京:科學(xué)出版社,2000:104-118.

[27] MAZZIA A,PUTTI M.Mixed-finite element and finite volume discretization for heavy brine simulations in groundwater[J].Journal of Computational and Applied Mathematics,2002,147(1):191-213.

[28] 王浩然,趙金熙.基于區(qū)域分解法的地下水有限元與邊界元耦合模型:淄博市王旺莊水源地地下水?dāng)?shù)值模擬[J].地質(zhì)論評(píng),2003,49(1):48-52.(WANG Haoran,ZHAO Jinxi.A groundwater flow domain decomposition model coupling the boundary and finite element methods:A groundwater numerical simulation of the Wangwangzhuang water source area of Zibo City[J].Geological Review,2003,49(1):48-52.(in Chinese))

[29] O’MALLEY D,VESSELINOV V V.Analytical solutions for anomalous dispersion transport[J].Advances in Water Resources,2014,68:13-23.

[30] 魏恒,肖洪良.地下水溶質(zhì)遷移模擬研究進(jìn)展[J].冰川凍土.2013 (6):1582-1589.(WEI Heng,XIAO Hongliang.Advances in groundwater solute transport simulating reasearch[J].Journal of Glaciology and Geocryology,2013 (6):1582-1589.(in Chinese))

[31] 程金茹,郭擇德.黃土包氣帶土壤水動(dòng)力彌散系數(shù)的測(cè)定研究[J].輻射防護(hù)通訊,2001,21(5):24-26.(CHENG Jinru,GUO Zede.Determination of hydrodynamic diffusion coefficientfor soils in loess aerated zone[J].Radiation Protection Bulletin,2001,21(5):24-26.(in Chinese))

[32] PORRO I,NEWMAN M E,DUNNIVANT F M.Comparison of batch and column methods for determining strontium distribution coefficients for unsaturated transport in basalt[J].Environmental Science & Technology,2000,34(9):1679-1686.

[33] 宋樹(shù)林,林泉.地下水彌散系數(shù)的測(cè)定[J].海岸工程,1998,17(3):61-65.(SONG Shulin,LIN Quan.Determination of the dispersive coefficientin groundwater[J].Coastal Engineering,1998,17(3):61-65.(in Chinese))

[34] QI Zhiming,FENG Shaoyuan,HELMERS M J.Modeling cadmium transport in neutral and alkaline soil columns at various depths[J].Pedosphere,2012,22(3):273-282.

[35] 劉勇,鄭軍芳,賈海紅.西北某粘土礦水動(dòng)力彌散系數(shù)的室內(nèi)測(cè)定[J].環(huán)境科學(xué)與技術(shù),2012(增刊2):319-321.(LIU Yong,ZHENG Junfang,JIA Honghai.In-lab determination of northwest clay mine hydrodynamic dispersion coefficient[J].Environmental Science & Technology,2012(Sup2):319-321.(in Chinese))

[36] 邵愛(ài)軍,劉廣明,楊勁松.土壤水動(dòng)力彌散系數(shù)的室內(nèi)測(cè)定[J].土壤學(xué)報(bào),2002 (2):184-189.(SHAO Aijun,LIU Guangming,YANG Jinsong.In-lab determination of soil hydrodynamic dispersion coefficient[J].Acta Pedologica Sinica,2002 (2):184-189.(in Chinese))

[37] BRESLER E,NAOR A.Estimating transport parameters in soils by a maximum likelihood approach[J].Soil Science Society of America Journal,1987,51(4):870-875.

[38] JIN Y,PRATT E,YATES M V.Effect of mineral colloids on virus transport through saturated sand columns[J].Journal of Environmental Quality,2000,29(2):532-539.

[39] IIDA T,UEKI K,TSUKAHARA H,et al.Point physical model of movement of ions through natural snow cover[J].Journal of Hydrology,2000,235(3):170-182.

[40] 郭建青,錢(qián)會(huì).分析一維砂柱彌散試驗(yàn)數(shù)據(jù)的反函數(shù)法[J].水利學(xué)報(bào),1999(2):43-48.(GUO Jianqing,QIAN Hui.The inverse function method for analyzing the data of one-dimension sand column dispersive test[J].Journal of Hydraulic Engineering,1999(2):43-48.(in Chinese))

[41] 李合蓮,陳家軍.阻滯系數(shù)的確定方法及其應(yīng)用[J].山東科技大學(xué)學(xué)報(bào):自然科學(xué)版,2002,21(2):103-106.(LI Helian,CHEN Jiajun.Determination and application of retardation factor[J].Journal of Shandong University of Science and Technology:Natural Science,2002,21(2):103-106.(in Chinese))

[42] MEDINA I.Determination of diffusion coefficients for supercritical fluids[J].Journal of Chromatography:A,2012,1250:124-140.

[43] DO N Y,LEE S R.Temperature effect on migration of Zn and Cd through natural clay[J].Environmental Monitoring and Assessment,2006,118(1/3):267-291.

[44] MORADKHANI D,ESKANDARI S,SEDAGHAT B,et al.A study on heavy metals mobility from Zinc plant residues in Iran[J].Physicochemical Problems of Mineral Processing,2013,49(2):567-574.

[45] CHOTPANTARAT S,ONG S K,SUTTHIRAT C,et al.Effect of pH on transport of Pb2+,Mn2+,Zn2+and Ni2+through lateritic soil:column experiments and transport modeling[J].Journal of Environmental Sciences,2011,23(4):640-648.

[46] 廖敏.鎘在土水系統(tǒng)中的遷移特征[J].土壤學(xué)報(bào),1998,35(2):179-185.(LIAO Min.Transport characteristics of cadmiumin soil-water system[J].Acta Pedologica Sinica,1998,35(2):179-185.(in Chinese))

[47] BOND W J,SMILES D.Influence of velocity on hydrodynamic dispersion during unsteady soil water flow[J].Soil Science Society of America Journal,1983,47(3):438-441.

[48] PANG L,CLOSE M,SCHNEIDER D,et al.Effect of pore-water velocity on chemical nonequilibrium transport of Cd,Zn,and Pb in alluvial gravel columns[J].Journal of Contaminant Hydrology,2002,57(3):241-258.

[49] 章明奎,鄭順安,王麗平.土壤中顆粒狀有機(jī)質(zhì)對(duì)重金屬的吸附作用[J].土壤通報(bào),2008,38(6):1100-1104.(ZHANG Mingkui,ZHENG Shun’an,WANG Liping.Adsorption of heavy metals by soil particulate organic matter[J].Chinese Journal of Soil Science,2008,38(6):1100-1104.(in Chinese))

[50] 余貴芬,蔣新,孫磊,等.有機(jī)物質(zhì)對(duì)土壤鎘有效性的影響研究綜述[J].生態(tài)學(xué)報(bào),2002,22(5):770-776.(YU Guifen,JIANG Xin,SUN Lei,et al.A review for effect of organic substances on the availability of cadmium in soils[J].Acta Ecologica Sinica,2002,22(5):770-776.(in Chinese))

[51] BEKHIT H M,HASSAN A E.Two-dimensional modeling of contaminant transport in porous media in the presence of colloids[J].Advances in Water Resources,2005,28(12):1320-1335.

[52] DE JONGE L W,KJAERGAARD C,MOLDRUP P.Colloids and colloid-facilitated transport of contaminants in soils[J].Vadose Zone J,2004,3(2):321-325.

[53] CORAPCIOGLU M Y,JIANG S.Colloid-facilitated groundwater contaminant transport[J].Water Resources Research,1993,29(7):2215-2226.

[54] SAIERS J E,HORNBERGER G M.The role of colloidal kaolinite in the transport of cesium through laboratory sand columns[J].Water Resources Research,1996,32(1):33-41.

[55] FLURY M,QIU H.Modeling colloid-facilitated contaminant transport in the vadose zone [J].Vadose Zone J,2008,7(2):682-697.

[56] SAIERS J E,HORNBERGER G M.The influence of ionic strength on the facilitated transport of cesium by kaolinite colloids[J].Water Resources Research,1999,35(6):1713-1727.

[57] BAUMANN T,FRUHSTORFER P,KLEIN T,et al.Colloid and heavy metal transport at landfill sites in direct contact with groundwater[J].Water Research,2006,40(14):2776-2786.

[58] 尚熳廷,馮杰,劉佩貴.大孔隙對(duì)土壤水動(dòng)力彌散系數(shù)影響的實(shí)驗(yàn)研究[J].灌溉排水學(xué)報(bào),2009 (5):52-54,66.(SHANG Manting FENG Jie,LIU Peigui.Effect of macroporous on soil hydrodynamic dispersion coefficient[J].Journal of Irrigation and Drainage,2009(5):52-54,66.(in Chinese))

[60] 李國(guó)敏,陳崇希.空隙介質(zhì)水動(dòng)力彌散尺度效應(yīng)的分形特征及彌散度[J].地球科學(xué):中國(guó)地質(zhì)大學(xué)學(xué)報(bào),1995,20(4):405-409.(LI Guomin,CHEN Chongxi.Fractal geometry and estimation of scale-dependent dispersivity in geologic media[J].Earth Science:Journal of China University of Geosciences,1995,20(4):405-409.(in Chinese))

[61] 王永森,陳建生,陳亮.考慮彌散系數(shù)尺度效應(yīng)的溶質(zhì)運(yùn)移模型研究[J].人民黃河,2008,30(11):60-62.(WANG Yongsen,CHEN Jiansheng,CHEN Liang.Research on solute transport model with scale-dependent dispersion coefficients[J].Yellow River,2008,30(11):60-62.(in Chinese))

Advances of dynamics model for heavy metal solute transport in aquifer

TAO An, LIU Shuguang, LOU Sha, DAI Chaomeng

(CollegeofCivilEngineering,TongjiUniversity,Shanghai200092,China)

On the basis of summarizing the study of the heavy metal pollution at home and abroad, the authors of this paper reviewed the development of dynamics model for heavy metal solute transport in aquifer, and summed up the calculation methods of the model.The authors also discussed the determination of model parameters and the influences of temperature, pH, pore-water velocity, organic matter and colloid, soil heterogeneity on the determination of parameters.At last, the authors pointed out the main problems existing in the study on the transport process of heavy metal solute in aquifer, and looked forward to the study that needs focusing on in the field.

groundwater; aquifer; heavy metal transport; numerical simulation; dynamics model; soil heterogeneity

10.3880/j.issn.1004-6933.2015.04.002

國(guó)家自然科學(xué)基金(41372240)

陶安(1990—),男,碩士研究生,研究方向?yàn)榈叵滤當(dāng)?shù)值模擬。E-mail:antao1990@163.com

劉曙光,教授。E-mail:liusgliu@#edu.cn

P641.69

A

1004-6933(2015)04-0008-07