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

    非飽和帶柴油入滲實驗研究及HSSM模擬

    2014-05-13 03:04:08宋興龍趙勇勝李璐璐秦傳玉杜鵬程吉林大學地下水資源與環(huán)境教育部重點實驗室吉林長春130021
    中國環(huán)境科學 2014年7期
    關(guān)鍵詞:中砂非飽和濕潤

    宋興龍,趙勇勝,李璐璐,秦傳玉,杜鵬程 (吉林大學地下水資源與環(huán)境教育部重點實驗室,吉林 長春 130021)

    非飽和帶柴油入滲實驗研究及HSSM模擬

    宋興龍,趙勇勝*,李璐璐,秦傳玉,杜鵬程 (吉林大學地下水資源與環(huán)境教育部重點實驗室,吉林 長春 130021)

    通過土柱模擬實驗研究了柴油在不同含水率、不同粒徑砂土中的入滲及殘留特征,并采用HSSM(Hydrocarbon Spill Screening Model)模擬了柴油在含水率為 6%中砂介質(zhì)中的入滲過程.結(jié)果表明:隨著介質(zhì)含水率增加,柴油在介質(zhì)中濕潤鋒推進速度先增大后減小,最快速度對應含水率處于相應介質(zhì)最大殘余含水率40%~50%范圍內(nèi),殘余柴油量隨介質(zhì)含水率增大而減小.在細砂、中砂和粗砂3種介質(zhì),柴油入滲平均濕潤鋒推進速度分別為0.42,0.52,0.73cm/min;平均殘余柴油量分別為98.10,68.70,48.79mL. 濕潤鋒推進速度及殘余柴油量均與介質(zhì)粒徑呈負相關(guān).HSSM 擬合柴油在含水率為 6%中砂介質(zhì)中的濕潤鋒推進速度為 0.5832cm/min,與實驗值 0.5689cm/min相比,相對誤差為2.51%.HSSM能較好的模擬柴油在非飽和帶入滲過程,對于土壤及地下水污染預報具有重要意義.

    柴油;非飽和帶;入滲;模擬

    在石油開采、煉制、運輸和使用過程中,難免會產(chǎn)生溢出和泄漏事故[1].石油類污染物泄漏后,會對土壤造成污染,嚴重時還會污染地下水,給人體健康帶來很大危害[2].研究石油類污染物在非飽和帶中的運移機制,選用合適、高效、便捷的模型來預測污染物在非飽和帶的運移速度、污染范圍及到達地下水的時間,對土壤、地下水污染預報及實際修復具有重要意義.

    柴油在非飽和帶中入滲屬于復雜的多孔介質(zhì)、多相流問題[3-5],國外研究者通過實驗及理論研究,先后建立了許多多相流運移模型.根據(jù)復雜程度,這些模型可以分為2類[2],一類是復雜模型,主要包括MISER[6]、STOMP[7-8]、TOUGH-2[9-10]、NAPL simulator[11]等.另一類是篩選模型,主要有API-LNAST[12]、OILSFSM[13]、HSSM[14-16]等.

    復雜模型往往需要大量的物理化學參數(shù),使用起來較為困難;而篩選模型將復雜的多相流問題簡化,需要的參數(shù)較少且能獲取污染物在環(huán)境中運移的主要信息[2].其中,HSSM 不僅能模擬LNAPL泄漏過程,還能模擬 DNAPL泄漏過程

    [14-16],且軟件成熟,使用方便,具有較好的應用前景.劉玉蘭[2]等基于HSSM建立了NAPL泄漏事故場地地下水污染風險快速評估與決策方法.Yoon[17]等曾將 HSSM 用于識別北美地區(qū)LNAPLs泄漏導致的地下水土污染,并取得了很好的效果.實驗室條件下,關(guān)于HSSM模擬效果的研究未見報道.

    雖然目前建立了許多多相流模型來研究柴油在非飽和帶中的運移,但是含水率對柴油在非飽和帶運移的影響研究較少,且研究結(jié)果存在差異[1,18-20],國內(nèi)對于篩選模型的研究和運用也鮮見報道[2].本文采用模擬柱實驗研究介質(zhì)含水率及粒徑對柴油在非飽和帶入滲和殘留特征的影響規(guī)律,并運用HSSM對其遷移過程進行模擬.

    1 材料與方法

    1.1 實驗材料

    實驗用柴油為-10#柴油,購自中國石油天然氣股份有限公司吉林長春銷售分公司.實驗所用非飽和帶介質(zhì)為石英砂,經(jīng)過篩分為細砂、中砂、粗砂3種,實驗用砂物理參數(shù)見表1.

    表1 實驗用砂物理參數(shù)Table 1 Properties of sands used in experiments

    1.2 實驗裝置

    實驗所用有機玻璃模擬柱內(nèi)徑為 5cm,長為60cm,實驗裝置如圖1所示.

    圖1 實驗裝置示意Fig.1 Experimental setup

    1.3 實驗方法

    為研究介質(zhì)含水率及粒徑對柴油在非飽和帶入滲的影響,分別設置了細砂含水率為 0%(烘干砂)、2%、4%、6%、8%、10%,中砂含水率為0%、2%、4%、6%、8%、10%和粗砂含水率為0%、1%、2%、3%、4%、5%、6%、7%共20組實驗.將配置好含水率的石英砂裝入模擬柱中,每裝填2cm夯實一次.石英砂裝填高度為54cm,利用蠕動泵向柱中連續(xù)滴入已用蘇丹紅(300mg/L)染色的柴油,柴油流量為2mL/min.

    在滲漏開始后的60min內(nèi),每隔5min記錄柴油濕潤鋒深度,滲漏60min后每隔10mim記錄濕潤鋒深度.當濕潤鋒面到達模擬柱底部時,停止供油.靜置24h,收集從模擬柱底部流出的柴油.模擬柱頂部進入的柴油量減去底部收集的柴油量即為殘余柴油量.比較不同介質(zhì)粒徑及含水率條件下柴油遷移速度及殘余量,分析介質(zhì)粒徑及含水率對柴油入滲特征的影響.

    為研究HSSM對柴油在非飽和帶遷移的模擬效果,HSSM模擬了柴油在含水率為6%中砂介質(zhì)入滲過程,并比較了模擬結(jié)果與實測值的差異,分析了可能的誤差來源.

    2 結(jié)果與分析

    2.1 含水率對柴油入滲特征的影響

    如圖2所示,隨著時間的推移,柴油在砂性介質(zhì)中的入滲深度不斷增加,入滲深度與時間基本呈線性關(guān)系,這與張博聞[18]、WANG[21]等的研究結(jié)果冪函數(shù)關(guān)系不同.這可能是由于當介質(zhì)含水率較低時,介質(zhì)只有表面和少數(shù)孔隙被水分占據(jù),剩余孔隙通道較多,柴油下滲驅(qū)替水分遇到的阻力較小,主要受到重力作用和介質(zhì)對柴油的吸附作用,下滲速度相對穩(wěn)定[19-20].

    圖2 柴油在不同含水率介質(zhì)中入滲過程Fig.2 The infiltration process of diesel in sandy soil with different soil moistures

    入滲深度與入滲時間擬合方程的斜率K代表單位時間內(nèi)濕潤鋒前進的距離,即濕潤鋒推進速度.濕潤鋒推進速度隨含水率的變化如圖3所示.

    表2 柴油入滲深度與入滲時間線性擬合方程Table 2 The functional relationship between infiltration depth and time

    圖3 濕潤鋒推進速度隨介質(zhì)含水率的變化Fig.3 The effect of soil moisture on migration velocity of soil wetting front

    從圖 3可以看出,柴油在粗砂中入滲時,介質(zhì)含水率從0%提高到4%,加快了濕潤鋒推進速度,介質(zhì)含水率從 4%提高到 7%,降低了濕潤鋒推進速度;柴油在中砂和細砂中入滲也具有相似規(guī)律,分別在介質(zhì)含水率為6%和8%時出現(xiàn)最大濕潤鋒推進速度.最大濕潤鋒推進速度對應含水率均處于相應介質(zhì)最大殘余含水率 40%~50%之間.這主要是由于,介質(zhì)中的水一方面占據(jù)了介質(zhì)中的部分孔隙,減小了介質(zhì)的孔隙度,對柴油入滲有阻礙作用;另一方面,水分子可以覆蓋介質(zhì)表面,和柴油競爭在介質(zhì)顆粒中的極性吸附點位,降低了介質(zhì)對柴油的吸附作用,對柴油入滲有促進作用[18-20].因此,含水率對柴油在非飽和帶濕潤鋒推進速度的最終影響是這兩種因素共同作用的結(jié)果.

    圖4 殘余柴油量隨介質(zhì)含水率變化Fig.4 The effect of soil moisture on the amount of residual diesel

    由圖4可以看出,柴油在細砂中入滲時,當含水率從 0%提高到 10%,殘余柴油量從 142.7mL下降到59.3mL,柴油在中砂和粗砂中入滲也具有相似規(guī)律.這是因為,含水率增加導致水分子逐漸在介質(zhì)表面形成水膜,使得柴油吸附位點減少,介質(zhì)對柴油的吸附作用減弱,殘留在介質(zhì)中的柴油量逐漸減少.

    2.2 介質(zhì)粒徑對柴油入滲特征的影響

    如圖3所示,柴油在砂性介質(zhì)中入滲時,隨著介質(zhì)粒徑增大,濕潤鋒推進速度增大.細砂、中砂和粗砂中平均濕潤鋒推進速率分別為 0.42,0.52, 0.73cm/min.這是因為介質(zhì)粒徑越大,孔隙度越大,柴油在下滲過程中的通道越大,所受阻力越小,濕潤鋒推進越快.李永濤[22]研究表明柴油在粗砂中遷移速度快于中砂,與本實驗結(jié)果相似.

    如圖4所示,柴油在砂性介質(zhì)中入滲時,隨著介質(zhì)粒徑增大,殘余柴油量下降. 細砂、中砂和粗砂中平均柴油殘余量分別為 98.10,68.70, 48.79mL.這是由于介質(zhì)粒徑越大,比表面積越小,可供柴油吸附的位點越少.這與李永霞[20]的研究結(jié)果相類似.

    3 HSSM對柴油在非飽和帶入滲過程的模擬

    3.1 HSSM原理介紹

    Charbeneau等[15-16]在1995年建立了HSSM (Hydrocarbon Spill Screening Model),用于預測LNAPL泄漏導致的土壤和地下水污染.Willson等[23]在2006年完善了HSSM,使得HSSM也能用于預測DNAPL泄漏導致的污染.

    HSSM將LNAPL在地下環(huán)境中的運移過程概化為3個階段,并分別由3個模塊實現(xiàn).第1階段為 LNAPL從地表面垂直下滲到毛細帶邊緣,由KOPT(Kinematic Oily Pollutant Transport)模塊模擬;第2階段為LNAPL在毛細帶邊緣擴展形成透鏡體并向地下水中溶解,由OILENS模塊模擬;第 3階段為溶出的污染組分隨地下水流遷移至特定受體,由 TSGPLUME(Transient Source Gaussian Plume Model)模塊模擬[16].

    HSSM 模擬柴油在非飽和入滲,只需用到KOPT模塊.KOPT模塊主要是基于達西定律、相對滲透率-飽和度-毛細壓力(K-S-P)關(guān)系及Green-Ampt方程,該模塊可以計算LNAPL在非飽和帶中的運移速度、到達含水層的時間,以及LNAPL在非飽和帶中的飽和度變化[15].KOPT模塊主要控制方程如下:

    式中:Vi為單位水力梯度下液相i的滲流速度,m/s; Kei為液相i的有效滲透系數(shù),m/s;Ksi為液相i的飽和滲透系數(shù),m/s;kri為液相 i的相對滲透率,無量綱;η為孔隙度,無量綱;So為NAPL飽和度,無量綱;t為時間,s;Keo為 NAPL的有效滲透系數(shù),m/s;Sw(avg)為年均補給量下的水飽和度,無量綱;Z為入滲深度,m; Hf為NAPL濕潤鋒面的毛細壓力,Ν;h為毛細水頭,m; Ks為NAPL飽和滲透系數(shù),m/s;Kro為NAPL相對滲透率,無量綱.

    柴油在非飽和帶入滲過程模擬不涉及到OILENS模塊和TSGPLUME模塊,在此不展開討論.

    3.2 HSSM模擬結(jié)果分析

    3.2.1 模型參數(shù) 以柴油在含水率為 6%中砂介質(zhì)入滲過程模擬為例,模型模擬需要的主要參數(shù)如表3所示.

    表3 HSSM主要參數(shù)Table 3 Main parameters of HSSM

    圖5 入滲深度與入滲時間關(guān)系HSSM模擬結(jié)果Fig.5 HSSM simulation results of infiltration depth with infiltration time

    3.2.2 模擬結(jié)果分析 入滲深度隨時間變化關(guān)系實驗結(jié)果及HSSM模擬結(jié)果如圖5所示.

    從圖 5可知,當入滲時間小于 70min時, HSSM 模擬結(jié)果與實驗值基本一致.當入滲時間大于70min后,HSSM模擬的入滲深度稍大于實驗值.HSSM 模擬的濕潤鋒推進速度0.5832cm/min,實驗值為 0.5689cm/min,相對誤差為 2.51%.HSSM 模擬的柴油濕潤鋒面到達模擬柱底部時間 92min,實驗值為 96min,相對誤差為4.17%.HSSM 模擬值與實驗值之間雖然有一定的誤差,但是誤差較小.這可能是由于模型對入滲過程簡化以及相關(guān)參數(shù)測量誤差引起的.

    4 結(jié)論

    4.1 柴油在砂性介質(zhì)中入滲時,入滲深度隨時間變化可以用線性方程擬合.隨著介質(zhì)含水率上升,濕潤鋒推進速率先增加后降低,最大濕潤鋒推進速率對應的含水率處于相應介質(zhì)最大殘余含水率的 40%~50%范圍內(nèi).隨著介質(zhì)含水率上升,殘余柴油量下降.

    4.2 柴油在砂性介質(zhì)中入滲時,隨著介質(zhì)粒徑增大,濕潤鋒推進速度增大,而殘余柴油量下降.

    4.3 HSSM 對柴油在非飽和帶中入滲過程模擬的結(jié)果較為理想,HSSM 擬合的柴油在含水率為6%中砂介質(zhì)中的濕潤鋒推進速度與實驗值相比,相對誤差僅為2.51%.

    [1] 李永霞,鄭西來,馬艷飛.土壤非飽和帶柴油入滲的試驗研究 [J].環(huán)境污染與防治, 2011,33(5):27-31.

    [2] 劉玉蘭,程莉蓉,丁愛中,等.NAPL泄漏事故場地地下水污染風險快速評估與決策 [J]. 中國環(huán)境科學, 2011,31(7):1219-1224.

    [3] 薛 強.石油污染物在地下環(huán)境系統(tǒng)中運移的多相流模型研究[D]. 沈陽:遼寧工程技術(shù)大學, 2003.

    [4] 王 麗.輕非水相液體(LNAPL)污染土壤多相流實驗研究 [D].西安:長安大學, 2009.

    [5] 邢巍巍.LNAPLs在土體中的運移特性研究 [D]. 北京:清華大學, 2005.

    [6] Abriola L M, Lang J, Rahtfelder K. Michigan soil vapor extraction remediation (MISER) Model-A computer program to model bioventing of organic chemicals in unsaturated geological material [R]. EPA, Cincinnati, OH, USA, EPA/600/R-97/099,1997.

    [7] White M D, Oostrom M. STOMP, subsurface transport over multiple phases, version 2.0 theory guide [R]. Richland, WA, USA, PNNL-12030Pacific North west National Lab, 2000.

    [8] White M D, Oostrom M. STOMP, subsurface transport over multiplephases, version3.1user’sguide [R]. Richland, WA, USA, PNNL-14478PacificNorthwestNationalLab, 2004.

    [9] Pruess K, Battistelli A. TMVOC, a numerical simulator for there-phase non-isothermal flows of multicomponent hydrocarbon mixtures in saturated-unsaturated heterogeneous media [R]. Berkeley, CA, USA, LBNL-49375Lawrence Berkeley National Lab, 2002.

    [10] Pruessl K, Oldenburg C, MoirdisG. TOUGH2user’s guide, version 2.0 [R]. Berkeley, CA, USA, LBNL-43134 Lawrence Brekeley NationalLab, 1999.

    [11] Guarnaccia J F, Pinder G F, Fishman M. NAPL simulator documentation final report [R]. USA, EPA Cooperative AgreementNo.CR-820499, 1997.

    [12] Huntley D, Beckett G. D. Persistence of LNAPL sources: relationship between risk reduction and LNAPL recovery [J]. Journal of Contaminant Hydrology, 2002,59:3-26.

    [13] Hussein M, Jin M, Weaver J W. Development and verification of a screening model for surface spreading of petroleum [J]. Journal of Contaminant Hydrology, 2002,57:281-302.

    [14] Weaver J W, Charbeneau R J, Lien B K. A screening model for nonaqueous phase liquid transport in the vadose zone using Green-Ampt and kinematic wave theory [R]. Water Resources Research, 1994,30:93-105.

    [15] Weaver J W, Chrabeneau R J, Tauxe J D, et al. The Hydrocarbon Spill Screening Model (HSSM) Volume 1: User's Guide [R]. EPA/600/R-94/039a, 1994.

    [16] Charbeneau R J, Weaver J W, Lien B K, et a1. The Hydrocarbon Spill Screening Model (HSSM) Volume 2: Theoretical Background and Source Codes [R]. EPA/600/R-94/039b, 1995.

    [17] Yoon H, werth C J, Barkan C P. L, et al. An environmental screening model to assess the consequences to soil and groundwater from railroad-tank-car spills of light non-aqueous phase liquids [J]. Journal of Hazardous Materials, 2009,165:332-344.

    [18] 張博聞,邵明安.初始含水率對土壤中原油入滲的影響 [J]. 農(nóng)業(yè)工程學報, 2010,26(3):9-13.

    [19] 童 玲.石油污染含水介質(zhì)的水理和水力學特征研究 [D]. 青島:中國海洋大學, 2008.

    [20] 李永霞.含水介質(zhì)中石油污染物遷移與殘留特征 [D]. 青島:中國海洋大學, 2011.

    [21] Wang Y, Shao M. Infiltration characteristics of non-aqueous phase liquids in undisturbed loessal soil cores [J]. Journal of Environmental Sciences, 2009,21(10):1424-1431.

    [22] 李永濤.LNAPLs在包氣帶中運移機理及模擬研究 [D]. 西安:長安大學, 2010.

    [23] Willson C S, Weaver J W, Charbeneau R J. A screening model for simulating DNAPL flow and transport in porous media: theoretical development [J]. Environmental Modeling and Software, 2006,21:l6-32.

    《中國環(huán)境科學》被Ei收錄

    根據(jù) Ei總部 2013年頒布的期刊收錄情況,《中國環(huán)境科學》已被 Ei數(shù)據(jù)庫作為源期刊收錄,詳見http://www.chinaeidata.com/periodical.htm

    《中國環(huán)境科學》編輯部

    2013-03-14

    Experimental study and HSSM simulation of diesel infiltration in unsaturated zone.

    SONG Xing-long, ZHAO

    Yong-sheng*, LI Lu-lu, QIN Chuan-yu, DU Peng-cheng (Key Laboratory of Groundwater Resources and Environment, Ministry of Education, Jilin University, Changchun 130021, China). China Environmental Science, 2014,34(7):1818~1823

    This paper presents the results of laboratory investigation and simulation of diesel infiltration in unsaturated zone. A series of one-dimensional column experiments were conducted to study the infiltration and residue characteristics of diesel in sandy soils under varied soil moistures and media sizes. Then, Hydrocarbon Spill Screening Model (HSSM) was used to simulate the infiltration processes of diesel in medium sand with soil moisture equals to 6%. The results indicated that as soil moisture grew, the migration velocity of soil wetting front first increased and then decreased, and the fastest migration velocity corresponded to 40%~50% of maximum residue moisture content. Negative correlation was found between the amount of residual diesel and soil moisture content. In fine, medium and coarse sand, the average migration velocity of soil wetting front were 0.42, 0.52, 0.73cm/min, and the average amount of residual diesel were 98.10, 68.70, 48.79mL, respectively. The migration velocity of soil wetting front and the amount of residual diesel were negatively correlated with media size. Fitted migration velocity by HSSM of soil wetting front in medium sand with moisture content equals to 6% was 0.5832cm/min, while the experimental value was 0.5689cm/min. The relative percentage error between fitted value and the experimental value was 2.51%. HSSM can accurately simulate the infiltration process of diesel oil in the unsaturated zone, and has great significance to soil and groundwater pollution forecast.

    diesel oil;unsaturated zone;infiltration;simulation

    X703.1

    A

    1000-6923(2014)07-1818-06

    宋興龍(1990-),男,湖南湘潭人,吉林大學碩士研究生,主要從事污染場地控制與修復研究.

    2013-08-13

    中國地質(zhì)調(diào)查局項目(1212011220987);北京市教育委員會市屬高校創(chuàng)新能力提升計劃項目----松散地層地下水污染控制與修復技術(shù)研究

    * 責任作者, 教授, zhaoyongsheng@jlu.edu.cn

    猜你喜歡
    中砂非飽和濕潤
    內(nèi)河航道整治護岸工程中砂枕護腳施工探討
    中國水運(2022年4期)2022-04-27 22:29:17
    利用建筑垃圾制備預拌抹灰砂漿及其性能研究
    散裝水泥(2021年6期)2021-12-30 09:47:56
    The Desert Problem
    溫度和含水率對中砂和粗砂導熱系數(shù)的影響
    非飽和原狀黃土結(jié)構(gòu)強度的試驗研究
    工程與建設(2019年1期)2019-09-03 01:12:24
    關(guān)于山石料廠廢料在混凝土中的應用
    商品混凝土(2019年7期)2019-08-06 02:38:38
    海邊的沙漠
    非飽和多孔介質(zhì)應力滲流耦合分析研究
    非飽和土基坑剛性擋墻抗傾覆設計與參數(shù)分析
    非飽和地基土蠕變特性試驗研究
    1024香蕉在线观看| 99国产综合亚洲精品| 三级毛片av免费| 久久精品国产99精品国产亚洲性色| www.精华液| 国产亚洲精品av在线| 麻豆成人av在线观看| 日韩精品青青久久久久久| 精品国产美女av久久久久小说| 亚洲国产欧美人成| 成人av在线播放网站| 亚洲欧美精品综合久久99| 中文字幕高清在线视频| 一级黄色大片毛片| 久久精品aⅴ一区二区三区四区| 19禁男女啪啪无遮挡网站| 国产精品av久久久久免费| 亚洲国产精品成人综合色| 欧美色欧美亚洲另类二区| 国内精品久久久久精免费| 国产精品一区二区精品视频观看| 精品国产三级普通话版| 欧美成人性av电影在线观看| 免费看日本二区| 国产成人精品久久二区二区免费| 又粗又爽又猛毛片免费看| 久久这里只有精品中国| 亚洲国产精品成人综合色| 欧美3d第一页| 国产午夜精品论理片| 午夜福利高清视频| 日韩三级视频一区二区三区| 非洲黑人性xxxx精品又粗又长| 亚洲 欧美一区二区三区| 亚洲精华国产精华精| 色视频www国产| 中文字幕av在线有码专区| 两个人的视频大全免费| 亚洲av日韩精品久久久久久密| 精品欧美国产一区二区三| 国内久久婷婷六月综合欲色啪| 国产毛片a区久久久久| 亚洲国产精品999在线| 国产精品精品国产色婷婷| x7x7x7水蜜桃| 91九色精品人成在线观看| 18禁观看日本| 最近视频中文字幕2019在线8| bbb黄色大片| 99久久精品一区二区三区| 成人av一区二区三区在线看| 午夜福利免费观看在线| 99久久精品国产亚洲精品| 国产1区2区3区精品| av在线蜜桃| 色吧在线观看| 久久亚洲真实| 99国产精品一区二区蜜桃av| 人妻久久中文字幕网| 啦啦啦观看免费观看视频高清| av黄色大香蕉| 欧美日韩瑟瑟在线播放| 中文字幕人成人乱码亚洲影| 我要搜黄色片| 老司机午夜十八禁免费视频| 啦啦啦观看免费观看视频高清| 国产精品久久久久久久电影 | 成年免费大片在线观看| 亚洲人成网站高清观看| 国产精品1区2区在线观看.| 精华霜和精华液先用哪个| 国产三级黄色录像| 十八禁人妻一区二区| 国产69精品久久久久777片 | 国产高清激情床上av| 国产高清有码在线观看视频| 亚洲欧美激情综合另类| 天堂av国产一区二区熟女人妻| 午夜两性在线视频| 免费av毛片视频| 精品国产乱码久久久久久男人| 亚洲色图 男人天堂 中文字幕| 成年女人毛片免费观看观看9| 欧美最黄视频在线播放免费| 一级黄色大片毛片| 欧美午夜高清在线| 男女视频在线观看网站免费| 级片在线观看| 欧美日韩综合久久久久久 | 噜噜噜噜噜久久久久久91| 久久久久国内视频| 伦理电影免费视频| 婷婷亚洲欧美| 国产一区二区在线观看日韩 | av视频在线观看入口| 成人三级做爰电影| h日本视频在线播放| 日韩精品青青久久久久久| 一二三四社区在线视频社区8| 亚洲av熟女| 国产成人av教育| 黑人欧美特级aaaaaa片| 国产视频一区二区在线看| 亚洲av美国av| 后天国语完整版免费观看| 在线看三级毛片| 欧美一级a爱片免费观看看| 亚洲色图av天堂| 久久人人精品亚洲av| 中文字幕人妻丝袜一区二区| 久久九九热精品免费| 国产成人一区二区三区免费视频网站| e午夜精品久久久久久久| 国产1区2区3区精品| 国产一区二区三区视频了| 国产又色又爽无遮挡免费看| 日本免费a在线| 国产成人av教育| 日韩大尺度精品在线看网址| 999精品在线视频| 国产免费av片在线观看野外av| 国产高清有码在线观看视频| 欧美中文综合在线视频| 免费在线观看影片大全网站| 九色成人免费人妻av| 亚洲精品久久国产高清桃花| 亚洲avbb在线观看| 又紧又爽又黄一区二区| 最新中文字幕久久久久 | 欧美日韩福利视频一区二区| 亚洲 欧美一区二区三区| 欧美又色又爽又黄视频| 国产一区二区三区在线臀色熟女| 亚洲人成电影免费在线| 亚洲人成网站高清观看| 国产99白浆流出| 一进一出好大好爽视频| 国产精华一区二区三区| 国产人伦9x9x在线观看| 三级男女做爰猛烈吃奶摸视频| 村上凉子中文字幕在线| 舔av片在线| 久久草成人影院| 香蕉久久夜色| 少妇人妻一区二区三区视频| 成人午夜高清在线视频| 色视频www国产| 亚洲成av人片在线播放无| 三级国产精品欧美在线观看 | 在线视频色国产色| 露出奶头的视频| 18禁黄网站禁片午夜丰满| 精品国内亚洲2022精品成人| 1000部很黄的大片| 亚洲中文av在线| 一级作爱视频免费观看| 国产激情欧美一区二区| 午夜福利在线观看免费完整高清在 | 国产精品野战在线观看| 又粗又爽又猛毛片免费看| 国产伦精品一区二区三区视频9 | 99热只有精品国产| 国产亚洲精品久久久com| 久久久久精品国产欧美久久久| 欧美激情久久久久久爽电影| 亚洲美女视频黄频| 啪啪无遮挡十八禁网站| 日本在线视频免费播放| 久久香蕉精品热| 精品国产美女av久久久久小说| 欧美黑人巨大hd| 欧美3d第一页| 精品福利观看| 日韩欧美一区二区三区在线观看| 真人一进一出gif抽搐免费| 亚洲欧美日韩高清在线视频| 12—13女人毛片做爰片一| 岛国视频午夜一区免费看| 男插女下体视频免费在线播放| 美女高潮的动态| 99久久精品一区二区三区| 婷婷精品国产亚洲av| avwww免费| 69av精品久久久久久| 久久久国产成人精品二区| 一本综合久久免费| 午夜福利欧美成人| 亚洲av熟女| 久久久久久久久免费视频了| 麻豆久久精品国产亚洲av| 亚洲熟妇熟女久久| 久久久久性生活片| 欧美日韩瑟瑟在线播放| 在线免费观看的www视频| 18禁裸乳无遮挡免费网站照片| 欧美大码av| 婷婷精品国产亚洲av| 99久久精品热视频| 成人av在线播放网站| 99久久精品一区二区三区| av国产免费在线观看| 日日干狠狠操夜夜爽| 老司机午夜福利在线观看视频| 国产精品野战在线观看| 久久国产乱子伦精品免费另类| 听说在线观看完整版免费高清| 在线看三级毛片| 国产精品99久久久久久久久| 国产免费av片在线观看野外av| 精品久久久久久久末码| 亚洲一区高清亚洲精品| 小说图片视频综合网站| 黄片小视频在线播放| 亚洲成av人片免费观看| 成人欧美大片| 欧美日本亚洲视频在线播放| 色播亚洲综合网| 国产v大片淫在线免费观看| 一区二区三区激情视频| 亚洲av日韩精品久久久久久密| 精品一区二区三区四区五区乱码| 欧美日韩国产亚洲二区| 欧美在线黄色| 亚洲五月婷婷丁香| 看免费av毛片| 日本免费一区二区三区高清不卡| 岛国在线免费视频观看| 久久人妻av系列| 动漫黄色视频在线观看| 麻豆国产97在线/欧美| xxxwww97欧美| 成年版毛片免费区| 波多野结衣高清作品| 国产淫片久久久久久久久 | 亚洲片人在线观看| avwww免费| 欧美日韩乱码在线| 脱女人内裤的视频| 日韩欧美精品v在线| 黄频高清免费视频| 精品一区二区三区av网在线观看| 午夜福利高清视频| 国产精品精品国产色婷婷| 偷拍熟女少妇极品色| 免费搜索国产男女视频| 久9热在线精品视频| 制服人妻中文乱码| av在线蜜桃| 成年女人永久免费观看视频| 国产97色在线日韩免费| 亚洲av五月六月丁香网| 国产成年人精品一区二区| 真人一进一出gif抽搐免费| 欧美最黄视频在线播放免费| 国产精品久久久av美女十八| 亚洲精品久久国产高清桃花| 亚洲av成人一区二区三| 亚洲国产精品sss在线观看| 久久亚洲真实| 一级毛片女人18水好多| www日本黄色视频网| 久久精品91无色码中文字幕| 日日夜夜操网爽| 午夜免费成人在线视频| 色精品久久人妻99蜜桃| 国产1区2区3区精品| 男女午夜视频在线观看| 国产成人影院久久av| 亚洲第一欧美日韩一区二区三区| 在线a可以看的网站| 五月玫瑰六月丁香| 国产成人av教育| 少妇的逼水好多| 午夜福利在线观看免费完整高清在 | 嫁个100分男人电影在线观看| 免费在线观看日本一区| 亚洲国产精品sss在线观看| 性色avwww在线观看| 窝窝影院91人妻| 色哟哟哟哟哟哟| 天天添夜夜摸| 国产精品永久免费网站| 午夜视频精品福利| 一个人观看的视频www高清免费观看 | 一级毛片女人18水好多| 久久精品人妻少妇| 久久亚洲真实| 色综合婷婷激情| 亚洲片人在线观看| 真人一进一出gif抽搐免费| 欧美乱码精品一区二区三区| 大型黄色视频在线免费观看| 99国产精品一区二区蜜桃av| 亚洲中文日韩欧美视频| 国产精品av视频在线免费观看| 十八禁人妻一区二区| 成人精品一区二区免费| 国产精品久久久av美女十八| 少妇裸体淫交视频免费看高清| 成人鲁丝片一二三区免费| 成人国产综合亚洲| 久久久久久久久久黄片| 桃色一区二区三区在线观看| 国产精品99久久99久久久不卡| 免费高清视频大片| 色综合欧美亚洲国产小说| 一区二区三区激情视频| 国产私拍福利视频在线观看| 成人性生交大片免费视频hd| 国产伦在线观看视频一区| 国产伦精品一区二区三区视频9 | 18美女黄网站色大片免费观看| 免费在线观看日本一区| 亚洲国产精品sss在线观看| 国产亚洲精品久久久com| 在线国产一区二区在线| 男人舔女人下体高潮全视频| 精品熟女少妇八av免费久了| 成年版毛片免费区| 九九热线精品视视频播放| 欧美极品一区二区三区四区| 精品国产亚洲在线| АⅤ资源中文在线天堂| 窝窝影院91人妻| 国产主播在线观看一区二区| 好男人在线观看高清免费视频| 成人永久免费在线观看视频| 免费高清视频大片| 国产久久久一区二区三区| 老汉色∧v一级毛片| 亚洲国产精品成人综合色| 国产激情偷乱视频一区二区| 少妇丰满av| 欧美色视频一区免费| 欧美乱码精品一区二区三区| 黄色视频,在线免费观看| 最近最新中文字幕大全电影3| 波多野结衣高清无吗| 99久久无色码亚洲精品果冻| 亚洲av中文字字幕乱码综合| 99热只有精品国产| 特大巨黑吊av在线直播| 女同久久另类99精品国产91| 亚洲av日韩精品久久久久久密| 午夜精品一区二区三区免费看| 在线视频色国产色| 国产真实乱freesex| 日韩精品中文字幕看吧| 淫秽高清视频在线观看| 制服人妻中文乱码| 国产一区在线观看成人免费| 成年女人毛片免费观看观看9| 亚洲中文字幕一区二区三区有码在线看 | 国产淫片久久久久久久久 | 麻豆国产av国片精品| 岛国在线观看网站| 国产一区二区三区在线臀色熟女| 欧美一级a爱片免费观看看| 美女免费视频网站| 日本撒尿小便嘘嘘汇集6| 伊人久久大香线蕉亚洲五| 性色avwww在线观看| 搞女人的毛片| 国产黄a三级三级三级人| 99热这里只有精品一区 | 日本免费a在线| 淫妇啪啪啪对白视频| 亚洲精品一区av在线观看| 91久久精品国产一区二区成人 | 日本一本二区三区精品| 少妇丰满av| 黄色片一级片一级黄色片| 免费av不卡在线播放| 亚洲av日韩精品久久久久久密| www.精华液| 国产男靠女视频免费网站| 最新中文字幕久久久久 | 国产aⅴ精品一区二区三区波| 大型黄色视频在线免费观看| 久久久色成人| 国产精品99久久99久久久不卡| 日韩三级视频一区二区三区| 最新在线观看一区二区三区| 欧美一区二区精品小视频在线| 欧美在线一区亚洲| 亚洲专区国产一区二区| 久久中文字幕一级| 免费一级毛片在线播放高清视频| 亚洲av美国av| 伦理电影免费视频| 婷婷精品国产亚洲av在线| 午夜激情欧美在线| 国产99白浆流出| 欧美zozozo另类| 国产1区2区3区精品| 香蕉av资源在线| e午夜精品久久久久久久| 亚洲熟女毛片儿| 日韩av在线大香蕉| 99久久99久久久精品蜜桃| 亚洲av美国av| www.熟女人妻精品国产| 欧美在线黄色| 麻豆国产av国片精品| 亚洲一区二区三区色噜噜| 欧美在线一区亚洲| 久久精品亚洲精品国产色婷小说| 天堂动漫精品| 国产人伦9x9x在线观看| 久久精品夜夜夜夜夜久久蜜豆| 可以在线观看的亚洲视频| 国产视频一区二区在线看| 神马国产精品三级电影在线观看| 搡老岳熟女国产| 日韩中文字幕欧美一区二区| 国产精品av视频在线免费观看| 变态另类丝袜制服| 非洲黑人性xxxx精品又粗又长| 麻豆一二三区av精品| 亚洲真实伦在线观看| 90打野战视频偷拍视频| 狂野欧美激情性xxxx| 婷婷精品国产亚洲av| 亚洲18禁久久av| 精品无人区乱码1区二区| 麻豆一二三区av精品| 午夜福利在线在线| а√天堂www在线а√下载| 国产淫片久久久久久久久 | 后天国语完整版免费观看| 国产av麻豆久久久久久久| 2021天堂中文幕一二区在线观| 波多野结衣高清作品| 久久中文字幕一级| 欧美黄色淫秽网站| www.精华液| 九九在线视频观看精品| 日韩欧美三级三区| 在线免费观看不下载黄p国产 | 中文字幕精品亚洲无线码一区| 岛国在线观看网站| 国产免费av片在线观看野外av| 亚洲色图 男人天堂 中文字幕| 国产精品一区二区三区四区免费观看 | 日韩欧美在线乱码| 国产欧美日韩精品亚洲av| 天堂影院成人在线观看| 美女免费视频网站| 黄色丝袜av网址大全| 亚洲美女视频黄频| 午夜成年电影在线免费观看| 久久久国产欧美日韩av| 国产成人av教育| 老熟妇仑乱视频hdxx| 51午夜福利影视在线观看| 国产午夜精品论理片| 淫秽高清视频在线观看| 精品久久久久久久久久免费视频| 久久久久久久午夜电影| 精品国产乱子伦一区二区三区| 美女午夜性视频免费| 日韩高清综合在线| 在线观看免费视频日本深夜| 午夜福利在线观看吧| 成人一区二区视频在线观看| 日韩欧美三级三区| 国产激情久久老熟女| 色吧在线观看| 观看免费一级毛片| 国产精品,欧美在线| 国产精品久久电影中文字幕| 国产一区在线观看成人免费| 亚洲欧美精品综合久久99| 美女大奶头视频| 美女高潮的动态| 一区福利在线观看| 色av中文字幕| 精品国产三级普通话版| 美女扒开内裤让男人捅视频| 亚洲美女视频黄频| 国产成人啪精品午夜网站| 亚洲成人精品中文字幕电影| aaaaa片日本免费| 国产精品久久久久久亚洲av鲁大| 在线免费观看不下载黄p国产 | 男女午夜视频在线观看| 日本a在线网址| 国产激情欧美一区二区| 搡老岳熟女国产| 男插女下体视频免费在线播放| 国产成人精品久久二区二区免费| 国产主播在线观看一区二区| 国产精品美女特级片免费视频播放器 | 巨乳人妻的诱惑在线观看| 成年免费大片在线观看| 神马国产精品三级电影在线观看| 国产真实乱freesex| 亚洲中文字幕一区二区三区有码在线看 | 久9热在线精品视频| 亚洲中文字幕一区二区三区有码在线看 | 九色成人免费人妻av| 三级毛片av免费| 国语自产精品视频在线第100页| 啪啪无遮挡十八禁网站| 99国产精品99久久久久| 黄色女人牲交| 99热6这里只有精品| 麻豆成人av在线观看| 露出奶头的视频| 丰满人妻熟妇乱又伦精品不卡| 动漫黄色视频在线观看| 亚洲av第一区精品v没综合| 日本在线视频免费播放| 亚洲精品一区av在线观看| 久久性视频一级片| 免费在线观看日本一区| 国内揄拍国产精品人妻在线| 国产美女午夜福利| 一级黄色大片毛片| 最近视频中文字幕2019在线8| 中文在线观看免费www的网站| 亚洲 国产 在线| 日日干狠狠操夜夜爽| 天堂av国产一区二区熟女人妻| 欧美午夜高清在线| 久久久国产成人精品二区| 亚洲精品在线美女| 国产欧美日韩精品亚洲av| 黑人欧美特级aaaaaa片| 九色成人免费人妻av| 看片在线看免费视频| 日韩三级视频一区二区三区| 日日摸夜夜添夜夜添小说| 亚洲国产欧美人成| 亚洲av片天天在线观看| 久久久久久大精品| 亚洲人成网站高清观看| 国产午夜精品久久久久久| 中文字幕久久专区| 麻豆国产97在线/欧美| 人人妻人人看人人澡| 国产一区二区三区在线臀色熟女| 91在线观看av| 欧美日韩乱码在线| av天堂中文字幕网| 一边摸一边抽搐一进一小说| 在线观看免费视频日本深夜| 真人一进一出gif抽搐免费| 国内精品美女久久久久久| 久久久国产成人免费| 精品电影一区二区在线| 在线免费观看不下载黄p国产 | 久久中文看片网| 久久国产精品人妻蜜桃| www.999成人在线观看| 窝窝影院91人妻| 成在线人永久免费视频| 久久久久免费精品人妻一区二区| 国产亚洲精品久久久com| 久久久久久九九精品二区国产| 国产真实乱freesex| 神马国产精品三级电影在线观看| 又爽又黄无遮挡网站| 午夜免费成人在线视频| 亚洲国产色片| 国产精品国产高清国产av| 久久久水蜜桃国产精品网| 在线观看午夜福利视频| 亚洲av片天天在线观看| 久久久水蜜桃国产精品网| 色av中文字幕| 在线观看日韩欧美| av福利片在线观看| 国产黄色小视频在线观看| 国产探花在线观看一区二区| 亚洲精品粉嫩美女一区| 国产aⅴ精品一区二区三区波| av欧美777| 99久国产av精品| 久久人人精品亚洲av| 国产成人一区二区三区免费视频网站| 国产一级毛片七仙女欲春2| 18美女黄网站色大片免费观看| 手机成人av网站| 午夜精品在线福利| 人人妻,人人澡人人爽秒播| 在线免费观看不下载黄p国产 | 欧美日本亚洲视频在线播放| svipshipincom国产片| 男人舔女人下体高潮全视频| 国产精品野战在线观看| 亚洲,欧美精品.| 亚洲va日本ⅴa欧美va伊人久久| 噜噜噜噜噜久久久久久91| 国产v大片淫在线免费观看| 日本五十路高清| 久久精品国产亚洲av香蕉五月| 久久热在线av| 亚洲av日韩精品久久久久久密| 久久午夜综合久久蜜桃| 成人国产综合亚洲| 90打野战视频偷拍视频| 欧美日韩国产亚洲二区| 色精品久久人妻99蜜桃| 在线视频色国产色| 蜜桃久久精品国产亚洲av| 亚洲熟妇熟女久久| 脱女人内裤的视频| 国产主播在线观看一区二区| 看免费av毛片| 99re在线观看精品视频| 久久热在线av| 在线国产一区二区在线| 欧美成人性av电影在线观看| 亚洲av片天天在线观看|