1 引言(Introduction)
地下水中非水相液體(Non-aqueous Phase Liquids，NAPLs)的污染和治理問題是國內外研究熱點(韓占濤等，2013;陳余道等，2004;劉漢樂等，2014;李卉等，2015;趙勇勝等，2015).重非水相液體(Dense Non-Aqueous Phase Liquids，DNAPLs)密度比水大，常見的DNAPL包括三氯乙烯(Trichloroethylene，TCE)、四氯乙烯(Perchloroethylene，PCE)等氯代烴，在世界各地地下水中被頻繁檢出(National Research Council, 1994; Harwell et al., 1999; Vogan et al., 1999; 馬長文等，2007).由于低水溶性、弱遷移性、難降解性并能穿透含水層而滯留在含水層底部，形成長期的污染源，其遷移和修復比一般的污染物甚至是LNAPLs(Light Non-Aqueous Phase Liquids)更加復雜.
地下水曝氣技術(Air Sparging, AS)被認為是去除飽和土壤和地下水中可揮發性有機化合物(Volatile Organic Compounds, VOCs)的最有效方法(Rabideau and Blayden, 2003; Lesson et al., 2005).它是一項將空氣注入到含水層飽和區中通過空氣流的吹脫作用去除VOCs并增強微生物降解效果的創新性原位修復技術(Kim and Annable, 2006)，具有低成本、高效率和原位操作的突出優勢.國內外眾多專家和學者的實驗和應用研究表明AS技術對于地下水NAPLs污染物的去除效果非常明顯，已成為地下水NAPLs修復技術的首選.
生物修復技術治理地下有機污染不僅費用低而且效果好，沒有或很少有二次污染.大量原位生物修復的研究工作證明，在微生物降解污染物的過程中會產生難溶于水的氣體，比如厭氧菌中的產甲烷菌在修復氯代烴所產生的甲烷氣體(Strack et al., 2005; Ye et al., 2009)，反硝化菌在生物脫氮時產生的氮氣(Soares, 1989, 1991; Soares, 2000)等.采用地下水曝氣技術或生物修復方法修復土壤和地下水中DNAPL污染時，會使含水層介質中出現大量人工注入的或生物作用產生的氣體.這些氣體的存在變化了原來的飽水地下水系統，形成了更加復雜的多相流地下水系統，氣體的存在必然對DNAPL在新的復雜多相流系統中的運移產生影響.
相關學者通過建立二維砂箱來定性或定量研究NAPL在非飽和區域的運移行為(Pantazidou and Sitar, 1993; Van Geel and Sykes, 1994; Illangasekare et al., 1995; Oostrom et al., 2003).其中光透法作為一種無損，非侵入的監測方法被廣泛用于二維砂箱的室內監測(Tidwell and Glass, 1994; Niemet and Selker, 2001; O'Carroll and Sleep, 2007; Mumford et al., 2009; Ye et al., 2009).CCD相機的應用更是極大地提高了監測時間的即時性和空間的高密度性，其中時間分辨率可達到1 s，空間分辨率小于1 mm2.基于以上研究基礎和科學問題，采用光透法電荷耦合裝置(Charge Coupled Device, CCD)監測技術，通過室內二維砂箱試驗，選擇TCE為目標污染物，以飽水條件下DNAPL的運移作為對照，研究人工注氣和生物產氣條件下DNAPL在孔隙介質中的運移，探討氣體存在對于DNAPL運移的影響，這一研究對于DNAPL運移及修復研究尤其是修復策略制定有重要的參考意義.然而，國內外卻未見相關研究報道.
2 光透法原理(The theory of light transmission method)
光透法原理主要基于比爾-朗伯定律，又稱比爾定律、布格-朗伯-比爾定律，是光吸收的基本定律，適用于所有的吸光物質，包括氣體、固體、液體、分子、原子和離子.當光照射于某一吸收介質的表面，在通過一定厚度的介質后，由于介質吸收了一部分光能，透射光的強度減弱，吸收介質的濃度愈大，介質的厚度愈大，則光強度的減弱愈顯著.比爾-朗伯定律可表達為：
式中，I是指穿過介質i后的光強，I0是入射光源的光強，μi是介質i的吸收系數，li是介質i的厚度.其中C是一個光學幾何參數，與發光點和觀察點的位置有關，對于準直光源或是光源(燈箱)和介質(砂箱)到接收器(CCD相機)的距離大致相同時，C可以作為一個常量忽略不計(Niemet and Selker, 2001).在本實驗中，CCD相機距砂箱的距離(1.8 m)遠大于燈箱與砂箱間的距離(12.5 cm), 符合上述假定條件，故C在此處忽略不計.
當光穿過一個多相系統時，結合菲涅耳定律(Fresnel's law)，比爾-朗伯定律亦可表達為：
式中，τj, k是指穿過相j及相k間界面的透射率，其它符號的意義與式(1)相同，此處不再贅述.根據菲涅耳定律，
式中，nj，nk分別是指兩個相鄰相j，k的折射率.本實驗中水、空氣、染色前TCE及石英砂的折射率分別為1.33299、1.00027、1.4782(程能林，2002)和1.547(袁繼祖，2007).染色后TCE的吸收系數不可忽略，透光性變弱，本試驗中其透光性介于水和氣體之間.根據比爾-朗伯定律、菲涅爾定律，結合本試驗具體條件可知，對于具體某一像素點，介質完全飽水時的光強最大，完全飽氣時的光強值最小，而充滿染色后TCE或處于中間狀態時的光強值則介于前二者之間.
3 材料與方法(Materials and methods)3.1 實驗裝置
在室內建立了透射光系統來監測DNAPL的運移形態，其中透射光監測系統主要包括3個部分：二維砂箱，燈箱和CCD相機.實驗系統的示意圖見圖 1.實驗裝置的具體細節和參數可以參考Ye等(2009)，楊靖等(2011)，簡單介紹如下：二維砂箱尺寸為：55 cm寬×45 cm高×1.28 cm厚，由兩塊厚度為10 mm的鋼化玻璃內夾一個厚度為12.8 mm的中心鋁框組成，并利用兩個鋁質邊框在外部進行固定，鋁質邊框的大小略大于玻璃.玻璃與中心框之間利用橡膠條和玻璃膠(GE Silicone Ⅱ)進行密封.砂箱具體示意圖見圖 2，孔14~16作為進水口，出水口(孔1~3)在實驗過程中固定在略高于砂箱頂部的高度，孔17和26連接水壓管測量砂箱內水壓，其它孔作為取樣口或注樣口使用.燈箱作為砂箱的唯一光源，位于砂箱的一側12.5 cm處，由六根平行的日光燈管(Panasonic，YZ18RR6500K)發光經過擴散板后保證光源的均勻性.CCD相機(AP2E，Apogee Instruments，Auburn，CA)連接到一臺計算機并位于砂箱另一側1.8 m處，置于與砂箱一體的木質暗箱內，用來接收透過砂箱的光照，并通過軟件Maxim DL(Ottawa，ON)記錄其光強值.孔隙介質為半透明石英砂(C190 Accusand，20/30 mesh，Unimin-Le Sueur，MN，USA)，平均粒徑為0.73 mm，顆粒均勻性指數為1.21，孔隙度約為0.338.采用分層填充的方式將石英砂濕裝填入砂箱內，約2 cm為一層，每層間充分攪拌均勻后并壓實.三氯乙烯(TCE)作為此次研究DNAPL代表污染物，實驗中利用油紅O(Oil Red O，購于Sigma公司)對TCE染色，染色濃度為100 mg·L-1.實驗中統一使用相同染色濃度的TCE，因此染色后的TCE在下文中亦簡稱為TCE.實驗過程中利用注射泵(LSP01-2A，保定蘭格恒流泵有限公司)或蠕動泵(L/S便攜式，Cole Parmer公司)將目標物質(空氣或TCE)注入砂箱內.
圖 1
圖 1透射光系統示意圖
圖 2
圖 2砂箱示意圖 (數字'1'等代表孔號)
3.2 二維砂箱試驗3.2.1 飽水條件下DNAPL注入實驗
首先，在飽水的砂箱(C1)中進行TCE注入試驗.在完全飽水砂箱中，通過孔8垂直插入的不銹鋼注樣針注入TCE，注入點位置位于砂箱頂部下部約7 cm處.通過蠕動泵向砂箱注入TCE，整個TCE注入過程持續34 min，0~5 min的注入速度為15 mL·min-1，6~34 min的注入速度為5 mL·min-1，TCE注入過程中只有出水口孔1處于打開狀態.TCE注入完成后，關閉注入口和出水口，砂箱靜置7 d.注入TCE前，利用CCD相機記錄飽和砂箱(C1)飽水時光強值(I01-w).TCE注入過程中利用CCD相機實時監測砂箱的光強值(I1)，在TCE注入過程中CCD相機的拍照間隔設為1 min，在砂箱靜置過程中則設為1 h.
3.2.2 人工注氣條件下DNAPL注入實驗
其次，在砂箱(C2)中分別先后進行了人工注氣和TCE注入試驗.在完全飽水的砂箱中，通過孔9垂直插入不銹鋼注樣針，在砂箱頂部往下約27 cm處注入空氣.為了不考慮氣體的溶解，保證在氣體的注入過程中氣體不會溶解在水中，注氣前利用經過充分曝氣(空氣，曝氣時間達24 h以上)的水充分驅替(約5個砂箱體積)砂箱中的水.本實驗中，采用了手動、注射泵和蠕動泵3種不同的注入方式，具體過程見表 1.氣體注入結束后，稱量排出水的體積為223.54 mL.注氣過程完成后，砂箱靜置24 h，箱內氣體分布未發現有明顯變化，氣體分布達到穩定.在人工注氣分布穩定的砂箱中，與飽水條件下類似，通過在孔8垂直插入不銹鋼注樣針，注入點位于砂箱頂部往下約7 cm處，利用蠕動泵以5 mL·min-1的速度往砂箱(C2)注入TCE，共計注入194.37 mL TCE，TCE注入完成后砂箱靜置6 d.在整個注入及靜置實驗過程中，CCD相機實時監測砂箱內的光強變化.在人工注入氣體(空氣，購于南京五十五研究所)前，CCD相機記錄砂箱(C2)飽水時的光強值(I02-w).在注入TCE前，CCD相機記錄砂箱(C2)人工注氣完成時的光強值(I02-g).TCE注入過程中利用CCD相機實時監測砂箱的光強值(I2)，在TCE注入和砂箱靜置過程中CCD相機的拍照間隔都設為3 min.
3.2.3 生物產氣條件下DNAPL注入實驗
最后，在砂箱(C3)內分別先后進行了生物產氣和TCE注入試驗.本試驗中采用的微生物是KB-1(SiREM公司免費提供)，是一種從三氯乙烯(TCE)甲醇污染場的土壤和地下水中提取的復合菌群，可用于降解所有的氯代乙烯，如四氯乙烯(PCE)，三氯乙烯(TCE)等.通過底部孔21、22、23接種KB-1菌液250 mL，通過砂箱底部各孔注入150 mL的500 mg·L-1甲醇溶液/3 d.微生物生長約3個月后，砂箱內有明顯的氣體分布，累計共有263 mL的水排出砂箱.此時通過孔9不銹鋼注樣針，注入點位于砂箱頂部以下約5 cm處，以1.3 mL·min-1的速度注入TCE，持續122 min，注入結束后稱量排出水的體積為173.38 mL，TCE注入結束后砂箱靜置7 d.砂箱(C3)接種KB-1前，利用CCD相機記錄砂箱(C3)飽水時的光強值(I03-w).在注入TCE前，CCD相機記錄砂箱(C3)生物產氣完成時的光強值(I03-g).TCE注入過程中利用CCD相機實時監測砂箱的光強值(I3)，在TCE注入過程中CCD拍攝間隔設為1 min，在砂箱靜置的過程中則設為1 h.
4 結果與討論(Results and discussions)4.1 系統誤差分析
由于透射光系統的穩定性和外界因素的干擾造成光強測量值的誤差，所以選取參照區域(area of reference，AOR)進行系統誤差分析.AOR選取標準為：實驗過程中該區域砂箱內介質狀態保持恒定，即光強信息理論上應保持不變，AOR位置見圖 3.此外，對AOR不同時刻光強變化進行統計分析：選擇飽水條件下的兩個不同時刻，簡稱為飽水狀態;選擇TCE注入前和注入后兩個時刻，簡稱為注入狀態，用以研究TCE注入時AOR光強情況.AOR光強變化統計結果見表 2和圖 4.據圖 4知：對于飽水狀態(飽水時兩個不同時刻)的光強變化量統計，不同試驗(飽水條件，人工注氣和生物產氣)系統誤差分布特點為中間集中，兩端較少，整體近似為正態分布;表 2中的統計參數卻不相同且均值都不為零，甚至與零值偏差較大.說明了透射光監測系統的系統誤差在不同時刻均服從正態分布(Bob et al., 2008, Tidwell and Glass, 1994, Niemet and Selker, 2001)，但參數具有一定的隨機性.對于試驗1-飽水條件下，飽水狀態和注入狀態理論上AOR介質都是完全飽水的，但是統計參數卻相差較大(均值分別為-2357和-244)，說明透射光監測系統的系統誤差與時間有關.相關研究亦表明：系統誤差與實驗過程中光源和CCD監測系統的開、關有關，同時燈光的輸出強度與室溫也相關(Niemet and Selker, 2001).
圖 3
圖 3不同條件下TCE注入完成后CCD圖片 (a.飽水條件，b.人工注氣和c.生物產氣, 其中虛線框指AOR：區域大小為10 cm×30 cm，與上邊界和左邊界對齊;白線框指AOI：區域大小均為10 cm×10 cm，均距下邊界2 cm，其中圖a和b中距右邊界25 cm，圖c中距右邊界20 cm)
圖 4
圖 4飽水狀態AOR(按表 2順序)光強變化量統計直方圖
為了進一步研究系統誤差并驗證透射光系統監測技術的有效性，選取一個特定區域，選取標準為：TCE注入完成時該區域發生了明顯的光強變化，稱之為影響區域(area of impact，AOI)，位置見圖 3，光強變化量統計參數見表 2.據表 2知，飽水狀態同一試驗中AOR和AOI的統計參數也表現出了差異(如試驗1-飽水條件下：AOR的均值為：-1059.3;AOI的均值為：-1547.3)，這說明系統誤差與空間位置相關.對于注入狀態，3組試驗中AOI光強變化量均值均遠遠大于AOR(以試驗2-人工注氣為例：AOI的均值為5192.3;AOR的均值為：117.8)，同時AOI的95%置信下限亦遠大于AOR區域的95%置信上限(以試驗2-人工注氣為例：AOI的95%置信下限為5183.7;AOR的置信上限為：121.3)，說明透射光監測系統雖然具有一定的誤差，但仍可有效地監測多相流系統中流體的運移.
為區分光強變化量是由系統誤差造成還是由TCE注入引起，此處引入一個閾值(Threshold Value，ThV)：系統誤差引起光強變化的95%置信上限.在TCE注入過程中，當光強變化量高于ThV時則認為是由TCE注入引起的，是有效的光強變化量.反之，當光強變化量低于或等于ThV時則認為是由于系統誤差造成的，是無效的光強變化量.
4.2 氣體分布及飽和度計算
3個砂箱中TCE注入前的光強分布如圖 5所示.其中圖 5a是砂箱(C1)飽水時的光強分布情況，顯示受填砂方式等人為因素的影響介質表現出一定的非均質性和成層性，實驗中難以保證絕對均質.圖 5b顯示了TCE注入砂箱(C2)前單點人工注氣產生的氣體分布情形.在人工注氣的過程中，氣體以垂直向上運動為主，最終在砂箱頂部形成約8~10 cm的人工氣體連續分布，圖 5b中實線勾勒了人工注氣形成氣體分布的外部輪廓.單點的注氣方式，使得連續氣體分布下方形成一個類似'倒三角'的遷移路徑，可見氣體遷移的'指狀'通道.圖 5c是菌群KB-1接種到砂箱(C3)約3個月以后，形成的最終氣體分布情況.氣體隨機不連續地分布于整個砂箱(C3)中，未能形成大范圍的連續氣體分布，在砂箱頂部形成了局部小范圍的連續氣體分布區域."指狀"通道現象與人工注氣(圖 5b)相比更加明顯，圖中紅點代表了被水包裹的封閉性氣泡.從圖 5c中亦可以看出八根長短不一的取樣針位置，這些取樣針用于實驗過程中采集水樣進行測試分析.
圖 5
圖 5光強分布圖
Niemet和Selker(2001)基于光透法原理建立了關于如何利用光強值求解流體飽和度的5個水/氣模型，根據本實驗條件，選擇水/氣模型(模型C具體表達式見文獻(Niemet and Selker, 2001)中的方程(11)，此處不詳細介紹.章艷紅等(2014)為了簡化模型，引入了新的參數C1(C1=Ires/Iw，其中Ires是指當砂箱只有殘余水時穿過的光強值，Iw是指砂箱完全飽水時穿過的光強值)，利用公式Sg=ln(I/Iw)/ln(C1)計算氣體飽和度，根據實驗具體得到參數取值為0.10，計算得到氣體飽和度分布如圖 6所示.
圖 6
圖 6砂箱氣體飽和度分布圖
4.3 DNAPL入滲和再分布
將TCE的運移分為入滲和再分布兩個階段，其中TCE注入階段為入滲期，TCE停止注入后砂箱進入再分布期.3組試驗中由TCE入滲和再分布引起光強值的變化如圖 7所示.相關研究表明：如果DNAPL的量足夠多，可以克服毛管壓力和驅逐孔隙中的水，DNAPL在重力作用下將繼續往下遷移直到到達弱透水層(Fetter，2011).在飽水條件下，TCE主要是在自身重力的影響下以垂向滲流為主，垂向滲流的同時存在一定的水平滲流(圖 7 d中污染羽對稱軸-紅色實線是一條曲線)，局部存在明顯的水平滲流(圖 7d中x方向約15 cm，depth方向約5 cm處).這可能與前面提到的介質局部非均質性有關，砂箱內出現局部相對較松散或密實的情形.5~15 min內，TCE污染羽前緣以蠕動泵5 mL·min-1的注入速度下10 min內在飽水砂箱(C1)中向下遷移的距離是20.47 cm，平均垂向下遷速度是2.05 cm·min-1.
圖 7
圖 7不同飽水條件下TCE入滲和重分布過程中光強變化量分布
與飽水條件下TCE遷移相比，人工注氣試驗中NAPL/水/氣三相系統中TCE的遷移仍是受自身重力影響以垂直向下運動為主，然而氣體的存在使得TCE的遷移更為復雜，污染羽整體形狀更加不規則，尤其是水/氣界面和氣體遷移通道處.實驗結果表明，在depth約5 cm，x方向10~20 cm處存在明顯的水平運動.在注入點附近，污染羽形狀呈"寬口"型，相對于飽水條件橫向擴散更加活躍，這可能與TCE的揮發有關，大范圍的連續分布氣體使得TCE能夠以氣態形式的存在.從圖 7b可以看出，TCE在飽水區域(左側)向下遷移的距離明顯大于非飽水區域(右側)，在本試驗中DNAPL在飽和區向下遷移速度更快，說明下遷過程優先驅替孔隙中的水.但是圖 7b中非飽和區域(右側)的橫向遷移比飽水區域(左側)更加明顯，說明TCE在進行橫向遷移時時更傾向于驅替孔隙中的空氣.DNAPL在橫向遷移時，當有足夠大的可利用的壓力使自由相的NAPL通過毛孔時，NAPL便會取代曾經占據的空氣或水.所需壓力的大小取決于作用在毛孔兩邊不同流體的毛細作用力的大小.毛細力作用在這兩個流體單元的方式在某種程度上可以解釋為可濕潤性的作用：即能進入毛孔的歸類為可濕潤流體.在空氣/NAPL/水的系統中，可濕潤性流體的濕潤性為：水 > NAPL > 空氣.因此，NAPL在進行橫向遷移時優先驅替孔隙中的空氣，占據空氣的孔隙空間.人工注氣條件下，在5 mL·min-1的注入速度下TCE污染羽前緣9 min內在砂箱(C2)中向下遷移的距離是15.57 cm，得到平均垂向下遷速度是1.72 cm·min-1.然而，在飽水條件TCE平均垂向下遷速度是2.05 cm·min-1，這說明氣體的存在對TCE的垂向向下遷移具有阻礙作用，與圖 7b中TCE在氣體內更易進行橫向遷移現象吻合.在TCE重分布的過程中連續氣體分布的界面處存在較明顯的光強變化(圖 7h和7k)，造成這一變化可能有以下原因：①TCE進入重分布時期由于連續氣體分布的存在使得TCE易以揮發態的形式存在，從而引起水/氣界面的下移;②TCE入滲過程中氣體被TCE驅替出來后由于出水口的設計不易排出砂箱而且TCE的注入阻礙了氣體的向上遷移，靜置過程中氣體通過原有的遷移通道開始向上遷移并在頂部累積造成水/氣界面下移;③試驗中通過注射泵或蠕動泵的外部壓力將TCE注入砂箱內而且注入點位于氣體分布內，TCE入滲過程中由于增加的外力影響使得氣體壓力變大，從而引起氣體體積壓縮;在重分布過程中這個外力的影響消失，氣體壓力變小，從而體積變大，造成水/氣界面下移.
在生物產氣的條件下，污染羽形狀整體呈"狹長"型，TCE的入滲過程還是受重力影響下的垂向下遷為主，遷移過程中橫向擴散卻弱于人工注氣和飽水的情況下，可能與大范圍的不連續氣體分布和TCE注入速度有關(圖 7c，7f).因為本實驗中的生物產氣在砂箱頂部未能形成大范圍的連續氣體分布，所以在注入點附近未發現人工注氣中出現的"寬口"型的污染羽形態，TCE在注入點附近的橫向遷移弱于人工注氣條件.對比圖 6b和7f，TCE注入結束后，污染池底部(Depth方向約0~5 cm處)光強值變化比較明顯的區域位于x方向約10~30 cm和40~50 cm處，然而生物產氣氣體飽和度分布圖(圖 6b)中這些區域的氣體飽和度卻相對較低.相反地，圖 7f中污染池底部(Depth方向約0~5 cm處)光強值變化較弱的區域(x方向約0~10 cm和30~40 cm處)生物產氣氣體飽和度卻相對較高.根據光透法原理可知，水/氣兩相系統中氣體飽和度和由氣體引起的光強變化量(ΔI1)存在正相關的關系(如圖 5c和6b).在污染池底部區域，由生物產氣引起的光強變化量(ΔI1)越大，氣體飽和度越大，但是由TCE注入引起的光強變化量(ΔI2)卻越小.因此說明TCE在水平滲流時，由TCE注入引起的光強變化量(ΔI2)和由生物產氣引起的光強變化量(ΔI1)存在負相關的關系.這是因為非飽和帶中TCE殘余飽和度更低從而引起的光強變化量較小，而且出水口的設計更利于砂箱系統排水導致氣體可能會滯留在TCE污染羽內部.在TCE重分布過程中污染羽形態整體沒有明顯變化，但上部TCE受重力影響繼續向下遷移，上部光強值有所降低，未發現人工注氣條件下的水/氣界面處明顯光強變化的現象，可見砂箱內若干離散光強變化點.這可能與重分布期內TCE的局部遷移引起氣體的重新分布有關.
為進一步研究入滲和重分布期間光強隨時間的變化，選擇了6個研究單元，每個研究單元大小為10×10個柵格(空間大小約為：0.5 cm×0.5 cm)，研究單元中心點具體坐標位置見表 3.其中：Up代表污染羽中心上部研究單元，Mi代表污染羽中心中部研究單元，Bo代表污染羽中心下部研究單元，In則是為了研究水/氣界面處光強變化而設定的研究單元，為了研究TCE底部污染羽而設定了分別代表污染羽左下，右下研究單元Bl和Br.同時，為了便于對比分析不同空間位置的光強變化，引入光強比例(Light ratio)這一統計參數，定義為目標物質引起的光強變化量(ΔI2)與對應飽水條件下光強值(I1)的比值.3組試驗中6個研究單元光強比例隨時間變化的關系見圖 8.水/氣界面研究單元(In)光強比例歷時曲線在飽水和生物產氣條件下近似為一條光強比例為0的直線，只有在人工注氣條件下該歷時曲線在TCE注入過程(入滲時期)中出現了較小程度的抬升，在靜置(TCE重分布)時期光強比例發生了較大的抬升，這一現象與圖 7h，7k現象一致，說明入滲過程中主要是以TCE的垂向下遷為主，下遷過程中優先驅替孔隙中的水，同時非飽和區域也發生少量橫向遷移從而驅替氣體;然而人工注氣條件下重分布期氣體的重分布現象十分明顯不可忽略.污染羽左下，右下的研究單元(Bl和Br)光強比例歷時曲線結果表明：在飽水條件和人工注氣條件下，TCE在底部堆積成污染池時左右橫向遷移沒明顯的選擇性(突破時間：tBl≈tBr)，兩者砂箱底部污染羽形態整體接近，左右對稱性良好.然而，在生物產氣條件下，TCE在底部堆積成污染池時優先向右橫向遷移然后再向左邊橫向遷移(突破時間：tBr= 41 min，tBl=92 min).在生物產氣條件下，圖 6b底部右側(x方向30~40 cm)處較之于左側(20~30 cm)聚集較明顯的氣體分布，TCE污染池優先向右(氣體分布明顯區域)遷移，這也說明了TCE橫向遷移時優先驅替氣體，這一現象與人工注氣條件下的圖 7b一致，而且與地下水流流向無關(注入過程中圖 2中孔1作為排水口，即地下水流向大致為底部向左).具體聯系污水寶或參見http://www.bnynw.com更多相關技術文檔。
圖 8
圖 8不同飽水條件下TCE入滲和再分布過程中光強比例變化歷時曲線
TCE注入完成后進入重新分布期，NAPL補給不足時，作用在自由DNAPL上的壓力就會消失，形成不定形的或者滲透到原來連續流體的NAPL，在毛細力作用下，NAPL會被截留在單獨的毛孔或成組的毛孔中(菲利普·B·貝迪恩特等，2010).在試驗①中，Up，Mi和Bo3個研究單元在TCE注入結束后，光強比例在重分布期先下降后穩定，表明這些研究單元重分布前期(1d)存在NAPL自由態的下遷，后期則是被截留在單獨的毛孔或成組的毛孔中不易移動.Bl和Br光強比例在重分布前期上升是由于接受來自上部的自由態NAPL的補給，后期下降則是由實驗誤差造成的，試驗中觀測到底部TCE污染池有部分進入下部取樣孔內.在試驗②中人工注氣條件下，Up，Mi和Bo(除Up外)也表現出和飽水條件下相似的變化，但是Bo下降幅度較之于Mi大，是因為Mi受到了氣體的影響(Mi偏大的原因)，Up主要是受氣體的影響導致TCE可能以揮發態形式存在.Bl和Br在重分布期光強比例有了小幅的抬升，然后保持穩定.在試驗3)中，Up，Mi和Bo光強比例在重分布期表現出下降趨勢，Br重分布期基本保持穩定，Bl卻表現出先升后降的變化，這是與研究單元對應的氣體飽和度有關，Bl研究單元的氣體飽和度較高(0.36)利于揮發，而Br則處于氣體低飽和度區域(0.18).
圖 9給出了3組試驗中光強比例在不同時刻的飽和度分布，在注入體積45 mL時，飽水條件下為單峰，而在人工注氣和生物產氣條件則表現為雙峰，存在一個光強比例接近零值的高峰，可能是揮發態的TCE引起的較小光強變化.在注入完成和靜置時，零值峰同時出現在3組試驗中，而且在人工注氣和生物產氣條件下的相對比例高于飽水條件.這是因為：除上面所述的揮發態TCE引起的較小光強變化外，隨著TCE污染羽和水接觸的表面積增大，此時以溶解態形式存在于砂箱中的TCE不容忽略，這一部分引起的光強變化表現出來.與注入體積(45 mL)相比，3組試驗中注入完成時光強比例峰值更大，這是由于TCE的不斷注入TCE并在砂箱內不斷累積引起的光強變化.除近零值峰外，光強比例直方圖隨時間表現均為由單峰到雙峰的變化規律.相對于飽水條件而言，人工注氣和生物產氣條件下兩峰的間距較大.雙峰中的光強比例較小峰可能對應遷移過程中截留TCE以殘余態形式存在時引起的光強變化，然而光強比例較大峰則代表底部污染池NAPL態TCE引起的光強變化.由于氣體的存在導致運移路徑上孔隙中TCE的截留量變小，造成下遷過程中光強比例變小;然而更多自由態TCE在底部聚積形成高飽和度的污染池，造成該區域光強比例變大，從而使兩峰的間距增大.飽和帶中的NAPL的剩余飽和度(LNAPL或DNAPL)比非飽和度帶的要高，有以下幾個原因：①NAPL與空氣密度比大于NAPL與水的密度比，有利于NAPL從非飽和區排出，②在飽和區，NAPL是不可濕潤流體，并被束縛在較大的孔隙中，③在有利的毛細管條件下，NAPL在包氣帶中有盡可能向外圍展開的趨勢.
圖 9
圖 9不同條件下光強比例不同時刻直方圖
5 結論(Conclusions)
1) 基于CCD相機的透射光系統是研究多相流遷移規律的一種有效的非侵入式監測方法，應用于實時監測包含水/NAPL, 水/氣，水/NAPL/氣多相系統中流體滲流過程，系統誤差服從不同的正態分布，參數值具有一定的隨機性.
2) 人工注氣在砂箱內注入點附近形成一個類似'倒三角'的'指狀'通道，最終在砂箱頂部形成了一個連續氣體分布區;而生物產氣在整個砂箱內形成多個隨機的指狀通道，同時很多氣體以隨機的不連續的氣泡形式分布在砂箱內，最終在砂箱頂部形成連續氣體分布區域.
3) DNAPL的入滲過程受重力影響以垂向向下遷移為主，垂向滲流時更易于驅替孔隙中的水份，然而橫向滲流時優先驅替孔隙中的空氣，與地下水流向無關.
4) 氣體的存在阻礙了TCE的垂向下遷：本次研究中飽和條件下，5~15 min內，TCE污染羽前緣以蠕動泵5 mL·min-1的注入速度下10 min內在飽水砂箱(C1)中向下遷移的距離是20.47 cm，平均垂向下遷速度是：2.05 cm·min-1;人工注氣條件下，在5 mL·min-1的注入速度下TCE污染羽前緣9 min內在砂箱(C2)中向下遷移的距離是15.57 cm，得到平均垂向下遷速度是1.72 cm·min-1.
5) 氣體的存在使得DNAPL污染羽的形狀更加不規則;氣體的存在也導致遷移路徑上孔隙中TCE的截留量變小.(來源：環境科學學報 作者：章艷紅)
