潘茜雯, 胡韜, 孫洪廣*

(1.河海大學(xué)力學(xué)與材料學(xué)院，水文水資源與水力工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室，南京 210098；2.南京市江南小化工集中整治工作現(xiàn)場指揮部，南京 210098)

隨著中國社會(huì)經(jīng)濟(jì)的發(fā)展，土壤污染已經(jīng)成為一個(gè)不容忽視的環(huán)境問題，對人們的生活和身體健康產(chǎn)生了嚴(yán)重的危害[1-4]。建國初期建造的化工廠因?yàn)榻?jīng)濟(jì)、效能、技術(shù)等問題，造成了不同程度的土壤污染，這些污染物隨著時(shí)間的推移不斷向土壤入滲[5]，污染范圍逐步擴(kuò)大，影響周圍的水土環(huán)境[6-7]。環(huán)境中廣泛存在著砷元素及其化合物，其中三氧化二砷是劇毒物[8]。因此對污染情況進(jìn)行準(zhǔn)確評(píng)估，對指導(dǎo)后期污染控制和修復(fù)工作具有重要的科學(xué)意義和現(xiàn)實(shí)意義。

污染物在土壤中的運(yùn)移是由對流、分子擴(kuò)散和機(jī)械彌散[9-12]3個(gè)物理過程以及溶質(zhì)在運(yùn)移過程中所發(fā)生的化學(xué)、物理及其他過程綜合作用的結(jié)果[13-14]。郭鐘群等[15]分析了各種模型的特點(diǎn)及其在離子型稀土溶質(zhì)運(yùn)移過程的適用性；探討了色譜塔板理論模型、地理信息系統(tǒng)應(yīng)用模型、人工神經(jīng)網(wǎng)絡(luò)分析離子型稀土溶質(zhì)運(yùn)移的可行性，為揭示離子型稀土浸取過程的溶質(zhì)遷移規(guī)律提供思路。張?zhí)飙^等[16]將降雨入滲與重金屬所在飽和非飽和土壤作為一個(gè)整體考慮，通過將入滲模型、對流-彌散模型和地下水稀釋模型進(jìn)行耦合，達(dá)到計(jì)算基于保護(hù)地下水的土壤風(fēng)險(xiǎn)值的目的。王先穩(wěn)[17]針對垃圾填埋場滲濾液入滲(點(diǎn)源)、地表河流入滲(線源)、農(nóng)田入滲(面源)3種污染情景分別建立了包氣帶土壤及淺層地下水中污染物遷移數(shù)值模型進(jìn)行多環(huán)芳烴(PAHs)遷移規(guī)律探究,由此建立了研究區(qū)土壤-地下水中PAHs的污染概念模型。

用于模擬飽和-非飽和多孔介質(zhì)中水分運(yùn)移和溶質(zhì)運(yùn)移的數(shù)值模型[18]，在土壤污染物運(yùn)移、地下水污染風(fēng)險(xiǎn)評(píng)價(jià)方面得到了廣泛運(yùn)用[19]，可對土壤污染物擴(kuò)散過程進(jìn)行準(zhǔn)確分析，在時(shí)間維度上可以預(yù)測未來土壤污染物的濃度和擴(kuò)散趨勢[20]。相比傳統(tǒng)的調(diào)查打孔和定點(diǎn)監(jiān)測更節(jié)約成本，能夠直接指導(dǎo)污染防治。劉玲等[21]應(yīng)用FEFLOW建立三維土壤-地下水模型，定量描述六價(jià)鉻在土壤-地下水中的遷移規(guī)律，并討論阻滯系數(shù)和反應(yīng)常數(shù)對溶質(zhì)運(yùn)移的影響。侯永莉等[22]利用Geostudio軟件模擬降雨、蒸發(fā)、蒸騰對重金屬污染物遷移的影響，為環(huán)境地質(zhì)評(píng)價(jià)和污染預(yù)測防治提供理論依據(jù)。柳山等[23]利用數(shù)值模擬方法進(jìn)行酸雨淋濾數(shù)值模擬，研究酸性降雨淋濾土壤的過程中，不同層位土壤中污染物的運(yùn)移特征及規(guī)律。董炎青等[24]利用有限元(computational fluid dynamics，CFD)軟件模擬分析甲苯在土壤中的泄漏過程，研究土壤孔隙度和含水率對甲苯擴(kuò)散的影響及運(yùn)移規(guī)律。因此為了保證土壤污染不影響到周圍居民的生產(chǎn)生活，利用地質(zhì)勘查、數(shù)值模擬等綜合手段建立一套土壤重金屬污染預(yù)測和評(píng)估模型，以期為土壤污染防控和修復(fù)提供科學(xué)指導(dǎo)。

1 研究區(qū)域概況

南京某化工廠占地超過27×104m2，場地位于長江南岸，靠近長江。地勢東高西低，南北高中間低。污染土壤深度2 m，除場地東南角高程15～18 m，其余地方高程10～15 m。場地周圍水系主要為場地西側(cè)的河道，其主要作用為防洪排澇，岸坡寬29.7 m，水面寬10 m，河底高程3.08～4.02 m，水面高程5.29～5.06 m，水體由南向北流動(dòng)。

2 研究方法

2.1 地質(zhì)勘查

為查明地塊范圍內(nèi)的地層結(jié)構(gòu)、地下水類型和地下水位，并提供各層土的含水率、容重、孔隙率、滲透系數(shù)(水平和垂直)和有機(jī)質(zhì)含量，土壤的物理力學(xué)指標(biāo)如表1所示。

表1 土壤的物理力學(xué)指標(biāo)

2.2 污染源調(diào)查

根據(jù)現(xiàn)有資料，首先對地塊土壤污染狀況進(jìn)行了初步分析。自20世紀(jì)50年代年以來，地塊一直為金陵石化化工一廠使用，結(jié)合工藝流程和“三廢”排放情況，主要產(chǎn)污環(huán)節(jié)有聚丙烯生產(chǎn)線、醫(yī)藥/凡士林生產(chǎn)線、白油生產(chǎn)線等工序。對土壤的潛在污染源為：各裝置區(qū)，罐區(qū)、鍋爐房，產(chǎn)生的主要污染物有石油烴類、重金屬等?；S和熱電廠臨近本地塊，為該地塊潛在污染源，通過資料分析，需要考慮的污染物有石油烴類、苯系物、重金屬、多環(huán)芳烴等[25]。

化工廠生產(chǎn)期間產(chǎn)生的污染物遷移途徑有：①廢水污水管道、罐組破損滲漏遷移進(jìn)入土壤和地下水；②物料在堆放區(qū)域通過地表裂縫滲漏進(jìn)入土壤和地下水；③物料在裝卸站區(qū)域發(fā)生泄漏進(jìn)入土壤和地下水。

根據(jù)資料分析、人員訪談、現(xiàn)場勘測結(jié)果，開展了該地塊的土壤采樣工作，對土壤樣品進(jìn)行檢測，鉛采樣篩選值為400 mg/kg，砷采樣數(shù)據(jù)篩選值為20 mg/kg，4-氯苯胺采樣數(shù)據(jù)篩選值為2.7 mg/kg。通過檢測，按照《土壤環(huán)境質(zhì)量 建設(shè)用地土壤環(huán)境風(fēng)險(xiǎn)管控標(biāo)準(zhǔn)(試行)》(GB36600—2018)，這片土壤重金屬含量已嚴(yán)重超過國家標(biāo)準(zhǔn)[26-27]。

2.3 數(shù)學(xué)模型

2.3.1 方法與目的

在自然過程作用下，土壤中的污染物會(huì)發(fā)生遷移和轉(zhuǎn)化，如降雨和地下水流動(dòng)會(huì)使污染物發(fā)生遷移和轉(zhuǎn)化等。特別是雨季降雨量大，易造成污染物進(jìn)一步轉(zhuǎn)移到下層土壤，造成地塊更大范圍的土壤污染[28]。通過數(shù)值模擬，模擬土壤重金屬在自然天氣中的遷移過程，有利于控制污染物的擴(kuò)散，減少和控制被污染土壤的分布范圍，降低修復(fù)成本。

2.3.2 水動(dòng)力學(xué)模型

土壤水采用Richards對流-彌散方程和伽遼金有限元法進(jìn)行數(shù)值求解[29]，基本方程為

(1)

式(1)中：θ為土壤體積含水量；t為時(shí)間；K(h)為非飽和土壤導(dǎo)水率；h為土壤水頭；S為吸附項(xiàng)。

2.3.3 土壤溶質(zhì)運(yùn)移模型

溶質(zhì)在土壤中的運(yùn)移受對流和水動(dòng)力彌散作用影響[30]，采用對流-彌散方程模型來描述溶質(zhì)的運(yùn)移，其表達(dá)式為[9]

(2)

式(2)中：c為濃度；qi為水流通量；DT為橫向彌散系數(shù)；DL為縱向彌散系數(shù)；ρ為土壤容重；S為溶質(zhì)吸附濃度。

2.4 初始條件和邊界條件

溶質(zhì)運(yùn)移初始條件為

c(x,y,z,t)=c0(x,y,z)，

x≥0；y≥0；z≥0；t=0

(3)

式(3)中:c0為初始溶質(zhì)濃度。

邊界條件如圖1所示，上邊界A′B′C′D′內(nèi)除區(qū)域R以外的區(qū)域設(shè)為大氣邊界條件，考慮土壤被蒸發(fā)情況。EF處有水流通過，設(shè)為可變水頭邊界條件，模擬河流通過左邊界A′D′DA。圓形區(qū)域R處設(shè)為固定水頭邊界條件，用于模擬上層土壤污染物有多天脈沖。其他剩余邊界A′ABB′、B′BCC′、C′CDD′、ABCD和A′ADD′除EF以外區(qū)域?yàn)榱阃窟吔鐥l件。

圖1 模擬區(qū)域形狀及污染源位置示意圖

2.5 模型參數(shù)

模型參數(shù)通過場地勘察、室內(nèi)土工試驗(yàn)以及參考文獻(xiàn)獲得。根據(jù)土壤顆粒大小分析確定土壤分布情況，結(jié)合土工試驗(yàn)和資料數(shù)據(jù)庫確定土壤水力學(xué)參數(shù)[31-32]，土壤水力特征系數(shù)及相關(guān)參數(shù)如表2所示，溶質(zhì)運(yùn)移參數(shù)如表3所示。

表2 土壤水力特征系數(shù)及相關(guān)參數(shù)

表3 溶質(zhì)運(yùn)移參數(shù)

3 結(jié)果與分析

為了研究不同因素對土壤重金屬運(yùn)移產(chǎn)生的影響，針對不同脈沖時(shí)間、不同溶質(zhì)濃度、不同土壤飽和含水率等因素進(jìn)行研究，將污染場地地層概化為4層，如表4所示。在0.5、1、1.5、2 m處設(shè)置4個(gè)觀察點(diǎn)。

表4 研究區(qū)地層概化結(jié)果

3.1 不同脈沖時(shí)間對溶質(zhì)運(yùn)移的影響

根據(jù)數(shù)學(xué)物理模型和相關(guān)初始邊界數(shù)值，可以模擬土壤在50、70、100 d脈沖時(shí)間下重金屬離子的運(yùn)移過程。

從圖2可以觀測到，污染源脈沖時(shí)間的變化對砂壤土、壤土影響最為明顯。比較0.5 m處砂壤土中污染物砷的濃度變化曲線[圖2(a)]，脈沖時(shí)間50、70、100 d時(shí)，均在第2天到達(dá)0.5 m處，在前62 d溶質(zhì)濃度擴(kuò)散行為相同，之后按照不同的脈沖時(shí)間繼續(xù)擴(kuò)散或下降；砂壤土中溶質(zhì)經(jīng)過峰值后通過土壤吸附等作用有下降的趨勢，脈沖時(shí)間越短，下降越早；脈沖時(shí)間100 d時(shí)，砂壤土中溶質(zhì)趨于飽和。比較1 m處壤土中污染物砷的濃度變化曲線[圖2(b)]，經(jīng)過12 d到達(dá)1 m處，在78 d之前溶質(zhì)擴(kuò)散相同，隨著脈沖時(shí)間的增加，土壤溶質(zhì)擴(kuò)散范圍更大，土壤的溶質(zhì)濃度更高，脈沖時(shí)間70 d中污染物的濃度比脈沖時(shí)間50 d高0.056 mg/cm3。比較1.5 m處淤泥質(zhì)壤土中污染物砷的濃度變化曲線[圖2(c)]，脈沖時(shí)間50 d時(shí)1.5 m深處淤泥質(zhì)壤土溶質(zhì)濃度略低為0.133 1 mg/cm3，脈沖時(shí)間70 d與100 d時(shí)，淤泥質(zhì)壤土中污染物濃度均為0.133 5 mg/cm3。比較2 m處黏壤土中污染物砷的濃度變化曲線[圖2(d)]，2 m深處黏壤土中污染物濃度相同，均為0.01 mg/cm3。因此脈沖時(shí)間越長，對深層土壤的影響越大，并推測隨著時(shí)間的增加及溶質(zhì)的不斷運(yùn)移，中層土壤、深層土壤也將受到影響。

圖2 不同脈沖時(shí)間下0.5、1、1.5、2 m處黏壤土中污染物砷的濃度變化曲線(線性吸附、污染區(qū)域?yàn)楣潭ㄋ^邊界、考慮土壤蒸發(fā))

3.2 不同土壤飽和含水率對溶質(zhì)運(yùn)移的影響

考察不同的土壤飽和含水率條件下污染物的運(yùn)移過程。主要模擬結(jié)果如圖3～圖5所示。

線性吸附、污染區(qū)域?yàn)楣潭ㄋ^邊界、考慮土壤蒸發(fā)

線性吸附、污染區(qū)域?yàn)楣潭ㄋ^邊界、考慮土壤蒸發(fā)

線性吸附、污染區(qū)域?yàn)楣潭ㄋ^邊界、考慮土壤蒸發(fā)

從圖3中可以看出，土壤中的溶質(zhì)到達(dá)各層土壤更慢，土壤深度越深，差距越明顯。在0.5 m的砂壤土中，具有較高飽和土壤含水率的土壤中溶質(zhì)濃度在第84天達(dá)到峰值，濃度為3.5 mg/cm3，具有較低飽和土壤含水率的土壤中溶質(zhì)濃度在第90天達(dá)到峰值，濃度為3 mg/cm3，較高土壤飽和含水率的土壤溶質(zhì)濃度更高且更快達(dá)到峰值。在1 m的壤土中，具有較高飽和土壤含水率的土壤中溶質(zhì)第4天到達(dá)1 m，溶質(zhì)在第100天濃度為1.33 mg/cm3，具有較低飽和土壤含水率的土壤中溶質(zhì)第8天到達(dá)1 m，溶質(zhì)在第100天濃度為0.75 mg/cm3。在1.5 m的淤泥質(zhì)壤土中，具有較高飽和土壤含水率的土壤中溶質(zhì)第22.5天到達(dá)1.5 m，溶質(zhì)在第100天濃度為0.133 5 mg/cm3，具有較低飽和土壤含水率的土壤中溶質(zhì)第26.5天到達(dá)1.5 m，溶質(zhì)在第100天濃度為0.024 5 mg/cm3。在2 m的黏壤土中，具有較高飽和土壤含水率的土壤中溶質(zhì)第36天到達(dá)2 m，溶質(zhì)在第100天濃度為0.01 mg/cm3，具有較低飽和土壤含水率的土壤中溶質(zhì)第40天到達(dá)2 m，溶質(zhì)在第100天濃度為0.000 3 mg/cm3。比較圖3可知，較高的土壤飽和含水率促進(jìn)污染物在土壤中的垂向運(yùn)移，提前4～6 d到達(dá)更深層土壤，污染物濃度提高2～6倍。

圖4為脈沖時(shí)間70 d時(shí)具有高飽和土壤含水率土壤不同深度的水平向溶質(zhì)運(yùn)移二維圖。圖5為脈沖時(shí)間70 d時(shí)具有低飽和土壤含水率土壤不同深度的水平向溶質(zhì)運(yùn)移二維圖。圖4中，當(dāng)土壤具有更高的飽和土壤含水率時(shí)，溶質(zhì)在橫向與縱向擴(kuò)散更加均勻。圖5中，當(dāng)土壤具有較低的飽和土壤含水率時(shí)，溶質(zhì)在縱向擴(kuò)散更快，溶質(zhì)運(yùn)移區(qū)域偏狹長。綜合考慮圖4、圖5可知，它們整體向有河流的方向偏移，在高飽和土壤含水率土壤中，溶質(zhì)擴(kuò)散更快，污染范圍更大。

4 結(jié)論

化工廠土壤重金屬的運(yùn)移受到很多不同因素的影響，通過建立脈沖時(shí)間50、70、100 d的土壤溶質(zhì)運(yùn)移模型以及不同飽和土壤含水率的模型，分析脈沖時(shí)間、土壤飽和含水率對重金屬運(yùn)移的影響，得出如下結(jié)論。

(1)脈沖時(shí)間對土壤溶質(zhì)運(yùn)移具有促進(jìn)作用，對淺層土壤作用明顯。隨著脈沖時(shí)間的不斷增加，對深層土壤的影響愈發(fā)明顯。對于受污染的土壤，必須采取措施減少污染物的滲漏時(shí)間，防止污染物在土壤中的進(jìn)一步擴(kuò)散，減少對更深層土壤的污染。

(2)具有較高土壤飽和含水率的土壤對溶質(zhì)的擴(kuò)散起促進(jìn)作用。由于土壤特性，較高土壤飽和含水率對橫向擴(kuò)散更加明顯，污染范圍更大，促進(jìn)污染物的下滲。降低土壤的飽和含水率對治理污染物的運(yùn)移有更好效果。