張 燁，施小清，鄧亞平，江思珉，吳劍鋒，吳吉春

(1.表生地球化學(xué)教育部重點(diǎn)實(shí)驗(yàn)室/南京大學(xué)地球科學(xué)與工程學(xué)院，江蘇南京 210023;2.同濟(jì)大學(xué)土木工程學(xué)院，上海 200092)

隨著石油行業(yè)的不斷發(fā)展，石化產(chǎn)品泄漏事件屢屢發(fā)生，其泄漏的燃料油及有機(jī)氯溶劑經(jīng)包氣帶土層進(jìn)入地下水，在含水層中以液態(tài)形式存在。難溶于水且密度大于水的，稱為重非水相液體(Dense Nonaqueous Phase Liquids，DNAPL)［1］。

目前針對 DNAPL 污染物的修復(fù)方法很多［2～3］，其中空氣曝氣(Air Sparging，AS)方法對于去除土壤和地下水中的揮發(fā)性有機(jī)物(Volatile Organic Compounds，VOCs)具有高效、低成本、修復(fù)周期短、原位操作等優(yōu)點(diǎn)。AS將空氣注入含水層，藉由揮發(fā)與生物降解等作用機(jī)制，通過空氣使VOCs揮發(fā)成氣相，隨空氣上升至非飽和帶，再配合土壤氣相抽提(Soil Vapor Extraction，SVE)將 VOCs抽出處理［4］。AS對于中、高亨利常數(shù)(低分子量、高揮發(fā)性)的有機(jī)物修復(fù)效果顯著，但對有機(jī)氯化物去除效率較差［5］。近年來，為解決有機(jī)氯化物的污染問題，部分學(xué)者對AS的應(yīng)用方式進(jìn)行了改進(jìn)，將注入冷空氣改為熱空氣或蒸汽［5～6］。因?yàn)楦邷氐目諝?或蒸汽)可使得土壤和地下水溫度升高，從而加快VOCs的揮發(fā)速率，且土壤顆粒在高溫下對有機(jī)污染物的吸附能力明顯降低，更有利于污染物的去除［6～7］。

已有的AS研究主要是基于現(xiàn)場修復(fù)經(jīng)驗(yàn)［8～9］以及室內(nèi)砂箱實(shí)驗(yàn)結(jié)合數(shù)值方法［7，10～11］，甚少有對NAPL有機(jī)污染修復(fù)的多相流模擬-優(yōu)化問題的研究，同時考慮各種影響因素的多目標(biāo)模擬優(yōu)化問題就更少。事實(shí)上，考慮多個影響因素的相互作用，根據(jù)多目標(biāo)優(yōu)化模擬結(jié)果綜合修復(fù)效率和實(shí)際需求設(shè)定現(xiàn)場操作方案，這樣能大大減少實(shí)驗(yàn)成本［12］。

近年來發(fā)展了一系列求解地下水模擬優(yōu)化問題的新方法，其中進(jìn)化算法(Evolution Algorithm，EA)［13］尤其引人注目，包括遺傳算法(Genetic algorithm，GA)［14～15］、粒 子 群 優(yōu) 化 方 法 (Particle Swarm Optimization，PSO)［16］、差 分 進(jìn) 化 算 法 (Differential Evolution，DE)［17］和 MH 方 法 (Metropolis-Hastings Algorithm)［18］等。相對于單相流的多目標(biāo)優(yōu)化，多相流多目標(biāo)優(yōu)化問題計(jì)算時間長，考慮變量多，具有非常大的研究價(jià)值［19］，EA在針對該類問題的研究分析較少。

由于被廣泛使用和不合理處置，三氯乙烯(Trichloroethylene，TCE)成為土壤和地下水中常見的有機(jī)污染物之一。本文選其作為典型有機(jī)氯化物污染物，利用T2VOC模擬注入蒸汽(和空氣)對DNAPL污染物的修復(fù)過程，比較同時注入蒸汽和空氣與只注入蒸汽對修復(fù)結(jié)果的影響。通過對蒸汽井(和空氣井)參數(shù)的敏感性分析，基于模擬-優(yōu)化方法建立修復(fù)過程的多目標(biāo)優(yōu)化模型，評估蒸汽焓和空氣井注入量參數(shù)對DNAPL污染物修復(fù)效率的影響。

1 研究方法

1.1 T2VOC

T2VOC是多相流數(shù)值模擬軟件TOUGH2的眾多模塊之一，可模擬三相(氣相、水相及NAPL相)，4個組分(溫度/熱量、水、空氣及VOCs)在三維非均質(zhì)孔隙介質(zhì)中的運(yùn)移［20～21］，目前已廣泛應(yīng)用于NAPL污染物遷移和修復(fù)問題的研究［6，11，22］，具體原理和數(shù)值求解方法可參見 Falta等［23］。

1.2 多目標(biāo)優(yōu)化算法

建立了多目標(biāo)優(yōu)化模型后，采用一種基于多種進(jìn)化算法的自適應(yīng)優(yōu)化方法AMALGAM(A Multialgorithm Genetically Adaptive Multi-objective Method)［24］進(jìn)行求解。該算法搭建了混合優(yōu)化框架，能夠同時使用多種進(jìn)化算法(常見的有GA、PSO、MH、DE)，通過遺傳適應(yīng)或自適應(yīng)的方式生成后代以發(fā)揮各種種群進(jìn)化算法的優(yōu)勢，算法求解過程見圖1。

圖1 AMALGAM優(yōu)化算法設(shè)計(jì)框圖(虛線部分為AMALGAM中調(diào)用的獨(dú)立算法)Fig.1 Flowchart describing the evolutionary process for AMALGAM(the dotted square symbols represent the individual algorithms used by AMALGAM)

該方法首先通過LHS(Latin Hypercube Sampling)方法生成一個隨機(jī)的初始種群Pt(種群大小為N)，然后對初始種群中的每個個體調(diào)用T2VOC并根據(jù)得到的觀測值計(jì)算每個個體的適應(yīng)度值。通過AMALGAM，各獨(dú)立算法各自生成子代種群合并為種群大小為N，每一個獨(dú)立算法貢獻(xiàn)的子代個數(shù)為由AMALGAM根據(jù)自適應(yīng)過程決定，其中t指示第幾代種群，q指示第幾個獨(dú)立算法。AMALGAM同時調(diào)用了多個獨(dú)立算法從相同的父代種群獨(dú)立生成子代種群，由于各個算法使用的進(jìn)化方法不一樣，各個算法得到的子代種群也不完全一樣，由此能夠保證種群的多樣性，使算法能更好地搜尋最優(yōu)解。合并Qt與Pt得到新的種群Rt(Rt=Pt∪Qt，種群大小為2N)，對 Rt采用快速非劣排序算法(Fast Nondominated Sorting，F(xiàn)NS)進(jìn)行非劣排序，搜選出最好的N個個體作為下一代父代種群Pt+1。通過這種精英搜索策略，父代種群Pt在多途徑搜索和自適應(yīng)過程中不斷進(jìn)化，確保最優(yōu)值始終存在于種群Rt中，并被賦給下一代父代種群Pt+1。

為了確保最好的算法能夠在下一次優(yōu)化中提供更多的子代個體，AMALGAM根據(jù)各獨(dú)立算法在每一次子代生成過程中的貢獻(xiàn)率Ψi更新算法所占權(quán)重以及N值:

AMALGAM能同時考慮多種不同的進(jìn)化搜索策略，自適應(yīng)地根據(jù)搜索過程中的貢獻(xiàn)率調(diào)整各個算法的權(quán)重，并考慮隨機(jī)性和確定性以降低優(yōu)化結(jié)果陷于局部最優(yōu)的風(fēng)險(xiǎn)［26］。與單獨(dú)算法相比，對于復(fù)雜的高維度多峰優(yōu)化問題，AMALGAM有明顯的優(yōu)勢。

2 蒸汽注入模擬參數(shù)敏感性分析

通過對一個理想范例進(jìn)行敏感性分析，為后文模擬-優(yōu)化模型的參數(shù)選取提供基礎(chǔ)。為了研究蒸汽(和空氣)在含水層中的運(yùn)移狀況，本試驗(yàn)在無污染物的情況下，向含水層持續(xù)注入蒸汽(和空氣)進(jìn)行敏感性數(shù)值分析。

2.1 蒸汽(空氣)注入模型概況

考慮二維研究區(qū)域:長、高分別為40 m、20 m，將其離散成20層40列共計(jì)800個網(wǎng)格單元。研究區(qū)頂部為大氣邊界，底部為隔水邊界，潛水面埋深為5 m，初始溫度為25℃。巖層顆粒密度為2.65 g/cm3，孔隙度為0.42，X、Y、Z方向的滲透率均為1×10-11m2。在左側(cè)邊界X=0 m，垂向深度Z=14～15 m處均勻布置2口蒸汽井，垂向深度Z=12～13 m處均勻布置2口空氣井。

2.2 蒸汽(空氣)注入?yún)?shù)敏感性的情景分析

對蒸汽井單井注入速率為1.2×10-1g/s，蒸汽溫度為100℃(2.677×106J/kg，純蒸汽)時進(jìn)行情景分析，圖2(a)為持續(xù)注入4 h時的溫度鋒面，靠近蒸汽井處，由于蒸汽液化放熱，井附近溫度高于100℃，蒸汽鋒面以井為圓心，成圓弧狀向外擴(kuò)散;圖2(b)為不同時刻60℃的溫度鋒面，從圖中可以看出蒸汽擴(kuò)散速率較穩(wěn)定，蒸汽隨時間呈近似勻速擴(kuò)散。

僅注入蒸汽時，對蒸汽焓和蒸汽井單井注入速率分別進(jìn)行敏感性分析，具體參數(shù)見表1。選取TCE作為DNAPL污染物，沸點(diǎn)為87.1℃，因此進(jìn)行敏感性分析時溫度鋒面選取85℃的溫度鋒面，蒸汽(和空氣)持續(xù)注入4 h。圖3(a)為蒸汽注入速率恒定，對蒸汽溫度(焓)進(jìn)行敏感性分析，溫度鋒面變化情況顯示，隨著蒸汽焓值的增大，蒸汽擴(kuò)散速率也隨之增大;當(dāng)蒸汽為非純蒸汽時，蒸汽擴(kuò)散速率變化明顯，當(dāng)蒸汽為純蒸汽時，蒸汽擴(kuò)散速率逐漸趨緩。圖3(b)為蒸汽溫度(焓)恒定，對蒸汽井單井注入速率進(jìn)行敏感性分析，溫度鋒面變化情況顯示，隨著蒸汽注入速率的增大，蒸汽擴(kuò)散速率也隨之明顯增大。

圖2 蒸汽井單井注入速率為1.2×10-1g/s、蒸汽溫度為100℃時的溫度鋒面Fig.2 Temperature fronts as the steam injection rate is 1.2 ×10-1g/s，steam temperature is 100℃

在同時注入蒸汽和空氣時，對空氣井單井注入速率進(jìn)行敏感性分析。蒸汽井和空氣井的參數(shù)見表2。圖4(a)為有無空氣注入情況下持續(xù)注入4 h時溫度鋒面的變化情況，其中空氣井單井注入速率為2.4×10-2kg/s。可以看出，有空氣注入時溫度鋒面形狀發(fā)生明顯變化，蒸汽被空氣壓迫貼近隔水底板運(yùn)移，且在水平方向上的運(yùn)移速率較只注入蒸汽時加快;圖4(b)為不同空氣井單井注入速率下溫度鋒面的變化情況，可以看出，空氣注入速率越大，空氣對蒸汽向下壓迫的作用越明顯，蒸汽在水平方向上的運(yùn)移速率越快。

表1 敏感性分析參數(shù)Table 1 Sensitive analysis of steam wells

圖3 對蒸汽焓和蒸汽井單井注入速率的敏感性分析Fig.3 Sensitivity analysis of the specific steam enthalpies and steam injection rates

圖4 對空氣井單井注入速率進(jìn)行敏感性分析Fig.4 Sensitivity analysis of the air injection rates

表2 敏感性分析實(shí)驗(yàn)蒸汽井和空氣井的參數(shù)設(shè)置Table 2 Parameters of steam wells and air wells for sensitive analysis

由以上幾組對比試驗(yàn)表明，蒸汽焓、蒸汽井和空氣井的注入速率以及蒸汽持續(xù)注入時間對于溫度鋒面的變化都有較大影響;空氣的注入能明顯改變蒸汽鋒面形式，抑制蒸汽向上運(yùn)移的趨勢，迫使蒸汽在含水層底部加速運(yùn)動，從而有效去除含水層底部的DNAPL。

3 對比算例

3.1 模型概況

建立一個理想二維非均質(zhì)模型，研究對比只注入蒸汽與同時注入蒸汽和空氣時AS對多孔介質(zhì)中DNAPL污染物的修復(fù)效果。研究區(qū)概念模型見圖5［6］，研究區(qū)長、高分別為 1.0 m、0.75 m，包括粗砂和細(xì)砂兩種介質(zhì)，實(shí)線框以內(nèi)為細(xì)砂介質(zhì)，以外為粗砂介質(zhì)。初始條件下，污染區(qū)域內(nèi)污染物處于殘余飽和度狀態(tài)且污染物呈均勻分布，即污染物殘留于死端孔隙且不可流動，污染區(qū)域以外無污染物分布，參數(shù)見表3。相對滲透率函數(shù)和毛管壓力函數(shù)分別采用Stone模型［27］和 Parker模型［28］;相對滲透率即多相流系統(tǒng)中兩相滲透率之比。模型被均勻離散為1200個網(wǎng)格，網(wǎng)格的大小為0.025 m×0.025 m;在研究區(qū)的左邊界垂向深度Z=-0.3～-0.2 m之間均勻布置4口空氣井，垂向深度Z=-0.4～-0.3 m之間均勻布置4口蒸汽井，在水平方向X=0.875 m、垂向深度Z=0～0.5 m之間均勻布置20口抽氣井，用于模擬水流運(yùn)動及SVE的抽提過程;參數(shù)見表4。除注入井外，頂部為大氣邊界，底部為隔水邊界。通過污染物飽和度的變化觀察污染物運(yùn)移情況并判定AS修復(fù)效果。

圖5 概念模型示意圖Fig.5 Problem definition and sketch

表3 模型的水文地質(zhì)及其它參數(shù)［6］Table 3 Hydrogeological parameters and other parameters of the simulation model

表4 蒸汽井和空氣井的參數(shù)設(shè)置Table 4 Parameters of steam wells and air wells

3.2 結(jié)果分析

3.2.1 只注入蒸汽

僅設(shè)置蒸汽井，持續(xù)注入5 h。隨著注入時間增長，DNAPL污染物覆蓋范圍逐漸減小，污染物主要聚集在細(xì)砂區(qū)，原本處于不可流動狀態(tài)下的污染物在蒸汽的推動下移出死端孔隙，向右運(yùn)移。隨著蒸汽的注入，污染物飽和度先增大再減小;在蒸汽推動和重力的共同作用下，少量DNAPL污染物運(yùn)移至隔水底板(圖6)，最終滯留于隔水底板。總體來看，注入蒸汽能有效推動DNAPL污染物在水平方向上的運(yùn)移，但不能阻止其在重力作用下向隔水底板的運(yùn)移，而滯留于隔水底板的污染物將降低污染物的抽提去除效率。

3.2.2 同時注入蒸汽和空氣

同時布置蒸汽井和空氣井，持續(xù)注入5 h。與只注入蒸汽時類似，隨著注入時間增長，DNAPL污染物覆蓋范圍逐漸減小，污染物主要聚集在細(xì)砂區(qū)，污染物受蒸汽驅(qū)動移出死端孔隙并向右運(yùn)移，污染物的飽和度先增大再減小，最大飽和度小于只注入蒸汽時的情況;少量DNAPL污染物運(yùn)移接近隔水底板，但并未發(fā)生聚集，說明空氣對DNAPL污染物存在一定的向上托舉的作用，減緩污染物向下運(yùn)移的趨勢;最終細(xì)砂區(qū)以外的DNAPL污染物幾乎都被去除(圖7)。這是由于空氣逼迫蒸汽貼近隔水底板運(yùn)移，可有效去除細(xì)砂區(qū)下端的污染物?？傮w來看，同時注入蒸汽和空氣能有效提高DNAPL污染物在水平方向上的運(yùn)移速率，減緩污染物向下運(yùn)移的趨勢，阻止污染物在隔水底板再冷凝聚集，壓迫蒸汽貼近隔水底板運(yùn)移以去除下部的污染物，從而提高污染物的抽提去除效率。

4 多目標(biāo)優(yōu)化模擬

4.1 優(yōu)化模型

實(shí)際生產(chǎn)中，選擇修復(fù)方案往往需要綜合考慮多方面因素，以選取較優(yōu)的實(shí)施方案。通過模型的模擬-優(yōu)化可得到一系列的Pareto解，以供實(shí)際需求選擇最佳方案。基于上文同時注入蒸汽和空氣的修復(fù)模型，對蒸汽焓和空氣井單井注入速率2個參數(shù)進(jìn)行多目標(biāo)函數(shù)優(yōu)化，包括殘余污染物質(zhì)量和修復(fù)成本2個目標(biāo)函數(shù)。殘余污染物質(zhì)量函數(shù)為修復(fù)2 h后含水層中殘余的DNAPL污染物的總質(zhì)量，修復(fù)成本函數(shù)為蒸汽成本與空氣成本之和，具體可表示為:

圖6 注入蒸汽2h(a)和4h(b)后DNAPL污染物的遷移分布情況Fig.6 Distribution of DNAPL contaminant after 2 hours(a)，4 hours(b)after injection of pure steam

圖7 同時注入蒸汽和空氣2 h(a)和4 h(b)后DNAPL污染物隨時間遷移分布情況Fig.7 Distribution of DNAPL contaminant after 2 hours(a)，4 hours(b)after co-injection of steam and air

其中:2 200≤Ts≤2 875，10-4≤Va≤10-5，a=2 ×10-4，b=105。

式中:MT——?dú)堄郉NAPL污染物的總質(zhì)量/kg;

Ts——蒸汽焓/(J·kg-1);

Va——空氣井單井注入速率/(kg·s-1);

a——加權(quán)系數(shù)/(kg·J-1);

b——加權(quán)系數(shù)/(s·kg-1);

J1——目標(biāo)函數(shù)/kg;

J2——目標(biāo)函數(shù)。

4.2 優(yōu)化結(jié)果及分析討論

為檢測AMALGAM算法，將其與常用的NSGA算法［26，29］進(jìn) 行 比 較，其 中 NSGA 設(shè) 置 種 群 參 數(shù) 與AMALGAM一致，種群大小均為50，種群代數(shù)為20，其它參數(shù)為默認(rèn)值。同時只調(diào)用一個獨(dú)立算法(GA、PSO、MH、DE)來觀察該獨(dú)立算法進(jìn)行多目標(biāo)優(yōu)化的結(jié)果。

圖8為AMALGAM和NSGA算法的結(jié)果比較。由圖8可以看出，計(jì)算種群帶代數(shù)為20代時，AMALGAM與 NSGA的 Pareto曲線基本一致，AMALGAM相較NSGA更為平滑;對比AMALGAM和NSGA計(jì)算第3代種群得到的Pareto鋒面可以看出，AMALGAM已收斂得到一個較平滑的Pareto鋒面，與第20代Pareto鋒面幾乎重合，NSGA計(jì)算至第3代種群仍未收斂得到Pareto鋒面。在進(jìn)行多項(xiàng)流多目標(biāo)優(yōu)化模擬中，多項(xiàng)流模型較單向流模型更為復(fù)雜，一次正問題的計(jì)算時間長，且需考慮多個目標(biāo)函數(shù)的相互作用關(guān)系，需進(jìn)行多次迭代才能得到較平滑完整的Pareto鋒面。由圖8可以得出，AMALGAM算法相對NSGA，收斂速度更快，體現(xiàn)出綜合了多算法的優(yōu)越性。

圖9為AMALGAM中各獨(dú)立算法求解多目標(biāo)模型的優(yōu)化結(jié)果對比?？梢钥闯觯珼E、GA和MH得到的Pareto曲線表現(xiàn)較好，與圖8中AMALGAM的Pareto曲線較為一致，但平滑程度稍差;而 PSO所得到的Pareto曲線曲線不完整。由圖10中每一代子代種群中各個算法貢獻(xiàn)的個體數(shù)可以看出，各算法在每一代中貢獻(xiàn)的個體數(shù)波動起伏變化，特別是MH算法的貢獻(xiàn)率，這正體現(xiàn)出AMALGAM計(jì)算過程中對多種算法的自適應(yīng)性?？傮w上，DE、GA所貢獻(xiàn)的個體數(shù)較多，而PSO所貢獻(xiàn)的個體數(shù)呈遞減趨勢，說明AMALGAM能夠根據(jù)獨(dú)立算法的表現(xiàn)自適應(yīng)調(diào)整各算法的貢獻(xiàn)率，保證子代種群最優(yōu)。

圖8 AMALGAM(a)和NSGA算法(b)求解多目標(biāo)模型的優(yōu)化結(jié)果Fig.8 Results of multi-objective optimization using AMALGAM(a)and NSGA(b)respectively

圖9 AMALGAM中各獨(dú)立算法求解多目標(biāo)模型的優(yōu)化結(jié)果Fig.9 Results of multi-objective optimization using individual algorithmsin AMALGAM

圖10 每一代子代種群中各獨(dú)立算法貢獻(xiàn)的個體數(shù)Fig.10 The number of individuals produced by each algorithm in AMALGAM

采用多目標(biāo)模型進(jìn)行地下水非水相污染物修復(fù)設(shè)計(jì)，由AMALGAM求解一次多目標(biāo)模型就可得到一系列較優(yōu)的 Pareto解，可根據(jù)實(shí)際情況選擇不同的Pareto解作為相應(yīng)的修復(fù)方案。這是多目標(biāo)優(yōu)化模型的最顯著優(yōu)勢。相較于獨(dú)立的優(yōu)化算法和其他多目標(biāo)優(yōu)化算法，AMALGAM收斂速度更快;在相同的計(jì)算成本條件下，能夠得到更光滑連續(xù)的Pareto鋒面，節(jié)約計(jì)算成本，減少得到最優(yōu)解集的時間。

5 結(jié)論

(1)無污染物情況下對蒸汽井參數(shù)的敏感性分析結(jié)果表明，蒸汽溫度(焓)、蒸汽井注入速率對于蒸汽的擴(kuò)散鋒面影響較大。

(2)注入蒸汽的同時注入空氣能改變蒸汽鋒面的形狀，通過控制空氣的注入量能控制蒸汽的擴(kuò)散方向;用同時注入蒸汽和空氣取代只注入蒸汽能夠抑制污染物在隔水底板的再冷凝現(xiàn)象，去除含水層底部的DNAPL污染物。

(3)利用混合優(yōu)化框架的AMALGAM優(yōu)化算法則能綜合多種獨(dú)立優(yōu)化算法的優(yōu)點(diǎn)，根據(jù)不同模型自適應(yīng)選擇獨(dú)立算法的貢獻(xiàn)率，收斂速度更快，能夠在計(jì)算種群代數(shù)盡可能少的情況下得出平滑完整的Pareto解。

(4)相較于單向流模型的多目標(biāo)優(yōu)化，多項(xiàng)流模型的多目標(biāo)優(yōu)化更為復(fù)雜，一次正問題的計(jì)算時間也更長，計(jì)算成本高，未來可考慮使用替代模型等進(jìn)行多項(xiàng)流模擬優(yōu)化，在保證計(jì)算精度的前提下可縮短計(jì)算時間。

［1］ 施小清，姜蓓蕾，吳吉春，等.非均質(zhì)介質(zhì)中重非水相污染物運(yùn)移受泄漏速率影響數(shù)值分析［J］.水科學(xué)進(jìn)展，2012，23(3):376-382.［SHI X Q，JIANG B L，WU J C，et al.Numerical analysis of the effect of leakage rate on dense non-aqueous phase liquid transport in heterogonous porous media［J］.Advances in Water Science，2012，23(3):376-382.(in Chinese)］

［2］ 鄭葦，趙勇勝，秦傳玉，等.表面活性劑強(qiáng)化空氣擾動修復(fù)污染地下水影響區(qū)域研究［J］.環(huán)境科學(xué)學(xué)報(bào)，2011，31(1):102-106.［ZHENG W，ZHAO Y S，QIN C Y，et al.Effects of srucfactant-enhanced air sparging on the influence zone in treatmentof contaminated groundwater［J］. Acta Scientiae Circumstantiae，2011， 31(1):102-106.(in Chinese)］

［3］ 白靜，孫超，趙勇勝.地下水循環(huán)井技術(shù)對含水層典型NAPL污染物的修復(fù)模擬［J］.環(huán)境科學(xué)研究，2014，27(1):78-85.［BAI J，SUN C，ZHAO Y S.Remediation simulation of groundwater circulation well of NAPL contaminated aquifer［J］.Research of Envirnmental Sciences，2014，27(1):78-85.(in Chinese)］

［4］ Johnson R H，Poeter E P.Insights into the use of time-lapse GPR data asobservationsforinverse multiphase flow simulations of DNAPL migration［J］.Journal of Contaminant Hydrology，2007，89(1/2):136-155.

［5］ 陳華清，李連義.地下水苯系物污染原位曝氣修復(fù)模擬研究［J］.中國環(huán)境科學(xué)，2010，30(1):46-51.［CHEN H Q，LI L Y.Simulation on in-situ air sparging remediation of BTEX-contaminated groundwater［J］.China Environmental Science，2010，30(1):46-51.(in Chinese)］

［6］ Finsterle S.Demonstration of optimization techniques for groundwater plume remediation using iTOUGH2［J］.Environmental Modelling ＆ Software，2006，21(5):665-680.

［7］ Kaslusky S F，Udell K S.A theoretical model of air and steam co-injection to preventthe downward migration of DNAPLs during steam-enhanced extraction［J］.Journal of Contaminant Hydrology，2002，55(3/4):213-232.

［8］ 王志強(qiáng)，武強(qiáng)，鄒祖光，等.地下水石油污染曝氣治理技術(shù)研究［J］.環(huán)境科學(xué)，2007，28(4):754-760.［WANG Z Q，WU Q，ZOU Z G，et al.Study on the Groundwater Petroleum Contaminant Remediation by Air Sparging［J］.Chinese Journal of Environmental Science，2007，28(4):754-760.(in Chinese)］

［9］ 王貞國，梁偉，楊詢昌.曝氣法對石油污染土壤的修復(fù)研討［J］.山東國土資源，2006，22(6/7):107-108.［WANG Z G，LIANG W，YANG X C.Soil recovery using air sparging method polluted by oil［J］.Land and Resources in Shandong Province，2006，22(6/7):107-108.(in Chinese)］

［10］ Peng S，Wang N，Chen J.Steam and air co-injection in removing residual TCE in unsaturated layered sandy porous media［J］.Jounrnal of Contaminant Hydrology，2013，153:24-36.

［11］ McCray J E，F(xiàn)alta R W.Numerical Simulation of Air Sparging for Remediation of NAPL Contaminant［J］.Ground Water，1997，35(1):99-110.

［12］ 駱乾坤，吳劍鋒，林錦，等.地下水污染監(jiān)測網(wǎng)多目標(biāo)優(yōu)化設(shè)計(jì)模型及進(jìn)化求解［J］.水文地質(zhì)工程地質(zhì)，2013，40(5):97-103.［LUO Q K，WU J F，LIN J，etal. An evolutionary-based multi-objective optimization model for groundwater monitoring network design［J］.Hydrogeology ＆ Engineering Geology，2013，40(5):97-103.(in Chinese)］

［13］ 駱乾坤.進(jìn)化算法及其在地下水污染修復(fù)系統(tǒng)優(yōu)化管理中的應(yīng)用研究［D］.南京:南京大學(xué)，2013.［LUO Q K. Developmentand application of evolutionary algorithms for optimal design of groundwaterremediation systems［D］. Nanjing:Nanjing University，2013.(in Chinese)］

［14］ 席裕庚，柴天佑，惲為民.遺傳算法綜述［J］.控制理論與應(yīng)用，1996，13(6):697-708.［XI Y G，CHAI T Y，YUN W M.Survey on genetic algorithm［J］.Control Theory and Application，1996，13(6):697-708.(in Chinese)］

［15］ 葛繼科，邱玉輝，吳春明，等.遺傳算法研究綜述［J］.計(jì)算機(jī)應(yīng)用研究，2008，25(10):2911-2916.［GE J K，QIU Y H，WU C M，et al.Summary of genetic algorithm research［J］. Application Research of Computers，2008，25(10):2911-2916.(in Chinese)］

［16］ 曾建潮，崔志華.一種保證全局收斂的PSO算法［J］. 計(jì)算機(jī)研究與發(fā)展，2004，41(8):1333-1338.［ZENG J C，CUI Z H.A guaranteed global convergence particle swarm optimizer［J］.Journal of Computer Research and Development，2004，41(8):1333-1338.(in Chinese)］

［17］ 劉波，王凌，金以慧.差分進(jìn)化算法研究進(jìn)展［J］.控制與決策，2007，22(7):721-729.［LIU B，WANG L，JIN Y H.Advances in differential evolution［J］.Control and Decision，2007，22(7):721-729.(in Chinese)］

［18］ 朱嵩，毛根海，劉國華，等.改進(jìn)的MCMC方法及其應(yīng)用［J］.水利學(xué)報(bào)，2009，40(8):1019-1023.［ZHU S，MAO G H，LIU G H，et al.Improved MCMC method and its application［J］.Journal of Hydraulic Engineering，2009，40(8):1019-1023.(in Chinese)］

［19］ Hubschwerlen N，Zhang K，Mayer G，et al.Using Tough2-MP on a cluster—optimization methodology and study of scalability［J］. Computers ＆Geosciences，2012，45:26-35.

［20］ Pruess K，Oldenburg C，Moridis G.TOUGH2 Users Guide ［R］. Berkeley:Lawrence Berkeley Laboratory，1999.

［21］ 施小清，張可霓，吳吉春.TOUGH2軟件的發(fā)展及應(yīng)用［J］.工程勘察，2009，10(10):29-39.［SHI X Q，ZHANG K N，WU J C.The history and application of TOUGH2 code［J］.Geotechnical Invesitigation ＆Survery，2009，10(10):29-39.(in Chinese)］

［22］ Kaslusky S F，Udell K S.Co-injection of air and steam forthe prevention ofthe downward migration of DNAPLs during steam enhanced extraction［J］.Journal of Contaminant Hydrology，2005，77:325-347.

［23］ Falta R W，Pruess K，F(xiàn)insterle S，et al.T2VOC User’s Guide［R］.Berkeley:Lawrence Berkeley Laboratory，1995.

［24］ Vrugt J A，Robinsion B A，Hyan J M.Self-Adaptive Multimethod Search for Global Optimization in Real-Parameter Spaces［J］.Transactions on Evolutionary Computation，2009，13(2):243-259.

［25］ Zhang X，Srinivasan R，Liew M V.On the use of multi-algorithm，genetically adaptive multi-objective method for multi-site calibration of the SWAT model［J］.Hydrological Processes，2010，24(8):955-969.

［26］ 吳劍鋒，彭偉，錢家忠，等.基于INPGA的地下水污染治理多多目標(biāo)優(yōu)化管理模型:Ⅰ理論方法與算例驗(yàn)證［J］.地質(zhì)論評，2011，57(2):277-284.［WU J F，PENG W，QIAN J Z，et al.INPGA-based multiobjective management model for optimal design of groundwater remediation system:Ⅰ.methodology［J］.Geological Review，2011，57(2):277-284.(in Chinese)］

［27］ Stone H L.Probability model for estimating threephase relative permeability［J］.Journal of Petroleum Technology，1970，22(2):214-218.

［28］ Parker J C，Lenhard R J，T Kuppusamy.A parametric model for constitutive properies governing multiphase flow in porous media［J］.Water Resources Research，1987，23(4):618-624.

［29］ Deb K，Pratap A，Agarwal S，et al.A Fast and Elitist Multiobjective Genetic Algorithm:NSGA-II［J］. Transactionson EvolutionaryComputation，2002，6(2):182-197.