土壤氣抽提多井方案的數值模擬研究

欒 旭,武曉峰,胡黎明 (清華大學水利水電工程系,水沙科學與水利水電工程國家重點實驗室,北京 100084)

土壤氣抽提多井方案的數值模擬研究

欒 旭,武曉峰*,胡黎明 (清華大學水利水電工程系,水沙科學與水利水電工程國家重點實驗室,北京 100084)

以存在石油烴污染的北京市某加油站場地為研究對象,設計了土壤氣抽提(SVE)原位修復系統,應用AIR3D軟件對SVE作用下土壤包氣帶的氣壓場進行數值模擬.劃定以土壤氣壓強低于或等于101.122kPa的區域為抽氣影響區域.模擬結果顯示,單井的抽氣影響區域不能滿足要求;多井抽氣時,井間區域土壤氣體真空度受多井共同作用,擴大了抽氣的影響區域.抽氣井不同布置情況下井間距大小對相互作用效果存在影響.研究發現,井間相互作用隨著井間距的增大而減小,井間距存在臨界值.井間距臨界值與抽氣井的布置形式相關.對于本研究中的2眼井、3眼井、4眼井情景,此臨界值分別在3～5m、5～8m、8～10m區間.超過此臨界值時井間相互作用將顯著減弱,并導致抽氣影響區域不能覆蓋整個井間范圍.

土壤氣抽提；數值模擬；AIR3D；多井方案

土壤氣抽提(SVE)是治理包氣帶內土壤石油烴污染的一種有效的原位物理修復方法[1-2].它通過抽氣井保持一定的真空度抽取非飽和土壤中的空氣,土壤中的有機污染物揮發轉化為氣態形式隨抽氣氣流離開土體,將抽出的氣體處理后再予排放,從而達到清除非飽和區的揮發性污染物,修復土壤包氣帶的目的.

在實施SVE前,首先需要針對場地設計現場抽氣方案.井數及井布置方案是影響工程的成本與治理效果的2個關鍵因素.數值模擬,尤其是三維數值模擬是研究這2項影響因素的重要手段.

20世紀80年代起,發達國家在土壤修復技術方面進行了大量研究.總體上,對于 SVE的作用機理已經有一定的認識[3],數學模型由二維向三維發展[4-6],并注重對地表覆蓋情況以及抽氣影響半徑等實際工程因素的討論[7-11].近幾年,隨著土壤多相流研究以及生物降解研究的不斷深入,SVE與曝氣、與生物降解相結合的治理手段正越來越受到重視[12-13].

我國對SVE的研究起步較晚,除在國外研究基礎上進行了數值模型研究外[14-15],近年來也開展了一些室內實驗[16]及現場試驗研究[17].目前國內外針對SVE的研究鮮見針對多井情況的專門討論,而單井往往不能滿足污染現場的治理要求,與工程實際應用仍有較遠距離.

本研究以北京市某加油站的石油烴污染場地作為研究對象,針對SVE系統在井群抽氣情況下的井布置方案及其抽氣效果進行數值模擬,以討論井群間各井的相互影響以及井布置方案、井間距等因素對抽氣效果的影響,并為現場治理提供指導.

1 材料與方法

1.1 SVE基本原理

揮發性有機污染物在土壤包氣帶中主要有4種存在形式:氣相、溶解相、吸附相和非水相流體(NAPL)相.其中氣相指污染物揮發成氣態的存在形式;溶解相指污染物溶于水的部分;吸附相指污染物吸附在土壤顆粒上的部分;NAPL相是指未溶于水的液態有機污染物.當NAPL進入土壤包氣帶后,將有部分揮發成氣態并同時有部分溶解于土壤水中.一般來說,在溶解相和氣相中的污染物均可與土壤顆粒發生吸附/解吸關系.對于SVE中的氣相—土壤吸附,通常情況下,由于土壤固相表面均有水膜覆蓋,氣、固相幾乎沒有接觸面存在,因而可忽略該過程.

SVE通過布置在污染土壤中的抽氣井抽出土壤孔隙中的空氣,形成負壓,使土壤氣中的揮發性有機污染物(VOCs)同時被抽出.氣相污染物濃度降低之后,溶解相、吸附相以及NAPL相的污染物將進一步向氣相轉化,從而以土壤氣為載體,通過不斷抽取氣相污染物的方式實現對土壤中污染物的去除.

土壤氣抽提的流場控制方程為:

式中: μ為氣體黏性系數; θ為土壤孔隙度; P為氣體壓強; t為時間; k~為土壤氣相滲透率張量.

1.2 AIR3D數值模型

運用AIR3D模型進行SVE的模擬.AIR3D為三維有限差分法計算軟件,可以模擬多孔介質中氣體滲透過程和氣壓場分布,由美國地質調查局(USGS)開發.穩態時其控制方程為:

式中:kxx, kyy, kzz分別為沿x,y,z方向的土壤氣相滲透率;φ =P2, P為壓強.

模型初始條件各差分網格氣壓均為大氣壓;邊界條件:側邊界處氣壓為大氣壓,底邊界(地下水水位處)為不透氣邊界,上邊界(地表)可與大氣連通或封閉,抽氣井花管處網格保持抽氣氣壓不變.

1.3 污染場地狀況

本研究對象為北京市某加油站場地,該加油站采用地下儲油罐儲油,在地表布有加油機.儲油罐下修筑了混凝土防滲槽以防止汽油泄漏污染儲油罐下部的土壤.儲油罐與加油機間通過輸油管相連.清華大學王博[18]對現場的土壤污染狀況進行了檢測.檢測階段在加油站范圍內9個位置布置取樣點,抽取加油站混凝土層下 1m深處氣樣,用氣相色譜聯用儀對樣本中總揮發性有機物(TVOC)濃度及苯、二甲苯、正己烷等11種污染物濃度進行測定,TVOC濃度結果如表1所列.并依據這9個采樣點的數據繪制了土層1m深處水平面上土壤氣體中TVOC濃度的等值線圖,如圖1所示.

表1 各檢測點位TVOC濃度[18]Table 1 TVOC concentration at each detection point[18]

圖1 土壤氣體中TVOC濃度等值線[18]Fig.1 TVOC concentration contour in soil gas[18]

2 數值模擬

2.1 多井聯合抽氣的疊加效應

根據研究區TVOC濃度的分布,確定將高濃度區10m×10m的正方形區域作為目標治理區域,如圖1所示.用AIR3D建立30m×30m×3m的三維模型,模擬 SVE穩態時的抽氣效果,場地地質參數及抽氣井設計參數如表2所列.

表2 模型相關參數Table 2 Parameters in the AIR3D model

為判斷抽氣效果的優劣,在模擬得到的土壤氣壓等值線圖中,劃定壓強低于或等于101.122kPa(0.998倍標準大氣壓)的區域為 SVE的有效影響區域.

首先模擬單井抽氣的情況,井布置在目標治理區域中心,模擬結果如圖 2所示.在單井情況下,SVE有效影響區域為圓形,影響半徑為2m左右,未能覆蓋整個目標治理區域.為了能夠覆蓋整個治理區域,需要設計多井抽氣的方案.

圖2 單井抽氣時SVE有效覆蓋范圍Fig.2 SVE effective area when single well is in use

對于多井抽氣的情況,井群中各井對土壤氣壓分布存在相互影響.圖3顯示了間隔5m布置2個相同抽氣井同時作用,2井花管位置連線上的土壤氣體真空度分布,井的設計參數如表2所示.真空度在抽氣井附近達到最大值,隨著距抽氣井距離的增加真空度降低.在2井中間的區域,真空度比單井作用時高,表明多井的抽氣效果疊加,使該范圍內的土壤氣壓進一步降低.

圖3 兩井同時抽氣情況下土壤中氣體真空度分布Fig.3 Soil vacuum distribution when two wells are extracting at the same time

多井聯合作用下抽氣效果的疊加同時體現在SVE影響范圍的變化上,多井聯合作用下影響區域的大小會大于多個抽氣井單獨作用時影響區域之和.以本研究北京市某加油站場地的模擬結果為例,當在同一SVE系統中布置相互間隔5m的 4個井聯合抽氣時,影響區域為一個近似的圓角四邊形,如圖4所示;而將4個單井SVE系統的影響區域疊加,則得到4個互不相連的圓形,如圖5所示.可以看出 4個抽氣井聯合作用提高了井間土壤區域的真空值,增大了SVE的影響區域.

圖4 四井聯合抽氣時的SVE有效覆蓋范圍Fig.4 SVE effective area for 4-well system

2.2 多井抽氣井間距對相互作用效果的影響

對于含有多個抽氣井的井群,在設計時應考慮各抽氣井的間距以及井的排列方式對壓強場的影響,這2項因素都會影響SVE的設計方案.本研究針對直線式、三角形式、矩形式3種不同的井群排布方式分別進行了3組模擬.3組模擬中,抽氣井間的區域分別受2眼井、3眼井、4眼井的共同影響.在每組模擬中再變化抽氣井的間距以討論井間距與抽氣影響區域的關系,模擬結果如圖6～圖8所示.

3組模擬結果表明,隨著井間距增大,多井的相互作用減弱,且抽氣井間距存在一臨界值,當井間距小于此值時,抽氣井井間區域土壤氣體真空度達到治理目標,此區域中的氣體污染物可被有效去除.當井間距大于此值時,井間區域不能完全達到治理所要求的真空度,有效治理范圍不能覆蓋全部井間區域.

圖5 四個單獨作用的抽氣井有效范圍的疊加Fig.5 Superimposition of the SVE effective area of 4 single functional well

對比3組模擬結果,對于不同的井布置方案,這一臨界值的大小不同:雙井抽氣時這一值在3～5m間;3井抽氣時,這一值在5～8m間;4井抽氣時,這一值在8～10m間.

3 結論

3.1 數值模擬的結果表明,SVE過程中多井聯合抽氣時井間區域土壤氣的壓強受到多個抽氣井的共同影響.多井聯合作用下影響區域的大小會大于多個抽氣井單獨作用時影響區域之和.

3.2 劃定101.122kPa土壤氣壓強等值線為抽氣的影響范圍邊界,改變抽氣井的井間距,抽氣井的井間相互影響隨著抽氣井井間距增大而減小.

3.3 對于2眼井、3眼井、4眼井的3種不同井布置,井間距分別存在一個臨界值,使多井的抽氣效果與影響范圍最佳.且這一值與井的不同布置有關,隨著井數增多而增大.在該案例中,2眼井時,這一值在3～5m間;3眼井時,這一值在5～8m間;4眼井時,這一值在8～10m間.

圖6 雙井聯合抽氣時抽氣井有效范圍Fig.6 SVE effective area for 2-well system

圖7 三井聯合抽氣時抽氣井的有效范圍Fig.7 SVE effective area for 3-well system

圖8 四井聯合抽氣時抽氣井的有效范圍Fig.8 SVE effective area for 4-well system

[1] 楊樂巍,黃國強,李鑫鋼.土壤氣相抽提(SVE)技術研究進展 [J].環境保護科學, 2006,32(6):62-65.

[2] 黃國強,李鑫鋼,徐世民.土壤氣相抽提作用機制探討和基本數學模型建立 [J]. 土壤學報, 2004,41(3):394-399.

[3] Pedersen T A, Curtis J T. Soil Vapor Extraction Technology[M]. New Jersey: Noyes Data Corpration, 1991.

[4] Massmann J W. Applying groundwater flow models in vapor extraction system design [J]. Journal of Environmental Engineering, 1989,115(1):129-149.

[5] Baehr A L, Hult M F. Evaluation of unsaturated zone air permeability through pneumatic tests [J]. Water Resour. Res., 1991, 27(10):2605-2617.

[6] Baehr A L, Joss C J. An updated model of induced airflow in the unsaturated zone [J]. Water Resour. Res., 1995,31(2):417-421.

[7] Kaleris V, Croise J. Estimation of cleanup time for continuous and pulsed soil vapor extraction [J]. Journal of Hydrology, 1997, 194(1-4):330-356.

[8] Wilson D J, Rodríguez-Maroto J M, Goamez-Lahoz C. Soil cleanup by in-situ aeration. XIX. Effects of spill age on soil vapor extraction remediation rates [J]. Separation Science and Technology, 1994, 29(13):1645-1671.

[9] Wilson D J, Gómez-Lahoz C, Rodríguez-Maroto J M. Soil Cleanup by in-situ Aeration. XVI. Solution and diffusion in mass-transport-limited operation and calculation of darcy's Constants [J]. Separation Science and Technology, 1994,29(9): 1133-1163.

[10] Shan C, Falta R W, Javandel I. Analytical solutions for steady state gas flow to a soil vapor extraction well [J]. Water Resour. Res., 1992,28(4):1105-1120.

[11] Mohr D H, Merz P H. Application of a 2D air flow model to soil vapor extraction and bioventing case studies [J]. Ground Water, 1995,33(3):433-444.

[12] Rahbeh M E, Mohtar R H. Application of multiphase transport models to field remediation by air sparging and soil vapor extraction [J]. Journal of Hazardous Materials, 2007,143(1/2): 156-170.

[13] Rathfelder K M, Lang J R, Abriola L M. A numerical model (MISER) for the simulation of coupled physical, chemical and biological processes in soil vapor extraction and bioventing systems [J]. Journal of Contaminant Hydrology, 2000,43(3/4): 239-270.

[14] 薛 強,梁 冰,王惠蕓.揮發性有機污染物去污過程的數值模型及應用 [J]. 化工學報, 2005,56(10):1962-1966.

[15] 黃國強,李鑫鋼,姜 斌,等.豎井SVE修復有機污染土壤的數值模擬分析 [J]. 化工學報, 2003,54(8):1134-1140.

[16] 李金惠,聶永豐,馬海斌,等.油污染土壤氣體抽排去污模型及影響因素 [J]. 環境科學, 2002,23(1):92-96

[17] 楊樂巍.土壤氣抽提(SVE)現場試驗研究 [D]. 天津: 天津大學, 2006.

[18] 王 博.加油站石油污染修復技術研究 [D]. 北京: 清華大學, 2010.

Simulation of multiple wells system in soil vapor extraction.

LUAN Xu, WU Xiao-feng*, HU Li-ming (State Key Laboratory of Hydroscience and Engineering, Department of Hydraulic Engineering, Tsinghua University, Beijing 100084, China). China Environmental Science, 2012,32(3)：535～540

A numerical simulation of the air pressure distribution in the contaminated vadose zone was conducted using the AIR3D model, based on the designed in-situ SVE remediation system for the contaminated gas-station site. An affected area was defined as the domain which has an air pressure lower than 101.122kPa. At first, a single well was used in the simulation, but the result turned out that the affected area was not sufficient to cover the contaminated zone. Then a multiple wells system was proposed, and the combined and enhanced effect was observed. Further, the influence of well interval and well distribution on the affected area was discussed in the multiple wells system. The numerical results showed that interaction between wells decreases along with the increase of the well interval. There existed a critical value representing the maximum well interval of the combined effect. The critical well interval varied with the change of the well distribution pattern. The critical well intervals in the 2-well, 3-well, 4-well systems discussed in this study are in intervals of 3-5m, 5-8m, 8-10m respectively. In well systems with larger well intervals than the critical values, interaction between wells would decrease significantly, and the affected area would not cover the entire inter-well area.

soil vapor extraction；numerical simulation；AIR3D；multiple wells system

X53

A

1000-6923(2012)03-0535-06

2011-06-12

北京市科技計劃項目(D07050601510000)

* 責任作者, 副教授, wuxiaofeng@tsinghua.edu.cn

欒 旭(1987-),男,北京人,清華大學碩士研究生,研究方向為地下水水環境.