重金屬在飽和-非飽和分層土壤中遷移過程的三維數(shù)值模擬

2022-10-20 00:57:46劉歡許諾孫洪廣胡韜

計算機輔助工程 2022年3期

劉歡，許諾，孫洪廣，胡韜

(1.河海大學力學與材料學院水文水資源與水利工程科學國家重點實驗室，南京 210098;2.南京市燕子磯片區(qū)整治開發(fā)有限責任公司，南京 211100)

0 引言

隨著經(jīng)濟的快速發(fā)展，冶金、采礦、農(nóng)業(yè)生產(chǎn)等各種人類活動，導(dǎo)致大量污染物進入土壤，使得土壤中污染物的含量遠高于自然含量。近年來，我國對生態(tài)環(huán)境問題越來越重視，老舊工業(yè)企業(yè)和廢棄場地污染修復(fù)工程越來越多。重金屬離子的污染與防治一直是土壤修復(fù)的重點與難點，重金屬污染物進入土壤不會被降解，長期富集在土壤中，具有吸附性強、易積累、隱蔽性強等特點。

目前，國內(nèi)外已有較多關(guān)于重金屬離子在土壤中遷移規(guī)律的研究，研究焦點主要包括重金屬離子在土壤中的吸附特性、遷移機制、分布規(guī)律等。王玉雪等通過土柱吸附實驗發(fā)現(xiàn)土壤初始pH值越低或溶液初始濃度越高，重金屬鉻Cr(Ⅵ B)穿透時間與達到吸附飽和的時間越早。林青等通過室內(nèi)置換實驗研究銅、鎘、鋅、鉛在石英砂中的遷移規(guī)律，認為靜態(tài)吸附實驗得到的吸附參數(shù)不能直接應(yīng)用于重金屬運移模擬。王國賢等測試土壤不同深度的重金屬元素,認為在垂直方向隨著深度的增加金屬含量減少。張茹吉等認為銅、鎘、鉛在土壤中垂直分布的規(guī)律隨土壤質(zhì)地不同而存在差異，黏壤土吸附重金屬的能力最強，砂壤土中重金屬的運移能力最強。JACQUES等分析土壤中重金屬元素的遷移，發(fā)現(xiàn)土壤含水量和水通量是影響重金屬元素遷移率的重要參數(shù)。DONG等對重金屬在垃圾填埋場的遷移行為進行數(shù)值模擬，結(jié)果認為在土壤和地下含水區(qū)最敏感的參數(shù)是阻滯因子和滲透系數(shù)。CHEN等發(fā)現(xiàn)重金屬溶解過程的轉(zhuǎn)化、重金屬離子隨地下水的遷移和重金屬在多孔介質(zhì)中的遷移是重金屬在土壤中遷移的3個重要方面。這些對層狀土中重金屬溶質(zhì)運移規(guī)律的研究主要集中于室內(nèi)的小型土柱，對野外較大場地一般簡化為均質(zhì)土壤。

本文以某實際污染場地為背景，建立飽和-非飽和層狀土壤溶質(zhì)運移模型，模擬分析含類重金屬砷(As)溶液連續(xù)通入情況下，不同土體構(gòu)型、邊界水位及污水入滲通量對土壤中As遷移距離及濃度分布的影響，從而更全面地了解重金屬的遷移規(guī)律，為實際污染場地土壤和地下水修復(fù)與防治提供參考。

1 數(shù)學模型的建立

1.1 飽和-非飽和水流運動數(shù)學模型

三維飽和-非飽和多孔介質(zhì)的水流運動基本方程為Richards方程

(1)

式中：為總水勢；為土體含水率；、和分別為、、主方向的非飽和滲透系數(shù)，是含水率或基質(zhì)勢的函數(shù)；為初始基質(zhì)勢；為水頭邊界條件；為通量邊界條件。

1.2 溶質(zhì)運移數(shù)學模型

離子在土壤中的遷移作用，可以分為物理、化學和生物作用。本文只考慮As在土壤中受對流、彌散和吸附作用，因此

,=,,;=;=;=

(2)

(3)

式中：為土壤溶質(zhì)含量；為土壤吸附量；為土壤干體積質(zhì)量即容重；為源匯項；為彌散系數(shù)；為達西流速；、分別為動力學吸附位點的1階吸附系數(shù)和1階解吸附系數(shù)；為分配系數(shù)。

2 幾何建模與邊界條件

本文數(shù)值實驗的土壤參數(shù)與土體構(gòu)型參照南京市某廢棄化工廠，設(shè)計幾何模型為100 m×65 m×12 m的三維層狀矩形地塊，見圖1。

圖1 土體三維幾何模型

該地塊垂向為砂土夾粉質(zhì)黏土的土體構(gòu)型，土壤縱剖面由上至下均勻分為3層。西側(cè)低水位為8 m，東側(cè)高水位為9 m，南、北兩側(cè)水頭沿程均勻坡降，在土壤飽和帶中形成一定的地下水流。底部設(shè)置為不透水邊界，大氣降水以2 mm/d、含As溶液以3 mm/d的入滲通量由表面流入。砂壤土與粉質(zhì)黏土的土壤參數(shù)見表1。

表 1 土壤參數(shù)表

3 計算結(jié)果及分析

滲透系數(shù)與彌散系數(shù)增加均能促進溶質(zhì)的運移，而土壤分配系數(shù)表征土壤介質(zhì)的吸附能力，分配系數(shù)越大，溶質(zhì)遷移過程中受到的阻滯作用越強。由2種土質(zhì)的土壤參數(shù)可知，與砂壤土相比，粉質(zhì)黏土在對流、彌散和吸附作用方面對As的運移均有極強的阻滯作用。

3.1 不同土體構(gòu)型對As遷移的影響

基于實際工程條件，設(shè)計土體構(gòu)型1和土體構(gòu)型2，見圖2。在水流條件不變的情況下，連續(xù)通入含As污水20 a，分析As在2種土體構(gòu)型中的遷移規(guī)律。

3.1.1 不同土體構(gòu)型土壤水中As運移規(guī)律

通量滲漏20 a后不同土體構(gòu)型土壤水中As的遷移結(jié)果和含量分布見圖3。含As污水定通量滲漏20 a，As沒有穿透第二層土壤，As在2個土體構(gòu)型中的運移區(qū)別很明顯，在土體構(gòu)型1中的水平運移距離遠大于在土體構(gòu)型2中，而構(gòu)型1中的垂向運移距離則相對較小。這主要是因為土壤的不同分層結(jié)構(gòu)影響溶質(zhì)運移擴散過程。

在圖3(a)土體構(gòu)型1中，由于表層砂壤土對As的阻滯作用較小，As在砂壤土中的水平遷移距離較大。而As穿透毛細水帶進入地下含水層后，在東、西水位差的作用下，As向下游遷移速度加快，且由于含水率的增加，土壤的有效滲透系數(shù)和彌散度也隨之增加，水平方向最遠遷移距離為42.5 m，下游遷移距離達25.0 m。由于粉質(zhì)黏土層對As的作用較強，As由砂壤土進入粉質(zhì)黏土后，水平污染范圍沿深度方向迅速縮小，垂向遷移距離為6.5 m。

(a)土體構(gòu)型1

(a)遷移結(jié)果

在圖3(a)土體構(gòu)型2中，由于表層粉質(zhì)黏土層對As的吸附性強，As在該土層中的運移速度緩慢。As在土壤表面水平最遠遷移距離僅為15.0 m，相較于在土體構(gòu)型1中表面水平最遠運移距離少21.2 m，溶質(zhì)在土體構(gòu)型2表層土壤中主要表現(xiàn)為垂向運移。As由土壤包氣帶穿過毛細管區(qū)進入含水層后，粉質(zhì)黏土層飽和滲透系數(shù)小，As在含水層中隨水流向下游運移速度變化不明顯。在土層交界處，As進入砂壤土后迅速往四周彌散，且在水力梯度作用下開始往下游遷移。As水平最遠運移距離為25.8 m，垂向運移距離為7.2 m，與土體構(gòu)型1相比，其水平方向運移距離小16.7 m，垂向多運移0.7 m。這是因為土體構(gòu)型1中粉質(zhì)黏土層阻滯As的繼續(xù)下滲，垂向運移距離僅為6.5 m。

圖3(b)中土壤水As含量在不同土體構(gòu)型中均沿深度方向減小，且土體構(gòu)型1中土壤水As含量普遍高于土體構(gòu)型2中。在毛細水帶中，土體構(gòu)型1的土壤水As含量沿深度方向減小速度略有波動，但在土體構(gòu)型2中沒有明顯變化。這是因為粉質(zhì)黏土中毛細水帶較厚，=1=0.5 m，土壤水中As含量在毛細管區(qū)得到很好的過渡。在兩土層交界處，土體構(gòu)型1中As穿過土層交界面后，含量在粉質(zhì)黏土層的2.5 m內(nèi)由2.33 mg/L迅速降低到0.02 mg/L，而在土體構(gòu)型2中溶質(zhì)穿透土層交界面，As在砂壤土層運移3.2 m才從0.19 mg/L降為0.03 mg/L。

3.1.2 不同土體構(gòu)型土壤As吸附量分布規(guī)律

不同土體構(gòu)型土壤吸附量分布見圖4。土體構(gòu)型2表層土壤As吸附量遠大于土體構(gòu)型1表層，土壤As的污染范圍集中在污染源附近，土壤大量地吸收污水滲漏中的As。土體構(gòu)型1中As在土壤表層砂壤土中大范圍污染，但是進入粉質(zhì)黏土層之后水平污染范圍迅速收縮，As繼續(xù)往水平方向與深層土壤的污染受阻。

(a)吸附量結(jié)果

據(jù)吸附平衡公式可知，土壤As吸附量直接與土壤水As含量及土壤質(zhì)地(分配系數(shù))相關(guān)，且在非飽和土壤中，含水率的增加也會導(dǎo)致土壤As吸附量增加。在圖4中，土體構(gòu)型1土壤As的吸附量分布可明顯分為非飽和砂壤土、飽和砂壤土與飽和粉質(zhì)黏土2個部分，均表現(xiàn)為距離污染源越遠土壤As吸附量越低，垂向則在毛細水帶及土層交界處發(fā)生轉(zhuǎn)折變化。毛細水帶含水率沿深度不斷增加，土壤As的吸附量也隨之增加，在毛細水帶中吸附量由10.9 mg/kg逐漸增加到27.7 mg/kg。在土層交界處，由于兩側(cè)土壤參數(shù)不同，但本文土層模型未設(shè)置類似飽和-非飽和土壤間存在的毛細水過渡帶，因此土壤As吸附量在該交界面發(fā)生含量突變，土壤中As吸附量在土層交界處由22.1 mg/kg驟增到66.2 mg/kg。

由圖4可知，土體構(gòu)型2中土壤As吸附量變化主要在土層交界面處，土壤As吸附量由粉質(zhì)黏土層的5.79 mg/kg迅速降為砂壤土層的1.70 mg/kg。在飽和-非飽和土壤過渡區(qū)，毛細管區(qū)較厚且粉質(zhì)黏土層的持水性能好，土壤含水率變化區(qū)間不大，為0.37～0.45，從而土壤As吸附量在該區(qū)域變化不明顯。

3.2 不同水位對As遷移的影響

改變模型東側(cè)(高)邊界水位，分別取值為8.5、9.0與9.5 m，分析As定通量通入10 a在土壤中的運移和濃度分布的變化。

3.2.1 不同水位土壤水中As運移規(guī)律

不同水位土壤水中As的遷移和含量變化見圖5。水位增加，As在土壤中的水平遷移距離增加，垂直方向遷移距離均為6.5 m，沒有顯著變化。污染源中心垂向土壤水As含量隨增加而降低。飽和土壤中水力梯度隨著水位的增加而增加，地下水流往下游運移速度加快，As隨水流向下游運移距離增加：當為9.5 m時，水平方向最遠運移距離為40.4 m；當=8.5 m時，水平運移僅24.5 m。在彌散作用下，非飽和土壤水中As水平遷移距離也增加，As在土壤表面水平遷移距離由20.9 m增加到32.5 m。

(a)遷移結(jié)果

由于第二層粉質(zhì)黏土層的阻滯作用，水位在8.5～9.5 m范圍內(nèi)變化時，地下水垂向速度變化較小，As在垂直方向運移距離大致保持在6.5 m，沒有明顯變化。隨著增大，As更多地在水平方向擴散，而溶質(zhì)通量不變，因此污染源處土壤水As含量在垂向相應(yīng)減少，圖5(b)中土壤表面污染源中心含量由8.45 mg/L降為7.35 mg/L。

3.2.2 不同水位時土壤As吸附量分布規(guī)律

不同水位時土壤As吸附量分布見圖6。東側(cè)(高)水位增加，土壤對As的吸附量在水平方向變化規(guī)律與土壤水中As含量變化一致，土壤As污染范圍增加，而土壤中As吸附量最大值減小，更多的As溶質(zhì)隨水流向更遠處運移。

(a)吸附量結(jié)果

由圖6(b)可知,邊界水位的變化對土壤As吸附量的影響表現(xiàn)為污染范圍的變化，及毛細水帶變動導(dǎo)致土壤As吸附量在垂向的變化。土壤As吸附量最大值隨的增加而減小，但邊界水位的增加導(dǎo)致污染源處地下水位抬高，毛細水帶也相應(yīng)上移，該區(qū)域土壤As吸附量增加。

3.3 不同污水通量對As遷移的影響

改變污染源的污水入滲通量，分別取值為2、3和4 mm/d，其他條件不變，含As溶液定通量連續(xù)通入10 a，分析污水入滲通量對As在土壤中遷移的影響。

3.3.1 不同污水通量土壤水中As運移規(guī)律

不同污水入滲通量土壤水As的遷移結(jié)果和含量變化見圖7。污水入滲通量增加，土壤水As污染范圍增加較小，主要為污染源處土壤水As濃度的增加。含As溶液入滲通量增加，流入土壤中As的總量增加，土壤水As含量增加。土壤表層為彌散度大的砂壤土，在彌散作用下水平遷移距離增加。污水入滲通量=4 mm/d時，污染源中心土壤表面含量為8.88 mg/L，水平方向遷移距離為35.0 m，與=2 mm/d時相比，污染源中心表面含量增加1.45 mg/L，水平方向遷移距離增加4.5 m。由于粉質(zhì)黏土層的阻滯作用，垂向遷移距離幾乎不變。