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

    鹽度對多孔介質(zhì)中DNAPL運移和分布的影響

    2017-08-16 08:43:37徐紅霞孫媛媛施小清廖朋輝吳吉春
    水文地質(zhì)工程地質(zhì) 2017年4期
    關(guān)鍵詞:砂箱均質(zhì)運移

    程 洲,徐紅霞,孫媛媛,施小清,廖朋輝,吳吉春

    (南京大學(xué)地球科學(xué)與工程學(xué)院水科學(xué)系/污染控制與資源化國家重點實驗室,江蘇 南京 210023)

    ?

    鹽度對多孔介質(zhì)中DNAPL運移和分布的影響

    程 洲,徐紅霞,孫媛媛,施小清,廖朋輝,吳吉春

    (南京大學(xué)地球科學(xué)與工程學(xué)院水科學(xué)系/污染控制與資源化國家重點實驗室,江蘇 南京 210023)

    選用四氯乙烯(PCE)作為典型DNAPL污染物,以NaCl作為地下水中溶解鹽代表,研究鹽度對DNAPL在飽和多孔介質(zhì)中運移和分布的影響。通過批次實驗測定NaCl水溶液/石英砂/PCE三相體系下的接觸角和界面張力,結(jié)果表明,PCE在石英砂表面的接觸角隨著水中NaCl濃度的增大而減小,而PCE和NaCl水溶液的界面張力隨著NaCl濃度的增大而增大,尤其當(dāng)氯化鈉濃度較高時(>0.1 mol/L),影響程度更為顯著。在此基礎(chǔ)上,采用透射光法監(jiān)測不同介質(zhì)情景下DNAPL在二維砂箱中的運移和分布,定量測定DNAPL在介質(zhì)中的飽和度。實驗結(jié)果表明,地下水鹽度的增加將促進DNAPL的垂向入滲,減少被截留在運移路徑上的DNAPL量,使得DNAPL運移路徑及累積形成的池狀DNAPL(pool)向水流方向偏移。在均質(zhì)多孔介質(zhì)和含有透鏡體的非均質(zhì)多孔介質(zhì)中,隨著鹽度的增加,DNAPL在橫向和垂向上的展布均呈現(xiàn)出增加趨勢,導(dǎo)致污染源區(qū)變大,且介質(zhì)中以離散狀存在的DNAPL量明顯增加。

    四氯乙烯; 多孔介質(zhì);透射光法; 鹽度;運移;重非水相液體;飽和度

    在多孔介質(zhì)中,毛細(xì)力、粘滯力和重力等主要作用力控制著重非水相液體(dense non- aqueous phase liquid,DNAPL)的運移和分布[1~3],這些作用力受DNAPL的理化性質(zhì)、地下水化學(xué)性質(zhì)、巖性等因素的影響。其中毛細(xì)力由界面張力和接觸角(含水層介質(zhì)潤濕性的衡量指標(biāo))確定。潤濕性是指一種流體在另一種不相溶流體存在時優(yōu)先鋪展于固體介質(zhì)表面的趨勢[4],通常用接觸角(θ)表示,當(dāng)θ在0°~70°時為水相潤濕,當(dāng)接觸角為70°~120°時為中間潤濕狀態(tài),當(dāng)接觸角大于120°時為油相潤濕狀態(tài)。一般情況下,含水層的介質(zhì)是水潤濕的,但是近些年的研究表明潤濕性可以隨著DNAPL化學(xué)性質(zhì)、孔隙結(jié)構(gòu)、含水層中礦物的性質(zhì)和地下水化學(xué)性質(zhì)等因素的改變而發(fā)生改變[5]。地下水化學(xué)性質(zhì)包括pH、離子強度等,這些因素由于自然和人類的影響而發(fā)生很大的改變,這些改變會影響粘滯力以及界面張力和潤濕性,進而影響DNAPL的運移和分布[6]。由于這些物理化學(xué)因素的復(fù)雜性,使得精確描述DNAPL在地下環(huán)境中的形態(tài)變得非常困難,然而,精確描述DNAPL在地下水環(huán)境中運移和分布又是極其重要的。

    溶解鹽在地下水中普遍存在,其中以氯化鈉最為常見。受人類和自然因素的影響,這些溶解鹽在地下水中的濃度變化范圍較大。關(guān)于化學(xué)物質(zhì)對DNAPL在多孔介質(zhì)中運移的影響,前人已經(jīng)有相關(guān)的研究[5,7~10]。Lord 等[7]探討了在NAPL/水/石英三相體系下化學(xué)條件的改變對界面張力和接觸角的影響。Ondovcin等[8]和Boufadel等[9]的理論研究表明水相中鹽濃度的升高會使地下水的密度和粘度升高。Molnar等[5]研究了在DNAPL存在條件下,介質(zhì)潤濕性的改變對DNAPL運移的影響。Chen等[10]通過一維砂柱實驗研究了離子強度對飽和多孔介質(zhì)中有機污染物運移的影響。這些研究大多數(shù)都是基于理論的研究以及在溶液、石英玻璃表面或者一維砂柱中進行的實驗。本文選用在地下環(huán)境中普遍存在的四氯乙烯(PCE)為典型DNAPL污染物,以NaCl為地下水中溶解鹽代表,通過接觸角和界面張力的測定,結(jié)合室內(nèi)二維砂箱實驗,采用透射光法定量研究地下水鹽度的變化對PCE在飽和多孔介質(zhì)中遷移和分布的影響。本文的研究結(jié)果可以拓展對DNAPL在地下水中遷移的認(rèn)識,尤其是為地下水鹽度較高地區(qū)的相關(guān)重非水相污染物研究提供一定的參考。

    1 材料與方法

    1.1 實驗材料

    本文選用純度為99.5%的四氯乙烯(PCE)(分析純)作為DNAPL的典型代表物,PCE購于Sigma- Aldrich化學(xué)品公司。PCE的密度為1.625 g/mL (25 ℃),分子量為165.82 g/mol,水中溶解度約224 mg/L (25 ℃),PCE和水的界面張力為44.4 mN/m[11]。用溶于有機物而不溶于水的染色劑油紅-O(Fisher Scientific,F(xiàn)air Lawn,NJ) 將無色的PCE染成紅色,染色濃度為0.05 g/L,以便于觀察和透射光法分析,染色后PCE的理化性質(zhì)無顯著改變[11]。氯化鈉為分析純,購于中東化玻有限公司。用超純水配制不同濃度的氯化鈉溶液備用,設(shè)置氯化鈉的濃度范圍為0.005 85 ~58.5 g/L,這與地下水中鹽度的變化范圍類似。選擇0.005 85 g/L和32 g/L的NaCl溶液作為四組砂箱實驗的背景溶液,研究鹽度對DNAPL遷移和分布的影響。本實驗所用流體的物理化學(xué)性質(zhì)見表1。

    表1 實驗中所用流體的物理化學(xué)性質(zhì)

    A=Air; l=W (water) or N (NAPL); n.a.=not applicable 。

    a來源于 Taylor 等[11],在 25 ℃條件下測定;b在 22 ℃條件下測定。

    砂箱中填充介質(zhì)為半透明石英砂,實驗選用六種規(guī)格的美國Vnimin公司生產(chǎn)的Accusand(表2)。Accusand具有較高的純度,有機物含量極低,并且具有較好的透光性,適用于在透射光法中使用。石英砂的清洗方法和裝填方法參見文獻(xiàn)[12]。

    1.2 界面張力和接觸角的測定

    接觸角測定采用懸滴法[7],利用2 mL的測微注射器將PCE滴到浸沒在水相溶液中的石英玻璃平板表面,待其穩(wěn)定后,利用接觸角測定儀(DSA 100, Kruss, Germany)測定,采用軸對稱形狀分析法(axisymmetric drop shape analysis,ADSA)測算接觸角。DNAPL和不同濃度NaCl溶液的界面張力用ADSA方法測定,儀器為旋轉(zhuǎn)滴界面張力儀(Model 500, Kino Industry,USA),樣品準(zhǔn)備與接觸角測定一樣。實驗在室溫條件下進行。

    表2 所用石英砂的相關(guān)性質(zhì)

    1.3 二維砂箱實驗

    采用一種玻璃與鋁合金材料構(gòu)建的二維砂箱對DNAPL在多孔介質(zhì)中的運移進行物理模擬(圖 1)。砂箱尺寸為60 cm (長)×45 cm (高)×1.4 cm (厚),砂箱由三個鋁制框和兩片鋼化玻璃組成,有利于DNAPL運移的可視化觀察。兩個U型井位于砂箱的兩側(cè),以便于實現(xiàn)均勻流流場。井在整個深度上都被不銹鋼鋼網(wǎng)包裹。砂箱頂端采用可拆卸的頂蓋實現(xiàn)承壓含水層狀態(tài),利用橡膠條密封鋁制框與鋼化玻璃之間的間隙。外部鋁框用螺絲固定,確保內(nèi)框與玻璃間的牢固,同時通過擠壓橡膠條達(dá)到隔水的目的。水流通過蠕動泵從左側(cè)三個均勻分布的流入口流入砂箱,出水從三個流出口流出進入廢液瓶。在砂箱中間位置的上端設(shè)置注入點模擬DNAPL從污染源區(qū)的泄漏。

    本實驗共設(shè)置四組砂箱。兩組砂箱裝填20~30目的石英砂模擬均質(zhì)背景介質(zhì)(圖1a),分別標(biāo)記為Fw- 1和Sw- 1;另兩組砂箱,以20~30 目的石英砂作為背景介質(zhì),設(shè)置不同粒徑組成的層狀透鏡體模擬非均質(zhì)介質(zhì)(圖1b),標(biāo)記為Fw- 2和Sw- 2。在所有四組砂箱中,底部和頂部均設(shè)置1 cm厚的細(xì)砂層(70~80目)以防止PCE從砂箱頂部或底部溢出。Fw- 1、Sw- 1和Fw- 2、Sw- 2四組砂箱的孔隙度分別為0.278、0.269、0.302和0.298。

    圖1 實驗砂箱示意圖Fig.1 Schematic diagram of the 2- D sandboxa—70~100目; b—70~80目; c—40~50目; d—50~70目; e—30~40目

    在每組實驗開始前,先從左端流入口泵入背景溶液對整個砂箱進行水力沖刷,流速恒定在1.5 m/d,該流速值與自然環(huán)境中常見的地下水流速相一致。其中,在Fw- 1和Fw- 2實驗中泵入的是0.005 85 g/L的 NaCl溶液,在Sw- 1和Sw- 2實驗中泵入的是32 g/L 的NaCl溶液。利用氣密注射泵(Cole Palmer 74900 series, Anjou, QC)將一定量的染色PCE以0.5 mL/min的恒定速度通過注入孔注入到砂箱中,PCE在砂箱中充分運移分布直至達(dá)到穩(wěn)定(24 h)。在這個過程中按照預(yù)設(shè)的時間間隔,利用透射光監(jiān)測系統(tǒng)記錄砂箱中PCE的運移和分布。

    1.4 分析方法

    本文使用一種改進的透射光法(light trans- mission method,LTM)分析染色PCE在砂箱中的遷移和再分布行為[13]。該透射光監(jiān)測系統(tǒng)主要包含二維砂箱、光源和一個CCD (AP2E, Apogee Instrumnets, Auburn, CA)相機以及相關(guān)的控制軟件(Maxim DL, Ottawa, ON)。砂箱的一側(cè)為由6根日光燈管組成的燈箱,作為光源,CCD相機放在砂箱的另一側(cè),相機鏡頭(Nikon D90)對準(zhǔn)砂箱中間位置,距離砂箱大約1.8 m。實驗基于染色PCE的光吸附峰值,選擇543 nm的帶通濾波器,濾波器安裝在相機鏡頭的前面。砂箱與相機之間的空間以遮光板封閉,以保證位于砂箱另一側(cè)的 CCD 相機所接受的光線皆來自燈箱發(fā)出穿過砂箱后的光線。在實驗期間,CCD相機設(shè)置為每隔一定時間拍照一次記錄光強變化,對DNAPL的入滲過程和再分布過程進行動態(tài)監(jiān)測。采集的圖像使用MATLAB R2012a 軟件分析圖像和計算飽和度。實驗過程的室內(nèi)溫度保持在22 ℃±1.0 ℃。

    砂箱中PCE的飽和度通過透射光光強值計算轉(zhuǎn)換而來,之前的研究已經(jīng)表明了兩者的函數(shù)關(guān)系[13]。本次實驗是模擬飽和含水層,只有水油兩相,DNAPL的飽和度(So)計算利用Bob等[13]建立的模型,兩相系統(tǒng)內(nèi)的油相飽和度計算公式為:

    (1)

    式中:Is與Ioil——在完全飽水和與完全飽油條件下的每一個計算單元的光強值;

    So——油相的飽和度;

    I——透射光光強。

    2 結(jié)果及討論

    2.1 界面張力和接觸角分析

    在飽和多孔介質(zhì)中,相與相之間的界面特性對油相(DNAPL)在地下水環(huán)境中的運移和最終形態(tài)有著重要影響。圖2顯示的是在不同氯化鈉濃度條件下,PCE/NaCl溶液界面張力以及PCE在石英表面上的接觸角(θ)。當(dāng)水相中無NaCl存在時,PCE/水的界面張力及接觸角分別為47.5 mN/m 和44.4°,這與Molnar等[5]報道的結(jié)果一致。從圖中可以看出,在由水相(NaCl水溶液)、油相(PCE)和固相(石英玻璃)組成的三相體系中,隨著水相中NaCl濃度的增大,水/油兩相之間的界面張力表現(xiàn)為增大的趨勢,PCE在石英玻璃表面的接觸角呈現(xiàn)出減小的趨勢。當(dāng)氯化鈉的濃度從1×10-4mol/L增加至1 mol/L時,PCE/NaCl溶液的界面張力從47.5 mN/m增加到55.5 mN/m(圖2a),接觸角從44.4°降低至36°(圖2b)。當(dāng)氯化鈉的濃度變化范圍為1×10-4~0.01 mol/L時,界面張力和接觸角的改變幅度均較??;當(dāng)氯化鈉的濃度超過0.1 mol/L時,界面張力急劇增加;當(dāng)氯化鈉的濃度超過0.01 mol/L時,接觸角的變化顯著。上述實驗結(jié)果表明,當(dāng)?shù)叵滤}度在較低范圍內(nèi)變化時(NaCl濃度< 0.1 mol/L),對DNAPL/水/石英砂三相之間的界面特性產(chǎn)生的影響較?。划?dāng)?shù)叵滤}度較高時(>0.1 mol/L),對上述三相體系的接觸角尤其是液液界面張力影響顯著,從而可能進一步影響到DNAPL在多孔介質(zhì)中的運移和分布。

    圖2 氯化鈉濃度對DNAPL/(NaCl溶液/石英砂體系中)界面張力和接觸角的影響Fig.2 Effects of NaCl concentration on interfacial tension and advancing contact angle in the PCE/(NaCl solution/quartz) system

    2.2 PCE在砂箱中的運移

    選擇0.005 85 g/L(1×10-4mol/L)和32 g/L(0.55 mol/L)的NaCl溶液作為背景溶液進行實驗,定量研究鹽度對DNAPL運移和分布的影響。圖3為不同時刻PCE的飽和度分布結(jié)果。PCE泄漏到砂箱中以后,PCE快速的垂向入滲,同時伴隨一定的側(cè)向遷移。實驗結(jié)束后,在運移路徑上存在一定量截留PCE,在由細(xì)砂構(gòu)成的透鏡體和砂箱底端細(xì)砂層上方有累積形成的DNAPL污染池(pool)。本次運移實驗的地下水流速為1.5 m/d,該側(cè)向水流對DNAPL污染羽向水流方向移動有一定的影響,但影響較弱,與課題組前期研究成果一致[14]。

    在均質(zhì)介質(zhì)中(圖3a,3b),注入的PCE在重力作用下快速向下運移,形成一個邊緣較尖銳的污染羽。在Fw- 1和Sw- 1兩組實驗中,均可以觀察到較為明顯的優(yōu)先流(finger flow)現(xiàn)象,這是由裝砂過程中人為因素導(dǎo)致的局部非均質(zhì)造成的,表明即使是很微小的非均質(zhì)環(huán)境也會對DNAPL的移動產(chǎn)生很強的影響[15]。在非均質(zhì)介質(zhì)中,地下水鹽度較低時(Fw- 2),PCE向下遷移到達(dá)了其中4個透鏡體(a,c,d,e)(圖3 c),在116 min時運移到底端形成三個pool,并在148 min時融為一個pool;而在Sw- 2實驗中,DNAPL運移到達(dá)所有的透鏡體(圖3 d),且實驗進行80 min后,PCE即到達(dá)砂箱底端形成兩個污染池,并在128 min時兩個pool融為一個pool。

    對DNAPL在砂箱底端細(xì)砂層上污染池(pool)的起始位置和長度進行分析,同時計算PCE污染羽的平均飽和度和底端pool的最大飽和度。在Fw- 1中(圖 3a),PCE污染羽幾乎均勻地分布在注入點的兩側(cè),當(dāng)PCE分布達(dá)到穩(wěn)定后(T= 24 h),pool的起始位置開始于X=15.27 cm,即注入點所在位置的上游14.73 cm,整個pool的長度為34.28 cm (表 3);而在Sw- 1中(圖 3b),DNAPL污染池開始于10.19 cm處,整個pool的長度達(dá)到48.75 cm,表明地下水鹽度的增大促進DNAPL污染池的長度增加。高鹽度地下水的粘滯性比低鹽度大,水流粘滯力的作用促進了DNAPL污染池的移動[16],因此在高鹽度組Sw- 1中所形成的pool池較長。在Fw- 1和Sw- 1兩組實驗中,運移達(dá)到穩(wěn)定后被截留在運移路徑上(計算區(qū)域為圖中的虛線框區(qū)域)的DNAPL量分別為15.72 mL和12.26 mL。界面張力的增大和接觸角的減小(圖2),將導(dǎo)致毛管阻力的增加[17],不利于DNAPL的向下運移。鹽度的增加會增加地下水粘滯性(表1),從而促進DNAPL遷移。綜合分析,鹽度變化造成的粘滯性改變在本研究中起主要作用,最終表現(xiàn)為促進DNAPL的向下運移并減小DNAPL在運移路徑上的截留量。非均質(zhì)情形下,當(dāng)DNAPL分布穩(wěn)定后,F(xiàn)w- 2和Sw- 2實驗中砂箱底部污染池的起始位置分別為12.76 cm和 21.7 cm,長度分別為40.85 和 37.02 cm,污染池長度縮短,這與Fw- 1和Sw- 1實驗中得到的結(jié)果不一致。這是因為在Sw- 2中PCE到達(dá)了每一個透鏡體,使得相較于Fw- 2情況下的到達(dá)底部的DNAPL的總量減小。Fw- 2和Sw- 2實驗中底部污染池的中心位置分別位于X=33.7 cm 和X=43.5 cm,表明非均質(zhì)介質(zhì)中地下水鹽度的增加也會促進DNAPL污染池向水流方向移動。

    圖3 運移過程中PCE的飽和度分布Fig.3 Infiltration and redistribution of PCE

    在所有4組砂箱實驗中,隨著地下水鹽度的增加,DNAPL污染羽的平均飽和度以及砂箱底部細(xì)砂層上的DNAPL污染池的最大飽和度均隨之降低(表3)。對砂箱中運移達(dá)到穩(wěn)定的PCE污染羽偏移角(注入點與質(zhì)量中心點所在直線與注入點所在垂直線形成的角度)進行分析計算。計算得到的數(shù)據(jù)(表3)表明,高鹽度水流促進DNAPL向水流方向偏移。均質(zhì)條件下,兩組實驗的偏移角分別為1.92°(Fw- 1)和 3.97°(Sw- 1);在非均質(zhì)條件下,兩組實驗的偏移角為18.7°(Fw- 2)和 29.5°(Sw- 2)。以上這些實驗結(jié)果均表明,地下水鹽度對飽和多孔介質(zhì)中DNAPL的遷移和形態(tài)分布有較大影響,地下水鹽度的增加將促進DNAPL的垂向入滲,并使得DNAPL運移路徑及累積形成的池狀DNAPL向水流方向偏移。

    表3 DNAPL污染羽在介質(zhì)中的分布特征

    為了進一步說明地下水鹽度變化對DNAPL在多孔介質(zhì)運移的影響,本文計算并分析了不同實驗條件下PCE入滲鋒面的運移速率,結(jié)果如圖4所示。鹽度對DNAPL入滲鋒面的運移有明顯影響,32 g/L NaCl水流條件下PCE入滲鋒面的運移速率明顯加快。在Fw- 1/Sw- 1和Fw- 2/Sw- 2四組實驗中,PCE到達(dá)砂箱底端細(xì)砂層的時間分別為30、14.0、92和56 min。結(jié)合不同濃度氯化鈉濃度條件下的接觸角和界面張力的測定,地下水鹽度增大時,界面張力隨之增大,同時接觸角減小使得石英砂潤濕性從水潤濕狀態(tài)變?yōu)閺娝疂櫇駹顟B(tài)。一般地,根據(jù)Young- Laplace方程,隨著界面張力的增大和接觸角的減小,PCE需要克服更大的毛管阻力向下運移,而本實驗結(jié)果顯示在高鹽度的地下水流條件下,PCE的下滲鋒面運移的更快。這進一步說明,本研究中地下水鹽度變化所導(dǎo)致的粘滯性改變主要影響著DNAPL在多孔介質(zhì)中的運移行為。

    圖4 PCE入滲鋒面平均運移速率Fig.4 Average advancing velocities of the PCE fronts

    2.3 DNAPL污染羽的形態(tài)特征

    2.3.1 展布范圍

    本文對DNAPL污染羽在水平和垂直方向上的展布范圍(二階矩)隨時間的變化規(guī)律進行了分析,如圖5所示。在前期注入過程中,由于存在注入壓力,污染羽快速擴展,表現(xiàn)為水平和垂直的展布范圍快速增大,隨著時間的延長,擴展逐漸減慢,最后達(dá)到穩(wěn)定。圖5a顯示的是在20~30目均質(zhì)條件下的展布范圍,在32 g/L NaCl的高鹽度地下水流條件下,由于增加的粘滯性起主要作用,無論是在垂直方向還是水平方向的展布范圍與0.005 85 g/L NaCl水流條件下相比均明顯增大。在水平方向上的展布隨著鹽度的增加,其增加幅度尤為顯著,主要是由地下水鹽度的增加促進砂箱底部的污染池向水流方向移動使得其長度和面積變大造成的。在非均質(zhì)條件下,PCE污染羽展布范圍的變化規(guī)律呈現(xiàn)出與均質(zhì)條件下相似的變化規(guī)律(圖5b),并且透鏡體的存在使得PCE橫向和縱向展布在注入期內(nèi)的增加趨勢表現(xiàn)得更加明顯。

    圖5 PCE在砂箱中的展布范圍Fig.5 Horizontal spread ) and vertical

    2.3.2 DNAPL污染羽的飽和度分布

    圖6顯示的是每一組砂箱實驗中不同時刻(與圖3中的時間對應(yīng))的PCE飽和度直方圖,橫坐標(biāo)表示飽和度(飽和度間隔0.01),縱坐標(biāo)為該飽和度區(qū)間PCE的體積分?jǐn)?shù)。從圖中可以看出,和0.005 85 g/L NaCl水流條件下相比,PCE飽和度峰值在高鹽度地下水流條件下均向左移動,即無論是在均質(zhì)還是在非均質(zhì)條件下,鹽度的增加都會導(dǎo)致砂箱中被截留PCE的飽和度降低,飽和度直方圖中峰值對應(yīng)的飽和度區(qū)間也相應(yīng)減小,還表現(xiàn)為DNAPL飽和度均值的降低(表3)。這進一步表明,即使是在較低地下水流速條件下(1.5 m/d),高鹽度地下水同樣會促進DNAPL遷移,使得DNAPL以較小的飽和度殘留。

    圖6 PCE四組實驗中PCE在不同時刻的飽和度直方圖Fig.6 Histograms of PCE saturation of the four experiments at different times

    2.3.3 “GTP”

    DNAPL飽和度低于殘余飽和度時稱為不連續(xù)的離散狀DNAPL(ganglia),相應(yīng)地,高于殘余飽和度的稱為pool?!癎TP”(ganglia to pool)即為離散狀DNAPL和池狀DNAPL的比值,本文采用該比值來表征PCE污染羽的結(jié)構(gòu)特征[1]。本文所用20~30目石英砂的殘留飽和度為0.11[18]。四組實驗(Fw- 1,Sw- 1,F(xiàn)w- 2,Sw- 2)均達(dá)到穩(wěn)定后,對應(yīng)的GTP值分別為0.046,0.083,0.076和0.13。在Fw- 1實驗中,GTP值為0.046即表明有1.41 mL(4.4%)的PCE以ganglia離散狀形式存在。分別比較均質(zhì)和非均質(zhì)介質(zhì)中在不同鹽度條件下的GTP值,可以看出地下水鹽度的變化對多孔介質(zhì)中PCE污染羽的結(jié)構(gòu)特征具有非常顯著的影響。地下水鹽度的增加不僅能夠減少殘留在運移路徑上的DNAPL總量,還使得污染區(qū)域中以離散狀(ganglia)存在的DNAPL量明顯增加。

    3 結(jié)論

    (1)地下水中溶解鹽(NaCl)的存在會對不同相之間的界面特性產(chǎn)生影響。在石英砂/水/PCE三相體系里,PCE和水之間的界面張力隨著水中NaCl濃度的增大而增大,PCE在石英砂表面的接觸角隨著NaCl濃度的增大而減小,尤其當(dāng)NaCl濃度高于0.1 mol/L時,影響顯著。

    (2)無論是在均質(zhì)還是在含有層狀透鏡體的非均質(zhì)介質(zhì)中,地下水鹽度的增加都會促進DNAPL入滲鋒面的垂直運移,使得截留在運移路徑上的殘留DNAPL總量減少。此外,地下水鹽度的增加還會使得DNAPL運移路徑及累積形成的連續(xù)池狀DNAPL(pool)向水流方向偏移。

    (3)在粗砂介質(zhì)(20~30目石英砂)中,地下水鹽度的增加直接影響DNAPL污染羽的形態(tài)特征。DNAPL在水平方向和垂直方向上的展布范圍均隨著地下水鹽度的增加而變大,污染羽的GTP值也隨之增加,介質(zhì)中以不連續(xù)離散狀(ganglia)形式存在的DNAPL量明顯增加。

    [ 1] 鄭菲,高燕維,徐紅霞,等.非均質(zhì)介質(zhì)性對DNAPL污染區(qū)結(jié)構(gòu)特征影響的實驗研究[J].水文地質(zhì)工程地質(zhì), 2016,43(5):140-148. [ZHENG F, GAO Y W, XU H X,etal. A experimental study of the influence heterogeneity on the source- zone architecture[J]. Hydrogeology & Engineering Geology, 2016, 43(5): 140-148.(in Chinese)]

    [ 2] 張燁,施小清,鄧亞平,等.結(jié)合蒸汽和空氣注入修復(fù)多孔介質(zhì)中的DNAPL污染物的多目標(biāo)多相流模擬優(yōu)化[J].水文地質(zhì)工程地質(zhì),2015,42(5):141-148. [ZHANG Y, SHI X Q, DENG Y P,etal. Multi- objective multi- phase optimization with steam and air co- injection for DNAPL contaminant remediation in porous media[J]. Hydrogeology & Engineering Geology, 2015,42(5):141-148.(in Chinese)]

    [ 3] LANHSU H, DEMOND A. Influence of organic acid and organic base interactions on interfacial properties in napl- water systems[J]. Environmental Science and Technology, 2007,41(3): 897-902.

    [ 4] CRAIG F F. The reservoir engineering aspects of water flooding. Monograph Series[J]. Society of Petroleum Engineers, Richardson, 1971,3.

    [ 5] MOLNAR I L, O’CARROL D M,etal. Impact of surfactant- induced wettability alterations on DNAPL invasion in quartz and iron oxidecoated sand systems [J]. Journal of Contaminant Hydrology , 2011,119(4): 1-12.

    [ 6] XING L, GUO H M, ZHAN Y H. Groundwater hydrochemical characteristics and processes alongow paths in the North China Plain [J]. Journal of Asian Earth Sciences, 2013, 70: 250-264.

    [ 7] LORD D, DEMOND A, HAYES K. Effects of organic base chemistry on interfacial tension, wettability, and capillary pressure in multiphase subsurface waste systems [J]. Transport in Porous Media , 2000,38:79-92.

    [ 8] ONDOVCIN T, MLS J. Relation Between Phase Density And Componene concentration in groundwater Flow Modeling [J]. Journal of Hydrology and Hydromechanics, 2007, 55(4):236-245.

    [ 9] BOUFADEL M C, SUIDAN M T. A numerical model for density and viscosity dependent flows in two- dimensional variably saturated porous media [J]. Journal of Contamiant Hydrology , 1999,40(2):285-296.

    [10] CHEN H, GAO B, LI H. Effects of pH and ionic strength on sulfamethoxazole and ciprooxacin transport in saturated porous media [J]. Journal of Contaminant Hydrology, 2011,126(2): 29-36.

    [11] TAYLOR T P, PENNELL K D. Surfactant enhanced recovery of tetrachloroethylene from a porous medium containing low permeability lenses 1. Experimental studies [J]. Journal of Contaminant Hydrology , 2001, 48(3): 325-350.

    [12] YE S J, SLEEP B E, CHIEN C. The impact of methanogenesis onow and transport in coarse sand [J]. Journal of Contaminant Hydrology , 2009,103(2): 48-57.

    [13] BOB M M, BROOKS M C, WOOD A L. A modified light transmission visualization method for DNAPL saturation measurements in 2- D models [J]. Advances in Water Resources, 2008, 31(5):727-742.

    [14] 高燕維, 徐紅霞,吳吉春,等. 基于透射光法探討水流流速對DNAPL運移分布的影響[J]. 環(huán)境科學(xué), 2015, 36(7): 209-216.[GAO Y W, XU H X, WU J C,etal. Laboratory investigation of DNAPL migration behavior and distribution at varying flow velocities based on light transmission method [J]. Environmental Science, 2015, 36(7):209-216.(in Chinese)]

    [15] ZHENG F, GAO B, SUN Y,etal. Removal of tetrachloroethylene from homogeneous and heterogeneous porous media: Combined effects of surfactant solubilization and oxidant degradation[J]. Chemical Engineering Journal, 2016,283(1):595-603.

    [16] PADGEET P K, HAYDEN N J. Mobilization of residual tetrachloroethylene during alcohol flusing of clay- containing porous media. [J]. Journal of Contaminant Hydrology , 1999, 40(3): 285-296.

    [17] DONG J, CHOWDHRY B, LEHARNE S. Investigation of the wetting behavior of coal tar in three phase systems and its modification by poloxamine block copolymeric surfactants[J]. Environmental Science and Technology, 2004, 38(8): 594-602.

    [18] SUCHOMEL E J, PENNELL K D. Reductions in contaminant mass discharge following partial mass removal from DNAPL source zones [J]. Environmental Science and Technology, 2006, 40(19): 6110-6116.

    責(zé)任編輯:張若琳

    Effect of salinity on DNAPL migration and distribution in saturated porous media

    CHENG Zhou, XU Hongxia, SUN Yuanyuan, SHI Xiaoqing, LIAO Penghui, WU Jichun

    (StateKeyLaboratoryofPollutionControlandResourceReuse,DepartmentofHydrosciences/SchoolofEarthSciencesandEngineering,NanjingUniversity,Nanjing,Jiangsu210023,China)

    In this work the effect of salinity on the migration and distribution of DNAPL in saturated porous media was investigated. Tetrachloroethylene (PCE) was selected as the representative DNAPL (Dense Non- aqueous Phase liquid), and sodium chloride (NaCl) was used as the typical dissolved salt in groundwater. Batch experiments were first conducted to measure the contact angles and interfacial tensions (IFTs) of the NaCl solutions/PCE/quartz system. The results show that the contact angles of PCE decrease and IFTs between PCE and NaCl solution increase with the increasing concentration of NaCl, and the effects are more obvious at high concentrations of NaCl (>0.1 mol/L). Four 2- D flow cell experiments were then conducted, and the light transmission method was used to monitor the migration process of PCE and quantitatively measure the PCE saturation distribution in the porous media. The experimental results show that the increasing salinity in groundwater promotes the vertical migration of DNAPL, decreases the volumes of DNAPL entrapped along the migration path, and leads to the inclination of the migration path and pooled DNAPL towards the water flow direction. Moreover, in both homogeneous and heterogeneous (lenses) porous media, the horizontal and vertical spreading of DNAPL increases with the increasing salinity, resulting in a larger source- zone and the increasing PCE ganglia.

    PCE; porous media; light transmission method; salinity; migration; DNAPL; Saturation distribution

    10.16030/j.cnki.issn.1000- 3665.2017.04.20

    2016- 08- 15;

    2016- 10- 17

    國家自然科學(xué)基金(U1503282、41030746和41172207);教育部博士點基金(20110091110039); 南京大學(xué)研究生跨學(xué)科科研創(chuàng)新基金(2014CL10)

    程洲(1986- ),男,博士研究生,主要從事土壤及地下水中有機污染物遷移及修復(fù)研究。E- mail:czhounju@ 126.com

    吳吉春(1968- ),男,教授,主要從事水文地質(zhì)學(xué)、土壤和地下水污染研究。E- mail: jcwu@nju.edu.cn

    P641.3

    A

    1000- 3665(2017)04- 0129- 08

    猜你喜歡
    砂箱均質(zhì)運移
    淺談砂箱斷面尺寸對砂箱剛性的影響
    基于Ansys workbench 的單壁砂箱有限元分析
    曲流河復(fù)合點壩砂體構(gòu)型表征及流體運移機理
    一種嵌裝式砂型
    玻璃(2020年2期)2020-12-11 22:20:04
    一種用于熔鑄磚的新型水玻璃砂型
    玻璃(2020年2期)2020-12-11 22:20:04
    東營凹陷北帶中淺層油氣運移通道組合類型及成藏作用
    Orlicz對偶混合均質(zhì)積分
    開采過程中上覆急傾斜巖層運移規(guī)律模擬與研究
    川西坳陷孝泉-新場地區(qū)陸相天然氣地球化學(xué)及運移特征
    非均質(zhì)巖心調(diào)堵結(jié)合技術(shù)室內(nèi)實驗
    久久这里有精品视频免费| 给我免费播放毛片高清在线观看| 老司机影院成人| 精品日产1卡2卡| 久久亚洲国产成人精品v| 亚洲一区二区三区色噜噜| 黄色日韩在线| 直男gayav资源| 免费黄网站久久成人精品| 午夜激情福利司机影院| 人妻系列 视频| 亚洲最大成人手机在线| 黄色视频,在线免费观看| 日韩成人av中文字幕在线观看| 麻豆国产av国片精品| 精品久久久久久久久久久久久| 国产精品99久久久久久久久| 日韩一本色道免费dvd| 国产精品一区二区三区四区久久| 男人和女人高潮做爰伦理| 中文字幕人妻熟人妻熟丝袜美| 少妇高潮的动态图| www.色视频.com| 丝袜美腿在线中文| 色尼玛亚洲综合影院| 99久久九九国产精品国产免费| 黄片无遮挡物在线观看| 69人妻影院| 51国产日韩欧美| 国产高清三级在线| 日韩欧美 国产精品| 国产精品一区二区三区四区久久| 国产精品美女特级片免费视频播放器| 国产精品人妻久久久久久| 国内少妇人妻偷人精品xxx网站| 九九久久精品国产亚洲av麻豆| 1000部很黄的大片| 久久久色成人| 麻豆乱淫一区二区| 国产精品一二三区在线看| 日本一二三区视频观看| 欧美日韩乱码在线| 国产高清有码在线观看视频| 国产真实伦视频高清在线观看| 高清在线视频一区二区三区 | 久久久久国产网址| 日本成人三级电影网站| 日韩欧美三级三区| 午夜福利在线观看吧| 亚洲美女搞黄在线观看| 白带黄色成豆腐渣| 日韩欧美一区二区三区在线观看| 插阴视频在线观看视频| 最近最新中文字幕大全电影3| 插逼视频在线观看| 精品久久久久久久久av| 秋霞在线观看毛片| 亚洲成av人片在线播放无| 亚洲精品乱码久久久久久按摩| 国产黄a三级三级三级人| 99在线视频只有这里精品首页| 国内精品久久久久精免费| 欧美三级亚洲精品| 在线免费观看的www视频| 日韩高清综合在线| 国产成人精品一,二区 | 久久精品国产99精品国产亚洲性色| 麻豆成人av视频| 国产三级在线视频| 69av精品久久久久久| 人妻久久中文字幕网| 色综合色国产| 久久久久久久久久黄片| 久久午夜福利片| 久久精品夜色国产| 中文字幕人妻熟人妻熟丝袜美| 国产成人一区二区在线| 中文亚洲av片在线观看爽| 最好的美女福利视频网| 在线观看一区二区三区| 久久精品综合一区二区三区| 激情 狠狠 欧美| 12—13女人毛片做爰片一| 精品一区二区免费观看| 日日摸夜夜添夜夜爱| 国产v大片淫在线免费观看| 免费一级毛片在线播放高清视频| 在线a可以看的网站| 日韩一本色道免费dvd| 91aial.com中文字幕在线观看| av免费在线看不卡| 嫩草影院精品99| 非洲黑人性xxxx精品又粗又长| 日韩成人伦理影院| 97人妻精品一区二区三区麻豆| av视频在线观看入口| 亚洲av电影不卡..在线观看| 不卡一级毛片| 国产一区二区激情短视频| 欧美xxxx性猛交bbbb| 在线国产一区二区在线| avwww免费| 国产精品久久视频播放| 免费一级毛片在线播放高清视频| 婷婷色综合大香蕉| 国产亚洲av嫩草精品影院| 大又大粗又爽又黄少妇毛片口| 一区二区三区高清视频在线| 天美传媒精品一区二区| 真实男女啪啪啪动态图| 蜜臀久久99精品久久宅男| 欧美日韩乱码在线| 免费黄网站久久成人精品| 欧美最黄视频在线播放免费| av卡一久久| 欧美在线一区亚洲| 美女内射精品一级片tv| 久久精品国产亚洲网站| 别揉我奶头 嗯啊视频| av在线老鸭窝| 国产伦在线观看视频一区| 日本色播在线视频| 国产亚洲av嫩草精品影院| 可以在线观看毛片的网站| 亚洲精品自拍成人| 成人综合一区亚洲| 91aial.com中文字幕在线观看| 69人妻影院| 国产av一区在线观看免费| 身体一侧抽搐| 两个人视频免费观看高清| 国产精品伦人一区二区| 欧美日韩综合久久久久久| АⅤ资源中文在线天堂| 人妻久久中文字幕网| 美女脱内裤让男人舔精品视频 | 简卡轻食公司| 91狼人影院| 日韩精品有码人妻一区| 一本久久精品| 搞女人的毛片| 国产一区亚洲一区在线观看| 韩国av在线不卡| 国产精品国产高清国产av| 99久久中文字幕三级久久日本| 精品欧美国产一区二区三| 免费av观看视频| 日本色播在线视频| 久久久国产成人免费| 色5月婷婷丁香| 少妇丰满av| 一本久久精品| 在线播放无遮挡| 久久精品夜夜夜夜夜久久蜜豆| 午夜免费激情av| 久久这里有精品视频免费| 黄片无遮挡物在线观看| 能在线免费看毛片的网站| 成人综合一区亚洲| 最新中文字幕久久久久| 久久九九热精品免费| 18+在线观看网站| 色噜噜av男人的天堂激情| 亚洲va在线va天堂va国产| 日韩一区二区视频免费看| 国产精品福利在线免费观看| 最近手机中文字幕大全| 国产探花极品一区二区| 久久久久久伊人网av| 日韩 亚洲 欧美在线| 国产精品乱码一区二三区的特点| 国产私拍福利视频在线观看| 一级黄片播放器| 黄色一级大片看看| 麻豆国产97在线/欧美| 夜夜夜夜夜久久久久| 人人妻人人澡人人爽人人夜夜 | 国产黄片美女视频| 久久久久久国产a免费观看| 成人二区视频| 22中文网久久字幕| 岛国在线免费视频观看| 99热精品在线国产| 婷婷六月久久综合丁香| 精品一区二区三区视频在线| www.av在线官网国产| 亚洲国产欧洲综合997久久,| 亚洲一区高清亚洲精品| 在线a可以看的网站| 有码 亚洲区| 女的被弄到高潮叫床怎么办| 国内精品美女久久久久久| 听说在线观看完整版免费高清| 欧美极品一区二区三区四区| 免费观看在线日韩| 一个人观看的视频www高清免费观看| 十八禁国产超污无遮挡网站| 中文欧美无线码| 一边摸一边抽搐一进一小说| 麻豆成人av视频| av免费在线看不卡| 亚洲精品456在线播放app| 女同久久另类99精品国产91| 免费大片18禁| 免费观看精品视频网站| 精品久久久久久久久亚洲| 日韩欧美精品免费久久| 在线播放无遮挡| 成人国产麻豆网| 亚洲欧美日韩高清专用| 波野结衣二区三区在线| 久久久精品94久久精品| 久久精品国产亚洲av天美| 成年女人看的毛片在线观看| 深爱激情五月婷婷| av在线老鸭窝| 91av网一区二区| 男人的好看免费观看在线视频| 国产欧美日韩精品一区二区| 亚洲综合色惰| 免费看日本二区| 久久精品91蜜桃| 淫秽高清视频在线观看| 乱码一卡2卡4卡精品| 国产三级中文精品| 久久九九热精品免费| 国产成人freesex在线| 色吧在线观看| 夜夜夜夜夜久久久久| av福利片在线观看| 精品免费久久久久久久清纯| 亚洲三级黄色毛片| 国产一区二区三区av在线 | 久久久久久九九精品二区国产| 久久久久久大精品| 人妻少妇偷人精品九色| 国产色爽女视频免费观看| 亚洲,欧美,日韩| 国产精品久久久久久久电影| 亚洲欧美成人综合另类久久久 | av在线天堂中文字幕| 一进一出抽搐动态| 精品国内亚洲2022精品成人| 日本黄大片高清| 亚州av有码| 国产中年淑女户外野战色| 亚洲av二区三区四区| 亚洲国产欧美人成| 国内精品久久久久精免费| 欧美日韩在线观看h| 性色avwww在线观看| 亚洲欧美成人精品一区二区| 又粗又硬又长又爽又黄的视频 | 一区二区三区高清视频在线| av卡一久久| 亚洲av中文字字幕乱码综合| 免费观看a级毛片全部| 国产精品嫩草影院av在线观看| 久久鲁丝午夜福利片| 欧美精品国产亚洲| 亚洲天堂国产精品一区在线| 成人国产麻豆网| 午夜免费激情av| 国产精品电影一区二区三区| 欧美不卡视频在线免费观看| 2021天堂中文幕一二区在线观| 久久精品国产亚洲av香蕉五月| 欧美性感艳星| 别揉我奶头 嗯啊视频| 一级黄色大片毛片| 国产乱人偷精品视频| 人妻制服诱惑在线中文字幕| 日产精品乱码卡一卡2卡三| 有码 亚洲区| 人妻少妇偷人精品九色| 国语自产精品视频在线第100页| 国产一区二区激情短视频| 亚洲久久久久久中文字幕| 午夜福利高清视频| 性色avwww在线观看| 成人午夜精彩视频在线观看| 日韩欧美 国产精品| 久久亚洲精品不卡| 欧美激情在线99| 大又大粗又爽又黄少妇毛片口| 久久精品夜夜夜夜夜久久蜜豆| 非洲黑人性xxxx精品又粗又长| 永久网站在线| 成人无遮挡网站| 亚洲性久久影院| 麻豆国产97在线/欧美| 国产乱人偷精品视频| 国产极品天堂在线| 久久久国产成人免费| 三级男女做爰猛烈吃奶摸视频| 综合色丁香网| 麻豆成人av视频| 国产一区二区在线av高清观看| 欧洲精品卡2卡3卡4卡5卡区| 深夜a级毛片| 99在线视频只有这里精品首页| 国产高清不卡午夜福利| 免费搜索国产男女视频| 国产在线男女| 国产精品人妻久久久影院| 91精品国产九色| 亚洲精品日韩在线中文字幕 | 一进一出抽搐gif免费好疼| 免费人成视频x8x8入口观看| 国产精品乱码一区二三区的特点| 精品免费久久久久久久清纯| 久久亚洲国产成人精品v| 国产高清视频在线观看网站| 九九爱精品视频在线观看| 国产午夜精品久久久久久一区二区三区| 国产精品美女特级片免费视频播放器| 亚洲第一区二区三区不卡| 97在线视频观看| 亚洲精品456在线播放app| 校园人妻丝袜中文字幕| 午夜爱爱视频在线播放| 99精品在免费线老司机午夜| 中国美白少妇内射xxxbb| 插逼视频在线观看| 久久久久久久午夜电影| 久久久国产成人精品二区| 在线观看av片永久免费下载| 精品一区二区三区视频在线| 亚洲中文字幕一区二区三区有码在线看| 99国产精品一区二区蜜桃av| 看十八女毛片水多多多| 卡戴珊不雅视频在线播放| 亚洲欧美日韩东京热| 国产一区二区在线av高清观看| 天堂网av新在线| 国产蜜桃级精品一区二区三区| 久久草成人影院| 日日啪夜夜撸| 少妇人妻精品综合一区二区 | 麻豆国产97在线/欧美| 久久久a久久爽久久v久久| 久久亚洲精品不卡| 九九爱精品视频在线观看| eeuss影院久久| 成年免费大片在线观看| 国产av在哪里看| 国产黄色视频一区二区在线观看 | 床上黄色一级片| 99热这里只有是精品50| av又黄又爽大尺度在线免费看 | av在线亚洲专区| 久久久精品94久久精品| 赤兔流量卡办理| 国国产精品蜜臀av免费| 国产欧美日韩精品一区二区| 国产一区二区在线av高清观看| 欧美最黄视频在线播放免费| 五月玫瑰六月丁香| 亚洲在线自拍视频| 爱豆传媒免费全集在线观看| 99热这里只有是精品50| 自拍偷自拍亚洲精品老妇| 亚洲欧美清纯卡通| 天堂av国产一区二区熟女人妻| av国产免费在线观看| 国产老妇伦熟女老妇高清| 亚洲中文字幕日韩| 亚洲欧美清纯卡通| 国产成人精品久久久久久| 99久国产av精品国产电影| 在线天堂最新版资源| 一边摸一边抽搐一进一小说| 插逼视频在线观看| 国产老妇伦熟女老妇高清| 精品熟女少妇av免费看| 插阴视频在线观看视频| 成年女人看的毛片在线观看| 亚洲国产色片| 人妻夜夜爽99麻豆av| 亚洲婷婷狠狠爱综合网| 人体艺术视频欧美日本| 亚洲va在线va天堂va国产| 女同久久另类99精品国产91| 一级毛片我不卡| 一级黄色大片毛片| 国内精品宾馆在线| 成人特级黄色片久久久久久久| 男女视频在线观看网站免费| 少妇高潮的动态图| 97热精品久久久久久| 夫妻性生交免费视频一级片| 欧美成人一区二区免费高清观看| 国产精品国产高清国产av| 一夜夜www| 国产免费男女视频| 国产精品久久久久久av不卡| 女同久久另类99精品国产91| 性色avwww在线观看| 欧美区成人在线视频| 午夜免费激情av| 在线天堂最新版资源| 午夜老司机福利剧场| 欧美区成人在线视频| 亚洲五月天丁香| 精品不卡国产一区二区三区| 亚洲精品久久国产高清桃花| .国产精品久久| 欧美zozozo另类| 国产精品.久久久| 亚洲欧美日韩卡通动漫| 黄色一级大片看看| 毛片一级片免费看久久久久| 一个人看的www免费观看视频| 欧美成人一区二区免费高清观看| 人体艺术视频欧美日本| 美女cb高潮喷水在线观看| 欧美日韩综合久久久久久| 插阴视频在线观看视频| 男的添女的下面高潮视频| 精品一区二区三区视频在线| av在线蜜桃| 国产老妇伦熟女老妇高清| 男的添女的下面高潮视频| 99热全是精品| 免费观看a级毛片全部| 91精品国产九色| 亚洲乱码一区二区免费版| 免费看日本二区| 看非洲黑人一级黄片| 精品日产1卡2卡| 国产不卡一卡二| 精品久久久久久成人av| 深夜精品福利| 国产在线男女| 18禁在线播放成人免费| 内射极品少妇av片p| 亚洲七黄色美女视频| 特级一级黄色大片| 亚洲精品影视一区二区三区av| 亚洲国产高清在线一区二区三| 亚洲四区av| 我的老师免费观看完整版| 在现免费观看毛片| 日韩一本色道免费dvd| 成人鲁丝片一二三区免费| 哪个播放器可以免费观看大片| 久久精品国产亚洲av天美| av在线亚洲专区| 国产黄色小视频在线观看| 91久久精品国产一区二区三区| 日韩一本色道免费dvd| 中国美女看黄片| 免费看光身美女| 久久久精品欧美日韩精品| 天堂影院成人在线观看| 国产一区二区三区av在线 | 久久精品91蜜桃| 亚洲欧美中文字幕日韩二区| 99国产极品粉嫩在线观看| 欧美丝袜亚洲另类| 日韩欧美精品v在线| 特大巨黑吊av在线直播| 91麻豆精品激情在线观看国产| 欧美一级a爱片免费观看看| 成年女人永久免费观看视频| 欧美高清成人免费视频www| 日本免费a在线| 在线观看66精品国产| 最后的刺客免费高清国语| 人妻夜夜爽99麻豆av| 国产视频内射| 久久久成人免费电影| 亚洲欧美日韩高清在线视频| 18禁在线播放成人免费| 人妻制服诱惑在线中文字幕| 99riav亚洲国产免费| 大型黄色视频在线免费观看| 性欧美人与动物交配| 成人国产麻豆网| 丝袜喷水一区| 久久精品综合一区二区三区| 内射极品少妇av片p| 好男人视频免费观看在线| 麻豆久久精品国产亚洲av| 51国产日韩欧美| 国产 一区精品| 亚洲不卡免费看| 亚洲国产欧洲综合997久久,| 亚洲国产日韩欧美精品在线观看| 欧美在线一区亚洲| 国产成人精品久久久久久| 99热这里只有是精品在线观看| 久久精品久久久久久噜噜老黄 | 麻豆久久精品国产亚洲av| 韩国av在线不卡| 国产高清激情床上av| 国产高清三级在线| 综合色av麻豆| 国语自产精品视频在线第100页| 亚洲最大成人手机在线| 国产午夜精品久久久久久一区二区三区| 免费黄网站久久成人精品| 国产成年人精品一区二区| 免费av观看视频| 日韩av在线大香蕉| 我的老师免费观看完整版| 99热这里只有是精品在线观看| 久久精品夜夜夜夜夜久久蜜豆| 日韩视频在线欧美| 欧美性感艳星| 欧美日韩一区二区视频在线观看视频在线 | 欧美成人a在线观看| 国产午夜福利久久久久久| 日本黄色片子视频| 亚洲av男天堂| 亚洲成av人片在线播放无| 久久韩国三级中文字幕| 大型黄色视频在线免费观看| 中文字幕免费在线视频6| 国产亚洲91精品色在线| 久久久久久久久久成人| 给我免费播放毛片高清在线观看| 国产三级在线视频| 精品人妻一区二区三区麻豆| 亚洲精品粉嫩美女一区| 岛国毛片在线播放| 乱系列少妇在线播放| 只有这里有精品99| 丝袜喷水一区| 欧美性感艳星| 日本熟妇午夜| 岛国在线免费视频观看| 黄色日韩在线| 一区二区三区免费毛片| 女同久久另类99精品国产91| av卡一久久| 国产精品乱码一区二三区的特点| 女人十人毛片免费观看3o分钟| 岛国毛片在线播放| 日日摸夜夜添夜夜爱| 夫妻性生交免费视频一级片| 中国美女看黄片| 你懂的网址亚洲精品在线观看 | 久久精品夜色国产| 国产三级中文精品| 午夜激情福利司机影院| 波多野结衣高清无吗| 久久99热6这里只有精品| 村上凉子中文字幕在线| 久久久国产成人精品二区| 免费观看a级毛片全部| 午夜精品国产一区二区电影 | 国国产精品蜜臀av免费| 成年版毛片免费区| 欧美xxxx性猛交bbbb| 成人国产麻豆网| 国产伦精品一区二区三区四那| 性插视频无遮挡在线免费观看| av在线蜜桃| 晚上一个人看的免费电影| 国产亚洲91精品色在线| 大香蕉久久网| 乱系列少妇在线播放| 久久精品夜夜夜夜夜久久蜜豆| 欧美变态另类bdsm刘玥| 亚洲av熟女| 欧美一区二区精品小视频在线| 亚洲国产精品久久男人天堂| 欧美性猛交黑人性爽| 网址你懂的国产日韩在线| 51国产日韩欧美| 日韩制服骚丝袜av| 在线国产一区二区在线| 最近手机中文字幕大全| 国产精品无大码| 亚洲在线观看片| 亚洲欧美日韩卡通动漫| ponron亚洲| 国产精品久久电影中文字幕| or卡值多少钱| 你懂的网址亚洲精品在线观看 | 黄色配什么色好看| 国产精品一区www在线观看| 精品国产三级普通话版| 午夜免费男女啪啪视频观看| 中国美女看黄片| 国产精品无大码| 国产色婷婷99| 亚洲aⅴ乱码一区二区在线播放| 久久热精品热| 久久精品国产亚洲av天美| 噜噜噜噜噜久久久久久91| 天堂中文最新版在线下载 | 国产三级在线视频| 亚洲av成人精品一区久久| 校园人妻丝袜中文字幕| 91久久精品电影网| 亚洲国产精品国产精品| 一级二级三级毛片免费看| 永久网站在线| 亚洲不卡免费看| 最近视频中文字幕2019在线8| 亚洲人成网站在线播放欧美日韩| 91精品国产九色| 国产精品三级大全| 一级黄色大片毛片| 啦啦啦啦在线视频资源| 我的女老师完整版在线观看| 精品少妇黑人巨大在线播放 | 日韩欧美在线乱码| 国产成人影院久久av| 久久99蜜桃精品久久| 天堂影院成人在线观看| 岛国在线免费视频观看| 亚洲中文字幕一区二区三区有码在线看|