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

    井內(nèi)嵌入同軸多孔與固體圓柱地下水滲流流型分析

    2021-07-23 06:13:08雷海燕戴傳山
    水文地質(zhì)工程地質(zhì) 2021年4期
    關(guān)鍵詞:水環(huán)流型滲透系數(shù)

    朱 琳,雷海燕,馬 非,戴傳山

    (天津大學(xué)中低溫?zé)崮芨咝Ю媒逃恐攸c(diǎn)實(shí)驗(yàn)室,天津 300350)

    垂直于圓管外的繞流現(xiàn)象是自然界較普遍存在的流體力學(xué)現(xiàn)象之一,如卡門渦街。在多孔介質(zhì)滲流層內(nèi)垂直于圓管外的繞流問題也是石油鉆探、地下水開采[1]、地?zé)崮荛_發(fā)、流化床等許多工業(yè)領(lǐng)域經(jīng)常遇到的物理現(xiàn)象,如地?zé)峋才c熱儲層的耦合傳熱問題[2],甚至在臨床醫(yī)學(xué)研究中也有相關(guān)的微血管外的間質(zhì)滲流研究[3]。在飽和多孔介質(zhì)內(nèi)嵌入物體,被嵌物體外的繞流、傳熱傳質(zhì)問題在很多工業(yè)應(yīng)用領(lǐng)域中有廣泛的應(yīng)用。針對多孔介質(zhì)內(nèi)嵌有簡單幾何形狀的物體結(jié)構(gòu),若能獲得其物體外部的繞流與熱質(zhì)擴(kuò)散的問題的簡單表達(dá)形式的理論解析解,對指導(dǎo)類似結(jié)構(gòu)下的工業(yè)設(shè)計(jì)有重要的作用。

    Pop等[4]首次獲得了考慮流體黏性耗散Brinkman模型圓柱外繞流的理論解析解。相對Pop等[4]研究的有滲流通過的多孔介質(zhì)內(nèi)嵌入圓柱的流場結(jié)構(gòu)而言,多孔介質(zhì)內(nèi)開挖有一圓形孔洞,且孔洞內(nèi)充滿純流體的空間結(jié)構(gòu)下的流場理論解析解獲得的比較早。Drost等[5]采用忽略慣性力僅考慮流體黏性力的Stokes流模型獲得了無限大多孔介質(zhì)內(nèi)有一圓形孔洞的流場理論解析解,得出孔內(nèi)平均流速與地層遠(yuǎn)端流速相比(即平均收斂因子)約等于2。Sano[6]進(jìn)一步分析了圓形孔洞的存在對周圍達(dá)西流的影響,得出在多孔介質(zhì)滲透系數(shù)較低情況下,孔洞圓心處速度約是遠(yuǎn)端地層滲流速度的3倍等重要結(jié)論。這一結(jié)論隨后得到了Momii等[7]在實(shí)驗(yàn)室內(nèi)基于激光多普勒方法進(jìn)行的實(shí)驗(yàn)驗(yàn)證。Raja等[8?10]進(jìn)一步推導(dǎo)了有均勻滲流通過的不同滲透系數(shù)下多孔介質(zhì)層內(nèi)有一規(guī)則或微變形的圓形孔洞或球形孔洞結(jié)構(gòu)下的流場內(nèi)流體速度(即收斂因子α)分布理論解。

    Li等[11]針對一種相對前人研究更為復(fù)雜的有滲流通過的無限大多孔介質(zhì)內(nèi)嵌一固體圓柱,且圓柱外包裹有水環(huán)的幾何結(jié)構(gòu)下的各區(qū)域流型結(jié)構(gòu)進(jìn)行了研究,其研究目的是獲取多孔介質(zhì)內(nèi)嵌入細(xì)長吸收碳纖維捕捉流體內(nèi)膠體顆粒物的機(jī)理。Sekhar[12]將流體環(huán)內(nèi)固體替換為多孔介質(zhì),得出該幾何結(jié)構(gòu)下的流動(dòng)問題理論解,分析了內(nèi)、外多孔介質(zhì)在不同滲流系數(shù)情況下,整體多孔介質(zhì)表現(xiàn)出的宏觀平均水文特征、局部速度分布以及流動(dòng)阻力特性。

    前人研究結(jié)果表明,有均勻滲流通過的多孔介質(zhì)內(nèi)嵌入固體圓柱、多孔層、有水環(huán)或無水環(huán)等復(fù)雜幾何結(jié)構(gòu)下流場的理論解析研究文獻(xiàn)并不多見,應(yīng)用領(lǐng)域或研究目的卻比較廣泛,如復(fù)合材料封井、井內(nèi)流速觀測[13]、控制污染物在地下水和土壤中的擴(kuò)散[14]等。值得注意的是,前人最多只對3層結(jié)構(gòu),即無限大多孔介質(zhì)-水環(huán)-環(huán)狀多孔層/圓形固體進(jìn)行了研究。本文借鑒以往的研究方法,對有滲流通過的無限大多孔介質(zhì)流場內(nèi)嵌入固體圓柱,且圓柱外部順序包裹有多孔介質(zhì)環(huán)、純流體環(huán)的復(fù)雜幾何結(jié)構(gòu)下的流場速度解析解進(jìn)行了理論分析,著重研究了不同滲透系數(shù)下的內(nèi)外多孔介質(zhì)、水環(huán)間隙情況等對圓柱外繞流流型以及縱橫界面上的速度分布的影響規(guī)律。由于地下水滲流速度較小,以往研究中常不考慮由彌散效應(yīng)引起的黏度系數(shù)的變化,即多孔介質(zhì)內(nèi)黏度系數(shù)與純流體內(nèi)的黏度系數(shù)視為相同[15]。此外,在井筒附近的地下水滲流可能會產(chǎn)生顆粒的沉積堵塞[16],但本研究忽略這一影響,認(rèn)為局部滲透系數(shù)不變。本文研究結(jié)果對類似結(jié)構(gòu)下的地埋管套管換熱器、地下水污染物吸收裝置設(shè)計(jì)、地下水測速裝置研發(fā)等有一定的理論指導(dǎo)意義。

    1 物理模型的建立與求解方法

    考慮無限大2D均質(zhì)多孔介質(zhì)(Pe)內(nèi)有單向滲流,垂直井筒內(nèi)另嵌入同軸均質(zhì)多孔介質(zhì)(Pi)與固體圓柱(S),在內(nèi)外多孔介質(zhì)之間有開放水環(huán)(Wen)的情況,假設(shè)多孔介質(zhì)均為各向同性材料:滲透系數(shù)為ke的無限大多孔介質(zhì)內(nèi)有黏性流體穩(wěn)定通過,繞過井中一包裹有滲透系數(shù)為ki的多孔層且半徑為r0的圓形固體(圖1)。水環(huán)區(qū)域(λ1r0

    圖1 物理模型示意圖Fig.1 Schematic diagram of the physical model.

    水環(huán)內(nèi)、外多孔介質(zhì)區(qū)由Brinkman方程[18]和連續(xù)性方程控制:

    式中:μeあ—多孔介質(zhì)內(nèi)流體有效黏度系數(shù);

    μen—水環(huán)內(nèi)純流體黏度系數(shù),μeあ=μen=μ;

    en—水環(huán);

    l=i、e—內(nèi)部、外無限大多孔介質(zhì)區(qū)。

    速度矢量(V)、橫縱坐標(biāo)(x、y)、各點(diǎn)到固體圓柱圓心半徑(r)和壓力(P),進(jìn)行無量綱化:

    式中:U∞—無限遠(yuǎn)處的流速大小。

    水環(huán)區(qū)控制式(1)(2)無量綱化:

    多孔區(qū)控制式(3)~(4)無量綱化:

    式中:χ2—滲透性能,χ2=r20/k。

    由速度分量和剪切應(yīng)力(T)的連續(xù)性[19?20],得到邊界條件(每個(gè)邊界條件都可以分解為切向和法向,切向剪切應(yīng)力包含壓力P):

    引入極坐標(biāo)系下的流函數(shù)(ψ):

    式中:θ—速度分量Vθ與橫坐標(biāo)軸x夾角。

    式(6)—(8)簡化為:

    結(jié)合無窮遠(yuǎn)處均勻流動(dòng)條件ψe∞→?rsinθ,各區(qū)域式(13)(14)的解為:

    式中:I1—第一類修正一階貝塞爾函數(shù);

    K1—第二類修正一階貝塞爾函數(shù)。

    其中系數(shù)A、B、C、D、E、F、G、H、M、N可由式(9)—(11)確定。

    值得注意的是,若r0→0,解的形式與Sekhar[12]一致:

    此外,若r0→0,λ1→0,解的形式為:

    笛卡爾坐標(biāo)系下的流速U:

    2 方法驗(yàn)證

    為了驗(yàn)證計(jì)算方法的正確性,本文以Raja等[8]的計(jì)算模型為例,進(jìn)行了程序正確性驗(yàn)證,該模型中包括純流體孔洞部分的Stokes流模型以及外部多孔介質(zhì)部分流動(dòng)的Brinkman模型。輸入?yún)?shù)與文獻(xiàn)[8]一致,結(jié)果見圖2。

    圖2 χe變化下縱向中心剖面(Y=0),橫向中心剖面(X=0)無量綱速度分布圖Fig.2 Dimensionless velocity distribution map of the longitudinal center profile(Y=0) and the lateral center profile(X=0) under the changes of χe

    該計(jì)算模型對應(yīng)于無限大多孔介質(zhì)內(nèi)有一孔洞的情況,即r0→0,λ1→0。其中,橫縱坐標(biāo)無量綱化:χe變化下經(jīng)過圓心的縱、橫剖面速度分布見圖2。值得注意的是,無量綱化特征長度為孔洞半徑r0c=λ2r0。隨著χe的增加(滲透系數(shù)變小),孔洞的存在感變強(qiáng),中心流速與遠(yuǎn)端流速比值αc變大到3,若χe為0(滲透系數(shù)極大),則αc→1。本文結(jié)果與Raja等[8]研究結(jié)果符合非常好。

    3 多層結(jié)構(gòu)計(jì)算結(jié)果

    在方法驗(yàn)證的基礎(chǔ)上,對圖1所示的多層結(jié)構(gòu)模型進(jìn)行理論求解。由于影響參數(shù)比較多,僅對不同半徑、以及內(nèi)外滲透系數(shù)變化情況下的流型結(jié)構(gòu)求解。計(jì)算工況見表1。值得注意的是,為了更清晰地表示外部尺寸固定(即λ1/λ2為定值)情況下,內(nèi)部圓柱尺寸的變化對于流場的影響,此處給出的r0是一組以r0=r02為標(biāo)準(zhǔn)的假設(shè)值(流場的變化只與λ1、λ2、χi、χe有關(guān))。同理,整個(gè)流場均如式(5)所述,計(jì)算過程中使用相應(yīng)的r0進(jìn)行參數(shù)的無量綱化為了更加直觀地看出不同參數(shù)的影響,在計(jì)算結(jié)果的整理中,橫縱中心剖面坐標(biāo)軸的無量綱化均以r02半徑為準(zhǔn):

    表1 不同工況的計(jì)算參數(shù)Table 1 Summary of calculation conditions under different parameters.

    3.1 固體圓柱半徑大小r0對流型結(jié)構(gòu)的影響

    圖3(a)(b)可以看出,固體半徑r0的變化對水環(huán)外的多孔區(qū)幾乎沒有影響。但隨著r0增大,水環(huán)內(nèi)流速峰值變高,這是由于流體經(jīng)過包裹有多孔層的固體結(jié)構(gòu)時(shí)流動(dòng)阻力變大(圖4)。此外,r0變大,固體圓柱在多孔介質(zhì)層中對流體的影響半徑增大,逐漸呈繞流趨勢,如圖5(a)(b)。更進(jìn)一步地,水環(huán)內(nèi)的多孔區(qū)縱向剖面的平均流速值變小,橫向剖面上的速度變大。這是因?yàn)樵搮^(qū)域內(nèi)部阻力變大,流入的流量越來越少,但繞流增加了橫向剖面上的流速。另外,在r0極大的情況下(工況3),多孔層內(nèi)流速急速下降到0,這是由于極薄多孔層對流場沒有太大影響,如圖6(a)。

    圖3 縱向、橫向中心剖面無量綱速度分布圖Fig.3 Dimensionless velocity distribution map of the longitudinal center profile and the lateral center profile

    圖4 典型工況下的壓力云圖(工況2)Fig.4 Pressure contours under typical working condition 2

    圖5 典型工況下的流線圖Fig.5 Streamline diagram under the typical working conditions

    以工況2為對比工況,分別研究參數(shù)λ1、λ2、χi、χe的變化對流型結(jié)構(gòu)的影響。此外,由于繞流,流體在水環(huán)內(nèi)橫向剖面上的流速變化更為明顯,因此水環(huán)內(nèi) 流型變化主要用該剖面上的流速峰值分析。

    3.2 λ1變化對流型結(jié)構(gòu)的影響

    圖3(c)(d)為λ1變化下的中心剖面無量綱速度分布圖。λ1的變化對水環(huán)外的多孔區(qū)(下文簡稱外多孔區(qū))幾乎沒有影響。但是,隨著λ1變小,水環(huán)內(nèi)流速峰值先變高后變低,這是由于縫隙變大,進(jìn)而剖面上速度分布展寬變大。即進(jìn)入水環(huán)內(nèi)多孔區(qū)(下文簡稱內(nèi)多孔區(qū))的流體變少,而水環(huán)內(nèi)流量增加流速變大,直至流體幾乎不通過內(nèi)多孔區(qū)(工況5),見圖6(b),且過程中隨著展寬變大,水環(huán)內(nèi)流量分配空間變大導(dǎo)致其流速峰值逐漸降低。類似地,由于縫隙變大,內(nèi)多孔區(qū)阻力相對水環(huán)阻力變大,該區(qū)域流量降低流速變小,因此流場也更容易平穩(wěn)。同時(shí),λ1的減小會使內(nèi)多孔層中的繞流效果相對增強(qiáng),因此該區(qū)域流型平穩(wěn)程度要耦合區(qū)域內(nèi)流量大小和圓柱外繞流強(qiáng)度兩種因素。

    圖6 典型工況下的速度云圖Fig.6 Speed contours under the typical working conditions

    3.3 λ2變化對流型結(jié)構(gòu)的影響

    圖3(e)(f)為λ2變化下的中心剖面無量綱速度分布圖。當(dāng)λ2很大時(shí),縫隙很大,極少流體流入內(nèi)多孔區(qū)(工況6)。隨著λ2的增大,水環(huán)內(nèi)流速峰值先升高后降低,這是由于縫隙變大流速分布展寬變大,流體越來越少地進(jìn)入內(nèi)多孔層且水環(huán)內(nèi)流量分配空間變大導(dǎo)致的。由于內(nèi)多孔區(qū)阻力相對水環(huán)阻力增大,該區(qū)域平均流速成倍降低且流型變化趨勢更平穩(wěn)。λ2減小,意味著整個(gè)障礙物尺寸相對減小,因此,外部多孔介質(zhì)區(qū)的流型會更加平緩。

    3.4 χi變化對流型結(jié)構(gòu)的影響

    圖3(g)(h)為χi變化下的中心剖面無量綱速度分布圖。χi的變化對外多孔區(qū)的流場幾乎沒有影響。但是,當(dāng)χi極大時(shí),流體幾乎不進(jìn)入內(nèi)多孔區(qū)。隨著χi的減小,水環(huán)內(nèi)流速峰值降低。這是因?yàn)閮?nèi)多孔區(qū)滲透性加強(qiáng),內(nèi)多孔層阻力相對水環(huán)阻力變小,內(nèi)多孔區(qū)流量變大水環(huán)內(nèi)流量減小。同理,內(nèi)多孔區(qū)橫向和縱向剖面流速均相應(yīng)變大,且變化明顯。當(dāng)內(nèi)多孔區(qū)滲透系數(shù)極大時(shí)(工況9),見圖6(c),水環(huán)與內(nèi)多孔層之間的速度階躍消失,出現(xiàn)“穿透”現(xiàn)象,這是因?yàn)閮?nèi)多孔層阻力與水環(huán)阻力同量級。

    3.5 χe變化對流型結(jié)構(gòu)的影響

    圖3(i)(j)為χe變化下的中心剖面無量綱速度分布圖。隨著χe變小水環(huán)峰值降低,內(nèi)多孔區(qū)流速越低流型越平緩,但外多孔區(qū)流速增大。這是由于內(nèi)外多孔區(qū)滲透系數(shù)差變小,整體結(jié)構(gòu)流型過渡更為平緩,水環(huán)內(nèi)的加速效應(yīng)相應(yīng)減弱,內(nèi)多孔區(qū)的流速相應(yīng)降低。χe≤10時(shí)外多孔區(qū)流速變化趨勢與χe較大時(shí)相反,這是由于與流體向水環(huán)內(nèi)的滲透作用相比,水環(huán)外的繞流作用占了主導(dǎo)地位,見圖5(c),流體極少進(jìn)入內(nèi)多孔區(qū)。

    4 結(jié)論

    (1)水環(huán)內(nèi)部參數(shù)如r0、λ1、χi的變化對外層多孔介質(zhì)區(qū)影響很??;水環(huán)直徑的減小會使外部流型更為平穩(wěn);水環(huán)外部流型主要受控于外部滲透系數(shù),但由于障礙物的存在,滲透系數(shù)對外部流型的影響是非單調(diào)的函數(shù)。

    (2)橫向剖面上水環(huán)內(nèi)的速度變化較大,固體圓柱r0減小,水環(huán)內(nèi)峰值增加。內(nèi)部與外部滲透系數(shù)同時(shí)增大,水環(huán)內(nèi)加速效應(yīng)減弱。水環(huán)間隙對水環(huán)內(nèi)流速峰值影響較大。

    (3)研究發(fā)現(xiàn)內(nèi)部滲透系數(shù)增加到某一臨界值情況下會出現(xiàn)“穿透”現(xiàn)象,橫向速度剖面從階梯形變?yōu)閽佄锞€形。

    猜你喜歡
    水環(huán)流型滲透系數(shù)
    參觀水環(huán)中心
    水平井油水兩相流型實(shí)驗(yàn)研究
    云南化工(2021年7期)2021-12-21 07:27:32
    基于Origin的滲透系數(shù)衰減方程在地?zé)崴毓嘀械膽?yīng)用
    多孔材料水滲透系數(shù)預(yù)測的隨機(jī)行走法
    輸水渠防滲墻及基巖滲透系數(shù)敏感性分析
    河北平原新近系熱儲層滲透系數(shù)規(guī)律性分析
    共流型轉(zhuǎn)子的有限元分析
    基于Taitel-Dukler方法的氣液兩相流型邊界計(jì)算軟件開發(fā)
    水環(huán)壓縮機(jī)組在乙炔生產(chǎn)中的應(yīng)用
    中國氯堿(2014年10期)2014-02-28 01:04:59
    水環(huán)真空泵故障分析及處理方案
    河南科技(2014年23期)2014-02-27 14:18:55
    香蕉丝袜av| 九色亚洲精品在线播放| 免费女性裸体啪啪无遮挡网站| 免费在线观看视频国产中文字幕亚洲 | 久久久久精品性色| 亚洲av欧美aⅴ国产| 国产精品嫩草影院av在线观看| 一级片免费观看大全| 国产无遮挡羞羞视频在线观看| av网站在线播放免费| 国产男女超爽视频在线观看| av电影中文网址| 日本午夜av视频| 欧美 日韩 精品 国产| 亚洲国产精品一区二区三区在线| 欧美成人精品欧美一级黄| 久久99蜜桃精品久久| 欧美 日韩 精品 国产| 最近最新中文字幕免费大全7| 日韩视频在线欧美| 色播在线永久视频| 夜夜骑夜夜射夜夜干| 一级片免费观看大全| 妹子高潮喷水视频| 男人操女人黄网站| 亚洲精品久久成人aⅴ小说| 免费av中文字幕在线| 国产在线视频一区二区| 久久久精品区二区三区| 久久这里只有精品19| 久久久国产精品麻豆| 狂野欧美激情性bbbbbb| 婷婷色综合大香蕉| 一级爰片在线观看| 男人爽女人下面视频在线观看| 爱豆传媒免费全集在线观看| 成年人免费黄色播放视频| 丝袜脚勾引网站| 99热国产这里只有精品6| 日韩,欧美,国产一区二区三区| 久久久久久人人人人人| 热99国产精品久久久久久7| 一区在线观看完整版| 看免费成人av毛片| a级毛片黄视频| 伦理电影免费视频| 精品国产露脸久久av麻豆| 国产乱人偷精品视频| 国产精品一区二区在线不卡| 肉色欧美久久久久久久蜜桃| 啦啦啦啦在线视频资源| 91在线精品国自产拍蜜月| 欧美精品一区二区大全| 久久久久久久亚洲中文字幕| 少妇人妻 视频| 欧美日韩亚洲高清精品| 久久久久视频综合| 国产av精品麻豆| 午夜福利乱码中文字幕| 国产精品久久久久久精品古装| 欧美精品av麻豆av| 国产成人91sexporn| 狠狠婷婷综合久久久久久88av| 亚洲少妇的诱惑av| www日本在线高清视频| 亚洲内射少妇av| 哪个播放器可以免费观看大片| 欧美激情极品国产一区二区三区| 成人二区视频| 超碰成人久久| 国产精品欧美亚洲77777| √禁漫天堂资源中文www| 日韩伦理黄色片| 亚洲精品一二三| 777久久人妻少妇嫩草av网站| 精品午夜福利在线看| 国产成人精品久久二区二区91 | 十分钟在线观看高清视频www| 国产av精品麻豆| 亚洲精品一二三| 青春草视频在线免费观看| 女的被弄到高潮叫床怎么办| 一二三四中文在线观看免费高清| 99久久中文字幕三级久久日本| 欧美精品一区二区免费开放| 午夜日本视频在线| 99精国产麻豆久久婷婷| 亚洲在久久综合| 日日啪夜夜爽| 亚洲欧美一区二区三区国产| 亚洲,一卡二卡三卡| 中文精品一卡2卡3卡4更新| 1024视频免费在线观看| 久久久国产精品麻豆| 免费看不卡的av| 国产亚洲精品第一综合不卡| 欧美 日韩 精品 国产| 下体分泌物呈黄色| 国产av国产精品国产| 熟女少妇亚洲综合色aaa.| 国产精品香港三级国产av潘金莲 | 日本免费在线观看一区| 十八禁网站网址无遮挡| 亚洲三级黄色毛片| 亚洲国产毛片av蜜桃av| 黑人巨大精品欧美一区二区蜜桃| 亚洲av综合色区一区| 免费看av在线观看网站| 国产 精品1| 黄片播放在线免费| 日韩中文字幕视频在线看片| 日韩av不卡免费在线播放| 亚洲欧美成人精品一区二区| 久久久久久久国产电影| 日本欧美国产在线视频| 精品第一国产精品| 日韩一卡2卡3卡4卡2021年| 亚洲精品日韩在线中文字幕| 一级毛片黄色毛片免费观看视频| 十八禁网站网址无遮挡| 成年女人在线观看亚洲视频| 精品99又大又爽又粗少妇毛片| 人妻系列 视频| 在线免费观看不下载黄p国产| 曰老女人黄片| 99九九在线精品视频| 亚洲色图综合在线观看| 亚洲内射少妇av| 美女国产视频在线观看| 美女大奶头黄色视频| 丝袜美腿诱惑在线| 亚洲中文av在线| 18禁国产床啪视频网站| 波多野结衣av一区二区av| 女的被弄到高潮叫床怎么办| 日韩中文字幕欧美一区二区 | 1024香蕉在线观看| 精品亚洲乱码少妇综合久久| 国产精品国产三级国产专区5o| 国产老妇伦熟女老妇高清| 久久精品久久精品一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 亚洲一码二码三码区别大吗| av在线app专区| av电影中文网址| 亚洲精品国产色婷婷电影| 男女边摸边吃奶| 两个人免费观看高清视频| 熟妇人妻不卡中文字幕| 在现免费观看毛片| 久久青草综合色| 久久久久久伊人网av| av网站免费在线观看视频| 欧美日韩综合久久久久久| 久久精品国产亚洲av涩爱| 黄色配什么色好看| 国产精品久久久av美女十八| 日韩中文字幕欧美一区二区 | 午夜免费鲁丝| 婷婷色综合大香蕉| 国产又爽黄色视频| 国产精品女同一区二区软件| 日韩精品有码人妻一区| 国产日韩欧美亚洲二区| 成人漫画全彩无遮挡| 热99久久久久精品小说推荐| 国产精品久久久久久久久免| 亚洲国产精品成人久久小说| 在线 av 中文字幕| 亚洲精品国产av成人精品| 久久99热这里只频精品6学生| 亚洲色图 男人天堂 中文字幕| 18在线观看网站| 男女边摸边吃奶| 久久久亚洲精品成人影院| 你懂的网址亚洲精品在线观看| 日韩,欧美,国产一区二区三区| a级片在线免费高清观看视频| 在线亚洲精品国产二区图片欧美| 亚洲欧美一区二区三区国产| 国产亚洲一区二区精品| 美女xxoo啪啪120秒动态图| 欧美日韩一级在线毛片| 18在线观看网站| 男女啪啪激烈高潮av片| 亚洲av电影在线进入| 在线天堂中文资源库| 宅男免费午夜| 在线观看免费视频网站a站| 丁香六月天网| 精品亚洲成国产av| 丝袜在线中文字幕| 秋霞在线观看毛片| 99热全是精品| 久久久久久免费高清国产稀缺| 亚洲成人一二三区av| 欧美精品亚洲一区二区| 久久久久网色| 色婷婷久久久亚洲欧美| 中文字幕精品免费在线观看视频| 久久精品国产综合久久久| 性色avwww在线观看| 咕卡用的链子| 啦啦啦在线免费观看视频4| 久久 成人 亚洲| 啦啦啦中文免费视频观看日本| 欧美精品av麻豆av| 欧美精品av麻豆av| 青青草视频在线视频观看| 最近中文字幕高清免费大全6| 色视频在线一区二区三区| 国产一区二区 视频在线| 欧美精品av麻豆av| 综合色丁香网| 男女边吃奶边做爰视频| 亚洲成国产人片在线观看| 可以免费在线观看a视频的电影网站 | 欧美国产精品一级二级三级| 2018国产大陆天天弄谢| 国产精品 欧美亚洲| 中文精品一卡2卡3卡4更新| 亚洲精品aⅴ在线观看| 国产亚洲av片在线观看秒播厂| 一区二区三区乱码不卡18| 亚洲欧洲国产日韩| 欧美av亚洲av综合av国产av | 久久人妻熟女aⅴ| 女人被躁到高潮嗷嗷叫费观| 丝袜美足系列| 亚洲精品乱久久久久久| 寂寞人妻少妇视频99o| 久久精品国产亚洲av天美| 亚洲精品av麻豆狂野| 街头女战士在线观看网站| av视频免费观看在线观看| 在线观看免费高清a一片| 2018国产大陆天天弄谢| 26uuu在线亚洲综合色| 十八禁网站网址无遮挡| 热re99久久精品国产66热6| 成人18禁高潮啪啪吃奶动态图| 亚洲av男天堂| 少妇人妻久久综合中文| 26uuu在线亚洲综合色| 国产97色在线日韩免费| 国产成人精品婷婷| 高清av免费在线| 久久久久人妻精品一区果冻| 观看美女的网站| xxxhd国产人妻xxx| 午夜91福利影院| 91精品伊人久久大香线蕉| 在线观看美女被高潮喷水网站| 久久综合国产亚洲精品| 午夜91福利影院| a 毛片基地| 超碰成人久久| 国产乱来视频区| 久久精品夜色国产| 香蕉精品网在线| 老女人水多毛片| 久久这里有精品视频免费| 伦精品一区二区三区| 国产在线免费精品| 亚洲成人av在线免费| 国产女主播在线喷水免费视频网站| 精品福利永久在线观看| 亚洲激情五月婷婷啪啪| 成人国产麻豆网| 亚洲综合色网址| 男女高潮啪啪啪动态图| 97在线视频观看| 亚洲综合色惰| 啦啦啦在线观看免费高清www| 国产精品欧美亚洲77777| 久久久久久久久免费视频了| 男男h啪啪无遮挡| 9191精品国产免费久久| 亚洲国产色片| 建设人人有责人人尽责人人享有的| 蜜桃国产av成人99| 天天躁狠狠躁夜夜躁狠狠躁| 十八禁高潮呻吟视频| 日韩在线高清观看一区二区三区| 九草在线视频观看| 电影成人av| 日本vs欧美在线观看视频| 熟女少妇亚洲综合色aaa.| 天美传媒精品一区二区| 国产片特级美女逼逼视频| 成年av动漫网址| 亚洲欧美成人精品一区二区| 永久网站在线| 又大又黄又爽视频免费| 一边亲一边摸免费视频| 亚洲av电影在线观看一区二区三区| 久久精品国产综合久久久| 韩国av在线不卡| 欧美在线黄色| 亚洲中文av在线| 亚洲人成电影观看| 国产一区二区在线观看av| 我的亚洲天堂| 9191精品国产免费久久| 韩国av在线不卡| 国产在线一区二区三区精| 国产精品国产av在线观看| 亚洲人成77777在线视频| 亚洲av中文av极速乱| 免费在线观看黄色视频的| 国产极品天堂在线| 丝袜美腿诱惑在线| 日本av手机在线免费观看| www.精华液| 亚洲精品自拍成人| 美女高潮到喷水免费观看| 男人添女人高潮全过程视频| 在线 av 中文字幕| 自拍欧美九色日韩亚洲蝌蚪91| 大话2 男鬼变身卡| 国产精品不卡视频一区二区| 日韩制服丝袜自拍偷拍| 国产成人a∨麻豆精品| 2018国产大陆天天弄谢| 免费观看a级毛片全部| 国产av精品麻豆| 老司机影院毛片| 狠狠婷婷综合久久久久久88av| 日韩av在线免费看完整版不卡| 国产精品国产三级专区第一集| 久久亚洲国产成人精品v| 欧美日本中文国产一区发布| www.av在线官网国产| av.在线天堂| 国产熟女欧美一区二区| 国产精品一国产av| 亚洲av男天堂| 我的亚洲天堂| 最近中文字幕高清免费大全6| 亚洲综合色网址| 色婷婷av一区二区三区视频| av片东京热男人的天堂| 国产一区二区三区综合在线观看| av线在线观看网站| 韩国精品一区二区三区| 激情视频va一区二区三区| 久久久久精品性色| 国产精品久久久久久av不卡| 中国国产av一级| 国产高清国产精品国产三级| 夫妻午夜视频| 欧美国产精品一级二级三级| 日韩av免费高清视频| 中国国产av一级| 看免费av毛片| av国产精品久久久久影院| 亚洲 欧美一区二区三区| 热re99久久国产66热| 日韩免费高清中文字幕av| 亚洲欧美精品综合一区二区三区 | 欧美激情 高清一区二区三区| 亚洲第一青青草原| 女性生殖器流出的白浆| 久久久亚洲精品成人影院| 国产熟女午夜一区二区三区| 美女国产视频在线观看| 欧美另类一区| 好男人视频免费观看在线| 水蜜桃什么品种好| 久久人人97超碰香蕉20202| 国产亚洲av片在线观看秒播厂| 1024视频免费在线观看| 少妇被粗大猛烈的视频| 日韩三级伦理在线观看| 国产精品免费大片| 精品少妇久久久久久888优播| 五月天丁香电影| 在线亚洲精品国产二区图片欧美| 国产色婷婷99| 亚洲国产日韩一区二区| 一级,二级,三级黄色视频| 亚洲美女视频黄频| 久久久久久伊人网av| 亚洲欧洲日产国产| 亚洲色图综合在线观看| 不卡av一区二区三区| 妹子高潮喷水视频| 少妇精品久久久久久久| 考比视频在线观看| 99国产精品免费福利视频| 肉色欧美久久久久久久蜜桃| 日韩伦理黄色片| 欧美日韩亚洲国产一区二区在线观看 | 在线 av 中文字幕| 国产精品久久久久久精品电影小说| 国产精品久久久久久精品古装| 午夜影院在线不卡| 久久亚洲国产成人精品v| 国产一区二区激情短视频 | 超碰成人久久| 国产精品久久久久久精品古装| 精品一区二区三卡| 岛国毛片在线播放| 麻豆av在线久日| 免费人妻精品一区二区三区视频| 大话2 男鬼变身卡| 免费观看在线日韩| 国产极品粉嫩免费观看在线| 日日撸夜夜添| 国产一区二区 视频在线| 2021少妇久久久久久久久久久| 久久久久久人人人人人| 亚洲视频免费观看视频| 在线天堂最新版资源| 大香蕉久久成人网| 国产精品免费大片| 久久 成人 亚洲| 亚洲欧美一区二区三区国产| 男女午夜视频在线观看| 深夜精品福利| 黄频高清免费视频| 亚洲美女视频黄频| 综合色丁香网| 久久久久精品性色| 精品亚洲成a人片在线观看| 一区二区日韩欧美中文字幕| 黑丝袜美女国产一区| 欧美精品一区二区大全| 亚洲第一区二区三区不卡| 成人亚洲欧美一区二区av| 亚洲精品自拍成人| 大片电影免费在线观看免费| 精品国产国语对白av| 在线观看一区二区三区激情| 你懂的网址亚洲精品在线观看| 波多野结衣一区麻豆| 久久久久精品久久久久真实原创| av片东京热男人的天堂| 国产亚洲最大av| 性色avwww在线观看| 黄色 视频免费看| 秋霞伦理黄片| 久久久国产精品麻豆| 亚洲激情五月婷婷啪啪| 久久久亚洲精品成人影院| 性少妇av在线| 观看美女的网站| 精品少妇一区二区三区视频日本电影 | 国产一区二区 视频在线| 99久国产av精品国产电影| 女人被躁到高潮嗷嗷叫费观| 人人妻人人爽人人添夜夜欢视频| 成人毛片60女人毛片免费| 亚洲欧美成人综合另类久久久| 午夜av观看不卡| 久久久久人妻精品一区果冻| 日韩中文字幕欧美一区二区 | 欧美中文综合在线视频| 最近的中文字幕免费完整| 成年人午夜在线观看视频| 又黄又粗又硬又大视频| 久久久久久人妻| 美女高潮到喷水免费观看| 亚洲色图 男人天堂 中文字幕| 在线观看国产h片| 国产精品国产三级国产专区5o| 久热这里只有精品99| 国产福利在线免费观看视频| 免费高清在线观看视频在线观看| 性少妇av在线| 午夜福利视频精品| 久久人人97超碰香蕉20202| 欧美日韩综合久久久久久| 波多野结衣一区麻豆| 成年女人毛片免费观看观看9 | 最近中文字幕2019免费版| 亚洲色图 男人天堂 中文字幕| 青春草亚洲视频在线观看| 777久久人妻少妇嫩草av网站| 久久婷婷青草| 伊人亚洲综合成人网| 精品午夜福利在线看| 亚洲综合精品二区| 成年av动漫网址| 黄色 视频免费看| 成年人免费黄色播放视频| 久久精品亚洲av国产电影网| 有码 亚洲区| 国产老妇伦熟女老妇高清| 国产又爽黄色视频| 国产有黄有色有爽视频| 成人国语在线视频| 亚洲精品国产av成人精品| av在线老鸭窝| 免费女性裸体啪啪无遮挡网站| 街头女战士在线观看网站| 久久综合国产亚洲精品| 国产精品麻豆人妻色哟哟久久| 欧美激情高清一区二区三区 | 日本免费在线观看一区| 久久久精品国产亚洲av高清涩受| 精品久久久精品久久久| 国产乱人偷精品视频| 成年人免费黄色播放视频| 亚洲综合色惰| 亚洲精品日本国产第一区| 少妇被粗大的猛进出69影院| 男女啪啪激烈高潮av片| 美女xxoo啪啪120秒动态图| 热re99久久精品国产66热6| 一级毛片黄色毛片免费观看视频| 精品久久蜜臀av无| 水蜜桃什么品种好| 黄片小视频在线播放| 国产精品偷伦视频观看了| 亚洲,一卡二卡三卡| 乱人伦中国视频| 黑丝袜美女国产一区| 可以免费在线观看a视频的电影网站 | 最近2019中文字幕mv第一页| 日韩一区二区视频免费看| 亚洲精品日韩在线中文字幕| 亚洲av电影在线观看一区二区三区| 亚洲精品在线美女| 国产探花极品一区二区| 街头女战士在线观看网站| 少妇的丰满在线观看| 丰满少妇做爰视频| 老熟女久久久| 欧美激情高清一区二区三区 | 成人二区视频| www.av在线官网国产| 各种免费的搞黄视频| www日本在线高清视频| 国产av一区二区精品久久| 国产一区二区三区av在线| 黄色毛片三级朝国网站| 如何舔出高潮| 国产男女超爽视频在线观看| 成年女人毛片免费观看观看9 | 免费观看a级毛片全部| 少妇猛男粗大的猛烈进出视频| 国产深夜福利视频在线观看| 成人亚洲欧美一区二区av| 香蕉国产在线看| 啦啦啦中文免费视频观看日本| 久久精品国产亚洲av天美| 亚洲一码二码三码区别大吗| 在线观看国产h片| 久久精品国产鲁丝片午夜精品| 天天操日日干夜夜撸| 色视频在线一区二区三区| 十八禁高潮呻吟视频| 天天躁夜夜躁狠狠久久av| 一二三四中文在线观看免费高清| 新久久久久国产一级毛片| 纵有疾风起免费观看全集完整版| 亚洲国产精品一区三区| 久久精品国产a三级三级三级| 成年人免费黄色播放视频| 亚洲综合色惰| 国产乱来视频区| 成人免费观看视频高清| xxx大片免费视频| 婷婷色av中文字幕| 黑丝袜美女国产一区| 久久久久人妻精品一区果冻| 国产一级毛片在线| 一级爰片在线观看| 99九九在线精品视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 热99久久久久精品小说推荐| 婷婷色综合大香蕉| 免费日韩欧美在线观看| 成人午夜精彩视频在线观看| 日韩中字成人| 免费日韩欧美在线观看| 熟女av电影| 看非洲黑人一级黄片| 夜夜骑夜夜射夜夜干| 18禁国产床啪视频网站| 亚洲综合色惰| 99热全是精品| 国产熟女欧美一区二区| 成人影院久久| 如何舔出高潮| 波野结衣二区三区在线| 亚洲内射少妇av| 国产日韩欧美视频二区| av有码第一页| 国产男女内射视频| 亚洲精品中文字幕在线视频| 午夜福利在线观看免费完整高清在| 我的亚洲天堂| 亚洲欧美色中文字幕在线| 波多野结衣一区麻豆| 女性被躁到高潮视频| 捣出白浆h1v1| 亚洲精品国产一区二区精华液| 香蕉精品网在线| 国产成人精品无人区| 久久 成人 亚洲| 国产又爽黄色视频| 欧美日韩一级在线毛片| 国产精品99久久99久久久不卡 | 午夜久久久在线观看| 婷婷色av中文字幕| 在线天堂最新版资源| 精品一区在线观看国产| 亚洲av在线观看美女高潮| 亚洲激情五月婷婷啪啪| 青青草视频在线视频观看| 在线亚洲精品国产二区图片欧美| 黑丝袜美女国产一区| 欧美精品国产亚洲|