曾 頡，顧正聰，黃 勇

黏土貫穿裂隙污染物遷移的實(shí)驗(yàn)和數(shù)值模擬研究

曾 頡1，顧正聰2，黃 勇1

（1.河海大學(xué) 地球科學(xué)與工程學(xué)院，江蘇 南京 210098；2.云南湖柏環(huán)?？萍加邢薰?，云南 昆明 650034）

為研究裂隙黏土中的溶質(zhì)運(yùn)移，研制一套適用于裂隙黏土中溶質(zhì)運(yùn)移的實(shí)驗(yàn)裝置，通過測定實(shí)驗(yàn)裝置內(nèi)示蹤劑溶液的電導(dǎo)率來監(jiān)測溶質(zhì)的運(yùn)移規(guī)律。實(shí)驗(yàn)結(jié)果發(fā)現(xiàn)模型底部監(jiān)測孔處的滲出液濃度在短時(shí)間內(nèi)就趨于穩(wěn)定達(dá)到最大，表明黏土中裂隙的存在會(huì)使溶質(zhì)順著裂隙快速穿過黏土層，形成優(yōu)先流現(xiàn)象，從而對黏土層下部的水體造成污染。采用數(shù)值方法對室內(nèi)試驗(yàn)進(jìn)行模擬分析，討論了提高水頭和改變巖土體兩種方法來進(jìn)行預(yù)測，結(jié)果顯示黏土層上部水頭變化對黏土中的溶質(zhì)濃度影響不大，且不同土體中如均出現(xiàn)貫穿型裂隙，那么污染物進(jìn)入含水層的時(shí)間基本一致，防滲失效，這時(shí)改變防滲材料沒有意義。

裂隙黏土；飽和滲流實(shí)驗(yàn)；溶質(zhì)運(yùn)移；優(yōu)先流；數(shù)值模擬

壓實(shí)的黏土在完整、無損、無裂隙等情況下滲透性差，防滲效果好，可以有效地防止污染物滲入到地下水中[1]。但是由于施工或管理的原因使黏土中產(chǎn)生了裂隙[2]，溶質(zhì)會(huì)沿著裂隙，以一種優(yōu)先流[3]的形式，快速到達(dá)相對于黏土層底部，會(huì)造成污染物的集中滲漏現(xiàn)象。在這種情況下，黏土防滲層被擊穿的時(shí)間縮短，防滲效果減弱甚至缺失。前人為了在室內(nèi)實(shí)驗(yàn)測試黏性土中污染物遷移的有效彌散系數(shù)和吸附參數(shù)，Barone和Rowe等[4]研究了揮發(fā)性有機(jī)化合物在自然黏土中的遷移過程。Edil等[5]獲得了壓實(shí)黏土中有機(jī)污染物運(yùn)移的彌散系數(shù)、吸附分布系數(shù)和降解系數(shù)。在室內(nèi)實(shí)驗(yàn)中污染物常采用示蹤劑來代替，模擬示蹤劑在裂隙和黏土基質(zhì)間的運(yùn)移轉(zhuǎn)化規(guī)律。黏土中污染物的運(yùn)移問題主要采用數(shù)值模擬的方法進(jìn)行研究。在研究土中污染物的運(yùn)移規(guī)律時(shí)，Rowe和Booker[6-7]采用有限層法進(jìn)行研究，Leo和Booker[8]采用邊界元法進(jìn)行研究，余開彪和胡斌等[9]采用有限單元法進(jìn)行靜態(tài)和動(dòng)態(tài)研究。MattaX和Kvte[6]把有蠟涂層的玻璃板用氫氟酸腐蝕成所需的微模型裂隙；后來Davis和Jonas[7]用抗蝕劑取代了蠟涂層。本文采用室內(nèi)實(shí)驗(yàn)的方法來研究水流和溶質(zhì)在裂隙黏土中的運(yùn)移規(guī)律，探討優(yōu)先流現(xiàn)象與黏土中裂隙的發(fā)育關(guān)系，并運(yùn)用數(shù)值法模擬實(shí)驗(yàn)，預(yù)測并分析了提高水頭和改變巖土體兩種情況時(shí)的溶質(zhì)運(yùn)移規(guī)律。

1 實(shí)驗(yàn)裝置與方案

1.1 實(shí)驗(yàn)裝置

實(shí)驗(yàn)裝置由四部分組成，包括圓柱形模型裝置、定/變壓力系統(tǒng)、聯(lián)通導(dǎo)管和電導(dǎo)率測量裝置（探頭已固定于模型內(nèi)）組成。實(shí)驗(yàn)?zāi)M一維條件下滲濾液滲漏穿過黏土層的情況，模型由高700 mm、直徑180 mm的圓柱體組成，裝置側(cè)面每隔50 mm設(shè)置一個(gè)直徑為5 mm的監(jiān)測孔。從圓柱形模型裝置底部開始，依次裝入厚70 mm的礫石，170 mm的細(xì)砂，300 mm的黏土土樣，余下的160 mm作為臨時(shí)貯存示蹤劑溶液的區(qū)域（圖1）。

圖1 實(shí)驗(yàn)裝置設(shè)計(jì)圖Fig.1 Design of experimental device

實(shí)驗(yàn)時(shí)，由蠕動(dòng)泵從原液池中不斷抽取示蹤劑注入頂部柱狀容器中，容器中示蹤劑保持固定的高度，多余的示蹤劑經(jīng)過導(dǎo)流管回流到原液池中。為了研究黏土中裂隙對水流和溶質(zhì)運(yùn)移的影響，在圓柱形模型的黏砂土中設(shè)置5 mm寬的裂隙，裂隙和黏砂土層一樣高，均為300 mm，為貫穿型裂隙。

1.2 實(shí)驗(yàn)步驟

（1）標(biāo)準(zhǔn)溶液的配置。在實(shí)驗(yàn)室配制了濃度分別為0.1、0.5、1、5、10、20、30、40、50、60 g/L的NaCl溶液，并測定了這些溶液的電導(dǎo)，通過一元線性回歸，獲得溶液濃度和電導(dǎo)的關(guān)系式：

式中，C為NaCl溶液濃度（mg/L）；σ25為25℃時(shí)的電導(dǎo)率值（ms/cm）。

（2）土樣裝填。先把黏土土樣放入盆中，加水充分飽和、沉淀。在圓柱形模型底部依次裝入礫石和細(xì)砂，兩層之間鋪一層濾網(wǎng)，防止細(xì)砂進(jìn)入礫石層。在砂層頂部鋪一層濾網(wǎng)，然后把制作好的裂隙放在裝置中間，與監(jiān)測孔軸線垂直。在裂隙四周充填飽和含砂黏土，當(dāng)裂隙與含砂黏土高度一樣時(shí)停止。土樣裝填完畢后，讓其靜止1～3天左右。

（3）加壓沉降。模型頂部用蓋子密封，往里加水，持續(xù)1～2周，待黏土沉降完全為止，然后把模型裝置里的水放干，再往里加入黏土，使其和裂隙一樣高，重復(fù)此步驟，直到黏土沉降不明顯為止。

（4）加入示蹤劑溶液進(jìn)行實(shí)驗(yàn)。

（5）測量監(jiān)測孔中示蹤劑溶液的電導(dǎo)率。

1.3 稀釋法測量示蹤劑溶液的電導(dǎo)率

實(shí)驗(yàn)過程中，由于監(jiān)測孔滲出液的體積比較小，一般在0.5～1 ml左右，因此須先收集監(jiān)測孔中的滲出液，用蒸餾水將其稀釋，測量稀釋液的電導(dǎo)率，根據(jù)式（2）換算出滲出液的電導(dǎo)率，其經(jīng)驗(yàn)計(jì)算公式為：

式中，σ1為監(jiān)測孔滲出液的電導(dǎo)率；v1為監(jiān)測孔滲出液的體積；σ2為蒸餾水的電導(dǎo)率；v2為蒸餾水的體積；σ為用蒸餾水稀釋后的稀釋液電導(dǎo)率。由于實(shí)驗(yàn)中用100 mL的容量瓶對滲出液進(jìn)行稀釋，因此監(jiān)測孔滲出液體積與蒸餾水體積之和為100 mL，即v1+v2=100 mL。

實(shí)驗(yàn)過程發(fā)現(xiàn)，將原溶液稀釋到100 mL，測其電導(dǎo)率，在利用式（2）計(jì)算出的原溶液電導(dǎo)率與原溶液的實(shí)際電導(dǎo)率有一定的偏差，一般偏大。本次實(shí)驗(yàn)中采用配制標(biāo)準(zhǔn)溶液測定溶液稀釋前后的電導(dǎo)率來進(jìn)行偏差校正。由此得出原溶液電導(dǎo)率與利用其稀釋液電導(dǎo)率換算出的原溶液電導(dǎo)率的關(guān)系式為：

式中，σ1為25℃時(shí)原溶液的電導(dǎo)率；σ為利用稀釋液的電導(dǎo)率換算出的原溶液電導(dǎo)率。

2 裂隙黏土中溶質(zhì)運(yùn)移數(shù)值模擬分析

2.1 模型離散

按照室內(nèi)實(shí)驗(yàn)圓柱形模型裝置的實(shí)際尺寸，采用Feflow軟件建立數(shù)值分析模型[10]，對模型采用三角形網(wǎng)格單元，將三維圓柱形模型剖分為59 360個(gè)節(jié)點(diǎn)，105 678個(gè)單元，網(wǎng)格剖分圖見圖2。

表1 模型參數(shù)反演值Tab.1 Inversion values of model parameters

圖2 模型網(wǎng)格剖分圖Fig.2 Meshing of the model

圖3 各監(jiān)測孔處滲出液濃度實(shí)驗(yàn)值與計(jì)算值對比Fig.3 Comparison of experimental and calculated values of exudate concentration at each monitoring hole

2.2 彌散參數(shù)計(jì)算

根據(jù)實(shí)測和計(jì)算的示蹤劑濃度，對模型進(jìn)行多次調(diào)試計(jì)算得到黏土、裂隙和礫石層的相關(guān)參數(shù)（表1）。為了驗(yàn)證模型的計(jì)算參數(shù)反演值的正確性，分別選取1#、4#和7#監(jiān)測孔處示蹤劑濃度實(shí)驗(yàn)值與計(jì)算濃度值進(jìn)行比較（圖3）。

從圖3中可看出，模型計(jì)算值先增大，最后均趨于穩(wěn)定，與實(shí)驗(yàn)值吻合較好，表明可以用驗(yàn)證的模型來預(yù)測示蹤劑在裂隙黏土中的遷移規(guī)律。

2.3 預(yù)測結(jié)果分析

（1）水頭變化對污染物遷移的影響

改變黏土層上部的水頭，分別為0.3、0.5、1、2 m四種工況，1#、3#和7#監(jiān)測孔處的溶質(zhì)濃度模擬計(jì)算結(jié)果見圖4。從圖中可以看出，隨著黏土層上部水頭的增大，各監(jiān)測孔處的溶質(zhì)濃度逐漸增加，但達(dá)到穩(wěn)定的時(shí)間差別較大。1#和3#監(jiān)測孔處的溶質(zhì)濃度隨著黏土層上部水頭的增大而增大，7#監(jiān)測孔處的溶質(zhì)濃度變化趨勢基本一致，相差不大，在短時(shí)間內(nèi)就趨近于原示蹤劑溶液濃度達(dá)到穩(wěn)定。說明在貫穿型裂隙黏土中，黏土層上部水頭的變化對黏土中的溶質(zhì)濃度分布有影響，對溶液順著裂隙穿過黏土層造成其下地下水的影響不大。

（2）土體滲透性對污染物遷移的影響

滲透系數(shù)分別取1×10-4、1×10-5、1×10-7、1×10-10m/s四種工況，1#、3#和7#監(jiān)測孔處的溶質(zhì)濃度模擬計(jì)算結(jié)果見圖5。從圖中可看出，隨著巖土體滲透系數(shù)的增大，相同時(shí)刻各監(jiān)測孔處溶質(zhì)的濃度也增大，但不同深度處的溶質(zhì)濃度增加速度不一樣。隨著土層滲透系數(shù)的增大，土層中的溶質(zhì)濃度趨近于原示蹤劑溶液濃度達(dá)到穩(wěn)定的時(shí)間縮短。黏土滲透系數(shù)為1×10-4m/s時(shí)1#和3#監(jiān)測孔處的溶質(zhì)濃度變化最大，說明隨著黏土滲透系數(shù)的增大，相同時(shí)刻監(jiān)測孔處的溶質(zhì)濃度會(huì)增大；黏土滲透系數(shù)為1×10-7m/s和1×10-10m/s時(shí)1#和3#監(jiān)測孔處的溶質(zhì)濃度變化趨勢基本一致，說明滲透系數(shù)較小時(shí)對黏土中溶質(zhì)的運(yùn)移規(guī)律影響較小，裂隙導(dǎo)致的優(yōu)先流對溶質(zhì)的運(yùn)移起著主導(dǎo)作用，黏土中的基質(zhì)流不明顯；7#監(jiān)測孔處的溶質(zhì)濃度變化趨勢基本一致，說明黏土滲透系數(shù)變化對溶質(zhì)運(yùn)移規(guī)律的影響不大。

圖4 黏土層上部水頭變化時(shí)的模擬結(jié)果對比曲線圖Fig.4 Curve of the simulation results for different upper heads of clay layer

綜上所述，在貫穿型裂隙土層中，滲透系數(shù)的變化對黏土中的溶質(zhì)濃度分布有影響，但對溶液順著裂隙穿過黏土層造成其地下水的影響不大。

3 結(jié)論

1）黏土中有裂隙存在時(shí)，溶質(zhì)會(huì)順著裂隙快速穿過黏土層，形成優(yōu)先流現(xiàn)象。

圖5 黏土滲透系數(shù)變化時(shí)的模擬結(jié)果對比曲線圖Fig.5 Curve of the simulation results for different clay permeability coefficients

2）預(yù)測結(jié)果顯示黏土層上部水頭變化對黏土中的溶質(zhì)濃度有影響，但對裂隙和黏土層下部地下水體中的溶質(zhì)濃度影響不大。由于貫穿裂隙的出現(xiàn)，污染物對不同巖土體內(nèi)部的滲入速度有差異，但對不同巖土層下方的含水層污染速度基本一致，此時(shí)如通過改變防滲材料來減少污染物入滲沒有作用。

[1]田素川. 裂 隙巖體注漿漿液擴(kuò)散規(guī)律研究[D].北京：中國礦業(yè)大學(xué)，2014.

[2]何 俊，萬 娟，王 宇.壓實(shí)黏土干燥裂隙及滲透性能研究[J].工程地質(zhì)學(xué)報(bào)，2012(3)：397-402.

[3]徐宗恒，徐則民，曹軍尉，等.土壤優(yōu)先流研究現(xiàn)狀與發(fā)展趨勢[J].土壤，2012（6）：905-916.

[4]BARONE F S，ROWE R K，QUIGLEY R M.A laboratory estimation of diffusion and adsorption coefficients for several volatile organics in a natural clay soil[J].Journal ofcontaminant hydrology，1992，10(3)：225-250.

[5]Edil T B, Wambold W S, Park J K. Partitioning of VOCs in Clay Liner Materials[C]// Geoenvironment 2000@ sCharacterization, Containment, Remediation, and Performance in Environmental Geotechnics. ASCE, 2014：775-790.

[6]ROWE R K, BOOKER J R. The analysis of pollutant migration in a non-homogenous soil[J]. Géotechnique, 1984, 34(4)：601-612.

[7]BOOKER J R, ROWE R K. A finite layer technique for calculating three-dimensional pollutant migration in soil[J]. Géotechnique, 1986, 36(2)：205-214.

[8]LEO C J, BOOKER J R. A boundary element method for analysis of contaminant transport in porous media I：homogeneous porous media[J]. International Journal for Numerical & Analytical Methods in Geomechanics, 1999, 23(14)：1701-1715.

[9]余開彪, 胡斌, 張挺. 垃圾滲濾液在土體中擴(kuò)散規(guī)律的模擬計(jì)算[J]. 巖石力學(xué)與工程學(xué)報(bào), 2004, 23(11)：1944-1948.

[10]賀國平，邵景力，崔亞莉，等.FEFLOW在地下水流模擬方面的應(yīng)用[J].成都理工大學(xué)學(xué)報(bào)，2003（4）：356-361.

（責(zé)任編輯 王利君）

Experimental and numerical simulation of movement of contaminants through cracks

ZENG Jie1,GU Zhengcong2,HUANG Yong1

(1.School of Earth Science and Engineering, Hohai University, Jiangsu Nanjing, 210098, China;2. Yunnan Hubai Environmental Protection Technology Co., Ltd. Yunnan Kunming, 650034)

In order to study the solute transport in fissured clay, experimental apparatus for solute transport in fissured clay are developed based on previous saturated seepage experiments. The solute transport is monitored by measuring the conductivity of the tracer solution in the experimental setup law. The results show that the concentration of exudate in the bottom of the model tends to be stable in a short time, which indicates that the existence of cracks in the clay will cause the solute to pass through the clay layer rapidly and form preferential flow phenomena. The upper head of the clay layer has little effect on the solute concentration in the clay, and the effect of the water head is similar in the different soils.

fractured clay；saturated seepage flow；solute transport；preferential flow；numerical simulation

P64

A

1673-9469（2017）01-0072-04

10.3969/j.issn.1673-9469.2017.01.016

2016-12-15

國家自然科學(xué)基金資助項(xiàng)目（41572209）

曾頡（1992-），男，江西贛州人，碩士，研究方向?yàn)樗牡刭|(zhì)。