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

    ERT法監(jiān)測(cè)高壓旋噴修復(fù)PAHs污染現(xiàn)場(chǎng)試驗(yàn)研究

    2023-11-28 12:53:44宋權(quán)威王颋軍杜顯元鄧永鋒
    中國環(huán)境科學(xué) 2023年11期
    關(guān)鍵詞:電阻率藥劑電阻

    柳 迪,宋權(quán)威,李 玲,胡 杰,王颋軍,杜顯元,鄧永鋒,柯 瀚

    ERT法監(jiān)測(cè)高壓旋噴修復(fù)PAHs污染現(xiàn)場(chǎng)試驗(yàn)研究

    柳 迪1,宋權(quán)威2,3,李 玲4,胡 杰1,王颋軍4,杜顯元2,3,鄧永鋒5,柯 瀚1

    (1.浙江大學(xué)巖土工程研究所,軟弱土與環(huán)境土工教育部重點(diǎn)試驗(yàn)室,浙江 杭州 310058;2.石油石化污染物控制與處理國家重點(diǎn)試驗(yàn)室,北京 102206;3.中國石油集團(tuán)安全環(huán)保技術(shù)研究院有限公司,北京 102206;4.中節(jié)能大地環(huán)境修復(fù)有限公司,北京 100085;5.東南大學(xué)交通學(xué)院,巖土工程研究所,江蘇 南京 210096)

    選用過硫酸鈉(PS)作為氧化劑,微米鐵/納米鐵/碳載鐵/NaOH作為激活劑,通過高壓旋噴工藝氧化修復(fù)多環(huán)芳烴(PAHs)污染土,利用高密度電阻率法(ERT)監(jiān)測(cè)旋噴過程中土壤的電阻值分布,分析了氧化藥劑在修復(fù)過程中的遷移轉(zhuǎn)化規(guī)律,比較了4組藥劑旋噴后48h的修復(fù)效果.結(jié)果表明:ERT準(zhǔn)確地反映了污染修復(fù)發(fā)生的位置,最大偏差小于0.5m;電阻變化值與萘修復(fù)值呈現(xiàn)正相關(guān)關(guān)系,通過擬合曲線估計(jì)的污染修復(fù)值與實(shí)際修復(fù)值之間存在誤差,誤差為10.3%~22.6%;碳載鐵粉顯示出良好的穩(wěn)定與緩釋性能,旋噴后穩(wěn)定氧化階段延長約8h;最終預(yù)設(shè)修復(fù)區(qū)域內(nèi)污染去除量:碳載鐵粉組>納米鐵粉組>微米鐵粉組≈NaOH組.

    高密度電阻率法;高壓旋噴;原位氧化;多環(huán)芳烴;現(xiàn)場(chǎng)試驗(yàn)

    由于環(huán)境監(jiān)察制度不嚴(yán)格和生產(chǎn)管理模式不健全,一些重污染企業(yè)在長期的歷史生產(chǎn)中,向土壤中排放了大量污染物,對(duì)當(dāng)?shù)厣鷳B(tài)環(huán)境造成了不利影響,危害人民群眾身體健康,嚴(yán)重影響土地的再次開發(fā)利用[1].多環(huán)芳烴(PAHs)作為常見的有機(jī)溶劑及生產(chǎn)原料的成分,在石油,煉焦,有機(jī)等化工行業(yè)應(yīng)用廣泛,可通過事故泄露,廢水排放,浸出液入滲等多種途徑進(jìn)入土壤[2].PAHs污染具有難降解,難遷移,易揮發(fā),毒性大等特點(diǎn),傳統(tǒng)物理修復(fù)技術(shù)難以對(duì)其進(jìn)行有效的修復(fù)治理.

    近些年來,高壓旋噴氧化技術(shù)成為國內(nèi)外處理PAHs污染場(chǎng)地的主要方式之一[3-4],其具有高效率,低成本的優(yōu)點(diǎn),能有效處理復(fù)雜地層中多種有機(jī)物形成的復(fù)合污染[5].對(duì)于PAHs污染,其主要污染區(qū)的穩(wěn)定表征和監(jiān)測(cè)結(jié)果對(duì)污染修復(fù)設(shè)計(jì),方案優(yōu)化和效果評(píng)估至關(guān)重要[6].然而旋噴試驗(yàn)過程中由于取樣困難,現(xiàn)有試驗(yàn)大多采取開挖的方式檢測(cè)旋噴效果[7],旋噴后的取樣結(jié)果無法準(zhǔn)確描述旋噴過程中氧化劑的遷移轉(zhuǎn)化過程[8-9];特別是當(dāng)污染土層表現(xiàn)出較為顯著的非均質(zhì)性時(shí),常規(guī)性取樣調(diào)查和地下水采樣難以達(dá)到預(yù)期監(jiān)測(cè)目標(biāo).高密度電阻率法(ERT)通過測(cè)試污染引起的土壤電阻率變化可以較好的解決監(jiān)測(cè)困難的問題,同時(shí)ERT監(jiān)測(cè)具有非侵入,可連續(xù),信息豐富等特點(diǎn),在石油勘探,污染泄露監(jiān)測(cè)等方面應(yīng)用廣泛[10],而針對(duì)較低濃度有機(jī)污染修復(fù)監(jiān)測(cè)的研究較少,并且現(xiàn)場(chǎng)監(jiān)測(cè)效果需進(jìn)一步驗(yàn)證[11-12].

    本文選取微米/納米/碳負(fù)載零價(jià)鐵,NaOH作為過硫酸鈉激活劑進(jìn)行高壓旋噴原位氧化試驗(yàn),通過ERT監(jiān)測(cè)旋噴過程中土壤的電阻值分布,以分析氧化藥劑在修復(fù)過程中的遷移轉(zhuǎn)化規(guī)律;結(jié)合取樣測(cè)試PAHs的濃度分布,對(duì)ERT測(cè)試結(jié)果進(jìn)行分析比較,解釋說明ERT在監(jiān)測(cè)PAHs污染修復(fù)過程的優(yōu)勢(shì)及局限性.

    1 材料與方法

    1.1 污染場(chǎng)地概況

    現(xiàn)場(chǎng)試驗(yàn)地點(diǎn)位于華北某煤氣化工廠遺留場(chǎng)地.該場(chǎng)地曾進(jìn)行煤化工產(chǎn)品生產(chǎn),煤矸石發(fā)電供熱等生產(chǎn)作業(yè),近40年的生產(chǎn)過程中含PAHs的生產(chǎn)原料浸出液滲入地基土中,其較強(qiáng)的揮發(fā)性造成較為嚴(yán)重的土壤,空氣污染,影響周邊居民生活及身體健康[11-12].

    圖1 試驗(yàn)場(chǎng)地區(qū)域位置

    1.1.1 試驗(yàn)區(qū)域條件 試驗(yàn)旋噴區(qū)域?yàn)?m寬,8m長的矩形區(qū)域,在試驗(yàn)區(qū)域自東向西間隔1m進(jìn)行取樣勘測(cè),場(chǎng)地9m的勘察深度內(nèi),自上而下為粉土層,粉質(zhì)砂土層,卵石層,粉質(zhì)黏土層;土壤理化性質(zhì)參照《土工試驗(yàn)方法標(biāo)準(zhǔn)GB/T 50123-2019》[13]及《土的工程分類標(biāo)準(zhǔn) GB/T 50145-2007》[14]進(jìn)行測(cè)試,取樣點(diǎn)位標(biāo)注在圖2土層分布圖中(水平坐標(biāo)0m為高密度電阻率測(cè)線起始點(diǎn)),土壤理化性質(zhì)見表1.

    圖2 試驗(yàn)區(qū)域地層分布

    地下水埋深約3.0~3.5m,屬孔隙潛水,主要賦存于粉質(zhì)砂土層及更深部的卵石層中,地下水主要受大氣降水補(bǔ)給,以蒸發(fā)和徑流方式排泄;地下水位受大氣降水影響顯著,旋噴試驗(yàn)后30~36h出現(xiàn)持續(xù)約6h的降雨天氣,地下水位上漲約為1.0m.試驗(yàn)區(qū)域地下水總體流向自西北向東南,水力梯度為2‰~4‰.

    表1 試驗(yàn)土層理化性質(zhì)

    1.1.2 試驗(yàn)場(chǎng)地污染分布 前期場(chǎng)地污染土取樣點(diǎn)位如圖2所示,取樣流程參照《土壤環(huán)境監(jiān)測(cè)技術(shù)規(guī)范HJ/T 166-2004》[15],《地塊土壤和地下水中揮發(fā)性有機(jī)物采樣技術(shù)導(dǎo)則 HJ 1019-2019》[16]進(jìn)行;有機(jī)物測(cè)試參照《土壤和沉積物揮發(fā)性有機(jī)物的測(cè)定吹掃捕集/氣象色譜-質(zhì)譜法 HJ 605-2011》[17]進(jìn)行;鐵離子測(cè)試參照《固體廢物金屬元素的測(cè)定電感耦合等離子體質(zhì)譜法HJ 766-2015》[18],《土壤 8種有效態(tài)元素的測(cè)定二乙烯三胺五乙酸浸提-電感耦合等離子體發(fā)射光譜法HJ 804-2016》[19]進(jìn)行.

    根據(jù)測(cè)試結(jié)果,場(chǎng)地內(nèi)主要超標(biāo)污染物為萘和苯,最大濃度分別達(dá)到586,22.2mg/kg,約為《土壤環(huán)境質(zhì)量建設(shè)用地土壤污染風(fēng)險(xiǎn)管控標(biāo)準(zhǔn)(試行) GB36600-2018》[20]中二類工業(yè)用地篩選值的8.15,6.85倍.

    由于萘的污染濃度遠(yuǎn)超苯,且在各點(diǎn)位均有檢出,因而本文選取萘作為特征污染物來描述場(chǎng)地整體污染及修復(fù)情況.萘主要分布于深度為3.0~9.0m含水層內(nèi),含水層內(nèi)萘濃度均大于200mg/kg,其中3.8~6.5m深度為污染濃度超過400mg/kg,此主污染區(qū)域也為本試驗(yàn)的目標(biāo)修復(fù)區(qū)域,詳細(xì)污染分布情況如圖3所示.

    1.2 現(xiàn)場(chǎng)高壓旋噴試驗(yàn)設(shè)計(jì)

    1.2.1 試驗(yàn)內(nèi)容 高壓旋噴工藝修復(fù)PAHs污染場(chǎng)地使用的氧化劑類型包括過氧化氫,Fenton試劑,高錳酸鉀,過硫酸鹽(PS)等,其中通過零價(jià)鐵活化過硫酸鈉的氧化體系由于原料易得,反應(yīng)條件寬松,反應(yīng)產(chǎn)物友好等優(yōu)點(diǎn)被廣泛應(yīng)用[21-22].Fe2+可活化過硫酸根(S2O82-)產(chǎn)生具有強(qiáng)氧化能力的硫酸根自由基(SO4-×),無選擇性的與有機(jī)物進(jìn)行氧化反應(yīng)[23].目前現(xiàn)場(chǎng)應(yīng)用零價(jià)鐵活化過硫酸鈉現(xiàn)場(chǎng)尺度的研究較少,其工程應(yīng)用效果有待進(jìn)一步驗(yàn)證.

    圖3 沿深度方向萘的分布

    為研究不同種類零價(jià)鐵作為激活劑在高壓旋噴過程中的遷移規(guī)律及修復(fù)效果,選取過硫酸鈉Na2S2O8(PS)作為氧化劑,微米鐵粉(Fe0[μm]),納米鐵粉(Fe0[m])及碳載零價(jià)鐵(BC-nZVI)3種不同種類鐵粉作為激活劑,另設(shè)置工程常見激活劑NaOH作為對(duì)照組,開展4類激活劑情況下的現(xiàn)場(chǎng)高壓旋噴試驗(yàn).

    為防止旋噴區(qū)域重疊,試驗(yàn)區(qū)域內(nèi)旋噴點(diǎn)位一字排開,與地下水流方向平行,相鄰孔間距約200cm,以高密度電阻率測(cè)線SY1起始點(diǎn)作為水平坐標(biāo)起點(diǎn),碳載零價(jià)鐵組,納米鐵粉組,微米鐵粉組,NaOH組旋噴中心分別位于水平坐標(biāo)15,17,19,21m處.旋噴作業(yè)從左到右依次進(jìn)行,旋噴深度為3.0~8.0m.

    旋噴后在各點(diǎn)位中心進(jìn)行取樣,分別測(cè)試深度3.0,4.5,6.5,8.5m處土壤中萘濃度,取樣流程及測(cè)試方法與節(jié)1.1.2中所述相同,試驗(yàn)平面布置及取樣點(diǎn)位見圖4.

    圖4 試驗(yàn)區(qū)域平面布置

    1.2.2 試驗(yàn)藥劑配比 參考現(xiàn)場(chǎng)常用旋噴藥劑使用情況[24],本試驗(yàn)對(duì)于萘含量為200~600mg/kg的污染區(qū)域,采用PS濃度為12%,NaOH:PS(物質(zhì)的量比)為2:1,旋噴水量/旋噴距離為300L/m的藥劑用量.

    根據(jù)室內(nèi)試驗(yàn)測(cè)試[25-26],過硫酸鈉與納米零價(jià)鐵的物質(zhì)量之比為 200:1,反應(yīng)溫度為60℃,其他條件保持最優(yōu)值時(shí),36h內(nèi)多環(huán)芳烴的降解率達(dá)到最高值;在過硫酸鈉與碳載零價(jià)鐵物質(zhì)的量之比為150:1,室溫條件下12h內(nèi)1,2-二氯乙烷(1,2-DCA)降解率最大可以達(dá)到78.45%.由于藥劑制備和旋噴過程中鐵粉氧化導(dǎo)致的激活劑損失.為保證修復(fù)效果,本試驗(yàn)零價(jià)鐵用量增加約20%,鐵基激活劑配比(物質(zhì)的量比)均為1:125.藥劑規(guī)格見表2,藥劑用量見表3.

    1.2.3 現(xiàn)場(chǎng)設(shè)備及工藝流程 高壓旋噴工藝的主要設(shè)備包括藥品攪拌裝置,高壓泵,高壓管道,控制臺(tái),鉆管,鉆頭,設(shè)備圖見圖5所示.每次注漿前60min,在配備有攪拌泵的藥劑桶中進(jìn)行過硫酸鈉,氫氧化鈉的混合,攪拌45min待藥劑充分溶解后加入Fe0粉末并繼續(xù)攪拌5min,攪拌期間攪拌桶處于密封狀態(tài)減少Fe0氧化;藥劑桶通過直徑90mm鋼絲軟管與XPB-90E型高壓注漿泵相連,并通過高壓管線輸送至ZP60型旋噴鉆機(jī)進(jìn)行旋噴施工.每次旋噴前于地面進(jìn)行試噴以確認(rèn)旋噴管道暢通,同時(shí)使管道內(nèi)殘余藥劑排出,此后按預(yù)定提升速度,鉆桿轉(zhuǎn)速,注漿壓力從上到下進(jìn)行旋噴作業(yè),在旋噴桿行至8.0m預(yù)定深度后按原旋噴參數(shù)從下至上進(jìn)行,提升至3.0m深度處結(jié)束旋噴.

    表2 試驗(yàn)藥劑規(guī)格

    表3 試驗(yàn)藥劑用量

    根據(jù)現(xiàn)場(chǎng)高壓旋噴修復(fù)工程中常用的旋噴參數(shù)范圍,設(shè)置注漿壓力為30MPa,提升速度為25cm/min,轉(zhuǎn)速為20r/min,單孔旋噴時(shí)長約為30min.

    1.2.4 高密度電阻率法(ERT)監(jiān)測(cè) 本試驗(yàn)使用儀器為賽盈地脈GD-20電法儀,采用Wenner- Schlumberger排列法進(jìn)行探測(cè),測(cè)量原理示意圖如圖6所示.測(cè)量時(shí),終端發(fā)射直流電以間距 a逐個(gè)測(cè)試電極間電流強(qiáng)度,測(cè)試結(jié)果的視電阻率為 Level 1層電阻率.待 Level 1層測(cè)試完成后,終端控制以3倍電極間距進(jìn)行測(cè)試,測(cè)試結(jié)果的視電阻率為 Level 2層電阻率,按照1,3,5倍類推,終端控制測(cè)試間距按照特定間隔逐漸增加,對(duì)應(yīng)的測(cè)量深度也逐漸加深,直至所需探測(cè)深度以及寬度后停止測(cè)試.

    圖6 Wenner-Schlumberger法測(cè)試原理示意

    當(dāng)供電主機(jī)向電極A及B供以一定的電流,電極M與N兩點(diǎn)間的電壓為D,則M與N兩點(diǎn)間的視電阻率為為:

    式中:s為視電阻率,W×;為裝置系數(shù),;DU為M與N兩點(diǎn)間電壓,mV;為A與B供電回路的電流,mA;AM,BM,AN及BN分別為A與M,B與M,A與N及B與N間的距離,m.

    本試驗(yàn)中測(cè)線長為35m,電極點(diǎn)間距50cm,預(yù)計(jì)探測(cè)深度約9m,測(cè)線SY1沿旋噴中心布設(shè),測(cè)線平面布置見圖4;測(cè)試序列分為3組,分別為背景值測(cè)試(T0),旋噴中射流階段測(cè)試(T1~T4),旋噴后遷移過程測(cè)試(T5~T9),監(jiān)測(cè)時(shí)長為72h,測(cè)試時(shí)間見表4.

    對(duì)旋噴中采集的數(shù)據(jù)使用Res2D軟件進(jìn)行電阻率反演,并通過Sufer軟件對(duì)反演圖像進(jìn)行進(jìn)一步圖形處理.將各時(shí)刻T電阻率與背景電阻率做差可得到電阻變化值r(式2).在3.2對(duì)旋噴過程中電阻率分析中,r可定量比較不同藥劑對(duì)土層電阻值分布的影響.

    表4 測(cè)試時(shí)間安排

    對(duì)旋噴后采集的數(shù)據(jù)通過Timelapse模塊進(jìn)行時(shí)序反演處理得到電阻變化百分比C(式3).時(shí)序反演每個(gè)監(jiān)測(cè)步驟中記錄的ERT數(shù)據(jù)獨(dú)立于任何其他信息進(jìn)行反演,在空間和時(shí)間域中引入的正則化進(jìn)一步穩(wěn)定結(jié)果,有效減少反演偽影[27],提高低濃度電阻率變化區(qū)域的靈敏度[28].

    計(jì)算公式如下:

    式中:rT時(shí)刻的電阻變化值;CT時(shí)刻的電阻變化百分比;si為T時(shí)刻的電阻率;so為背景電阻率.

    2 結(jié)果與分析

    2.1 背景電阻率分布

    圖7顯示了T0時(shí)測(cè)得的場(chǎng)地背景電阻率.0~ 3.0m毛細(xì)帶電阻率低于20W×m,對(duì)照萘的背景濃度分布,推測(cè)僅存在少量VOCs污染;深度3.0~9.0m可見均勻中電阻層,電阻率為20~35W×m,污染物主要富集于此深度的砂質(zhì)粉土及礫石層空隙中;深部8.5m粉黏土層滲透系數(shù)小,污染物滲透困難,電阻率偏低,電阻率值在10~22W×m.

    圖7 旋噴前土壤反演電阻率

    2.2 旋噴過程中電阻率分布

    測(cè)試序列T1~T4顯示了各組旋噴試驗(yàn)結(jié)束時(shí)刻的電阻率測(cè)試計(jì)算結(jié)果.圖8(a)~(d)為電阻率分布,垂直方向上可以明顯看到柱狀的低阻(小于10W×m)異常區(qū)域,圖9(c)中較寬的射流半徑可能受高滲透性地層及此前旋噴影響,總體上四組藥劑射流半徑相近,約為110cm,根據(jù)可溶性試劑旋噴試驗(yàn)[29],射流半徑約為可溶性鹽的遷移最大遷移距離,(a)~(d)中低阻區(qū)域即為PS的分布范圍.

    根據(jù)圖8(e)~(h)為電阻變化值分布,總體上四組電阻變化規(guī)律相近,旋噴注漿過程的電阻率主要受地下水稀釋藥劑影響,在藥劑驅(qū)替過程中射流前端藥劑濃度逐漸降低,徑向的電阻變化值從旋噴中心到射流前端逐漸下降;垂直方向上電阻變化最大值出現(xiàn)于深度4~5m處,這是由于旋噴過程中地層深部藥劑收到較大壓力向淺部遷移,聚集于中部砂質(zhì)粉土層區(qū)域,電阻變化值<-20W×m的區(qū)域面積NaOH組

    圖9顯示了不同區(qū)間電阻變化值在反演斷面中覆蓋的面積.A-1~A-4孔位電阻最大變化值在-22.33~-35.74W×m范圍內(nèi),平均電阻變化值分別為-8.88,-10.9,-15.4,-12.8W×m,與NaOH組相比,其余三組平均電阻變化值分別增加了22.74%,73.42%, 44.14%;對(duì)比碳載鐵,微米鐵兩組可見,在鐵粉粒徑相近的情況下,碳載鐵組平均電阻變化值比微米鐵組增加了17.4%,有效反應(yīng)半徑增加8.05%,表明相同粒徑下碳負(fù)載對(duì)零價(jià)鐵激活效果有著良好的增強(qiáng)作用.

    圖8 旋噴中土壤反演電阻率

    圖9 電阻變化值區(qū)間面積

    此外,由于電阻異常區(qū)域邊緣與旋噴中心位置藥劑濃度差異懸殊,前端低濃度區(qū)域難以達(dá)到修復(fù)目標(biāo)[30-31],因而電阻率<0W×m的區(qū)域A0在工程設(shè)計(jì)中無法作為有效修復(fù)區(qū)域.本文中以NaOH組作為對(duì)照組,取其平均電阻變化值為有效反應(yīng)閾值來判定其余各組有效反應(yīng)區(qū)域/半徑.

    根據(jù)圖9區(qū)域面積統(tǒng)計(jì),各組電阻率<-8.8W×m區(qū)域-8.8的面積分別為10.92,12.18,14.14,13.16m2,有效反應(yīng)半徑分別為78,87,101cm,94cm.與NaOH組相比,Fe0[mm]-,Fe0[nm]-,BC-nZVI分別增加了11.54%,29.49%,20.51%,零價(jià)鐵激活過硫酸鈉的有效反應(yīng)半徑隨鐵粉粒徑減小而增加.

    2.3 旋噴后電阻率分布

    圖10展示了T5~T9時(shí)的電阻率測(cè)試計(jì)算結(jié)果,相鄰測(cè)試間隔為12h.T6~T9時(shí)旋噴區(qū)域電阻率僅為10~20W×m,(a)~(e)電阻率圖像中難以準(zhǔn)確識(shí)別藥劑遷移轉(zhuǎn)化行為;為此對(duì)數(shù)據(jù)進(jìn)行時(shí)序反演處理得到電阻變化百分比分布(如圖(f)~(j)),以更加清晰展現(xiàn)旋噴區(qū)域整體電阻變化情況.

    從時(shí)序反演圖像中可知,T5~T9電阻變化百分比分布情況與2.1中射流階段相近,電阻變化百分比從旋噴中心向四周逐漸遞減,深度5~8m處為電阻變化的中心區(qū)域,變化百分比低于-30%的區(qū)域在納米鐵,零價(jià)鐵旋噴區(qū)域分布面積更大;低阻區(qū)域范圍通過擴(kuò)散作用及地下水流作用逐漸擴(kuò)大,其中變化百分比較小(-20%~0%)的區(qū)域在順?biāo)鞣较蜻w移距離更長.

    根據(jù)圖11中平均電阻變化百分比時(shí)間曲線,可將旋噴后藥劑遷移轉(zhuǎn)化過程分為3個(gè)階段:T5~T6為污染回流階段,旋噴區(qū)域中受污染的孔隙水在被氧化藥劑驅(qū)替后,受對(duì)流擴(kuò)散影響重新向旋噴區(qū)域.平均電阻變化百分比因污染濃度上升,氧化劑濃度降低,電阻變化幅度迅速下降;T6~T8為穩(wěn)定反應(yīng)階段,此階段整體反應(yīng)快速進(jìn)行,修復(fù)區(qū)域穩(wěn)步擴(kuò)張,平均電阻變化值快速下降;T8~T9為拖尾階段;平均電阻變化百分比趨近穩(wěn)定值.

    圖12顯示了電阻變化百分比各區(qū)間在反演截面中的覆蓋面積,通過面積的時(shí)間變化曲線可以進(jìn)一步分析各階段藥劑的遷移轉(zhuǎn)化過程.

    T5~T6時(shí)間段污染回流導(dǎo)致整體影響區(qū)域面積下降,其中NaOH組受影響最為嚴(yán)重;同時(shí)強(qiáng)氧化物迅速與污染物結(jié)合,地下水中導(dǎo)電粒子濃度迅速下降,導(dǎo)致低阻中心區(qū)域電阻變化百分比降低,電阻變化百分比超過30%的區(qū)域降至-20~-30%,此階段氧化反應(yīng)發(fā)生區(qū)域也主要位于該范圍內(nèi).T6時(shí)刻的修復(fù)效果 NaOH組

    圖10 旋噴后土壤反演電阻率

    T6~T8穩(wěn)定氧化階段中,氧化反應(yīng)發(fā)生區(qū)域主要位于A-2,A-3點(diǎn)位下方A-10%~-30%范圍內(nèi),污染降解引起中心反應(yīng)區(qū)域電阻逐步降低,A<-30%比T6時(shí)刻增加82.12%;整體旋噴影響區(qū)域約比T6時(shí)刻增加34.71%.值得一提的是,T7~T8時(shí)間段A-1,A-2,A-3點(diǎn)位旋噴范圍內(nèi)電阻變化接近停止,而A-4點(diǎn)位仍保持穩(wěn)定擴(kuò)張趨勢(shì).結(jié)合圖10(g)~(i)可判定, T7時(shí)刻A-1,A-2,A-3點(diǎn)位穩(wěn)定氧化階段完成,T8時(shí)刻A-4點(diǎn)位穩(wěn)定氧化階段完成,此時(shí)刻的修復(fù)效果NaOH組

    圖11 旋噴后平均電阻變化百分比

    圖12 旋噴后各變化區(qū)間面積

    T8~T9拖尾階段中,由于藥劑濃度已降至低值,濃度梯度驅(qū)動(dòng)的遷移行為逐漸減弱,受到地下水流,土層分布等因素影響,氧化藥劑分布范圍向西側(cè)緩慢擴(kuò)張,擴(kuò)散過程修復(fù)了A-3,A-4點(diǎn)位下方區(qū)域的少量殘留污染,因而圖10中A-1,A-2點(diǎn)位下方3~6m污染區(qū)域的電阻率小幅度降低;同時(shí),中心反應(yīng)區(qū)域<-30%面積僅從30.96m2增加至31.91m2,面積擴(kuò)張接近停滯;外緣0~-10%區(qū)域緩慢收縮,-10%~-20%區(qū)域緩慢擴(kuò)張,兩區(qū)域收縮擴(kuò)張面積接近互補(bǔ),整體擴(kuò)張接近停止,平均電阻變化百分比變化率接近于零,表明A-1~A-4修復(fù)區(qū)域氧化反應(yīng)均接近完成,整體平均電阻變化百分比穩(wěn)定至-23%左右,此時(shí)刻的修復(fù)效果NaOH組≈Fe0[μm]組

    2.4 不同激活劑修復(fù)效果討論

    圖13顯示了修復(fù)后48h萘的實(shí)測(cè)濃度分布,Fe0[nm]組,BC-nZVI組修復(fù)效果明顯優(yōu)于其他兩組,沿深度方向修復(fù)較為均勻,深度4m下萘濃度均降至100mg/kg以下,圖13(b)萘修復(fù)值分布中顯示出相同F(xiàn)e0摻量的碳載鐵組污染去除量更大;NaOH組在污染回流階段受對(duì)流影響嚴(yán)重,更易發(fā)生水平向的遷移擴(kuò)散,整體修復(fù)效果較差;Fe0[μm]組僅在深度4.5~8.0m深度展示出較好的修復(fù)效果.

    圖14 萘-電阻率相關(guān)性

    由圖14可見,電阻變化值與萘修復(fù)量呈現(xiàn)較強(qiáng)的正相關(guān)關(guān)系,擬合曲線斜率約為(-53.6±6.0),2為0.97.對(duì)比圖14(b) (c)萘實(shí)測(cè)值的等值線分布,電阻變化值等值線基本準(zhǔn)確的反映了污染修復(fù)發(fā)生的位置,即便在污染較輕的區(qū)域(此區(qū)域ERT監(jiān)測(cè)相對(duì)不靈敏[34-35]),其位置誤差也未超過0.5m;不可忽視的是,通過擬合曲線預(yù)估的修復(fù)值與實(shí)際修復(fù)值間存在誤差,對(duì)比圖15中萘實(shí)測(cè)修復(fù)值與對(duì)應(yīng)的估計(jì)值,整體誤差為10.3%~22.6%,引起誤差的原因是ERT監(jiān)測(cè)對(duì)局部電阻過高/過低區(qū)域監(jiān)測(cè)較為不靈敏.相似調(diào)查[36]中發(fā)現(xiàn)在DNAPLs污染嚴(yán)重的區(qū)域誤差甚至可以達(dá)到50%.

    3 結(jié)論

    3.1 ERT監(jiān)測(cè)結(jié)果基本準(zhǔn)確的反映了污染修復(fù)發(fā)生的位置,最大偏差小于0.5m;電阻變化值與萘修復(fù)值呈現(xiàn)正相關(guān)關(guān)系,通過擬合曲線預(yù)估的修復(fù)值與實(shí)際修復(fù)值間存在誤差,誤差為10.3%~22.6%.

    3.2 根據(jù)電阻率變化可將旋噴修復(fù)分為4個(gè)階段:射流階段,污染回流階段,穩(wěn)定氧化階段,拖尾階段.射流階段,地層出現(xiàn)明顯的柱狀低阻區(qū)域,藥劑在深度3~6m區(qū)間內(nèi)富集;徑向電阻變化值沿射流方向迅速減小,有效影響半徑NaOH組

    3.3 本試驗(yàn)在保持較高納米鐵粉添加量的情況下,納米鐵組、碳載鐵組的土壤修復(fù)效果較為理想;碳載鐵的緩釋性能使Fe2+的穩(wěn)定釋放,在保持較高的反應(yīng)速率的同時(shí)提高了反應(yīng)的穩(wěn)定性,延長穩(wěn)定反應(yīng)時(shí)間約8h.旋噴后48h,碳載鐵組預(yù)設(shè)旋噴區(qū)域平均萘修復(fù)值為358.6mg/kg,約為納米鐵修復(fù)量的1.25倍,各組污染修復(fù)量BC-nZVI組>Fe0[nm]組>Fe0[μm]組≈NaOH組.

    [1] 杜延軍,金 飛,劉松玉,等.重金屬工業(yè)污染場(chǎng)地固化/穩(wěn)定處理研究進(jìn)展 [J]. 巖土力學(xué), 2011,32(1):116-124. Du Y J, Jin F, Liu S Y, et al. Review of stabilization/solidification technique for remediation of heavy metals contaminated lands [J]. Rock and Soil Mechanics, 2011,32(1):116-124.

    [2] 王 穎.土壤氣相抽提修復(fù)苯系物污染場(chǎng)地模擬及效果預(yù)測(cè)研究 [D]. 保定:河北農(nóng)業(yè)大學(xué), 2020. Wang Y. Study on simulation and prediction of BTET contaminated site remediation effect by soil vapor extraction [D]. Baoding: Hebei Agricultural University, 2020.

    [3] Zhang W, Mi Y, Jiao W. Study on the migration mechanisms of water-soluble agents in high-pressure rotary jetting remediation [J]. Environ Sci Pollut Res Int, 2022,29:74038-74050.

    [4] Ma B, Wang Z, Yuan X, et al. In situ stabilization of heavy metals in a tailing pond with a new method for the addition of mineral stabilizers—high-pressure rotary jet technology [J]. Environmental Science and Pollution Research, 2020,27:15388-15400.

    [5] Ranc B, Faure P, Croze V, et al. Selection of oxidant doses for in situ, chemical oxidation of soils contaminated by polycyclic aromatic hydrocarbons (PAHs): A review [J]. Journal of Hazardous Materials, 2016,312:280-297.

    [6] Brusseau M L, Nelson N T, Zhang Z, et al. Source-zone characteri- zation of a chlorinated-solvent contaminated Superfund site in Tucson, AZ [J]. Journal of Contaminant Hydrology, 2015,90(1/2):21-40.

    [7] Xia W Y, Du Y J, Li F S, et al. In-situ solidification/stabilization of heavy metals contaminated site soil using a dry jet mixing method and new hydroxyapatite based binder [J]. Journal of Hazardous Materials, 2019,369:353-361.

    [8] 宛 召.高壓旋噴工藝在上海某污染場(chǎng)地修復(fù)中的應(yīng)用研究 [D]. 長春:吉林大學(xué), 2017. Wan Z, Study on the application of high pressure jet grouting technology in the remediation of a contaminated Site in Shanghai [D]. Changchun: Jilin University, 2017.

    [9] 董萬濤.新型復(fù)合氧化劑在石油污染場(chǎng)地修復(fù)的應(yīng)用研究 [D]. 蘭州:蘭州理工大學(xué), 2021. Dong Wan-tao. Research on application of new complex oxidants in petroleum contaminated site [D]. Lanzhou: Lanzhou University of Technology, 2021.

    [10] Site Characterization technologies for DNAPL Investigations [Z]. EPA 542-R-04-017, 2004.

    [11] Lane J W , Day L F D, Casey C C. Geophysical monitoring of a field-scale bio-stimulation pilot project [J]. Ground Water, 2010, 44(3).

    [12] Peale J G D, Mueller J, Molin J.Successful ISCR enhanced bioremediation of a TCE DNAPL source utilizing EHC and KB-1 [J]. Remediation Journal, 2010,20(3):63-81.

    [13] GB/T 50123-2019 土工試驗(yàn)方法標(biāo)準(zhǔn) [S]. GB/T 50123-2019 Standard for geotechnical testing method [S].

    [14] GB/T 50145-2007 土的工程分類標(biāo)準(zhǔn) [S]. GB/T 50145-2007 Standard for classification of engineering soils [S].

    [15] HJ/T 166-2004 土壤環(huán)境監(jiān)測(cè)技術(shù)規(guī)范[S]. HJ/T 166-2004 The technical specification for soil environmental monitoring [S].

    [16] HJ 1019-2019 地塊土壤和地下水中揮發(fā)性有機(jī)物采樣技術(shù)導(dǎo)則: [S]. HJ 1019-2019 Technical guideline for site soil and groundwater sampling of volatile organic compounds [S].

    [17] HJ 605-2011 土壤和沉積物揮發(fā)性有機(jī)物的測(cè)定吹掃捕集/氣象色譜-質(zhì)譜法[S] HJ 605-2011 Soil and sediment-determination of volatile organic compounds -purge and trap gas chromatography/mass spectrometry method [S].

    [18] HJ 766-2015 固體廢物金屬元素的測(cè)定電感耦合等離子體質(zhì)譜法[S]. HJ 766-2015 Solid waste - determination of metals inductively coupled plasma mass spectrometry (ICP-MS) [S].

    [19] HJ 804-2016 土壤8種有效態(tài)元素的測(cè)定二乙烯三胺五乙酸浸提-電感耦合等離子體發(fā)射光譜法[S]. HJ 804-2016 Soil - determination of bioavailable form of eight elements - Extraction withbuffered DTPA solution/Inductively coupled plasma optical emission spectrometry [S].

    [20] GB36600-2018 土壤環(huán)境質(zhì)量建設(shè)用地土壤污染風(fēng)險(xiǎn)管控標(biāo)準(zhǔn)(試行) [S]. GB36600-2018 Soil environmental quality risk control standard for soil contamination of development land [S].

    [21] 潘棟宇,侯梅芳,劉超男,等.多環(huán)芳烴污染土壤化學(xué)修復(fù)技術(shù)的研究進(jìn)展 [J]. 安全與環(huán)境工程, 2018,25(3):54-66. Pan D Y, Hou M F, Liu C N, et al. Review of chemical remediation technology of polycyclic aromatic hydrocarbons contaminated soil. [J]. Safety and Environmental Engineering, 2018, 25(3):54-66.

    [22] Killian P F, Bruell C J, Marley M C, et al. Iron (II) activated persulfate oxidation of MGP contaminated soil [J]. Journal of Soil Contamination, 2007,16(6):523-537.

    [23] Do S H, Kwon Y J, Kong S H. Effect of metal oxides on the reactivity of persulfate/Fe(II) in the remediation of diesel-contaminated soil and sand [J]. Journal of Hazardous Materials, 2010,182(1):933-936.

    [24] Pac T, Baldock J, Brodie B, et al. In situ chemical oxidation: Lessons learned at multiple sites [J]. Remediation Journal, 2019,29:75-91.

    [25] 張 羽,高春陽,陳昌照,等.零價(jià)鐵活化過硫酸鈉體系降解污染土壤中的多環(huán)芳烴 [J]. 環(huán)境工程學(xué)報(bào), 2019,13(4):955-962. Zhang Y, Gao C Y, Chen C Z, et al. Degradation of PAHs in contaminated soil by zero valent iron activatedsodium persulfate system [J]. Chinese Journal of Environmental Engineering, 2019,13(4): 955-962.

    [26] Zhang B W, Guo Y, Huo J Y. Combining chemical oxidation and bioremediation for petroleum polluted soil remediation by BC-nZVI activated persulfate [J]. Chemical Engineering Journal, 2021,382(1): 1230 55,2-9.

    [27] Kim J H, Yi M J, Park S G, et al. 4-D inversion of DC resistivity monitoring data acquired over a dynamically changing earth model [J]. Journal of Applied Geophysics, 2009,68(4):522-532.

    [28] Karaoulis M C, Kim J H, Tsourlos P I. 4D active time constrained resistivity inversion [J]. Journal of Applied Geophysics, 2011,73(1): 25-34.

    [29] 張文杰,秘永寶.高壓旋噴修復(fù)中水溶性藥劑徑向遷移規(guī)律研究 [J]. 巖土工程學(xué)報(bào), 2023,45(5):1017-1023. Zhang Wen-jie, Mi Yong-bao. Radial migration of water-soluble agents in high-pressure rotary jetting remediation [J]. Chinese Journal of Geotechnical Engineering, 2023,45(5):1017-1023.

    [30] 韓瑞瑞.碳載雙金屬活化過硫酸鈉修復(fù)石油污染土壤 [D]. 北京:北京化工大學(xué), 2022.DOI:10.26939/d.cnki.gbhgu.20 22.000769. Han Rui-rui. Remediation of petroleum-contaminated soil by sodium persulfate activated by carbon-suppprted bimetallic [D]. Beijing University of Chemical Technology, 2022.DOI:10.26939/d.cnki. gbhgu, 2022.000769.

    [31] Thao T T, Kim C, Hwang I. Application of nanosized zerovalent iron activated persulfate for treating groundwater contaminated with phenol. [J] Journal of Soil and Groundwater Environment, 2017,22(1):41–48.

    [32] Przemyslaw D, Leonidas P, Allaa, et al. Impact of peroxydisulfate in the presence of zero valent iron on the oxidation of cyclohexanoic acid and naphthenic acids from oil sands process-affected water [J]. Environmental Science & Technology, 2012,46(16):89-91.

    [33] Kim C, Ahn J Y, Kim T Y, et al. Activation of persulfate by nanosized zero-valent iron (NZVI): Mechanisms and transformation products of NZVI [J]. Environmental Science & Technology, 2018,52(6):3625– 3633.

    [34] Furukawa Y, KIM J W, Watkins J, et al. Formation of ferrihydrite and associated iron corrosion products in permeable reactive barriers of zero-valent iron [J]. Environ -mental Science & Technology, 2002, 36(24):5469-5475.

    [35] Goi A, Trapido M, Kulik N, et al. Ozonation and Fenton treatment for remediation of diesel fuel contaminated soil [J]. Ozone Science & Engineering, 2006,28(1):37-46.

    [36] Christopher Power, Jason I. Gerhard, MariosKaraoulis, Panagiotis Tsourlos, Antonios Giannopoulos, Evaluating four-dimensional time- lapse electrical resistivity tomogram -phy for monitoring DNAPL source zone remediation [J]. Journal of Contaminant Hydrology, 2014, s162–163:27-46.

    Field study on ERT monitoring high-pressure rotary jet repairing PAHs pollution.

    LIU Di1, SONG Quan-wei2,3, LI Ling4, HU Jie1, WANG Ting-jun4, DU Xian-yuan2,3, DENG Yong-feng5, KE Han1

    (1.Key Laboratory of Soft Soils and Geoenvironmental Engineering, Ministry of Education, Institute of Geotechnical Engineering, Zhejiang University, Hangzhou 310058, China;2.State Key Laboratory of Petroleum Pollution Control, Beijing 102206, China;3.CNPC Research Institute of Safety and Environment Technology, Beijing102206, China;4.CECEP DADI Environmental Remediation Co., Ltd., Beijing 100085, China;5.Institute of Geotechnical Engineering, Transportation College, Southeast University, Nanjing 210096, China)., 2023,43(11):5933~5943

    In this study, sodium persulfate (PS) was selected as the oxidant, Fe0[μm], Fe0[nm, BC-nZVI and NaOH were chosen as activators, and repair PAHs-contaminated soil through high-pressure rotary spraying process by mixture of PS and activators. The high-density electrical resistivity tomography (ERT) technique was employed to monitor the distribution of soil resistivity during the rotary jet process, enabling an analysis of the migration and transformation patterns of the oxidation agents during the remediation process. Additionally, the remediation effects of the 4 agent groups were compared after a 48-hour jet period.The experimental results demonstrated that ERT accurately reflected the locations of pollution remediation, with a maximum deviation of less than 0.5m. The changes in resistivity showed a positive correlation with naphthalene remediation. However, discrepancies existed between the estimated remediation values obtained from fitting curves and the actual remediation values, with errors ranging from 10.3% to 22.6%. Among the tested reagents, carbon-supported iron powder exhibited excellent stability and sustained release performance, leading to an extension of the stable oxidation stage by approximately 8hours after jet. Ultimately, the predetermined remediation area achieved the highest level of pollutant removal in the following order: BC-nZVI>Fe0[nm]>Fe0[μm]≈NaOH.

    multi-electrode resistivity method (ERT);high pressure rotary jet;in-situ oxidation;PAHs;field test

    X53

    A

    1000-6923(2023)11-5933-11

    柳 迪(1997-),男,山東濱州人,浙江大學(xué)碩士研究生,主要從事環(huán)境土工研究.d_liu@zju.edu.cn.

    柳 迪,宋權(quán)威,李 玲,等.ERT法監(jiān)測(cè)高壓旋噴修復(fù)PAHs污染現(xiàn)場(chǎng)試驗(yàn)研究 [J]. 中國環(huán)境科學(xué), 2023,43(11):5933-5943.

    Liu D, Song Q W, Li L, et al. Field study on ERT monitoring high-pressure rotary jet repairing PAHs pollution [J]. China Environmental Science, 2023,43(11):5933-5943.

    2023-05-23

    國家重點(diǎn)研發(fā)計(jì)劃(2019YFC1806000,2019YFC1806004);中國石油集團(tuán)科技創(chuàng)新基金資助項(xiàng)目(2022DJ6906);浙江省領(lǐng)雁研發(fā)攻關(guān)計(jì)劃(2022C03095)

    * 責(zé)任作者, 教授, boske@126.com

    猜你喜歡
    電阻率藥劑電阻
    小麥拌種防效好 藥劑選擇要當(dāng)心
    電阻大小誰做主
    巧測(cè)電阻
    選用對(duì)口藥劑 適期防治病蟲
    油田污水處理藥劑的選用
    AF-701藥劑在黃瓜上的應(yīng)用試驗(yàn)
    電阻焊
    三維電阻率成像與高聚物注漿在水閘加固中的應(yīng)用
    基于multisim13.0負(fù)電阻在電源中作用的仿真
    隨鉆電阻率測(cè)井的固定探測(cè)深度合成方法
    亚洲国产欧美人成| 韩国av一区二区三区四区| 免费一级毛片在线播放高清视频| 三级男女做爰猛烈吃奶摸视频| 色在线成人网| 中文字幕熟女人妻在线| 热99re8久久精品国产| 国产又黄又爽又无遮挡在线| 三级男女做爰猛烈吃奶摸视频| 3wmmmm亚洲av在线观看| 我要搜黄色片| aaaaa片日本免费| 国产在线精品亚洲第一网站| 九九在线视频观看精品| 免费av毛片视频| 亚洲av熟女| 非洲黑人性xxxx精品又粗又长| 亚洲在线观看片| av福利片在线观看| 国产蜜桃级精品一区二区三区| 欧美成人性av电影在线观看| 欧美3d第一页| 国产成人一区二区在线| 国产男靠女视频免费网站| 国产精品日韩av在线免费观看| 久久午夜亚洲精品久久| 精品一区二区三区av网在线观看| 亚洲第一区二区三区不卡| 久久久成人免费电影| 欧美在线一区亚洲| 亚洲久久久久久中文字幕| 一个人看视频在线观看www免费| 此物有八面人人有两片| 亚洲狠狠婷婷综合久久图片| 免费在线观看成人毛片| 在线国产一区二区在线| 国产乱人视频| a级毛片免费高清观看在线播放| 久久午夜亚洲精品久久| 大又大粗又爽又黄少妇毛片口| 亚洲电影在线观看av| 91麻豆精品激情在线观看国产| 男女视频在线观看网站免费| 尤物成人国产欧美一区二区三区| 欧美日韩综合久久久久久 | 欧美性猛交╳xxx乱大交人| 亚洲成人久久爱视频| 精品国内亚洲2022精品成人| 午夜福利高清视频| 国产精品爽爽va在线观看网站| 国产美女午夜福利| 赤兔流量卡办理| 中文字幕av成人在线电影| 国产av不卡久久| 男人和女人高潮做爰伦理| 国产av一区在线观看免费| 网址你懂的国产日韩在线| 国产高清激情床上av| 天天躁日日操中文字幕| 久久久久精品国产欧美久久久| 国产黄片美女视频| 一区福利在线观看| 色精品久久人妻99蜜桃| 亚洲第一电影网av| 熟女电影av网| 天堂动漫精品| 日韩 亚洲 欧美在线| 美女被艹到高潮喷水动态| 国产高清视频在线播放一区| 国产亚洲91精品色在线| 99精品久久久久人妻精品| 午夜福利18| 少妇的逼好多水| 男人和女人高潮做爰伦理| 啦啦啦观看免费观看视频高清| 女同久久另类99精品国产91| 国产成年人精品一区二区| 一本一本综合久久| 亚洲三级黄色毛片| 又粗又爽又猛毛片免费看| 成人特级av手机在线观看| 久久久久性生活片| 久久精品国产亚洲av涩爱 | 一进一出抽搐gif免费好疼| 亚洲熟妇中文字幕五十中出| 国产精品一区www在线观看 | 最近最新免费中文字幕在线| 国产综合懂色| 两个人的视频大全免费| 淫秽高清视频在线观看| 久久久国产成人精品二区| 国产精品乱码一区二三区的特点| 日本 av在线| 一区二区三区免费毛片| 亚洲av成人精品一区久久| 999久久久精品免费观看国产| 久久人人爽人人爽人人片va| 很黄的视频免费| 亚洲无线观看免费| 一a级毛片在线观看| 毛片一级片免费看久久久久 | 亚洲国产精品sss在线观看| 国产黄色小视频在线观看| 亚洲av不卡在线观看| 国产精品一区二区三区四区免费观看 | 精品久久久噜噜| 国产高潮美女av| 久久精品综合一区二区三区| 很黄的视频免费| 国产精品99久久久久久久久| 久久精品国产亚洲av香蕉五月| 一本久久中文字幕| 亚洲在线观看片| 亚洲 国产 在线| 最近视频中文字幕2019在线8| 麻豆av噜噜一区二区三区| 不卡视频在线观看欧美| 日韩欧美一区二区三区在线观看| 亚洲av免费高清在线观看| 国产欧美日韩精品一区二区| 日韩欧美精品v在线| 91久久精品国产一区二区成人| 日本在线视频免费播放| 一区二区三区高清视频在线| 在线免费十八禁| 在线播放国产精品三级| 一区二区三区高清视频在线| 少妇裸体淫交视频免费看高清| 欧美成人免费av一区二区三区| 一级黄色大片毛片| 久久人人爽人人爽人人片va| 美女高潮的动态| 老熟妇乱子伦视频在线观看| 天堂动漫精品| 免费看av在线观看网站| 三级男女做爰猛烈吃奶摸视频| 国产男人的电影天堂91| xxxwww97欧美| 一个人看的www免费观看视频| 国产精品一及| 欧美xxxx性猛交bbbb| 搡老熟女国产l中国老女人| 日韩中字成人| 久久久色成人| 91狼人影院| 色噜噜av男人的天堂激情| 桃红色精品国产亚洲av| 中文资源天堂在线| 国产高清激情床上av| 久久午夜亚洲精品久久| 久久久久久久久久成人| 日日撸夜夜添| 热99re8久久精品国产| 国产精品乱码一区二三区的特点| 88av欧美| 久久99热6这里只有精品| 黄色配什么色好看| 久久精品影院6| 夜夜夜夜夜久久久久| 天天一区二区日本电影三级| 亚州av有码| 国产美女午夜福利| 看十八女毛片水多多多| 国产成人福利小说| 如何舔出高潮| 亚洲成人久久性| 免费看美女性在线毛片视频| 国产在视频线在精品| 12—13女人毛片做爰片一| 国产精品精品国产色婷婷| 蜜桃久久精品国产亚洲av| 国产极品精品免费视频能看的| a级毛片免费高清观看在线播放| 婷婷丁香在线五月| 国产精品一区二区性色av| 黄色视频,在线免费观看| 久久精品国产亚洲av天美| 我的老师免费观看完整版| 亚洲性夜色夜夜综合| 九九爱精品视频在线观看| 国内精品美女久久久久久| 我要看日韩黄色一级片| 久久人妻av系列| 欧美+亚洲+日韩+国产| 狠狠狠狠99中文字幕| 最近最新中文字幕大全电影3| 神马国产精品三级电影在线观看| 日本 欧美在线| av专区在线播放| 亚洲精品成人久久久久久| 免费观看的影片在线观看| 亚洲精华国产精华液的使用体验 | 欧美一区二区精品小视频在线| 欧美日韩乱码在线| 精品久久久久久久久久免费视频| 亚洲成人精品中文字幕电影| 亚洲av免费高清在线观看| 欧美+日韩+精品| 久久中文看片网| 很黄的视频免费| 在线a可以看的网站| 变态另类丝袜制服| 精品久久久久久成人av| 亚洲av中文av极速乱 | av在线蜜桃| 欧美日韩乱码在线| 免费人成在线观看视频色| 久久国产乱子免费精品| 99在线视频只有这里精品首页| 久久精品国产亚洲av天美| 国产精品综合久久久久久久免费| 久久亚洲真实| 日日摸夜夜添夜夜添av毛片 | 国产av麻豆久久久久久久| 亚洲精品成人久久久久久| 1024手机看黄色片| 色综合色国产| 国产精品爽爽va在线观看网站| 尤物成人国产欧美一区二区三区| 亚洲自拍偷在线| 久久国内精品自在自线图片| 成人国产综合亚洲| 丰满乱子伦码专区| 老女人水多毛片| 午夜福利18| 国产精品嫩草影院av在线观看 | 九九热线精品视视频播放| 亚洲avbb在线观看| 国产av在哪里看| 亚洲欧美清纯卡通| 成人国产综合亚洲| 亚洲性夜色夜夜综合| 久久国产精品人妻蜜桃| 色在线成人网| 18禁黄网站禁片午夜丰满| 免费无遮挡裸体视频| 国产成人福利小说| 国产aⅴ精品一区二区三区波| 亚洲真实伦在线观看| 高清毛片免费观看视频网站| 亚洲va日本ⅴa欧美va伊人久久| 国产av不卡久久| 亚洲美女搞黄在线观看 | 亚洲无线在线观看| 老师上课跳d突然被开到最大视频| 日韩欧美 国产精品| 99久久无色码亚洲精品果冻| 国产精品亚洲美女久久久| 亚洲自偷自拍三级| 欧美国产日韩亚洲一区| 午夜福利18| 亚洲国产色片| 熟女电影av网| 男女视频在线观看网站免费| 国产免费男女视频| 久久久久久久午夜电影| 久久久精品大字幕| 久久精品久久久久久噜噜老黄 | www日本黄色视频网| 国产人妻一区二区三区在| 国产乱人伦免费视频| 国产精品自产拍在线观看55亚洲| 可以在线观看毛片的网站| 男女做爰动态图高潮gif福利片| xxxwww97欧美| 国产乱人视频| 1024手机看黄色片| 国产精品久久久久久久电影| 丝袜美腿在线中文| 伦精品一区二区三区| 久久人人爽人人爽人人片va| 国产 一区 欧美 日韩| 精品久久久久久久久亚洲 | 久久精品夜夜夜夜夜久久蜜豆| 美女 人体艺术 gogo| av.在线天堂| 91在线精品国自产拍蜜月| 国产激情偷乱视频一区二区| 日韩欧美三级三区| 又爽又黄a免费视频| 91av网一区二区| 热99re8久久精品国产| 亚洲电影在线观看av| 69av精品久久久久久| 亚洲国产日韩欧美精品在线观看| 床上黄色一级片| 日本在线视频免费播放| 久久精品国产自在天天线| 亚洲人成伊人成综合网2020| ponron亚洲| www.www免费av| 国产精品乱码一区二三区的特点| 欧美一区二区亚洲| 午夜免费男女啪啪视频观看 | 亚洲四区av| .国产精品久久| 69人妻影院| 日韩欧美在线乱码| 久久精品综合一区二区三区| 成人一区二区视频在线观看| 97超视频在线观看视频| 又爽又黄a免费视频| 黄片wwwwww| 成年免费大片在线观看| 午夜免费成人在线视频| 51国产日韩欧美| av在线蜜桃| 俺也久久电影网| 黄色一级大片看看| 日本免费a在线| 欧美高清成人免费视频www| 久久精品影院6| 观看免费一级毛片| 成人国产麻豆网| 日本欧美国产在线视频| 最近最新中文字幕大全电影3| 桃红色精品国产亚洲av| 直男gayav资源| 夜夜看夜夜爽夜夜摸| 婷婷六月久久综合丁香| 日本 欧美在线| 国产aⅴ精品一区二区三区波| 22中文网久久字幕| 日韩欧美免费精品| 精品久久久久久成人av| 999久久久精品免费观看国产| 我要看日韩黄色一级片| 久久精品国产亚洲av天美| 一进一出抽搐gif免费好疼| 国产高清有码在线观看视频| 真实男女啪啪啪动态图| 婷婷精品国产亚洲av在线| 99久久无色码亚洲精品果冻| 免费观看精品视频网站| 少妇人妻一区二区三区视频| 午夜福利视频1000在线观看| 国产亚洲精品av在线| 国产精品久久久久久精品电影| 三级毛片av免费| 亚洲经典国产精华液单| 成年女人永久免费观看视频| 国内久久婷婷六月综合欲色啪| 亚洲精品成人久久久久久| 啦啦啦啦在线视频资源| 成人三级黄色视频| 91精品国产九色| 亚洲成人精品中文字幕电影| 极品教师在线视频| avwww免费| 色av中文字幕| 日本 av在线| 亚洲久久久久久中文字幕| 3wmmmm亚洲av在线观看| 欧洲精品卡2卡3卡4卡5卡区| 久久精品夜夜夜夜夜久久蜜豆| 午夜福利成人在线免费观看| 国产精品av视频在线免费观看| 亚洲男人的天堂狠狠| 精华霜和精华液先用哪个| 国产69精品久久久久777片| 午夜福利在线在线| 国产成人av教育| 国产乱人视频| 国产视频内射| 97人妻精品一区二区三区麻豆| 桃色一区二区三区在线观看| 村上凉子中文字幕在线| 99热这里只有是精品在线观看| av天堂在线播放| 美女 人体艺术 gogo| 欧美+亚洲+日韩+国产| 热99re8久久精品国产| 又粗又爽又猛毛片免费看| 国产美女午夜福利| 成年人黄色毛片网站| 日韩中文字幕欧美一区二区| 真人一进一出gif抽搐免费| 亚洲第一区二区三区不卡| 久久国内精品自在自线图片| 国产精品综合久久久久久久免费| 在线免费观看不下载黄p国产 | 国产人妻一区二区三区在| 一区二区三区高清视频在线| 午夜激情福利司机影院| 91午夜精品亚洲一区二区三区 | 少妇丰满av| 中文字幕免费在线视频6| 欧美日韩精品成人综合77777| 老司机福利观看| 国内揄拍国产精品人妻在线| 中文在线观看免费www的网站| 国产免费av片在线观看野外av| 亚洲第一电影网av| 国产成年人精品一区二区| 成人综合一区亚洲| 露出奶头的视频| netflix在线观看网站| 乱人视频在线观看| 黄色配什么色好看| 久久久久国内视频| 免费在线观看影片大全网站| 热99re8久久精品国产| 国产av在哪里看| 性欧美人与动物交配| 美女大奶头视频| 91麻豆av在线| 色综合亚洲欧美另类图片| 观看免费一级毛片| 精品久久国产蜜桃| 尾随美女入室| 99国产精品一区二区蜜桃av| 欧美日韩综合久久久久久 | 精品福利观看| 午夜免费激情av| 身体一侧抽搐| 少妇人妻精品综合一区二区 | 国产精品99久久久久久久久| 99久久无色码亚洲精品果冻| av.在线天堂| 久久久国产成人免费| 97碰自拍视频| 99久国产av精品| 久久亚洲精品不卡| 国产成人a区在线观看| 久久人人爽人人爽人人片va| 在线播放无遮挡| 最近中文字幕高清免费大全6 | 观看美女的网站| 欧美不卡视频在线免费观看| 亚洲va在线va天堂va国产| 国产午夜福利久久久久久| 日韩亚洲欧美综合| 日韩欧美国产一区二区入口| 精品午夜福利视频在线观看一区| 99热网站在线观看| 国产探花在线观看一区二区| 99九九线精品视频在线观看视频| 成年版毛片免费区| 午夜精品在线福利| 日本五十路高清| 婷婷丁香在线五月| 亚洲一级一片aⅴ在线观看| 久久久精品大字幕| 听说在线观看完整版免费高清| 国产精品99久久久久久久久| 一区二区三区免费毛片| 日韩欧美 国产精品| 日本撒尿小便嘘嘘汇集6| 亚洲欧美日韩无卡精品| 尤物成人国产欧美一区二区三区| 亚洲精品粉嫩美女一区| 国产免费男女视频| 俺也久久电影网| 国产爱豆传媒在线观看| 久久草成人影院| 精品一区二区免费观看| 免费观看精品视频网站| 我要搜黄色片| 国产老妇女一区| av女优亚洲男人天堂| 亚洲自拍偷在线| 久久九九热精品免费| 变态另类成人亚洲欧美熟女| 一本一本综合久久| 99国产极品粉嫩在线观看| 99精品久久久久人妻精品| 日本-黄色视频高清免费观看| 简卡轻食公司| 97人妻精品一区二区三区麻豆| 久久久久精品国产欧美久久久| 国产av不卡久久| 亚洲无线观看免费| 国模一区二区三区四区视频| 22中文网久久字幕| 婷婷丁香在线五月| 美女高潮喷水抽搐中文字幕| 午夜日韩欧美国产| 99久久无色码亚洲精品果冻| 国产女主播在线喷水免费视频网站 | 国产美女午夜福利| 精品久久久久久久久亚洲 | 我要看日韩黄色一级片| h日本视频在线播放| 欧美色视频一区免费| 亚洲乱码一区二区免费版| 亚洲久久久久久中文字幕| 国产高清视频在线观看网站| 欧美日韩国产亚洲二区| 18禁裸乳无遮挡免费网站照片| ponron亚洲| 91精品国产九色| 国产精品国产三级国产av玫瑰| 中文字幕高清在线视频| 久久亚洲真实| 日本免费一区二区三区高清不卡| 亚洲精品456在线播放app | 亚洲中文日韩欧美视频| 少妇猛男粗大的猛烈进出视频 | 国产精品98久久久久久宅男小说| 成人av在线播放网站| 亚洲自拍偷在线| 午夜久久久久精精品| 久久精品综合一区二区三区| 老熟妇乱子伦视频在线观看| 午夜爱爱视频在线播放| 国产不卡一卡二| 日韩大尺度精品在线看网址| 免费看美女性在线毛片视频| 在线天堂最新版资源| 少妇人妻精品综合一区二区 | av福利片在线观看| 51国产日韩欧美| 男女下面进入的视频免费午夜| 99热这里只有是精品50| 美女xxoo啪啪120秒动态图| 日本成人三级电影网站| 免费在线观看成人毛片| 久久久久久久久大av| 日韩欧美在线乱码| 一边摸一边抽搐一进一小说| 人人妻,人人澡人人爽秒播| 日韩欧美在线二视频| 国产成人福利小说| 婷婷亚洲欧美| 国产成人aa在线观看| 日韩中字成人| 欧美潮喷喷水| 日韩av在线大香蕉| 在线观看av片永久免费下载| 亚洲精品乱码久久久v下载方式| 国产成人福利小说| 看十八女毛片水多多多| 日韩 亚洲 欧美在线| 国产主播在线观看一区二区| 99久久精品一区二区三区| 岛国在线免费视频观看| 国产精品一区二区性色av| 极品教师在线视频| 变态另类丝袜制服| 免费无遮挡裸体视频| 18禁黄网站禁片午夜丰满| 国产精品av视频在线免费观看| a级一级毛片免费在线观看| 国产精品久久久久久av不卡| 亚洲精品亚洲一区二区| 直男gayav资源| 狂野欧美激情性xxxx在线观看| 三级男女做爰猛烈吃奶摸视频| 韩国av在线不卡| 99久久中文字幕三级久久日本| 国产亚洲精品av在线| 午夜免费激情av| 国产一区二区三区av在线 | 久久久久久久久久久丰满 | 99精品在免费线老司机午夜| 欧美在线一区亚洲| 中文字幕av成人在线电影| 国产精品永久免费网站| 国产精品爽爽va在线观看网站| 女同久久另类99精品国产91| 亚洲一区高清亚洲精品| 欧美精品国产亚洲| 人妻制服诱惑在线中文字幕| 麻豆成人午夜福利视频| 久久这里只有精品中国| netflix在线观看网站| 午夜福利视频1000在线观看| 97人妻精品一区二区三区麻豆| 大又大粗又爽又黄少妇毛片口| 丝袜美腿在线中文| 日韩精品青青久久久久久| 又爽又黄a免费视频| 超碰av人人做人人爽久久| 免费观看在线日韩| 午夜福利成人在线免费观看| 亚洲av日韩精品久久久久久密| 麻豆成人av在线观看| 干丝袜人妻中文字幕| 免费高清视频大片| 少妇高潮的动态图| 国内精品久久久久久久电影| 美女黄网站色视频| 亚洲图色成人| 极品教师在线免费播放| av在线蜜桃| 欧洲精品卡2卡3卡4卡5卡区| 最近视频中文字幕2019在线8| 亚洲第一区二区三区不卡| 色av中文字幕| 亚洲最大成人av| 亚洲电影在线观看av| 欧美成人免费av一区二区三区| 色综合站精品国产| 亚洲中文字幕日韩| 成熟少妇高潮喷水视频| a级毛片免费高清观看在线播放| 国产亚洲精品久久久久久毛片| 国产精品久久久久久久久免| 免费不卡的大黄色大毛片视频在线观看 | 欧美精品啪啪一区二区三区| 一进一出抽搐动态| 久久久久久伊人网av| 亚洲va日本ⅴa欧美va伊人久久| 国产美女午夜福利| 亚洲一级一片aⅴ在线观看| 欧美精品啪啪一区二区三区| 午夜激情欧美在线| 美女高潮喷水抽搐中文字幕| 亚洲色图av天堂| 麻豆av噜噜一区二区三区| 九色成人免费人妻av| av福利片在线观看| 成年女人毛片免费观看观看9| 久久精品综合一区二区三区| 亚洲电影在线观看av|