• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    表面活性劑強(qiáng)化含水層修復(fù)模擬的尺度提升研究

    2019-12-26 04:32:32陳夢佳吳劍鋒孫曉敏吳吉春
    中國環(huán)境科學(xué) 2019年12期
    關(guān)鍵詞:拉普拉斯滲透系數(shù)外殼

    陳夢佳,吳劍鋒*,宋 健,孫曉敏,林 錦,吳吉春

    表面活性劑強(qiáng)化含水層修復(fù)模擬的尺度提升研究

    陳夢佳1,吳劍鋒1*,宋 健1,孫曉敏2,林 錦2,吳吉春1

    (1.南京大學(xué)地球科學(xué)與工程學(xué)院,水科學(xué)系,表生地球化學(xué)教育部重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210023;2.南京水利科學(xué)研究院,江蘇 南京 210029)

    為了尋找到能夠在保證模擬精度的前提下,大大提高表面活性劑強(qiáng)化含水層修復(fù)(Surfactant-enhanced aquifer remediation, SEAR)模型計(jì)算效率的合適的滲透系數(shù)尺度提升方法,在概略介紹拉普拉斯-外殼法的基礎(chǔ)上,利用算術(shù)平均尺度提升法與拉普拉斯-外殼法建立大尺度模型,與小尺度模型進(jìn)行對比.結(jié)果表明,拉普拉斯-外殼法所建模型對含水層污染物殘余質(zhì)量的最大計(jì)算誤差在所有情況下均優(yōu)于算術(shù)平均尺度提升法,含水層非均質(zhì)性越強(qiáng),拉普拉斯-外殼法優(yōu)越性越明顯;并且拉普拉斯-外殼法對于污染羽質(zhì)心位置與形狀的刻畫效果也更好.大尺度模型能大幅度減小SEAR的計(jì)算成本,應(yīng)用算術(shù)平均尺度提升法可減少至原計(jì)算時(shí)間的6.5%左右,應(yīng)用拉普拉斯-外殼法可減少至4.5%左右.

    數(shù)值模擬;含水層修復(fù);尺度提升;算術(shù)平均尺度提升法;拉普拉斯-外殼法

    隨著工業(yè)的迅速發(fā)展,非水相液體(Non- Aqueous Phase Liquid,NAPL)污染成為最具挑戰(zhàn)性的地下水污染問題之一[1-3].由于大部分NAPL具有低水溶性與高界面張力的性質(zhì),抽取-處理方法(pump-and-treat, PAT)等許多傳統(tǒng)的地下水治理方法難以取得較好的修復(fù)效果[4-6].表面活性劑強(qiáng)化含水層修復(fù)(Surfactant-enhanced aquifer remediation, SEAR)利用表面活性劑增加NAPL的水溶性或是減小NPAL的界面張力,以提高NAPL在地下水中的溶解度或增加NAPL在孔隙介質(zhì)中的流動性,再采用PAT方法治理,顯著提高了污染物去除效率[7-9].

    UTCHEM是一種考慮多組分、多相流三維運(yùn)移過程的化學(xué)驅(qū)數(shù)值模擬軟件,可模擬復(fù)雜的相行為、化學(xué)與物理變化及過程[10].Qi等[11]和Luo等[12]利用UTCHEM模擬SEAR過程尋找表面活性劑強(qiáng)化含水層修復(fù)的最優(yōu)方案.尋找最優(yōu)方案需要進(jìn)行多次運(yùn)算,耗費(fèi)極大計(jì)算成本,國內(nèi)外學(xué)者利用智能優(yōu)化算法、替代模型等方法提高計(jì)算效率[13-15],但較少考慮提升滲透系數(shù)尺度的方法來減少計(jì)算成本.

    滲透系數(shù)尺度提升是指在保持模擬精度的情況下,通過公式轉(zhuǎn)換、經(jīng)驗(yàn)確定等方法將小尺度下的滲透系數(shù)轉(zhuǎn)換為大尺度下的等效滲透系數(shù)[16].由于對計(jì)算精度與計(jì)算效率的追求,國內(nèi)外眾多學(xué)者針對各類情形提出了有效的滲透系數(shù)尺度提升法[17-21].拉普拉斯-外殼法(Laplacian with skin)針對簡單拉普拉斯法(simple Laplacian)的不足進(jìn)行了改進(jìn),由于并非基于經(jīng)驗(yàn),此法更能適用于各類場地[22-24].本文在對NAPL運(yùn)移模擬尺度提升研究[25]的基礎(chǔ)上,進(jìn)一步對SEAR過程進(jìn)行尺度提升研究,以求提高數(shù)值模型計(jì)算效率.

    本文主要目的是:(1)在滲透率對數(shù)方差不同的非均質(zhì)孔隙介質(zhì)中,與小尺度模型模擬結(jié)果對比,比較算術(shù)平均尺度提升法和拉普拉斯-外殼法的提升效果;(2)比較小尺度剖分和大尺度剖分所建立模型的模擬運(yùn)算時(shí)間.

    1 研究方法

    1.1 拉普拉斯-外殼法

    簡單拉普拉斯法是一種基于求解流動方程的滲透系數(shù)尺度提升方法,它假設(shè)所得的等效滲透系數(shù)張量的方向平行于大尺度網(wǎng)格邊界[17],所以當(dāng)大尺度網(wǎng)格內(nèi)小尺度滲透系數(shù)產(chǎn)生的總流量不平行于邊界時(shí),大尺度模型無法準(zhǔn)確再現(xiàn)小尺度模型的運(yùn)移行為;并且區(qū)域邊界水頭由人為給定[17]使得區(qū)域邊界的模擬不準(zhǔn)確.因此,針對這些缺陷,Gómez- Hernández假設(shè)等效滲透系數(shù)為全張量,并劃定區(qū)域外殼,提出拉普拉斯-外殼法[22].

    拉普拉斯-外殼法基于求解流動方程,為了保持模擬精度,需要滿足以下兩個(gè)條件:

    (1)大尺度網(wǎng)格(block)中的水頭與小尺度網(wǎng)格(cell)中的平均水頭相等,即:

    式中:表示大尺度網(wǎng)格范圍;n表示大尺度網(wǎng)格內(nèi)小尺度網(wǎng)格的數(shù)目;h(,)表示大尺度網(wǎng)格內(nèi)水頭,m;h(,)表示小尺度網(wǎng)格內(nèi)水頭,m.

    (2)通過大尺度網(wǎng)格中的流量與小尺度網(wǎng)格中的平均流量相等,即:

    式中:q()表示大尺度網(wǎng)格內(nèi)流量,m3/s;q()表示小尺度網(wǎng)格內(nèi)流量,m3/s.

    拉普拉斯-外殼法利用外殼區(qū)域估計(jì)大尺度網(wǎng)格的邊界條件,使得區(qū)域邊界的水流模擬比簡單拉普拉斯法中的人為給定更加符合實(shí)際,因此建立模型時(shí),選取的區(qū)域要比實(shí)際需計(jì)算區(qū)域稍大,如圖1,圖中同時(shí)給出了大尺度網(wǎng)格內(nèi)水頭和流量的定義方式.

    對于小尺度網(wǎng)格的多套邊界條件,計(jì)算小尺度網(wǎng)格中的流動:

    式中:表示流體密度,kg/m3;表示重力加速度,m/s2;表示小尺度網(wǎng)格中孔隙介質(zhì)的滲透率,m2;?h表示小尺度網(wǎng)格內(nèi)壓力水頭梯度;表示動力粘滯系數(shù), Pa·s.

    將等效滲透系數(shù)定義如下:

    運(yùn)用計(jì)算而得的大尺度網(wǎng)格中的參數(shù),建立超定系統(tǒng)方程如下:

    式中:kBxx,kBxy,kByy表示未知的等效滲透系數(shù)全張量的成分.

    最后,上述方程組可由最小二乘法[26]求解,得到通過拉普拉斯-外殼法尺度提升后的大尺度網(wǎng)格等效滲透系數(shù),Li等[27]將此方法在程序中實(shí)現(xiàn).通過UTCHEM程序進(jìn)一步計(jì)算大尺度模型中SEAR過程地下水中NAPL污染羽隨時(shí)間變化的空間分布.

    1.2 SEAR過程

    表面活性劑具有兩親性[28],在地下環(huán)境中可以起到增加有機(jī)物溶解度和降低油/水界面張力的作用.當(dāng)水相中的表面活性劑濃度超過臨界膠束濃度,即會形成疏水端朝內(nèi)、親水端向外的膠束,該膠束將不溶、難溶于水的有機(jī)物溶于其中,而達(dá)到增溶的效果;表面活性劑溶于水后在表面富集,疏水基豎在水面,親水基插入水中,即將水的表面變成了疏水基緊密排列的表面,從而降低界面張力[29].SEAR過程即是先將表面活性劑注入地下,然后利用各種PAT技術(shù)治理有機(jī)污染物.

    2 算例應(yīng)用

    2.1 算例概述

    表1 含水層參數(shù)與液相參數(shù)

    研究對象是150m×85m×5m的二維承壓含水層,地下水流向?yàn)樽宰笙蛴?由于某種人為原因,造成四氯乙烯(perchloroethylene, PCE)泄漏,并穿過弱透水層進(jìn)入含水層,泄漏量為1m3/d,泄漏了730d. UTCHEM模型中涉及的含水層參數(shù)與液相參數(shù)見表1,含水層中PCE的初始飽和度分布如圖2所示.

    圖2 含水層中PCE的初始飽和度與井的布置方案

    2.2 模型建立及尺度提升方案

    先建立小尺度模型作為大尺度模型的對照,根據(jù)參數(shù)與PCE原始分布將含水層平面上剖分為150×85個(gè)離散單元,每個(gè)單元為1m×1m,垂向?yàn)橐粚?則小尺度模型離散為12750個(gè)單元格.模型中場地的上下邊界為隔水邊界,左右邊界為定水頭邊界,左邊界水頭為100.075m,右為100m.在污染源上游布設(shè)兩口表面活性劑注入井,注入速率為100m3/d,表面活性劑濃度為1g/cm3,在下游布設(shè)兩口抽水井,抽水速率為100m3/d,井位如圖2.

    保持模型各參數(shù)不變,根據(jù)注入井與抽水井位置建立非均勻離散的大尺度模型,利用算術(shù)平均尺度提升法建立大尺度模型,與基于拉普拉斯-外殼法建立的模型對比,分析兩種提升方法的模擬效果.算術(shù)平均尺度提升法將原始的離散單元150×85(單元格大小1m×1m)提升到34×21,單元格數(shù)減少到714個(gè),如圖3a.拉普拉斯-外殼法建立的模型中,取區(qū)域邊緣為外殼,外殼上下邊界為隔水邊界,左右邊界為定水頭邊界,軸方向外殼的大小取小尺度的20網(wǎng)格,軸方向取小尺度的10網(wǎng)格,外殼不進(jìn)行計(jì)算,因此將離散單元提升到26×17,單元格數(shù)減少到442個(gè),如圖3b.

    圖3 算術(shù)平均尺度提升法模型(a)與拉普拉斯-外殼法模型(b)的非均勻剖分情況

    為了比較算術(shù)平均尺度提升法與拉普拉斯-外殼法在各種滲透率場中對污染物修復(fù)結(jié)果的刻畫效果,本文分別建立含水層滲透率對數(shù)方差為0.5、1、5的模型進(jìn)行對比,采用序貫高斯模擬生成均值為10-10m2的二維非均質(zhì)滲透率場[30-32],3種情況的滲透率分布如圖4所示.

    圖4 不同非均質(zhì)條件下的含水層滲透率分布

    3 結(jié)果與分析

    3.1 均質(zhì)情況

    表2 均質(zhì)滲透率場中污染羽在空間上的分布特征

    結(jié)合表2與圖5可知,大尺度模型對于場地內(nèi)PCE殘余質(zhì)量的估算在可接受范圍內(nèi),絕對誤差小于0.5′105kg,隨著修復(fù)時(shí)間的增加,絕對誤差呈逐漸增大趨勢,于修復(fù)后期趨于穩(wěn)定.由圖6與表2可知,大尺度模型可以大致刻畫PCE污染羽在含水層中的空間分布,刻畫精度隨著修復(fù)時(shí)間的增長而降低.以此可作為非均質(zhì)情況的對照.

    圖5 大尺度模型計(jì)算PCE殘余質(zhì)量的絕對誤差

    3.2 非均質(zhì)情況

    由圖7可得,3種情況下,對PCE殘余質(zhì)量的估計(jì)在整個(gè)修復(fù)期間內(nèi)均是拉普拉斯-外殼法較為準(zhǔn)確,并且含水層非均質(zhì)性越強(qiáng),拉普拉斯-外殼法的優(yōu)越性越明顯.從圖8中可知,對于抽水井單元中PCE質(zhì)量的估算,拉普拉斯-外殼法比算術(shù)平均法更優(yōu)越,更接近小尺度模型的計(jì)算結(jié)果.

    為了對比兩種尺度提升方法對于修復(fù)過程中污染羽空間分布的刻畫效果,圖9~10和表3~4分別給出了對應(yīng)3種滲透率場在300,600d時(shí)由不同模型計(jì)算得到的污染羽分布和污染羽空間矩.

    結(jié)合圖9和表3可得,SEAR過程經(jīng)過300d后,含水層內(nèi)PCE殘余量約為原來的1/2,兩種大尺度模型對于3類滲透率場中PCE殘余量的估算中均是拉普拉斯-外殼法更為準(zhǔn)確;對于質(zhì)心位置的估算,水流方向上,拉普拉斯-外殼法3種情況下都優(yōu)于算術(shù)平均尺度提升法,垂直于水流方向上,當(dāng)滲透率對數(shù)方差較大時(shí),拉普拉斯-外殼法的估計(jì)也更優(yōu)越;對于污染羽范圍的估計(jì),從圖中可以看出拉普拉斯-外殼法模型計(jì)算得的污染羽形狀與小尺度模型計(jì)算結(jié)果更為接近,二階矩的計(jì)算結(jié)果也證明了這一結(jié)論.

    圖8 大尺度模型計(jì)算抽水井單元中PCE質(zhì)量的絕對誤差

    表3 修復(fù)300d后非均質(zhì)滲透率場中污染羽在空間上的分布特征

    圖10與表4是SEAR過程進(jìn)行了600d之后的PCE殘留飽和度分布與空間矩計(jì)算結(jié)果,與修復(fù)300d之后的結(jié)果類似,拉普拉斯-外殼法對于PCE殘余質(zhì)量與污染羽質(zhì)心位置的估算比算術(shù)平均尺度提升法更準(zhǔn)確,對污染羽范圍的估計(jì)也更具優(yōu)越性.

    因此,結(jié)合前面的分析結(jié)果可以得到,相較于算術(shù)平均尺度提升法,利用拉普拉斯-外殼法建立模型可以更有效地刻畫SEAR過程中NAPL污染羽的變化,且含水層非均質(zhì)性越強(qiáng),此方法優(yōu)越性越明顯.

    表4 修復(fù)600d后非均質(zhì)滲透率場中污染羽在空間上的分布特征

    圖10 非均質(zhì)滲透率場中修復(fù)600d后PCE殘留飽和度分布

    3.3 計(jì)算效率

    為了測試大、小尺度模型的計(jì)算成本,本文分別記錄了各模型調(diào)用UTCHEM程序的計(jì)算時(shí)間,滲透系數(shù)尺度提升的計(jì)算由于耗時(shí)極短忽略不計(jì).本算例采用配置為Intel(R) Core(TM) i5處理器、3.2?GHz、3.4?G內(nèi)存的計(jì)算機(jī),具體運(yùn)行時(shí)間見表5.

    由表5可得,算術(shù)平均尺度提升后的模型運(yùn)算時(shí)間只占小尺度模型的6.5%左右,拉普拉斯-外殼法所建模型由于外殼部分不進(jìn)行計(jì)算,運(yùn)算效率提升更明顯,運(yùn)算時(shí)間只占小尺度模型的4.5%左右,因此,利用尺度提升方法來減少SEAR過程的計(jì)算成本是可行的.

    表5 不同尺度提升方法的模型運(yùn)行時(shí)間對比

    注:a時(shí)間比:尺度提升后模型計(jì)算時(shí)間與尺度提升前模型計(jì)算時(shí)間的百分比.

    4 結(jié)論

    4.1 含水層均質(zhì)情況下,隨著修復(fù)時(shí)間的增加,大尺度模型計(jì)算結(jié)果的絕對誤差逐漸增大,但誤差仍然在可接受范圍內(nèi),并且空間矩的計(jì)算表明大尺度模型對質(zhì)心位置與污染羽范圍的計(jì)算較準(zhǔn)確,因此,含水層均質(zhì)時(shí)大尺度模型可以基本重現(xiàn)SEAR過程.

    4.2 含水層非均質(zhì)情況下,拉普拉斯-外殼法比算術(shù)平均尺度提升法更能重現(xiàn)小尺度滲透率場,其所建模型在對污染物殘余質(zhì)量、污染羽質(zhì)心位置與污染羽范圍的估算更具優(yōu)越性,且非均質(zhì)性越強(qiáng),優(yōu)越性越明顯.拉普拉斯-外殼法對污染羽在水流方向上更準(zhǔn)確的刻畫使得這一方法值得在實(shí)踐中推廣.

    4.3 對各尺度模型計(jì)算成本的對比表明,算術(shù)平均尺度提升法與拉普拉斯-外殼法對滲透系數(shù)進(jìn)行尺度提升后的模型具有更高的計(jì)算效率,算術(shù)平均尺度提升法所建模型運(yùn)算時(shí)間約為小尺度模型的6.5%,拉普拉斯-外殼法所建模型約為4.5%.

    [1] Delshad M, Pope G A, Sepehrnoori K. A compositional simulator for modeling surfactant enhanced aquifer remediation, 1formulation [J]. Journal of Contaminant Hydrology, 1996,23(4):303-327.

    [2] Schaerlaekens J, Mertens J, Linden J V, et al. A multi-objective optimization framework for surfactant-enhanced remediation of DNAPL contaminations [J]. Journal of Contaminant Hydrology, 2006,86(3/4):176-194.

    [3] Wu M, Cheng Z, Wu J F, et al. Precise simulation of long-term DNAPL migration in heterogeneous porous media based on light transmission micro-tomography [J]. Journal of Environmental Chemical Engineering, 2017,5(1):725-734.

    [4] 李璐璐,趙勇勝,王賀飛,等.溶劑萃取分離地下強(qiáng)化抽出處理液中的污染物和表面活性劑 [J]. 中國環(huán)境科學(xué), 2014,34(4):912-916. Li L L, Zhao Y S, Wang H F, et al. Solvent extraction separation of surfactant and pollutants during SEAR [J]. China Environmental Science, 2014,34(4):912-916.

    [5] 宋 健,吳劍鋒,楊 蘊(yùn),等.基于含水層DNAPL污染修復(fù)替代模型的多目標(biāo)優(yōu)化研究 [J]. 中國環(huán)境科學(xué), 2016,36(11):3390-3396. Song J, Wu J F, Yang Y, et al. A Kriging-based surrogate model for multi-objective optimization of DNAPL-contaminated aquifer remediation [J]. China Environmental Science, 2016,36(11):3390- 3396.

    [6] Haley J L, Hanson B, Enfield C, et al. Evaluating the effectiveness of groundwater extraction systems [J]. Groundwater Monitoring & Remediation, 2010,11(1):19-124.

    [7] Jiang X, Lu W X, Hou Z Y, et al. Ensemble of surrogates-based optimization for identifying an optimal surfactant-enhanced aquifer remediation strategy at heterogeneous DNAPL-contaminated sites [J]. Computers & Geosciences, 2015,84(2015):37-45.

    [8] Pennell K D, Jin M, Abriola L M, et al. Surfactant enhanced remediation of soil columns contaminated by residual tetrachloroethylene [J]. Journal of Contaminant Hydrology, 1994,16(1): 35-53.

    [9] 程 洲,徐紅霞,吳吉春,等.地下水中Tween80對DNAPL運(yùn)移和分布的影響 [J]. 中國環(huán)境科學(xué), 2019,39(3):1068-1077. Cheng Z, Xu H X, Wu J C, et al. Effects of Tween 80 in groundwater on DNAPL migration and distribution [J]. China Environmental Science, 2019,39(3):1068-1077.

    [10] Austin U O T. Volume Ⅱ: Technical documentation for UTCHEM- 9.0 a three-dimensional chemical flood simulator [M]. University of Texas, USA, 2000:1-1.

    [11] Qi Q Y, Lu W, Miao T, et al. Application of ensemble surrogates and adaptive sequential sampling to optimal groundwater remediation design at DNAPLs-contaminated sites [J]. Journal of Contaminant Hydrology, 2017,207:31-38.

    [12] Luo J N, Lu W X. Sobol’ sensitivity analysis of NAPL-contaminated aquifer remediation process based on multiple surrogates [J]. Computers & Geosciences, 2014,67(2014):110-116.

    [13] Ouyang Q, Lu W, Hou Z, et al. Chance-constrained multi-objective optimization of groundwater remediation design at DNAPLs- contaminated sites using a multi-algorithm genetically adaptive method [J]. Journal of Contaminant Hydrology, 2017,200(2017): 15-23.

    [14] Qin X S, Huang G H, Chakma A, et al. Simulation-based process optimization for surfactant-enhanced aquifer remediation at heterogeneous DNAPL-contaminated sites [J]. Science of the Total Environment, 2007,381(1-3):17-37.

    [15] Luo J N, Lu W X. Comparison of surrogate models with different methods in groundwater remediation process [J]. Journal of Earth System Science, 2014,123(7):1579-1589.

    [16] Farmer C L. Upscaling: a review [J]. International Journal for Numerical Methods in Fluids, 2002,40(1/2):63-78.

    [17] Wen X H, Gómez-Hernández J J. Upscaling hydraulic conductivities in heterogeneous media: An overview [J]. Journal of Hydrology, 1996,183(1/2):ix-xxxii.

    [18] Mantoglou A, Gelhar L W. Effective hydraulic conductivities of transient unsaturated flow in stratified soils [J]. Water Resources Research, 1987,23(1):57–67.

    [19] Russo D. Upscaling of hydraulic conductivity in partially saturated heterogeneous porous formation [J]. Water Resources Research, 1992,28(28):397-409.

    [20] Haldorsen H H, Lake L W. A new approach to shale management in field scale simulation models [J]. Society of Petroleum Engineers, 1984,24(4):447-457.

    [21] Sanchez-Vila X, Carrera J, Girardi J P. Scale effects in transmissivity [J]. Journal of Hydrology, 1996,183(1):1-22.

    [22] Gómez-Hernández J J. A stochastic approach to the simulation of block conductivity values conditioned upon data measured at a smaller scale [D]. California: Stanford University, 1991.

    [23] Rubin Y, Gómez-Hernández J J. A stochastic approach to the problem of upscaling of conductivity in disordered media: Theory and unconditional numerical simulation [J]. Water Resources Research, 1990,26(4):691–701.

    [24] Fernandez-Garcia D, Gómez-Hernández J J. Impact of upscaling on solute transport: Traveltimes, scale dependence of dispersivity, and propagation of uncertainty [J]. Water Resources Research, 2007,43(2): W02423.

    [25] 陳夢佳,吳劍鋒,孫曉敏,等.地下水典型非水相液體污染運(yùn)移模擬的尺度提升研究 [J]. 水文地質(zhì)工程地質(zhì), 2019. Chen M J, Wu J F, Sun X M, et al. Upscaling of PCE transport modeling based on UTCHEM in heterogeneous porous media. Hydrogeology & Engineering Geology, 2019.

    [26] Press W H, Flannery B P, Teukolsky S A, et al. Numerical recipes in C [M]. Cambridge: Cambridge University Press, 1988:735.

    [27] Li L P, Zhou H Y, Gómez-HernándezJ J. Transport upscaling using multi-rate mass transfer in three-dimensional highly heterogeneous porous media [J]. Advances in Water Resources, 2011,34(4):478-489.

    [28] Mulligan C N, Yong R N, Gibbs B F. Surfactant-enhanced remediation of contaminated soil: areview [J]. Engineering Geology, 2001,60(1-4): 371-380.

    [29] 白 靜.表面活性劑強(qiáng)化地下水循環(huán)井技術(shù)修復(fù)NAPL污染含水層研究 [D]. 長春:吉林大學(xué), 2013. Bai J. Remediation of NAPL contaminated aquifer with surfactant- enhanced groundwater circulation well [D]. Changchun Jilin University, 2013.

    [30] West M R, Kueper B H. Numerical simulation of DNAPL source zone remediation with in situ chemical oxidation (ISCO) [J]. Advances in Water Resources, 2012,44(2012):126-139.

    [31] Schaerlaekens J, Carmeliet J, Feyen J. Multi-objective optimization of the setup of a surfactant-enhanced DNAPL remediation [J]. Environmental Science & Technology, 2005,39(7):2327-2333.

    [32] Liu L. Modeling for surfactant-enhanced groundwaterremediation processes at DNAPLs-contaminated sites [J]. Journalof Environmental Informatics, 2005,5(2):42-52.

    Upscaling of reactive transport modeling for surfactant-enhanced aquifer remediation.

    CHEN Meng-jia1, WU Jian-feng1*, SONG Jian1, SUN Xiao-min2, LIN Jin2, WU Ji-chun1

    (1.Key Laboratory of Surficial Geochemistry, Ministry of Education, Department of Hydrosciences, School of Earth Sciences and Engineering, Nanjing University, Nanjing 210023, China;2.Nanjing Hydraulic Research Institute, Nanjing 210029, China)., 2019,39(12):5040~5048

    In order to find an appropriate permeability upscaling method which can greatly improve the computational efficiency of surfactant-enhanced aquifer remediation (SEAR) model while ensuring the accuracy of the simulation, two coarse-scale models based on Laplacian with skin method and the arithmetic mean upscaling method were applied in this study, and their performance was also compared with fine-scale models. The results indicated that the maximum calculation error of the residual mass of aquifer pollutants by the model based on Laplacian with skin was better than the model built by the arithmetic mean upscaling method in all cases. The superiority of Laplacian with skin became more significant when the aquifer heterogeneity was stronger. In addition, Laplacian with skin achieved a better simulation effect on the centroid location and the shape of the pollutant. Particularly, the use of coarse-scale model can greatly reduce the computational cost of SEAR, leading to the significant computational cost-savings, e.g., about 6.5% of the original runtime by using the arithmetic mean upscaling method and 4.5% by Laplacian with skin.

    numerical simulation;aquifer remediation;upscaling;arithmetic mean upscaling method;Laplacian with skin

    X523;P641

    A

    1000-6923(2019)12-5040-09

    陳夢佳(1995-),女,南京大學(xué)碩士研究生,主要從事地下水?dāng)?shù)值模擬研究.

    2019-05-21

    國家自然科學(xué)基金資助項(xiàng)目(41772254);國家重點(diǎn)研發(fā)計(jì)劃(2016YFC0402800)

    * 責(zé)任作者, 教授, jfwu@nju.edu.cn

    猜你喜歡
    拉普拉斯滲透系數(shù)外殼
    基于Origin的滲透系數(shù)衰減方程在地?zé)崴毓嘀械膽?yīng)用
    U盾外殼組件注塑模具設(shè)計(jì)
    塑料外殼注射模設(shè)計(jì)
    模具制造(2019年7期)2019-09-25 07:30:00
    多孔材料水滲透系數(shù)預(yù)測的隨機(jī)行走法
    輸水渠防滲墻及基巖滲透系數(shù)敏感性分析
    隔爆外殼水壓試驗(yàn)工藝探討
    河北平原新近系熱儲層滲透系數(shù)規(guī)律性分析
    基于超拉普拉斯分布的磁化率重建算法
    位移性在拉普拉斯變換中的應(yīng)用
    含有一個(gè)參數(shù)的p-拉普拉斯方程正解的存在性
    国内揄拍国产精品人妻在线| 中文在线观看免费www的网站| 精品免费久久久久久久清纯| www日本黄色视频网| 亚洲中文字幕日韩| 中文资源天堂在线| 老师上课跳d突然被开到最大视频| 欧美成人精品欧美一级黄| 久久精品夜色国产| 亚洲av免费高清在线观看| 青青草视频在线视频观看| 日韩强制内射视频| 亚洲av免费在线观看| 亚洲国产精品sss在线观看| 午夜福利在线在线| 我的老师免费观看完整版| 看十八女毛片水多多多| 看十八女毛片水多多多| 国产毛片a区久久久久| 精品久久久噜噜| 乱码一卡2卡4卡精品| 亚洲国产精品久久男人天堂| 久久久久久久国产电影| 18禁动态无遮挡网站| 国产麻豆成人av免费视频| 能在线免费看毛片的网站| 一二三四中文在线观看免费高清| 亚洲最大成人av| 亚洲国产欧洲综合997久久,| 国产亚洲91精品色在线| 亚洲精品乱久久久久久| 97人妻精品一区二区三区麻豆| 国产三级中文精品| 成人毛片60女人毛片免费| 亚洲va在线va天堂va国产| 黄片wwwwww| 波多野结衣高清无吗| 边亲边吃奶的免费视频| 免费观看在线日韩| 你懂的网址亚洲精品在线观看 | 国产私拍福利视频在线观看| 久久精品国产自在天天线| 神马国产精品三级电影在线观看| 在线观看av片永久免费下载| 午夜福利网站1000一区二区三区| av福利片在线观看| 精品少妇黑人巨大在线播放 | 麻豆国产97在线/欧美| 成人特级av手机在线观看| 国产一区二区在线av高清观看| 一边亲一边摸免费视频| 亚洲自偷自拍三级| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 99九九线精品视频在线观看视频| 久久久久久久午夜电影| 国国产精品蜜臀av免费| 免费不卡的大黄色大毛片视频在线观看 | 真实男女啪啪啪动态图| 国产黄色视频一区二区在线观看 | 国产白丝娇喘喷水9色精品| 欧美日本视频| 国产伦精品一区二区三区四那| 99久久精品一区二区三区| 国产欧美日韩精品一区二区| 国产高清三级在线| 真实男女啪啪啪动态图| 久久久a久久爽久久v久久| 国产乱人偷精品视频| 啦啦啦啦在线视频资源| 欧美成人精品欧美一级黄| 超碰97精品在线观看| 女人被狂操c到高潮| 蜜桃亚洲精品一区二区三区| 成人亚洲精品av一区二区| 色视频www国产| 午夜激情福利司机影院| 国产美女午夜福利| 亚洲成人久久爱视频| 国产人妻一区二区三区在| 亚洲国产色片| 久久久国产成人免费| 午夜爱爱视频在线播放| 欧美色视频一区免费| 校园人妻丝袜中文字幕| 91久久精品国产一区二区成人| 国产成人a∨麻豆精品| 69人妻影院| 狠狠狠狠99中文字幕| 亚洲av熟女| 日韩高清综合在线| 久久午夜福利片| 毛片女人毛片| 最近2019中文字幕mv第一页| 国产精品国产三级国产专区5o | 只有这里有精品99| 精华霜和精华液先用哪个| 婷婷六月久久综合丁香| 国产av不卡久久| 成年女人永久免费观看视频| 久久人人爽人人爽人人片va| 最近视频中文字幕2019在线8| 深爱激情五月婷婷| 日本爱情动作片www.在线观看| av在线播放精品| 久久韩国三级中文字幕| 久久久久网色| 麻豆av噜噜一区二区三区| 久久久久久久国产电影| 又爽又黄a免费视频| 亚洲在线观看片| 久久这里有精品视频免费| 欧美日韩国产亚洲二区| 欧美激情久久久久久爽电影| 精品国产露脸久久av麻豆 | 亚洲三级黄色毛片| 特级一级黄色大片| 欧美最新免费一区二区三区| 夜夜看夜夜爽夜夜摸| 欧美又色又爽又黄视频| 国产精品av视频在线免费观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 春色校园在线视频观看| 乱系列少妇在线播放| 干丝袜人妻中文字幕| 亚洲国产精品国产精品| 午夜久久久久精精品| 可以在线观看毛片的网站| 伦精品一区二区三区| 国产精华一区二区三区| 51国产日韩欧美| 久久99热6这里只有精品| 别揉我奶头 嗯啊视频| 午夜精品一区二区三区免费看| 小蜜桃在线观看免费完整版高清| 精品少妇黑人巨大在线播放 | 99久久无色码亚洲精品果冻| 99久久精品国产国产毛片| 狂野欧美白嫩少妇大欣赏| 中文字幕免费在线视频6| 国产成人精品一,二区| 日日摸夜夜添夜夜添av毛片| 天天一区二区日本电影三级| 国产精品一区二区在线观看99 | 午夜精品一区二区三区免费看| 高清视频免费观看一区二区 | 亚洲欧美精品自产自拍| 22中文网久久字幕| 深夜a级毛片| 亚洲怡红院男人天堂| 噜噜噜噜噜久久久久久91| 91精品伊人久久大香线蕉| 国产精品美女特级片免费视频播放器| 亚洲国产精品久久男人天堂| 国产成人午夜福利电影在线观看| 欧美不卡视频在线免费观看| 国产午夜福利久久久久久| 亚洲自拍偷在线| 亚洲av免费在线观看| 69人妻影院| 精品人妻熟女av久视频| 亚洲精品乱码久久久v下载方式| 国产精品人妻久久久久久| 国产av码专区亚洲av| 少妇人妻一区二区三区视频| АⅤ资源中文在线天堂| 卡戴珊不雅视频在线播放| 久99久视频精品免费| 午夜精品国产一区二区电影 | 日韩欧美精品v在线| 久久久久久伊人网av| a级毛片免费高清观看在线播放| 亚洲五月天丁香| 久久精品国产亚洲网站| 99热这里只有是精品50| 伦理电影大哥的女人| av在线蜜桃| 婷婷六月久久综合丁香| 超碰av人人做人人爽久久| 毛片女人毛片| 蜜桃亚洲精品一区二区三区| 欧美97在线视频| 久久亚洲国产成人精品v| 人体艺术视频欧美日本| 国产一区二区三区av在线| 日日撸夜夜添| 免费看美女性在线毛片视频| 大香蕉97超碰在线| 中文字幕免费在线视频6| 中国美白少妇内射xxxbb| 亚洲成人av在线免费| 2021天堂中文幕一二区在线观| 高清日韩中文字幕在线| av视频在线观看入口| 美女高潮的动态| 99热这里只有是精品50| 亚洲国产精品合色在线| 中文字幕av成人在线电影| av国产免费在线观看| 成人特级av手机在线观看| 99国产精品一区二区蜜桃av| 简卡轻食公司| 久久草成人影院| 国产色爽女视频免费观看| 国产伦在线观看视频一区| 亚洲av二区三区四区| 国产精品1区2区在线观看.| 日韩av在线免费看完整版不卡| 久久精品国产自在天天线| 噜噜噜噜噜久久久久久91| 大香蕉97超碰在线| 特大巨黑吊av在线直播| a级一级毛片免费在线观看| 国产熟女欧美一区二区| 亚洲精品日韩在线中文字幕| 日日撸夜夜添| 我要搜黄色片| 亚洲va在线va天堂va国产| 又爽又黄无遮挡网站| 国产欧美日韩精品一区二区| 国产精品1区2区在线观看.| 国产三级中文精品| 99热网站在线观看| 欧美性感艳星| 亚洲成人久久爱视频| 观看免费一级毛片| 国内精品一区二区在线观看| 国内少妇人妻偷人精品xxx网站| 亚洲综合色惰| 精品久久久久久久末码| 超碰97精品在线观看| 夜夜爽夜夜爽视频| 免费无遮挡裸体视频| 91av网一区二区| 国产高清国产精品国产三级 | 一边摸一边抽搐一进一小说| 精品久久久久久成人av| av在线蜜桃| 亚洲精品乱久久久久久| 成人午夜高清在线视频| 免费黄色在线免费观看| 丝袜喷水一区| 久久久久免费精品人妻一区二区| 插逼视频在线观看| 七月丁香在线播放| 免费观看性生交大片5| 国产麻豆成人av免费视频| 欧美潮喷喷水| 国产熟女欧美一区二区| 久久精品夜夜夜夜夜久久蜜豆| 观看美女的网站| 三级毛片av免费| 久久这里有精品视频免费| 国产极品精品免费视频能看的| 久久亚洲精品不卡| 综合色丁香网| 热99re8久久精品国产| 精品久久久久久久人妻蜜臀av| 亚洲欧美精品自产自拍| 日韩欧美三级三区| 免费观看性生交大片5| 国产精品久久久久久久久免| 一个人看视频在线观看www免费| 亚洲精品久久久久久婷婷小说 | 麻豆成人午夜福利视频| 老女人水多毛片| 欧美日韩综合久久久久久| 国产精品久久久久久精品电影| 卡戴珊不雅视频在线播放| 欧美3d第一页| 91狼人影院| 国产精品无大码| 国产精品女同一区二区软件| 久久人人爽人人片av| a级一级毛片免费在线观看| 身体一侧抽搐| 国产精品国产三级国产av玫瑰| 亚洲综合精品二区| 色尼玛亚洲综合影院| 国产高清不卡午夜福利| 精品一区二区免费观看| 丰满少妇做爰视频| 一卡2卡三卡四卡精品乱码亚洲| 精品酒店卫生间| 五月伊人婷婷丁香| 免费观看人在逋| 免费看av在线观看网站| 日韩在线高清观看一区二区三区| 在线观看av片永久免费下载| 国产成人精品婷婷| 亚洲欧美日韩高清专用| 青春草视频在线免费观看| 精品一区二区三区人妻视频| 国产精品一区二区三区四区久久| 2021少妇久久久久久久久久久| 中文字幕久久专区| 亚洲精品亚洲一区二区| 男女国产视频网站| 国产成人aa在线观看| 国产精品不卡视频一区二区| 亚洲自拍偷在线| 日韩av在线大香蕉| 亚洲欧美精品综合久久99| 老女人水多毛片| 91精品国产九色| videossex国产| 一级毛片电影观看 | 免费看a级黄色片| 国产黄片美女视频| 精品人妻偷拍中文字幕| 国产一区二区在线av高清观看| 能在线免费观看的黄片| 国产欧美日韩精品一区二区| 人妻夜夜爽99麻豆av| 成年女人永久免费观看视频| 免费黄色在线免费观看| 国产午夜精品一二区理论片| 欧美日本亚洲视频在线播放| 国产黄片视频在线免费观看| 亚洲欧美精品综合久久99| 99视频精品全部免费 在线| 亚洲av.av天堂| 国产精品嫩草影院av在线观看| 亚洲三级黄色毛片| 波多野结衣高清无吗| 国产精品女同一区二区软件| av线在线观看网站| 亚洲经典国产精华液单| 国产精品一及| 午夜爱爱视频在线播放| 欧美日本视频| 久久亚洲精品不卡| 午夜老司机福利剧场| 久久久a久久爽久久v久久| 日韩人妻高清精品专区| av在线天堂中文字幕| 中文字幕亚洲精品专区| av天堂中文字幕网| 久久人人爽人人爽人人片va| 青春草国产在线视频| 黄色一级大片看看| 日日撸夜夜添| 中文乱码字字幕精品一区二区三区 | 国产视频首页在线观看| 国产69精品久久久久777片| 亚洲av二区三区四区| 亚洲av福利一区| 久久精品国产亚洲网站| 亚洲在久久综合| 久久这里只有精品中国| 亚洲四区av| 国产精品野战在线观看| 黄片无遮挡物在线观看| 午夜福利在线观看吧| 国产欧美日韩精品一区二区| 久久久亚洲精品成人影院| 综合色丁香网| 久久久久久久久久久丰满| 身体一侧抽搐| 美女xxoo啪啪120秒动态图| 亚洲精品乱码久久久v下载方式| 亚洲中文字幕一区二区三区有码在线看| 国产高清国产精品国产三级 | 两个人的视频大全免费| 九九热线精品视视频播放| 欧美成人精品欧美一级黄| 中文亚洲av片在线观看爽| 成年女人永久免费观看视频| 久久精品国产99精品国产亚洲性色| 欧美不卡视频在线免费观看| 亚洲熟妇中文字幕五十中出| 成人鲁丝片一二三区免费| 欧美区成人在线视频| 亚洲av成人av| 国产精品国产三级国产专区5o | 女人十人毛片免费观看3o分钟| 国产一区有黄有色的免费视频 | 99久久九九国产精品国产免费| 亚洲av.av天堂| 免费黄色在线免费观看| 欧美不卡视频在线免费观看| 精品午夜福利在线看| 欧美丝袜亚洲另类| 在线a可以看的网站| 欧美性猛交黑人性爽| 欧美高清性xxxxhd video| 久久热精品热| 一二三四中文在线观看免费高清| 啦啦啦观看免费观看视频高清| 国产精品福利在线免费观看| 日日摸夜夜添夜夜爱| 国产日韩欧美在线精品| 亚洲av中文av极速乱| 国产高清有码在线观看视频| 又黄又爽又刺激的免费视频.| 久久99精品国语久久久| 久久午夜福利片| 赤兔流量卡办理| 精品久久久久久久久亚洲| 日本一二三区视频观看| 日韩在线高清观看一区二区三区| 日韩制服骚丝袜av| 欧美日韩在线观看h| 美女大奶头视频| 免费观看的影片在线观看| 卡戴珊不雅视频在线播放| 中文字幕精品亚洲无线码一区| 伊人久久精品亚洲午夜| 亚洲熟妇中文字幕五十中出| 蜜桃亚洲精品一区二区三区| 人妻夜夜爽99麻豆av| 亚洲av熟女| 亚洲成av人片在线播放无| www.色视频.com| 国产乱来视频区| 夜夜看夜夜爽夜夜摸| 听说在线观看完整版免费高清| 亚洲天堂国产精品一区在线| 亚洲精品,欧美精品| 午夜福利视频1000在线观看| 国产精品乱码一区二三区的特点| 国产精品,欧美在线| 亚洲一区高清亚洲精品| 夜夜看夜夜爽夜夜摸| 日韩成人av中文字幕在线观看| 国产淫语在线视频| 婷婷六月久久综合丁香| 久久精品国产亚洲网站| 熟女电影av网| 天美传媒精品一区二区| 国产色婷婷99| 精华霜和精华液先用哪个| 国产片特级美女逼逼视频| www.av在线官网国产| 欧美精品国产亚洲| 老司机福利观看| 亚洲精品自拍成人| 毛片一级片免费看久久久久| 国产午夜精品一二区理论片| 高清毛片免费看| 91精品一卡2卡3卡4卡| 亚洲国产精品专区欧美| 国产亚洲一区二区精品| 免费观看在线日韩| 69人妻影院| 美女被艹到高潮喷水动态| 久久精品夜夜夜夜夜久久蜜豆| 日韩 亚洲 欧美在线| 三级经典国产精品| 深爱激情五月婷婷| 国产一区有黄有色的免费视频 | 国产精品野战在线观看| 久久久成人免费电影| 韩国高清视频一区二区三区| 国产一级毛片七仙女欲春2| 偷拍熟女少妇极品色| 美女被艹到高潮喷水动态| 久久亚洲精品不卡| 超碰97精品在线观看| 91久久精品电影网| 亚洲婷婷狠狠爱综合网| av在线亚洲专区| 久久精品综合一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 91aial.com中文字幕在线观看| 中文乱码字字幕精品一区二区三区 | 精品无人区乱码1区二区| 精品国产一区二区三区久久久樱花 | 成年女人看的毛片在线观看| 综合色丁香网| 秋霞伦理黄片| av免费在线看不卡| 精品99又大又爽又粗少妇毛片| 中文字幕av成人在线电影| 亚洲最大成人手机在线| 麻豆成人av视频| 国产亚洲91精品色在线| 亚洲天堂国产精品一区在线| 少妇丰满av| 久久99热6这里只有精品| av在线蜜桃| 亚洲综合精品二区| 国产精品国产高清国产av| 深爱激情五月婷婷| 午夜视频国产福利| 美女黄网站色视频| 一卡2卡三卡四卡精品乱码亚洲| 日韩人妻高清精品专区| 日韩高清综合在线| 成人高潮视频无遮挡免费网站| 国产精品精品国产色婷婷| 99久久精品一区二区三区| 亚洲精品色激情综合| 欧美成人免费av一区二区三区| 国产av在哪里看| 日韩亚洲欧美综合| 精品国产一区二区三区久久久樱花 | 你懂的网址亚洲精品在线观看 | h日本视频在线播放| 一级毛片电影观看 | av在线蜜桃| 最后的刺客免费高清国语| 一个人看的www免费观看视频| 欧美激情在线99| 欧美激情国产日韩精品一区| 国产视频首页在线观看| 亚洲丝袜综合中文字幕| 长腿黑丝高跟| 国产成人a区在线观看| 97超视频在线观看视频| av线在线观看网站| 少妇被粗大猛烈的视频| 国产单亲对白刺激| 国产精品国产三级专区第一集| 久久国内精品自在自线图片| 精品不卡国产一区二区三区| 亚洲国产精品成人综合色| 看免费成人av毛片| 噜噜噜噜噜久久久久久91| 久久这里有精品视频免费| 亚洲18禁久久av| 久久久精品大字幕| 亚洲一级一片aⅴ在线观看| 亚洲国产精品国产精品| 精品国产三级普通话版| 亚洲美女搞黄在线观看| av免费观看日本| 晚上一个人看的免费电影| 97人妻精品一区二区三区麻豆| 女人被狂操c到高潮| 国产 一区精品| 麻豆精品久久久久久蜜桃| 只有这里有精品99| 建设人人有责人人尽责人人享有的 | 国产乱人偷精品视频| 亚洲怡红院男人天堂| 国产精品永久免费网站| 久久亚洲精品不卡| 国产精品久久电影中文字幕| 国内揄拍国产精品人妻在线| 亚洲av不卡在线观看| 久久精品久久久久久噜噜老黄 | 搡女人真爽免费视频火全软件| 免费黄网站久久成人精品| 久久精品久久久久久噜噜老黄 | 久久精品久久久久久噜噜老黄 | 精品久久久久久久久亚洲| 国产一级毛片在线| 五月伊人婷婷丁香| 18禁动态无遮挡网站| 久久久久精品久久久久真实原创| 特大巨黑吊av在线直播| eeuss影院久久| 免费黄网站久久成人精品| 国产伦精品一区二区三区视频9| 国产久久久一区二区三区| 日本wwww免费看| 又爽又黄a免费视频| 97超视频在线观看视频| 午夜精品一区二区三区免费看| 国产成人精品婷婷| 日韩一区二区视频免费看| 久久久精品94久久精品| 欧美日本视频| 天美传媒精品一区二区| 亚洲最大成人中文| 毛片一级片免费看久久久久| 国产av码专区亚洲av| 日本欧美国产在线视频| 国产免费福利视频在线观看| 国产亚洲av嫩草精品影院| 大话2 男鬼变身卡| 51国产日韩欧美| 亚洲精品乱码久久久久久按摩| 青春草国产在线视频| 亚洲精品成人久久久久久| 内地一区二区视频在线| 免费av毛片视频| 老女人水多毛片| 变态另类丝袜制服| 97超碰精品成人国产| 成人漫画全彩无遮挡| 精品少妇黑人巨大在线播放 | 亚洲国产精品合色在线| 嫩草影院入口| 精品国产露脸久久av麻豆 | 久久久精品欧美日韩精品| 尾随美女入室| 亚洲av熟女| 99视频精品全部免费 在线| 国产老妇伦熟女老妇高清| 日韩一本色道免费dvd| 啦啦啦啦在线视频资源| 精品国产三级普通话版| 国产综合懂色| АⅤ资源中文在线天堂| a级一级毛片免费在线观看| 99九九线精品视频在线观看视频| 黄色配什么色好看| 日韩欧美精品免费久久| 身体一侧抽搐| 国产视频内射| 3wmmmm亚洲av在线观看| 国产精品永久免费网站| 国产高清视频在线观看网站| 亚洲人成网站在线观看播放| 日韩成人av中文字幕在线观看| 国产精品1区2区在线观看.| 99久国产av精品国产电影| 亚洲精品亚洲一区二区| 国产精品伦人一区二区| 热99在线观看视频| 久久久久久久亚洲中文字幕| 有码 亚洲区| 高清视频免费观看一区二区 | videossex国产|