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

    土壤氣相抽提過(guò)程中多孔介質(zhì)擾動(dòng)的數(shù)值分析

    2013-01-18 07:00:58章光新楊?lèi)傛i董李勤
    中國(guó)環(huán)境科學(xué) 2013年3期
    關(guān)鍵詞:真空度氣相增幅

    范 偉,章光新*,路 瑩,楊?lèi)傛i,3,董李勤

    (1.中國(guó)科學(xué)院東北地理與農(nóng)業(yè)生態(tài)研究所,中國(guó)科學(xué)院濕地生態(tài)與環(huán)境重點(diǎn)實(shí)驗(yàn)室,吉林 長(zhǎng)春 130012;2.吉林大學(xué)環(huán)境與資源學(xué)院,吉林 長(zhǎng)春 130026;3.卡迪夫大學(xué)地球與海洋科學(xué)學(xué)院,英國(guó) 卡迪夫 CF 103YE)

    全球工業(yè)化進(jìn)程中的場(chǎng)地污染問(wèn)題近年不斷凸顯、引發(fā)了各種環(huán)境風(fēng)險(xiǎn)與安全隱患,場(chǎng)地污染修復(fù)因而成為全社會(huì)關(guān)注的熱點(diǎn),污染修復(fù)理論、方法與技術(shù)的研究作為污染防治工作的科技支撐,是當(dāng)前環(huán)境領(lǐng)域重點(diǎn)關(guān)注的課題[1-2].與以往的修復(fù)研究聚焦于地下環(huán)境條件、污染物性質(zhì)與系統(tǒng)運(yùn)行特征對(duì)修復(fù)過(guò)程及機(jī)制的影響相比[3-4],近年來(lái)一些研究開(kāi)始關(guān)注場(chǎng)地人為修復(fù)活動(dòng)對(duì)外部環(huán)境的“反作用”,即探討修復(fù)過(guò)程可能引起的地下環(huán)境物理、化學(xué)與生物變化,例如注氣修復(fù)帶來(lái)的地下水水丘與側(cè)向遷移[5]、強(qiáng)化生物修復(fù)引起的水化學(xué)[6]及微生物群落響應(yīng)[7]等,這些變化會(huì)對(duì)場(chǎng)地修復(fù)過(guò)程產(chǎn)生反饋.

    土壤氣相抽提(SVE)是當(dāng)前全球廣泛使用的土壤與地下水污染高效修復(fù)技術(shù)之一[8],主要利用真空設(shè)備產(chǎn)生負(fù)壓梯度、驅(qū)使一定地下區(qū)域內(nèi)(抽提影響帶)的空氣在孔隙通道中運(yùn)移,攜帶揮發(fā)性污染組分流向抽提系統(tǒng),最后在地面進(jìn)行收集處理[9].已有研究集中于對(duì)比分析土壤介質(zhì)、污染組分、系統(tǒng)運(yùn)行參數(shù)等因素對(duì)SVE過(guò)程及其修復(fù)效果的影響,揭示了環(huán)境特征與系統(tǒng)調(diào)控影響下的SVE特征和機(jī)理[4,10].作者前期研究發(fā)現(xiàn),在室內(nèi)及實(shí)際污染場(chǎng)地的Air Sparging-SVE聯(lián)合修復(fù)中,還可觀(guān)測(cè)到地下多孔介質(zhì)顆粒在氣流擾動(dòng)驅(qū)替作用下發(fā)生遷移,多孔介質(zhì)的孔隙結(jié)構(gòu)與滲透性均發(fā)生改變[11],這一現(xiàn)象與 Kong等[12]的室內(nèi)觀(guān)測(cè)一致.考慮到多孔介質(zhì)的顆粒配級(jí)、滲透性及非均勻性等與抽提影響帶的氣體流型、氣-液相間接觸以及污染修復(fù)過(guò)程密切相關(guān)[4,13],刻畫(huà)抽提過(guò)程對(duì)多孔介質(zhì)的擾動(dòng)成為SVE機(jī)理研究及工程設(shè)計(jì)不可忽視的關(guān)鍵問(wèn)題.基于此,本文從理論機(jī)理方面開(kāi)展土壤氣相抽提過(guò)程中多孔介質(zhì)擾動(dòng)的數(shù)值分析,為上述有關(guān)問(wèn)題的定量化研究奠定基礎(chǔ).

    1 材料與方法

    1.1 SVE流動(dòng)過(guò)程理論模型

    1.1.1 SVE水氣二相流動(dòng)的連續(xù)性方程 假定流體不可壓縮,流動(dòng)過(guò)程符合達(dá)西定律.在 SVE影響帶內(nèi)土壤水-抽提氣流二相流動(dòng)區(qū)選取典型單元體,則根據(jù)質(zhì)量守恒定律在Δt內(nèi)流入與流出單元體的流體質(zhì)量差應(yīng)與其內(nèi)流體的變化量相等:

    式中:mx、my、mz分別為通過(guò)各向單位面積的質(zhì)量;vi為滲流速度;iρ為i相流體密度;φ為單元體孔隙度;Si為i相飽和度.

    根據(jù)連續(xù)介質(zhì)流體動(dòng)力學(xué)理論,i相的孔隙流速可用Darcy定律表示為:

    式中:k為介質(zhì)固有滲透率;kri為i相的相對(duì)滲透率;iμ為i相動(dòng)力黏滯系數(shù);pi為i相的壓力;g為重力加速度;?Hi為i相壓力頭差.

    將滲流速度的 Darcy表達(dá)式代入式(2)得到任意i相流體(水/氣)的滲流連續(xù)性方程:

    1.1.2 模型變量及其耦合關(guān)系 對(duì)式(4)中的變量進(jìn)行定義,并構(gòu)建不同變量之間的耦合關(guān)系[14-15]如下:

    式中:Sw,Sg分別為水相和氣相的飽和度;Pg,Pw分別為氣相和水相的壓力;m為 van Genuchten-Mualem 模型常量.由式(4)~式(8)共同組成 SVE二相流動(dòng)過(guò)程模型,聯(lián)立可實(shí)現(xiàn)SVE過(guò)程中飽和度、壓力場(chǎng)及速度場(chǎng)的求解.

    1.2 多孔介質(zhì)擾動(dòng)及其定量表征

    對(duì)于 SVE影響帶內(nèi)體積為Ve的典型單元體,n為孔隙度,SVE運(yùn)行之前單元體內(nèi)初始?xì)怏w及土壤水占據(jù)的體積為Vg_0、Vw_0,則初始?xì)?、水相飽和度Sg_0、Sw_0為:

    式中:tφ與0φ分別為t時(shí)刻與初始狀態(tài)的空隙度.由式(17)可定量識(shí)別Δt時(shí)段內(nèi)的滲透率變化.

    1.3 SVE過(guò)程及多孔介質(zhì)擾動(dòng)的概念模型

    1.3.1 物理模型及參數(shù)賦值 如圖1所示,構(gòu)建單井SVE軸對(duì)稱(chēng)物理模型.在(x,y)典型垂向剖面上,包氣帶厚度為M,單側(cè)寬度為L(zhǎng),地下水面設(shè)置為y=0平面,SVE井抽提段在(0,l1)~(0,l2)之間.

    圖1 場(chǎng)地SVE示意Fig.1 Definition sketch of the SVE physical system

    場(chǎng)地相關(guān)的模型參數(shù)如下表所示:

    1.3.2 模型定解條件 模型的定解條件包括初始條件及邊界條件兩部分:

    初始條件

    t=0時(shí),空間各處均為二相等勢(shì)體,水相壓力為 :pw_0=pa+ 5[m]·ρg· g-ρw·g ·y;氣 相 壓 力為:pg_0=pa+ρg·(5 [ m]-y)·g .

    表1 模型常量參數(shù)賦值表[11,19]Table 1 Summary of physical characteristics in field SVE experiment

    邊界條件

    其中,Γ1為一類(lèi)邊界,Γ2為二類(lèi)邊界.

    2 數(shù)值模擬及結(jié)果分析

    上述數(shù)學(xué)模型的有限元三角剖分、邊界條件和初始條件輸入以及偏微分方程的求解均通過(guò)Femlab實(shí)現(xiàn).考慮到實(shí)際應(yīng)用中一般SVE抽提真空度 Δp=pa-ps在 10~30kPa之間,故本次分別模擬SVE抽提段壓力為0.7×105Pa和0.9×105Pa的兩種壓力情景(相應(yīng)的真空度Δp約為 31kPa和 11kPa).參數(shù)C與介質(zhì)特性有關(guān),且由(16)式可知函數(shù) Δφ=f(C)的一階導(dǎo)數(shù)大于 0,本次選取C=0.1、0.3、0.5、0.8不同取值時(shí)的4種模擬情景進(jìn)行對(duì)比分析.

    2.1 氣相壓力場(chǎng)及速度場(chǎng)分布

    通過(guò)二相流動(dòng)模擬分析可獲取穩(wěn)定狀態(tài)下模擬區(qū)壓力場(chǎng)及速度場(chǎng)的空間分布(圖 2).結(jié)果顯示:越接近 SVE抽提段,空間各處氣相的真空度、壓力梯度及氣流速度均越大,以有效真空度為 0.1kPa 為界線(xiàn),ps=0.7×105Pa(Δp=31kPa)時(shí)抽提有效覆蓋半徑達(dá)到9m,ps=0.9×105Pa(Δp=11kPa)時(shí)則減小至 8.5m,可知抽提真空度大的條件下SVE影響半徑也較大,在半徑小于6m范圍內(nèi),氣流速度可達(dá)到 0.0001m/s.進(jìn)一步模擬計(jì)算穩(wěn)態(tài)條件下的氣相飽和度,如圖中短劃曲線(xiàn)所示,包氣帶部分區(qū)域出現(xiàn)Sg=0,即抽提作用下地下水位出現(xiàn)一定的抬升,但在兩種壓力情景下抬升的幅度差異并不明顯,短劃曲線(xiàn)Sg=0幾乎重合.

    圖2 穩(wěn)態(tài)不同情境下的氣相真空度與速度場(chǎng)分布Fig.2 Spatial distribution of air vaccum and velocity at steady stage of SVE

    2.2 特定壓力及C值情景下孔隙度及滲透率隨時(shí)間的動(dòng)態(tài)變化

    在ps=0.7×105Pa、C=0.8情景下,對(duì)上述數(shù)學(xué)模型進(jìn)行瞬態(tài)模擬分析.選取模擬區(qū)內(nèi)與抽提段距離遞增的4個(gè)典型點(diǎn)P1(1,1),P2(1.5,1.5), P3(2,2),P4(3,3),則模型求解可分別得到任意時(shí)刻各點(diǎn)的孔隙度變幅,并根據(jù)孔隙度~滲透率關(guān)系得到滲透率隨時(shí)間的動(dòng)態(tài)變化.結(jié)果如圖 3所示:SVE過(guò)程中孔隙度及滲透率隨時(shí)間呈現(xiàn)先增加后穩(wěn)定的動(dòng)態(tài)變化,包氣帶介質(zhì)整體達(dá)到穩(wěn)定所需的時(shí)長(zhǎng)約為 60min;與抽提段距離越近,孔隙度及滲透率變化響應(yīng)越快、達(dá)到平穩(wěn)所需的時(shí)間也越短,增幅越大,如P1點(diǎn)在約40min后達(dá)到穩(wěn)定、孔隙度增幅為0.0387,P4點(diǎn)約在60min后達(dá)到穩(wěn)定、孔隙度增幅為0.0031,相應(yīng)地,P1和P4點(diǎn)介質(zhì)滲透率分別增加至 2.22×10-11m2和 1.25×10-11m2.

    圖3 SVE過(guò)程典型點(diǎn)介質(zhì)孔隙度及滲透率隨時(shí)間動(dòng)態(tài)變化Fig.3 Temporal trends of porosity and permeability at four spatial points during SVE

    2.3 不同壓力情景參數(shù)C不同取值條件下的孔隙度變化

    基于上述不同壓力情景的模擬分析,結(jié)合多孔介質(zhì)擾動(dòng)的定量表征,得到參數(shù)C不同取值條件下穩(wěn)定后的孔隙度變化,如圖4所示.離SVE抽提段越近,介質(zhì)孔隙度變化越大.相同抽提壓力下,伴隨著C值的增加,空間同一位置的孔隙度變幅逐漸增大,ps=0.9×105Pa時(shí),Δφ=0.005分布界線(xiàn)的橫向最大距離在C=0.1,0.3,0.5,0.8時(shí)分別為1.36,4.04,5.34,6.75m,且C=0.1和 0.8時(shí)的最大Δφ分別約為 0.009和 0.055,即影響范圍及程度隨著C的增加而增大;對(duì)比相同C值條件下(C=0.8),不同 SVE 抽提壓力情景(ps=0.9,0.7×105Pa),可知形成真空度越大的條件下(抽提壓力0.7×105Pa)孔隙度增幅越大,Δφ=0.005分布界線(xiàn)的橫向最大距離增加至 7.24m,最大Δφ則達(dá)到0.066.綜上所述,孔隙度增幅與抽提真空度和參數(shù)C取值呈正相關(guān).

    3 結(jié)論

    3.1 對(duì)于特定的場(chǎng)地,SVE影響帶的空間形態(tài)與抽提真空度密切相關(guān),抽提真空度越大,影響半徑越大.就本模擬場(chǎng)地而言,SVE抽提段真空度為11kPa和31kPa的2種壓力情景下,其影響半徑分別約為8.5和9.0m;

    3.2 SVE過(guò)程中多孔介質(zhì)的孔隙度及滲透率呈現(xiàn)顯著的時(shí)空動(dòng)態(tài),二者隨時(shí)間呈現(xiàn)先增加后穩(wěn)定的變化規(guī)律,達(dá)到穩(wěn)定所需的時(shí)長(zhǎng)及其變化幅度則與離抽提段的空間距離成反相關(guān),距離抽提段較近的 P1點(diǎn)在約40min后達(dá)到穩(wěn)定、孔隙度增幅為 0.0387,較遠(yuǎn)的P4點(diǎn)約在60min后達(dá)到穩(wěn)定、孔隙度增幅為 0.0031,相應(yīng)地,P1和 P4點(diǎn)介質(zhì)滲透率分別從 1.18×10-11m2增加至 2.22×10-11、1.25×10-11m2;

    3.3 在相同抽提壓力下,孔隙度增幅與關(guān)鍵參數(shù)C值成正相關(guān),抽提壓力為 0.9×105Pa條件下,C=0.1和 0.8時(shí)的孔隙度最大增幅分別約為0.009和 0.055;相同參數(shù)C條件下,孔隙度增幅與抽提壓力成正相關(guān),C=0.8條件下,抽提壓力為 0.7×105Pa時(shí)的孔隙度最大增幅則達(dá)到0.066.

    圖4 不同抽提壓力及C取值條件下的孔隙度增幅空間分布Fig.4 Comparisons of changes in porosity under different SVE settings

    [1] Reddy K R. Technical Challenges to in-situ remediation of polluted sites [J]. Geotechnical and Geological Engineering,2010,28(3):211-221.

    [2] 駱永明.中國(guó)污染場(chǎng)地修復(fù)的研究進(jìn)展、問(wèn)題與展望 [J]. 環(huán)境監(jiān)測(cè)管理與技術(shù), 2011,23(3):1-6.

    [3] Rogers S W, Ong S K. Inflence of porous media, airflow rate, and air channel spacing on benzene NAPL removal during air sparging [J]. Environmental Science and Technology, 2000,34(5):764-770.

    [4] Qin C Y, Zhao Y S, Zheng W, et al. Study on influencing factors on removal of chlorobenzene from unsaturated zone by soil vapor extraction [J]. Journal of Hazardous Materials, 2010,176(1-3):294-299.

    [5] 范 偉,楊?lèi)傛i,路 瑩,等.層間地下水污染曝氣修復(fù)的影響帶[J]. 化工學(xué)報(bào), 2011,62(9):2600-2607.

    [6] Davis G B, Patterson B M, Johnston C D. Aerobic bioremediation of 1, 2 dichloroethane and vinyl chloride at field scale [J]. Journal of Contaminant Hydrology, 2009, 107(1-2):91-100.

    [7] Kao C M, Chen C Y, Chen S C, et al. Application of insitu biosparging to remediate a petroleum-hydrocarbon spill site: Field and microbial evaluation [J]. Chemosphere, 2008,70(8):1492-1499.

    [8] Soares A A, Pinho M T, Albergaria J T, et al. Sequential application of soil vapor extraction and bioremediation processes for the remediation of ethylbenzene-contaminated soils [J]. Water,Air, and Soil Pollution, 2012, DOI: 10.1007/s11270-011-1051-y.

    [9] H?ier C K, Sonnenborg T O, Jensen K H, et al. Model analysis of mechanisms controlling pneumatic soil vapor extraction [J].Journal of Contaminant Hydrology, 2009,103(3-4):82-98.

    [10] Albergaria J T, Alvim-Ferrza M da C, Delerue-Matos C. Soil vapor extraction in sandy soils: Influence of airflow rate [J].Chemosphere, 2008,73(9):1557-1561.

    [11] 范 偉.石油類(lèi)污染場(chǎng)地地下水原位空氣擾動(dòng)修復(fù)研究 [D].長(zhǎng)春:吉林大學(xué), 2011.

    [12] Kong X Z, Kinzelbach W, Stauffer F. Migration of air channels:an instability of air flow in mobile saturated porous media [J].Chemical Engineering Science, 2009,64(7):1528-1535.

    [13] Gidda T, Stiver W H, Zytner R G. Measuring aqueous-air and sorbed-aqueous mass transfer coefficients for application in soil vapor extraction [J]. Journal of Environmental Engineering,2011,137(10):880-888.

    [14] Hoeg S, Sch?ler H F, Warnatz J. Assessment of interfacial mass transfer in water-unsaturated soils during vapor extraction [J].Journal of Contaminant Hydrology, 2004,74(1-4):163-195.

    [15] 薛 強(qiáng),馮夏庭,梁 冰,等.水氣二相流系統(tǒng)的 K-S-P模型參數(shù)反演的最優(yōu)估計(jì) [J]. 水科學(xué)進(jìn)展, 2005,16(4):488- 493.

    [16] El-Farhan Y H, Denovio N M, Herman J S, et al. Mobilization and transport of soil particles during infiltration experiments in an agricultural field, Shenandoah Valley, Virginia [J]. Environmental Science and Technology, 2000,34(17):3555-3559.

    [17] Tsai Y J. Air distribution and size changes in the remediated zone after air sparging for soil particle movement [J]. Journal of Hazardous Materials, 2008,158(2/3):438-444.

    [18] 路 瑩.北京平谷地區(qū)雨洪水地下回灌堵塞機(jī)理分析與模擬研究 [D]. 長(zhǎng)春:吉林大學(xué), 2009.

    [19] Leij F J, Russell W B, Lesch S M. Closed-form expressions for water retention and conductivity data [J]. Ground Water, 1997,35(5):848-858.

    猜你喜歡
    真空度氣相增幅
    2023年全省一般公共預(yù)算收支情況一覽表
    上半年全省一般公共預(yù)算收支情況一覽表
    氣相過(guò)渡金屬鈦-碳鏈團(tuán)簇的研究
    直接空冷機(jī)組提高真空度的分析與改進(jìn)
    能源工程(2020年6期)2021-01-26 00:55:18
    擠奶機(jī)集乳器真空度評(píng)估研究
    2016年汽車(chē)產(chǎn)銷(xiāo)雙雙超過(guò)2800萬(wàn)輛增幅比上年較快提升
    新型釩基催化劑催化降解氣相二噁英
    預(yù)縮聚反應(yīng)器氣相管“鼓泡”的成因探討
    氣相防銹技術(shù)在電器設(shè)備防腐中的應(yīng)用
    儀器儀表業(yè)今年產(chǎn)銷(xiāo)增幅或達(dá)15%
    好看av亚洲va欧美ⅴa在| x7x7x7水蜜桃| 亚洲精品乱码久久久v下载方式| 99热这里只有是精品50| 国产精品1区2区在线观看.| 嫩草影视91久久| 88av欧美| 欧美xxxx黑人xx丫x性爽| 久久久国产成人精品二区| 99热精品在线国产| 亚洲中文字幕一区二区三区有码在线看| 国产高清激情床上av| 噜噜噜噜噜久久久久久91| 日韩亚洲欧美综合| 国产成人a区在线观看| 免费av不卡在线播放| 深夜a级毛片| 成人特级av手机在线观看| 国产亚洲精品久久久com| 成人亚洲精品av一区二区| 亚洲无线观看免费| 99国产极品粉嫩在线观看| АⅤ资源中文在线天堂| 精品人妻熟女av久视频| 天堂av国产一区二区熟女人妻| 三级男女做爰猛烈吃奶摸视频| 男女下面进入的视频免费午夜| 国内毛片毛片毛片毛片毛片| 国产亚洲欧美98| 免费av观看视频| 国产高清视频在线播放一区| 国产精品一区二区三区四区久久| 亚洲欧美清纯卡通| 亚洲天堂国产精品一区在线| 精品国产亚洲在线| 欧美精品国产亚洲| 国产精品久久久久久人妻精品电影| 国产单亲对白刺激| 欧美绝顶高潮抽搐喷水| 精品久久久久久,| 国产精品爽爽va在线观看网站| 一本一本综合久久| 国产欧美日韩精品一区二区| 两个人的视频大全免费| 麻豆国产97在线/欧美| 亚洲av免费高清在线观看| 国产精品一区二区免费欧美| 国产精品久久久久久久电影| 亚洲精品色激情综合| 欧美成人一区二区免费高清观看| 国产精品久久电影中文字幕| 亚洲经典国产精华液单 | 在现免费观看毛片| 免费人成在线观看视频色| 在线观看免费视频日本深夜| 在线免费观看的www视频| 色精品久久人妻99蜜桃| 又爽又黄a免费视频| 亚洲色图av天堂| 亚洲精品在线观看二区| 亚洲av一区综合| 成年免费大片在线观看| 99久久九九国产精品国产免费| 97人妻精品一区二区三区麻豆| 一个人免费在线观看电影| 最新在线观看一区二区三区| 成年版毛片免费区| 他把我摸到了高潮在线观看| 色播亚洲综合网| 色综合欧美亚洲国产小说| 成人高潮视频无遮挡免费网站| 精品久久久久久,| 免费大片18禁| 亚洲av不卡在线观看| 757午夜福利合集在线观看| 日韩欧美精品v在线| av在线天堂中文字幕| 两人在一起打扑克的视频| 韩国av一区二区三区四区| 国产精品影院久久| 黄色视频,在线免费观看| 无人区码免费观看不卡| 少妇熟女aⅴ在线视频| 丁香欧美五月| 最近中文字幕高清免费大全6 | 欧美高清性xxxxhd video| 国产精品一区二区三区四区免费观看 | 99久久精品一区二区三区| 国产精品综合久久久久久久免费| 久久人人爽人人爽人人片va | 精品福利观看| 国内揄拍国产精品人妻在线| 99久久无色码亚洲精品果冻| 亚洲av不卡在线观看| 99久久九九国产精品国产免费| 男人舔女人下体高潮全视频| 中国美女看黄片| 级片在线观看| 亚洲中文字幕一区二区三区有码在线看| 亚洲熟妇熟女久久| 真人一进一出gif抽搐免费| 怎么达到女性高潮| 国产精品久久久久久亚洲av鲁大| 国产精品女同一区二区软件 | 757午夜福利合集在线观看| 国内揄拍国产精品人妻在线| 国产淫片久久久久久久久 | 欧美黑人欧美精品刺激| 午夜精品久久久久久毛片777| 国产激情偷乱视频一区二区| 亚洲国产精品sss在线观看| 岛国在线免费视频观看| 久久精品国产亚洲av天美| 国产黄片美女视频| 我的女老师完整版在线观看| 久99久视频精品免费| 丁香六月欧美| 国产精品久久久久久亚洲av鲁大| 日本a在线网址| 精品日产1卡2卡| 久久久久国产精品人妻aⅴ院| 久久亚洲精品不卡| 真人一进一出gif抽搐免费| 免费人成在线观看视频色| 亚洲专区中文字幕在线| 十八禁网站免费在线| 久久久国产成人免费| 丁香六月欧美| 精品久久久久久久久亚洲 | 如何舔出高潮| 国产亚洲精品久久久久久毛片| 香蕉av资源在线| 一本久久中文字幕| 脱女人内裤的视频| ponron亚洲| 琪琪午夜伦伦电影理论片6080| 日韩 亚洲 欧美在线| 欧美一级a爱片免费观看看| 成年女人看的毛片在线观看| 九九热线精品视视频播放| 国产麻豆成人av免费视频| 国产激情偷乱视频一区二区| 3wmmmm亚洲av在线观看| 亚洲五月天丁香| 久久久久精品国产欧美久久久| 精品不卡国产一区二区三区| 国产一区二区激情短视频| 老司机午夜十八禁免费视频| 成人国产综合亚洲| 88av欧美| 日韩有码中文字幕| 两性午夜刺激爽爽歪歪视频在线观看| 热99在线观看视频| 精品熟女少妇八av免费久了| 熟妇人妻久久中文字幕3abv| 亚洲av.av天堂| 欧美乱妇无乱码| 脱女人内裤的视频| 欧美又色又爽又黄视频| 两个人视频免费观看高清| 色av中文字幕| 亚洲成人久久性| 久久久久久大精品| 中文字幕免费在线视频6| 欧美最新免费一区二区三区 | 日韩欧美在线二视频| 中文亚洲av片在线观看爽| 五月玫瑰六月丁香| 久久99热6这里只有精品| 在线十欧美十亚洲十日本专区| avwww免费| 欧美成人一区二区免费高清观看| 久久久久久久久中文| 久久国产乱子免费精品| 一本久久中文字幕| 欧美+日韩+精品| 婷婷色综合大香蕉| 亚洲av免费高清在线观看| 成人精品一区二区免费| 日韩成人在线观看一区二区三区| 青草久久国产| 色综合婷婷激情| 婷婷精品国产亚洲av在线| 美女xxoo啪啪120秒动态图 | 国产av不卡久久| 成人性生交大片免费视频hd| 久久久久久久久久成人| 日本精品一区二区三区蜜桃| 成年女人毛片免费观看观看9| 成人高潮视频无遮挡免费网站| 亚洲综合色惰| 亚洲人与动物交配视频| 内射极品少妇av片p| 国产免费av片在线观看野外av| 青草久久国产| 熟妇人妻久久中文字幕3abv| 精品午夜福利在线看| 在线免费观看不下载黄p国产 | 搡老岳熟女国产| 免费观看的影片在线观看| 天堂av国产一区二区熟女人妻| 日韩欧美精品v在线| 美女cb高潮喷水在线观看| 9191精品国产免费久久| 宅男免费午夜| 男人狂女人下面高潮的视频| 午夜福利在线观看免费完整高清在 | 国产精品伦人一区二区| 日韩中文字幕欧美一区二区| 观看美女的网站| 别揉我奶头 嗯啊视频| 亚洲电影在线观看av| 午夜福利免费观看在线| 乱码一卡2卡4卡精品| 嫩草影视91久久| 欧美成人a在线观看| 天堂网av新在线| 深夜a级毛片| 18禁裸乳无遮挡免费网站照片| 国产精品电影一区二区三区| 国产一区二区在线观看日韩| 国产真实乱freesex| 人妻制服诱惑在线中文字幕| 一二三四社区在线视频社区8| 色视频www国产| 国语自产精品视频在线第100页| 男女床上黄色一级片免费看| 成年女人毛片免费观看观看9| 在线天堂最新版资源| 丰满乱子伦码专区| www.熟女人妻精品国产| 91字幕亚洲| 狠狠狠狠99中文字幕| 国产又黄又爽又无遮挡在线| 精品人妻1区二区| 啪啪无遮挡十八禁网站| 波多野结衣巨乳人妻| av欧美777| 99精品久久久久人妻精品| 欧美日韩亚洲国产一区二区在线观看| av在线蜜桃| 深爱激情五月婷婷| 国产在视频线在精品| 欧美一区二区精品小视频在线| 亚洲七黄色美女视频| 亚洲在线观看片| 美女大奶头视频| 久久久久九九精品影院| 欧美高清成人免费视频www| av欧美777| 能在线免费观看的黄片| 在线天堂最新版资源| 一级毛片久久久久久久久女| 看免费av毛片| 国产免费av片在线观看野外av| 免费在线观看亚洲国产| 夜夜看夜夜爽夜夜摸| 亚洲性夜色夜夜综合| 欧美日韩亚洲国产一区二区在线观看| 国产极品精品免费视频能看的| 久久久久久大精品| 最近最新中文字幕大全电影3| 国产乱人视频| 我要搜黄色片| 欧洲精品卡2卡3卡4卡5卡区| 精品一区二区三区人妻视频| 两个人视频免费观看高清| 十八禁国产超污无遮挡网站| 成人三级黄色视频| 最近在线观看免费完整版| 亚洲一区二区三区不卡视频| 丰满人妻一区二区三区视频av| 中文字幕精品亚洲无线码一区| 国产探花极品一区二区| 哪里可以看免费的av片| 99久国产av精品| 成人特级黄色片久久久久久久| 国产精品自产拍在线观看55亚洲| 天堂影院成人在线观看| 婷婷精品国产亚洲av在线| 亚洲国产精品成人综合色| 十八禁国产超污无遮挡网站| av国产免费在线观看| 久久久久九九精品影院| 欧美日韩瑟瑟在线播放| 久久久久久国产a免费观看| 人妻久久中文字幕网| 亚洲 欧美 日韩 在线 免费| 精品免费久久久久久久清纯| 欧美潮喷喷水| 丰满的人妻完整版| 欧美中文日本在线观看视频| 国内精品美女久久久久久| 黄色视频,在线免费观看| 亚洲国产精品久久男人天堂| 久久欧美精品欧美久久欧美| 日本 欧美在线| 免费在线观看日本一区| 久久精品国产亚洲av香蕉五月| 欧美乱妇无乱码| 亚洲无线观看免费| 美女黄网站色视频| 中文字幕免费在线视频6| 亚洲成人久久性| 亚洲欧美日韩卡通动漫| 国产av麻豆久久久久久久| 国产不卡一卡二| 欧美xxxx性猛交bbbb| 757午夜福利合集在线观看| 99在线人妻在线中文字幕| 精品99又大又爽又粗少妇毛片 | 亚洲乱码一区二区免费版| www.www免费av| 国产成+人综合+亚洲专区| 精品乱码久久久久久99久播| 国产精品久久电影中文字幕| 一级a爱片免费观看的视频| 好男人在线观看高清免费视频| 亚洲国产精品合色在线| 久久精品综合一区二区三区| 又爽又黄a免费视频| 国产av在哪里看| 毛片女人毛片| 伦理电影大哥的女人| 国产欧美日韩精品亚洲av| 91麻豆av在线| 每晚都被弄得嗷嗷叫到高潮| 99riav亚洲国产免费| 免费一级毛片在线播放高清视频| 精品久久久久久,| 午夜精品久久久久久毛片777| 成人国产综合亚洲| 亚洲性夜色夜夜综合| 久久国产精品人妻蜜桃| 一a级毛片在线观看| 午夜精品久久久久久毛片777| 亚洲经典国产精华液单 | 在线a可以看的网站| 长腿黑丝高跟| 欧美色欧美亚洲另类二区| 1000部很黄的大片| 一夜夜www| 欧洲精品卡2卡3卡4卡5卡区| 欧美在线黄色| 性欧美人与动物交配| 神马国产精品三级电影在线观看| 午夜福利视频1000在线观看| 在线a可以看的网站| 亚洲av熟女| 午夜视频国产福利| 淫秽高清视频在线观看| 成人一区二区视频在线观看| АⅤ资源中文在线天堂| 天堂网av新在线| 久久久久久久久大av| 国产v大片淫在线免费观看| 久久欧美精品欧美久久欧美| 一进一出抽搐动态| 悠悠久久av| 亚洲成av人片免费观看| 国产亚洲精品久久久久久毛片| 色播亚洲综合网| 欧美绝顶高潮抽搐喷水| 欧美性猛交黑人性爽| 国内精品久久久久精免费| 高清日韩中文字幕在线| 日本三级黄在线观看| 亚洲国产精品999在线| 免费在线观看亚洲国产| 免费在线观看影片大全网站| 国产亚洲av嫩草精品影院| 久久国产精品影院| 性色avwww在线观看| 成年版毛片免费区| 一区二区三区免费毛片| 亚洲第一区二区三区不卡| 中亚洲国语对白在线视频| 久久久久国产精品人妻aⅴ院| 男女之事视频高清在线观看| 日本熟妇午夜| 国产高清有码在线观看视频| 色5月婷婷丁香| 亚洲最大成人av| 天堂av国产一区二区熟女人妻| 有码 亚洲区| 99热精品在线国产| 色综合站精品国产| 91在线观看av| 搡老妇女老女人老熟妇| 日韩精品中文字幕看吧| 欧美成人一区二区免费高清观看| 国产精品久久久久久久电影| 在现免费观看毛片| 国产爱豆传媒在线观看| 国产成人aa在线观看| 一区二区三区高清视频在线| 三级毛片av免费| 好男人电影高清在线观看| 国产美女午夜福利| 国产精品久久久久久亚洲av鲁大| 日本黄色片子视频| 深夜a级毛片| 国产精品久久久久久久久免 | 非洲黑人性xxxx精品又粗又长| 日本免费一区二区三区高清不卡| 精品熟女少妇八av免费久了| 人人妻,人人澡人人爽秒播| xxxwww97欧美| 午夜福利18| 老司机午夜福利在线观看视频| 无遮挡黄片免费观看| 一本久久中文字幕| 国产午夜福利久久久久久| 亚洲国产色片| 午夜福利在线在线| 欧美3d第一页| 国内揄拍国产精品人妻在线| 国产大屁股一区二区在线视频| 色综合亚洲欧美另类图片| 99久久成人亚洲精品观看| 亚州av有码| 亚洲三级黄色毛片| 欧美日韩中文字幕国产精品一区二区三区| 在线国产一区二区在线| 国产伦人伦偷精品视频| 日韩精品中文字幕看吧| 亚洲精品久久国产高清桃花| 三级国产精品欧美在线观看| 可以在线观看的亚洲视频| 久久欧美精品欧美久久欧美| 男人狂女人下面高潮的视频| 久久亚洲精品不卡| 真实男女啪啪啪动态图| АⅤ资源中文在线天堂| 国产亚洲精品综合一区在线观看| 亚洲av电影在线进入| 成人永久免费在线观看视频| 最近中文字幕高清免费大全6 | 1000部很黄的大片| 国内精品久久久久久久电影| 美女免费视频网站| 丰满人妻熟妇乱又伦精品不卡| 美女黄网站色视频| 免费人成在线观看视频色| 亚洲无线在线观看| 日本与韩国留学比较| 久久99热6这里只有精品| 亚洲国产精品999在线| 亚洲avbb在线观看| 色吧在线观看| 国产v大片淫在线免费观看| 亚洲精品日韩av片在线观看| 欧美高清性xxxxhd video| 高清日韩中文字幕在线| 波多野结衣巨乳人妻| 天美传媒精品一区二区| av天堂中文字幕网| 简卡轻食公司| 亚洲精品一区av在线观看| 欧美另类亚洲清纯唯美| 久久精品国产清高在天天线| 日日摸夜夜添夜夜添小说| 我要搜黄色片| 免费av不卡在线播放| 丝袜脚勾引网站| 黄片无遮挡物在线观看| 国产爽快片一区二区三区| 爱豆传媒免费全集在线观看| 国产午夜精品一二区理论片| 免费观看性生交大片5| 国产老妇女一区| 成年人午夜在线观看视频| 2021天堂中文幕一二区在线观| 狠狠精品人妻久久久久久综合| 亚州av有码| 国产成人午夜福利电影在线观看| av国产精品久久久久影院| 国产有黄有色有爽视频| 久久人人爽人人片av| 又大又黄又爽视频免费| 青春草视频在线免费观看| 国产综合精华液| 大陆偷拍与自拍| 国产淫片久久久久久久久| 黄片无遮挡物在线观看| 国产亚洲av片在线观看秒播厂| 欧美另类一区| 欧美激情久久久久久爽电影| 男人和女人高潮做爰伦理| 只有这里有精品99| 尤物成人国产欧美一区二区三区| 一级a做视频免费观看| 久久精品国产亚洲av涩爱| 久久久久九九精品影院| 免费高清在线观看视频在线观看| 亚洲av在线观看美女高潮| 精品久久国产蜜桃| 欧美激情久久久久久爽电影| 国产片特级美女逼逼视频| 国产精品无大码| 精品亚洲乱码少妇综合久久| 精品一区在线观看国产| 黄色日韩在线| 水蜜桃什么品种好| 视频区图区小说| 午夜精品一区二区三区免费看| 在线观看国产h片| 白带黄色成豆腐渣| 水蜜桃什么品种好| 国内揄拍国产精品人妻在线| 久久这里有精品视频免费| 亚洲色图综合在线观看| 午夜福利网站1000一区二区三区| 欧美成人精品欧美一级黄| 国产精品人妻久久久影院| 亚洲,欧美,日韩| 丰满少妇做爰视频| 免费高清在线观看视频在线观看| 新久久久久国产一级毛片| 狂野欧美激情性xxxx在线观看| 亚洲经典国产精华液单| 最近的中文字幕免费完整| 最新中文字幕久久久久| 久热这里只有精品99| 色播亚洲综合网| 国产爱豆传媒在线观看| 欧美高清性xxxxhd video| 国产一区二区亚洲精品在线观看| 黄色一级大片看看| 亚洲成人久久爱视频| 欧美日韩精品成人综合77777| 国产色婷婷99| 国产黄片视频在线免费观看| 两个人的视频大全免费| 狠狠精品人妻久久久久久综合| 中文精品一卡2卡3卡4更新| 国产精品99久久久久久久久| 久久久色成人| 3wmmmm亚洲av在线观看| 久久这里有精品视频免费| 卡戴珊不雅视频在线播放| 国产女主播在线喷水免费视频网站| 麻豆乱淫一区二区| 成人国产麻豆网| 少妇的逼水好多| 18禁在线无遮挡免费观看视频| 22中文网久久字幕| 干丝袜人妻中文字幕| 国产黄频视频在线观看| 高清在线视频一区二区三区| 九色成人免费人妻av| 国产精品久久久久久久久免| 亚洲精品日韩av片在线观看| 观看免费一级毛片| 国产综合精华液| 亚洲欧美一区二区三区黑人 | 日韩av在线免费看完整版不卡| 亚洲人成网站在线播| 午夜福利在线观看免费完整高清在| 免费看日本二区| 久久人人爽人人爽人人片va| 赤兔流量卡办理| 日韩欧美精品免费久久| 欧美最新免费一区二区三区| 建设人人有责人人尽责人人享有的 | 久久久久九九精品影院| 国产黄频视频在线观看| 中文字幕亚洲精品专区| a级毛片免费高清观看在线播放| 中文资源天堂在线| 国产精品熟女久久久久浪| 制服丝袜香蕉在线| 国产91av在线免费观看| 噜噜噜噜噜久久久久久91| 18禁动态无遮挡网站| 免费大片黄手机在线观看| 欧美日韩精品成人综合77777| 蜜桃亚洲精品一区二区三区| 中文精品一卡2卡3卡4更新| 久久久精品免费免费高清| av卡一久久| 亚洲电影在线观看av| 如何舔出高潮| a级毛片免费高清观看在线播放| 亚洲性久久影院| 亚洲国产色片| 午夜福利高清视频| 别揉我奶头 嗯啊视频| 爱豆传媒免费全集在线观看| 九草在线视频观看| av又黄又爽大尺度在线免费看| 欧美老熟妇乱子伦牲交| 美女cb高潮喷水在线观看| 国产免费福利视频在线观看| 欧美高清性xxxxhd video| 欧美性猛交╳xxx乱大交人| 人妻一区二区av| 老司机影院成人| 午夜福利高清视频| 国产亚洲av片在线观看秒播厂| 自拍欧美九色日韩亚洲蝌蚪91 | 中文字幕制服av| 久久久成人免费电影| 亚洲无线观看免费| 日本爱情动作片www.在线观看| 日韩av不卡免费在线播放| 久久久a久久爽久久v久久| 777米奇影视久久| 少妇高潮的动态图| av卡一久久| 国产一区二区三区av在线| 久热这里只有精品99| 中国国产av一级| 两个人的视频大全免费| 最后的刺客免费高清国语|