• <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%
    亚洲国产中文字幕在线视频| 桃色一区二区三区在线观看| 少妇人妻精品综合一区二区 | 波多野结衣高清无吗| 99国产综合亚洲精品| 有码 亚洲区| 久久精品国产清高在天天线| 亚洲 国产 在线| 在线国产一区二区在线| 亚洲一区高清亚洲精品| 欧美黑人巨大hd| 亚洲激情在线av| 亚洲精品在线美女| av专区在线播放| 国产97色在线日韩免费| 久久6这里有精品| 黄片小视频在线播放| 大型黄色视频在线免费观看| 一个人看视频在线观看www免费 | 黑人欧美特级aaaaaa片| 无人区码免费观看不卡| 一夜夜www| 丰满人妻一区二区三区视频av | 啦啦啦免费观看视频1| 国产视频一区二区在线看| 波多野结衣巨乳人妻| 欧美另类亚洲清纯唯美| tocl精华| 男女之事视频高清在线观看| 欧美黄色片欧美黄色片| 日韩欧美三级三区| 午夜免费成人在线视频| 免费看光身美女| 亚洲av成人精品一区久久| 99久久精品国产亚洲精品| 久久精品91无色码中文字幕| 久久香蕉国产精品| 国产一区在线观看成人免费| 在线天堂最新版资源| ponron亚洲| 老汉色∧v一级毛片| 午夜免费成人在线视频| 日韩欧美 国产精品| 男人舔奶头视频| 精品久久久久久,| 久久99热这里只有精品18| 亚洲精品国产精品久久久不卡| 日韩 欧美 亚洲 中文字幕| 草草在线视频免费看| 五月伊人婷婷丁香| 亚洲精品粉嫩美女一区| 90打野战视频偷拍视频| 少妇裸体淫交视频免费看高清| 国产精品一及| 丰满乱子伦码专区| 热99re8久久精品国产| 三级毛片av免费| 老司机福利观看| 丁香六月欧美| 国语自产精品视频在线第100页| 亚洲精品在线美女| 欧美在线黄色| 有码 亚洲区| 一区二区三区免费毛片| 久久久色成人| 在线播放国产精品三级| 亚洲国产日韩欧美精品在线观看 | 床上黄色一级片| 18禁裸乳无遮挡免费网站照片| 非洲黑人性xxxx精品又粗又长| 日韩欧美在线二视频| 99国产精品一区二区三区| 亚洲成人久久性| 欧美av亚洲av综合av国产av| 国产精品影院久久| 亚洲av成人不卡在线观看播放网| 好男人电影高清在线观看| 免费看a级黄色片| 91麻豆精品激情在线观看国产| 国产精品,欧美在线| 听说在线观看完整版免费高清| 成人特级av手机在线观看| 欧美黑人巨大hd| 日韩亚洲欧美综合| 人妻久久中文字幕网| 午夜精品一区二区三区免费看| 少妇的逼水好多| 亚洲电影在线观看av| 啦啦啦免费观看视频1| 日日干狠狠操夜夜爽| 精品久久久久久成人av| 久久香蕉精品热| 亚洲精品色激情综合| a在线观看视频网站| 亚洲第一欧美日韩一区二区三区| 成人亚洲精品av一区二区| 最近最新中文字幕大全电影3| svipshipincom国产片| 色尼玛亚洲综合影院| 国产在线精品亚洲第一网站| 久久久久免费精品人妻一区二区| 成人特级黄色片久久久久久久| 婷婷精品国产亚洲av在线| 综合色av麻豆| 国产精品自产拍在线观看55亚洲| 欧美激情在线99| 人妻久久中文字幕网| 搞女人的毛片| 变态另类成人亚洲欧美熟女| 午夜影院日韩av| 国产精品一区二区三区四区久久| 一区二区三区免费毛片| 可以在线观看的亚洲视频| 欧洲精品卡2卡3卡4卡5卡区| 色老头精品视频在线观看| 精品人妻1区二区| 亚洲av电影在线进入| 色精品久久人妻99蜜桃| 黄色丝袜av网址大全| 在线观看舔阴道视频| 国产精品久久久久久亚洲av鲁大| 亚洲成人中文字幕在线播放| 在线观看美女被高潮喷水网站 | 99久久精品热视频| 男女视频在线观看网站免费| 日本三级黄在线观看| 少妇高潮的动态图| 美女大奶头视频| 男女床上黄色一级片免费看| 精品免费久久久久久久清纯| 中文字幕精品亚洲无线码一区| 尤物成人国产欧美一区二区三区| av专区在线播放| 亚洲精品一卡2卡三卡4卡5卡| 97人妻精品一区二区三区麻豆| 19禁男女啪啪无遮挡网站| 性色avwww在线观看| 高清在线国产一区| 亚洲人成网站在线播放欧美日韩| 精品一区二区三区人妻视频| 99在线人妻在线中文字幕| 色播亚洲综合网| 国内久久婷婷六月综合欲色啪| 免费无遮挡裸体视频| 国产精品久久久人人做人人爽| 99久久无色码亚洲精品果冻| 国产v大片淫在线免费观看| 啦啦啦观看免费观看视频高清| 看免费av毛片| 久久精品人妻少妇| 舔av片在线| 久久久久性生活片| 少妇的丰满在线观看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 欧美成狂野欧美在线观看| 日本黄色片子视频| 成人欧美大片| 欧美日韩综合久久久久久 | 亚洲最大成人手机在线| 国产单亲对白刺激| 国产伦精品一区二区三区四那| 757午夜福利合集在线观看| 国内少妇人妻偷人精品xxx网站| 国产主播在线观看一区二区| 舔av片在线| 精品99又大又爽又粗少妇毛片 | 女人高潮潮喷娇喘18禁视频| 色综合亚洲欧美另类图片| 最近最新中文字幕大全免费视频| 女警被强在线播放| 亚洲天堂国产精品一区在线| 欧美乱码精品一区二区三区| 亚洲中文字幕一区二区三区有码在线看| 久久精品人妻少妇| 国产av不卡久久| 国产精品久久视频播放| 日本撒尿小便嘘嘘汇集6| 国产精品电影一区二区三区| 成人特级黄色片久久久久久久| 午夜影院日韩av| 九九热线精品视视频播放| 国产精品嫩草影院av在线观看 | 精品久久久久久久毛片微露脸| 日本撒尿小便嘘嘘汇集6| 少妇高潮的动态图| 19禁男女啪啪无遮挡网站| 99久国产av精品| 成人特级黄色片久久久久久久| 国产三级黄色录像| 亚洲av电影不卡..在线观看| 伊人久久大香线蕉亚洲五| 岛国在线免费视频观看| 怎么达到女性高潮| 国产精品香港三级国产av潘金莲| 一级作爱视频免费观看| 变态另类成人亚洲欧美熟女| 日本免费a在线| 真实男女啪啪啪动态图| 久久国产乱子伦精品免费另类| 免费人成视频x8x8入口观看| 日日摸夜夜添夜夜添小说| 国产成人系列免费观看| 欧美成狂野欧美在线观看| 成人永久免费在线观看视频| 美女高潮的动态| 在线a可以看的网站| 九色成人免费人妻av| 午夜福利成人在线免费观看| 国产精品影院久久| 毛片女人毛片| av专区在线播放| 波多野结衣高清无吗| 天美传媒精品一区二区| 免费看日本二区| 在线观看免费视频日本深夜| 岛国在线免费视频观看| 在线观看舔阴道视频| 宅男免费午夜| 久久精品国产99精品国产亚洲性色| 日本成人三级电影网站| а√天堂www在线а√下载| 欧美一级a爱片免费观看看| 精品久久久久久,| 男人和女人高潮做爰伦理| 内射极品少妇av片p| 欧美黄色片欧美黄色片| 亚洲不卡免费看| 又粗又爽又猛毛片免费看| 久久亚洲真实| 搞女人的毛片| 精品久久久久久久毛片微露脸| 内地一区二区视频在线| 午夜免费成人在线视频| 男人的好看免费观看在线视频| 成人精品一区二区免费| 亚洲国产欧洲综合997久久,| 国产成人啪精品午夜网站| 伊人久久大香线蕉亚洲五| 成人无遮挡网站| 在线观看舔阴道视频| 亚洲av美国av| 国产高清videossex| av中文乱码字幕在线| 青草久久国产| 久久久久国产精品人妻aⅴ院| 五月玫瑰六月丁香| 国产精品永久免费网站| 国产一区在线观看成人免费| 欧美黄色淫秽网站| 国内少妇人妻偷人精品xxx网站| 欧美zozozo另类| av专区在线播放| 久久久成人免费电影| 国产97色在线日韩免费| 成人国产一区最新在线观看| 桃红色精品国产亚洲av| 国产黄片美女视频| 亚洲av熟女| 真人做人爱边吃奶动态| 欧美黄色片欧美黄色片| 国产精品 欧美亚洲| 亚洲美女黄片视频| 欧美色欧美亚洲另类二区| 欧美日韩综合久久久久久 | av专区在线播放| 亚洲国产高清在线一区二区三| eeuss影院久久| 久久久久九九精品影院| 免费在线观看日本一区| 女同久久另类99精品国产91| 国产国拍精品亚洲av在线观看 | 欧美日韩国产亚洲二区| 欧美午夜高清在线| 18禁黄网站禁片午夜丰满| 欧美乱色亚洲激情| 午夜免费观看网址| 麻豆国产av国片精品| 亚洲国产日韩欧美精品在线观看 | 男女视频在线观看网站免费| 色综合欧美亚洲国产小说| 欧美+日韩+精品| av欧美777| 亚洲精华国产精华精| 真实男女啪啪啪动态图| 午夜福利高清视频| 99国产综合亚洲精品| 一进一出抽搐动态| 亚洲国产日韩欧美精品在线观看 | 热99re8久久精品国产| 一个人看视频在线观看www免费 | 特大巨黑吊av在线直播| 国产精品亚洲一级av第二区| 欧美午夜高清在线| 一级作爱视频免费观看| 亚洲午夜理论影院| 搡老岳熟女国产| 国产麻豆成人av免费视频| 国产视频一区二区在线看| 亚洲aⅴ乱码一区二区在线播放| 国产av麻豆久久久久久久| 97超视频在线观看视频| 在线播放国产精品三级| 丰满的人妻完整版| 亚洲第一电影网av| 色噜噜av男人的天堂激情| 亚洲人成伊人成综合网2020| 日韩欧美国产一区二区入口| 九九热线精品视视频播放| 亚洲在线自拍视频| 波多野结衣高清作品| 国产av一区在线观看免费| 午夜福利欧美成人| 一卡2卡三卡四卡精品乱码亚洲| 日韩免费av在线播放| 青草久久国产| 精品一区二区三区视频在线观看免费| 最近最新免费中文字幕在线| 成人国产综合亚洲| 免费人成视频x8x8入口观看| 激情在线观看视频在线高清| 久久久久久久精品吃奶| 成年版毛片免费区| 欧美在线一区亚洲| 色老头精品视频在线观看| 日本免费a在线| 91在线观看av| 黄色视频,在线免费观看| 成人精品一区二区免费| 露出奶头的视频| 欧美最黄视频在线播放免费| 色吧在线观看| 精品久久久久久久久久久久久| 91av网一区二区| 最新在线观看一区二区三区| 2021天堂中文幕一二区在线观| 国产淫片久久久久久久久 | 在线观看66精品国产| 狠狠狠狠99中文字幕| 搡老熟女国产l中国老女人| 欧美性猛交黑人性爽| 脱女人内裤的视频| 国产亚洲欧美在线一区二区| 看黄色毛片网站| 久久香蕉国产精品| 看黄色毛片网站| 久久精品综合一区二区三区| 欧美最黄视频在线播放免费| 国产精品日韩av在线免费观看| 淫妇啪啪啪对白视频| 91麻豆精品激情在线观看国产| 色av中文字幕| 国产精品av视频在线免费观看| 国产69精品久久久久777片| 国产男靠女视频免费网站| 三级毛片av免费| 精品电影一区二区在线| 在线免费观看不下载黄p国产 | 小说图片视频综合网站| 亚洲五月婷婷丁香| 动漫黄色视频在线观看| 麻豆成人午夜福利视频| 午夜福利高清视频| 很黄的视频免费| 亚洲精品在线观看二区| 欧美黄色淫秽网站| 亚洲成人精品中文字幕电影| 国产精品久久久久久久久免 | 亚洲激情在线av| 日本在线视频免费播放| 亚洲,欧美精品.| 两个人看的免费小视频| 国产av不卡久久| 免费在线观看日本一区| 久久久久久国产a免费观看| 国产激情偷乱视频一区二区| 中文字幕熟女人妻在线| 国产探花极品一区二区| 日韩av在线大香蕉| 国产精品亚洲一级av第二区| tocl精华| 好男人电影高清在线观看| 中文字幕精品亚洲无线码一区| 日韩中文字幕欧美一区二区| 国产欧美日韩一区二区精品| 91字幕亚洲| 亚洲人成伊人成综合网2020| 少妇人妻一区二区三区视频| 欧美日韩亚洲国产一区二区在线观看| 中文字幕人成人乱码亚洲影| 欧洲精品卡2卡3卡4卡5卡区| 国产精品,欧美在线| 99热这里只有是精品50| 有码 亚洲区| 一级毛片高清免费大全| 免费看美女性在线毛片视频| 欧美乱色亚洲激情| 91久久精品国产一区二区成人 | 中文字幕熟女人妻在线| 久久精品国产亚洲av涩爱 | 99久久综合精品五月天人人| 久久婷婷人人爽人人干人人爱| 老司机福利观看| av欧美777| www.熟女人妻精品国产| 99精品欧美一区二区三区四区| 国产高清videossex| 国产精品,欧美在线| 国产成人a区在线观看| 免费高清视频大片| 亚洲午夜理论影院| 亚洲av熟女| 亚洲精品一卡2卡三卡4卡5卡| 老鸭窝网址在线观看| 又紧又爽又黄一区二区| 真人一进一出gif抽搐免费| 久久久久久九九精品二区国产| 成人av在线播放网站| 精品人妻偷拍中文字幕| 国产69精品久久久久777片| www日本黄色视频网| 欧美另类亚洲清纯唯美| 免费电影在线观看免费观看| 美女被艹到高潮喷水动态| 噜噜噜噜噜久久久久久91| 国产熟女xx| 国产高清有码在线观看视频| 在线观看66精品国产| 毛片女人毛片| 国产精品影院久久| 欧美日韩亚洲国产一区二区在线观看| 搡女人真爽免费视频火全软件 | av在线天堂中文字幕| 精品久久久久久久人妻蜜臀av| www.999成人在线观看| 日日摸夜夜添夜夜添小说| 乱人视频在线观看| 亚洲无线观看免费| 男女视频在线观看网站免费| 熟女电影av网| 国产美女午夜福利| 97人妻精品一区二区三区麻豆| 欧美成人免费av一区二区三区| 3wmmmm亚洲av在线观看| 18禁黄网站禁片午夜丰满| 婷婷精品国产亚洲av在线| 蜜桃亚洲精品一区二区三区| 男女下面进入的视频免费午夜| 精品久久久久久成人av| 久久99热这里只有精品18| 亚洲在线自拍视频| 嫩草影院精品99| 国产探花极品一区二区| 性色avwww在线观看| 国产高清视频在线播放一区| 九色国产91popny在线| 亚洲国产欧洲综合997久久,| xxx96com| 亚洲激情在线av| 国产av一区在线观看免费| 网址你懂的国产日韩在线| 中文亚洲av片在线观看爽| 51午夜福利影视在线观看| 国产成人av激情在线播放| 美女黄网站色视频| 国产aⅴ精品一区二区三区波| 99国产极品粉嫩在线观看| tocl精华| 国产精品一区二区三区四区免费观看 | 欧美国产日韩亚洲一区| 日本黄色视频三级网站网址| 国产三级黄色录像| 一本精品99久久精品77| www日本黄色视频网| 色视频www国产| 亚洲 国产 在线| 真人一进一出gif抽搐免费| 99视频精品全部免费 在线| 天天躁日日操中文字幕| 在线十欧美十亚洲十日本专区| 国产三级黄色录像| 啦啦啦观看免费观看视频高清| 国产免费av片在线观看野外av| aaaaa片日本免费| 好男人电影高清在线观看| 精品欧美国产一区二区三| 精品不卡国产一区二区三区| 亚洲国产精品成人综合色| 国产精品乱码一区二三区的特点| 丰满乱子伦码专区| 日日夜夜操网爽| 一二三四社区在线视频社区8| 国产精品 欧美亚洲| 国产一区二区在线观看日韩 | 国产视频一区二区在线看| 可以在线观看的亚洲视频| 五月伊人婷婷丁香| 日韩大尺度精品在线看网址| 免费大片18禁| 日本黄色视频三级网站网址| 成年免费大片在线观看| 91麻豆精品激情在线观看国产| 精品国内亚洲2022精品成人| 一个人免费在线观看电影| 国产视频一区二区在线看| 狂野欧美白嫩少妇大欣赏| www日本黄色视频网| 午夜精品一区二区三区免费看| 亚洲av电影在线进入| 亚洲av二区三区四区| 午夜精品久久久久久毛片777| 国产av在哪里看| 国产av一区在线观看免费| 欧美一区二区国产精品久久精品| 免费看光身美女| 午夜精品在线福利| 亚洲成人中文字幕在线播放| 午夜精品一区二区三区免费看| 一个人免费在线观看的高清视频| 国产探花极品一区二区| 男人和女人高潮做爰伦理| 成人三级黄色视频| 日韩大尺度精品在线看网址| 两个人视频免费观看高清| av在线蜜桃| 亚洲欧美精品综合久久99| 亚洲熟妇中文字幕五十中出| 在线十欧美十亚洲十日本专区| 99riav亚洲国产免费| 精品人妻1区二区| 在线播放国产精品三级| 最近视频中文字幕2019在线8| 手机成人av网站| 亚洲国产色片| 一a级毛片在线观看| av专区在线播放| 国产精品久久久久久人妻精品电影| 全区人妻精品视频| 欧美丝袜亚洲另类 | 波多野结衣高清无吗| 高清在线国产一区| 嫩草影院精品99| 亚洲人成网站在线播放欧美日韩| 午夜激情福利司机影院| 午夜免费观看网址| 成人永久免费在线观看视频| 久久午夜亚洲精品久久| www日本黄色视频网| 日日夜夜操网爽| 日韩av在线大香蕉| 97碰自拍视频| 国产综合懂色| 一夜夜www| 成人av一区二区三区在线看| 亚洲人成网站在线播| 国产一区二区激情短视频| 亚洲av二区三区四区| 免费在线观看成人毛片| 日韩欧美三级三区| 欧美精品啪啪一区二区三区| 久99久视频精品免费| 搡女人真爽免费视频火全软件 | 亚洲,欧美精品.| 老司机深夜福利视频在线观看| 黄片大片在线免费观看| 亚洲精品一区av在线观看| 国产高清三级在线| 99国产综合亚洲精品| 亚洲中文字幕一区二区三区有码在线看| 色视频www国产| 十八禁网站免费在线| 久久久久国产精品人妻aⅴ院| 1024手机看黄色片| avwww免费| 在线国产一区二区在线| 亚洲五月婷婷丁香| 日韩 欧美 亚洲 中文字幕| 欧美最黄视频在线播放免费| 欧美成人a在线观看| 亚洲精品一区av在线观看| 久久久精品大字幕| 欧美在线一区亚洲| 人人妻,人人澡人人爽秒播| 国产在线精品亚洲第一网站| av视频在线观看入口| 97碰自拍视频| 99热这里只有是精品50| 欧美在线一区亚洲| 国产av在哪里看| 国产色婷婷99| 精品久久久久久久末码| 久久国产乱子伦精品免费另类| xxx96com| 亚洲黑人精品在线| av天堂在线播放| 国产色婷婷99| 久久天躁狠狠躁夜夜2o2o| 日本五十路高清| 亚洲欧美日韩高清在线视频| 欧美在线一区亚洲| 国产男靠女视频免费网站| 人妻丰满熟妇av一区二区三区| 欧美成人性av电影在线观看| 伊人久久大香线蕉亚洲五| 国产又黄又爽又无遮挡在线| 91在线精品国自产拍蜜月 | 欧美一级毛片孕妇| 色综合站精品国产| 全区人妻精品视频| 一卡2卡三卡四卡精品乱码亚洲| 成年女人永久免费观看视频| 亚洲中文字幕日韩| 欧美极品一区二区三区四区| 黄色视频,在线免费观看| 12—13女人毛片做爰片一| 天堂影院成人在线观看|