• <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
    99久久人妻综合| 国产免费一区二区三区四区乱码| 亚洲人成77777在线视频| 天天躁夜夜躁狠狠躁躁| 日韩一区二区视频免费看| 亚洲欧美成人精品一区二区| av女优亚洲男人天堂| 一本—道久久a久久精品蜜桃钙片| 国产日韩一区二区三区精品不卡| 日本欧美国产在线视频| 人人妻人人澡人人看| 亚洲精品日韩在线中文字幕| 蜜臀久久99精品久久宅男| 成人毛片60女人毛片免费| 国产成人精品福利久久| 日韩av免费高清视频| 狂野欧美激情性bbbbbb| 精品久久久久久电影网| 亚洲精品久久久久久婷婷小说| 国产探花极品一区二区| 久久99一区二区三区| 亚洲经典国产精华液单| 777米奇影视久久| 久久久久久人妻| 黄网站色视频无遮挡免费观看| 国产精品国产三级国产av玫瑰| 成人国语在线视频| 免费久久久久久久精品成人欧美视频 | 晚上一个人看的免费电影| 欧美日韩一区二区视频在线观看视频在线| 十八禁高潮呻吟视频| 亚洲国产成人一精品久久久| 午夜视频国产福利| 99精国产麻豆久久婷婷| 精品久久久久久电影网| 亚洲熟女精品中文字幕| 亚洲精品久久午夜乱码| 国产精品.久久久| 久久久久国产精品人妻一区二区| 久久久国产一区二区| 在线观看三级黄色| 国产一区二区在线观看av| 亚洲第一av免费看| 久久久久久久精品精品| 亚洲美女视频黄频| 欧美日韩一区二区视频在线观看视频在线| 亚洲欧美一区二区三区黑人 | 永久网站在线| 久久久久网色| 夫妻午夜视频| 中文乱码字字幕精品一区二区三区| 国产av国产精品国产| 久久久久国产网址| 丰满迷人的少妇在线观看| 成人国产av品久久久| 伊人久久国产一区二区| 国产男女超爽视频在线观看| 久久午夜福利片| 91午夜精品亚洲一区二区三区| 岛国毛片在线播放| 又黄又粗又硬又大视频| 欧美日韩成人在线一区二区| 欧美+日韩+精品| 亚洲av在线观看美女高潮| 欧美+日韩+精品| 一级毛片黄色毛片免费观看视频| 中国美白少妇内射xxxbb| 日韩av免费高清视频| 日本黄色日本黄色录像| 国产乱人偷精品视频| 人妻 亚洲 视频| 国产免费一级a男人的天堂| 久久久亚洲精品成人影院| 国产高清不卡午夜福利| av电影中文网址| 激情五月婷婷亚洲| 日韩中文字幕视频在线看片| av片东京热男人的天堂| 一本久久精品| 久久久久久人人人人人| 国内精品宾馆在线| 日韩不卡一区二区三区视频在线| 熟女av电影| 飞空精品影院首页| 久久久久久久久久久免费av| 女人被躁到高潮嗷嗷叫费观| 国产精品久久久av美女十八| 免费少妇av软件| 国产亚洲最大av| 熟女电影av网| 婷婷成人精品国产| 免费黄色在线免费观看| 在线观看免费视频网站a站| 天天影视国产精品| 免费观看无遮挡的男女| 超色免费av| 亚洲一码二码三码区别大吗| 免费大片黄手机在线观看| 两性夫妻黄色片 | 中文字幕另类日韩欧美亚洲嫩草| 青春草亚洲视频在线观看| 边亲边吃奶的免费视频| 国产精品国产av在线观看| 少妇精品久久久久久久| 妹子高潮喷水视频| 青春草亚洲视频在线观看| 亚洲av.av天堂| 国产成人精品无人区| 成人免费观看视频高清| 中文字幕另类日韩欧美亚洲嫩草| 国产爽快片一区二区三区| 久久影院123| 美女主播在线视频| 美女脱内裤让男人舔精品视频| 亚洲,欧美,日韩| 丝袜在线中文字幕| 国产麻豆69| 中文字幕另类日韩欧美亚洲嫩草| 久久精品国产鲁丝片午夜精品| 国产69精品久久久久777片| av卡一久久| 成人手机av| 老司机影院毛片| 成年动漫av网址| 国产亚洲av片在线观看秒播厂| 精品亚洲成a人片在线观看| 天堂俺去俺来也www色官网| 亚洲成国产人片在线观看| 欧美日韩一区二区视频在线观看视频在线| 欧美日韩亚洲高清精品| 少妇精品久久久久久久| 国产成人午夜福利电影在线观看| 一二三四中文在线观看免费高清| 日韩av在线免费看完整版不卡| 久久久久久久国产电影| 人妻少妇偷人精品九色| 国产亚洲精品久久久com| 欧美日韩成人在线一区二区| 亚洲少妇的诱惑av| 视频中文字幕在线观看| 国产视频首页在线观看| 欧美国产精品va在线观看不卡| 一二三四中文在线观看免费高清| 岛国毛片在线播放| 亚洲第一av免费看| 久久久久久人妻| 国产亚洲av片在线观看秒播厂| 又黄又爽又刺激的免费视频.| 欧美xxxx性猛交bbbb| 黑丝袜美女国产一区| 女的被弄到高潮叫床怎么办| 亚洲四区av| 晚上一个人看的免费电影| av女优亚洲男人天堂| 午夜激情久久久久久久| 成人国产av品久久久| 在线天堂中文资源库| 高清欧美精品videossex| 成人18禁高潮啪啪吃奶动态图| 亚洲美女视频黄频| 一本色道久久久久久精品综合| 男女国产视频网站| 激情五月婷婷亚洲| xxx大片免费视频| 午夜激情久久久久久久| av福利片在线| 丝袜人妻中文字幕| 久久99蜜桃精品久久| 免费高清在线观看视频在线观看| 美女脱内裤让男人舔精品视频| 国产精品麻豆人妻色哟哟久久| 夜夜骑夜夜射夜夜干| 又黄又粗又硬又大视频| 在线 av 中文字幕| 亚洲国产精品国产精品| 国语对白做爰xxxⅹ性视频网站| 久久这里只有精品19| 免费人妻精品一区二区三区视频| 亚洲av电影在线进入| 22中文网久久字幕| av播播在线观看一区| 99精国产麻豆久久婷婷| 91精品三级在线观看| 欧美激情 高清一区二区三区| 久热久热在线精品观看| 国产成人精品一,二区| 国产免费现黄频在线看| 黑丝袜美女国产一区| 免费少妇av软件| 国产欧美日韩一区二区三区在线| 亚洲欧美日韩另类电影网站| 欧美精品亚洲一区二区| 欧美日韩亚洲高清精品| 伦精品一区二区三区| 国产av一区二区精品久久| av在线app专区| 免费看av在线观看网站| 亚洲精品中文字幕在线视频| av线在线观看网站| 99国产精品免费福利视频| 欧美精品一区二区大全| 最近最新中文字幕大全免费视频 | 欧美少妇被猛烈插入视频| 亚洲av电影在线观看一区二区三区| 日本爱情动作片www.在线观看| 婷婷色麻豆天堂久久| 亚洲成色77777| 狠狠婷婷综合久久久久久88av| 精品酒店卫生间| 亚洲av日韩在线播放| 日韩精品有码人妻一区| 男女下面插进去视频免费观看 | 亚洲天堂av无毛| 国产熟女午夜一区二区三区| 成人国产av品久久久| 国产成人一区二区在线| 精品视频人人做人人爽| 日韩一区二区视频免费看| 国产在线免费精品| 欧美人与性动交α欧美精品济南到 | 国产免费视频播放在线视频| 国产精品熟女久久久久浪| 欧美xxxx性猛交bbbb| av一本久久久久| 日韩大片免费观看网站| 人人澡人人妻人| 中文字幕av电影在线播放| 亚洲国产色片| 午夜影院在线不卡| 久久久a久久爽久久v久久| 午夜福利视频精品| 国产精品久久久久成人av| 最近2019中文字幕mv第一页| 久久久久精品性色| 国产日韩欧美亚洲二区| 女人久久www免费人成看片| 国产伦理片在线播放av一区| 国产xxxxx性猛交| 午夜视频国产福利| 久久毛片免费看一区二区三区| 精品酒店卫生间| 一级毛片 在线播放| 最后的刺客免费高清国语| 多毛熟女@视频| 国产精品三级大全| 精品少妇内射三级| 国产精品蜜桃在线观看| 深夜精品福利| 制服丝袜香蕉在线| 国产成人精品久久久久久| 大香蕉久久网| 欧美xxⅹ黑人| 亚洲国产av新网站| av国产精品久久久久影院| 色婷婷久久久亚洲欧美| 欧美日韩视频精品一区| 少妇人妻 视频| 黑人欧美特级aaaaaa片| 欧美日韩视频高清一区二区三区二| xxx大片免费视频| 亚洲国产精品一区三区| a级毛片黄视频| 精品人妻在线不人妻| 久久久亚洲精品成人影院| 成人手机av| 国产一区二区在线观看av| 亚洲成国产人片在线观看| 久久97久久精品| 国产精品嫩草影院av在线观看| 国产乱来视频区| 精品人妻熟女毛片av久久网站| 国产精品国产三级专区第一集| 亚洲欧美精品自产自拍| 波野结衣二区三区在线| 国产精品久久久久久久电影| 国产女主播在线喷水免费视频网站| 五月伊人婷婷丁香| 久久97久久精品| 丝袜在线中文字幕| 精品人妻一区二区三区麻豆| 黑人巨大精品欧美一区二区蜜桃 | 波野结衣二区三区在线| 精品福利永久在线观看| 亚洲精品色激情综合| 老司机亚洲免费影院| 亚洲国产看品久久| 欧美+日韩+精品| 亚洲欧美成人精品一区二区| 亚洲精品一二三| 人体艺术视频欧美日本| 各种免费的搞黄视频| 韩国高清视频一区二区三区| 女人久久www免费人成看片| 亚洲色图 男人天堂 中文字幕 | 久久av网站| 99国产综合亚洲精品| 丝袜人妻中文字幕| 黄色怎么调成土黄色| 夜夜爽夜夜爽视频| 免费播放大片免费观看视频在线观看| 久久久久久久精品精品| a级片在线免费高清观看视频| 亚洲国产看品久久| a级毛片在线看网站| 亚洲情色 制服丝袜| 国产成人a∨麻豆精品| 性色avwww在线观看| 熟女人妻精品中文字幕| 精品一区二区三卡| 国产成人精品无人区| av天堂久久9| 亚洲av在线观看美女高潮| 亚洲精品视频女| 亚洲精品久久午夜乱码| 一本大道久久a久久精品| 欧美成人午夜精品| 国产精品久久久久久精品古装| av黄色大香蕉| 亚洲精品美女久久久久99蜜臀 | 日韩电影二区| 精品人妻在线不人妻| 女性生殖器流出的白浆| 国产国语露脸激情在线看| 国产成人免费无遮挡视频| 亚洲精品一区蜜桃| av网站免费在线观看视频| 男女高潮啪啪啪动态图| 一本久久精品| 少妇人妻精品综合一区二区| 黄色 视频免费看| 欧美+日韩+精品| 你懂的网址亚洲精品在线观看| 国产欧美日韩综合在线一区二区| 亚洲成人手机| 国产高清三级在线| 久久毛片免费看一区二区三区| 国产在线视频一区二区| 日韩人妻精品一区2区三区| www.色视频.com| 亚洲精品久久午夜乱码| 久久久久久人妻| 国产成人精品无人区| 一级毛片电影观看| 大香蕉97超碰在线| 97精品久久久久久久久久精品| 人人澡人人妻人| 亚洲精品av麻豆狂野| 9191精品国产免费久久| 亚洲国产欧美日韩在线播放| 亚洲欧美成人精品一区二区| av有码第一页| 国产免费现黄频在线看| 色94色欧美一区二区| 欧美激情国产日韩精品一区| 成人综合一区亚洲| 老女人水多毛片| 亚洲人成网站在线观看播放| 色网站视频免费| 亚洲av国产av综合av卡| 国产免费福利视频在线观看| 黑丝袜美女国产一区| 天堂中文最新版在线下载| 精品人妻偷拍中文字幕| 99re6热这里在线精品视频| 久久精品人人爽人人爽视色| 日产精品乱码卡一卡2卡三| 国产精品欧美亚洲77777| 另类精品久久| 亚洲欧美精品自产自拍| av一本久久久久| 国产老妇伦熟女老妇高清| 中文字幕亚洲精品专区| 久久久精品免费免费高清| 大陆偷拍与自拍| av免费在线看不卡| 欧美少妇被猛烈插入视频| 久热久热在线精品观看| 大香蕉久久成人网| 免费黄色在线免费观看| 久久免费观看电影| 毛片一级片免费看久久久久| 国产一级毛片在线| 日日啪夜夜爽| 99re6热这里在线精品视频| 久久99一区二区三区| 我的女老师完整版在线观看| 亚洲婷婷狠狠爱综合网| 精品人妻一区二区三区麻豆| 免费观看在线日韩| 女人精品久久久久毛片| 精品国产国语对白av| 亚洲欧美日韩另类电影网站| 制服人妻中文乱码| 国产精品久久久久成人av| kizo精华| 欧美人与性动交α欧美精品济南到 | 成人国语在线视频| av电影中文网址| 精品国产一区二区三区四区第35| 国产精品国产三级专区第一集| 成人无遮挡网站| 久久av网站| 日日撸夜夜添| 少妇的逼好多水| 日本av手机在线免费观看| 亚洲五月色婷婷综合| 午夜免费观看性视频| 亚洲熟女精品中文字幕| 亚洲欧美中文字幕日韩二区| 少妇高潮的动态图| 免费不卡的大黄色大毛片视频在线观看| 精品久久久精品久久久| 另类精品久久| 日韩精品免费视频一区二区三区 | 午夜精品国产一区二区电影| 久久99一区二区三区| 中文字幕av电影在线播放| 丝瓜视频免费看黄片| 日韩电影二区| 伊人亚洲综合成人网| 啦啦啦视频在线资源免费观看| 狂野欧美激情性bbbbbb| 天堂中文最新版在线下载| 亚洲 欧美一区二区三区| av在线老鸭窝| 免费看光身美女| 波多野结衣一区麻豆| 久久国产亚洲av麻豆专区| 中文字幕亚洲精品专区| www.av在线官网国产| 啦啦啦在线观看免费高清www| 99久国产av精品国产电影| 精品人妻熟女毛片av久久网站| 久久精品国产自在天天线| 午夜福利视频精品| 婷婷成人精品国产| 深夜精品福利| 亚洲精品久久久久久婷婷小说| 日韩大片免费观看网站| 国产黄色视频一区二区在线观看| 欧美97在线视频| 搡女人真爽免费视频火全软件| 国产成人免费无遮挡视频| 亚洲婷婷狠狠爱综合网| 男女边吃奶边做爰视频| 成年人午夜在线观看视频| 蜜桃在线观看..| 国产成人91sexporn| 丝袜美足系列| 欧美97在线视频| 久久精品国产自在天天线| 一级毛片我不卡| 老女人水多毛片| 精品人妻偷拍中文字幕| 日韩不卡一区二区三区视频在线| 黄色毛片三级朝国网站| 久久精品熟女亚洲av麻豆精品| 午夜日本视频在线| 人妻人人澡人人爽人人| 国产精品.久久久| 成人国语在线视频| a级毛色黄片| 国产一区二区激情短视频 | 中文字幕精品免费在线观看视频 | 菩萨蛮人人尽说江南好唐韦庄| 嫩草影院入口| 丁香六月天网| 久久99精品国语久久久| 美女内射精品一级片tv| 成年动漫av网址| 日韩制服骚丝袜av| 欧美人与性动交α欧美软件 | 日韩中字成人| 久久国产亚洲av麻豆专区| 美女中出高潮动态图| 欧美精品国产亚洲| 久久毛片免费看一区二区三区| 国产精品女同一区二区软件| 伦精品一区二区三区| 国产国拍精品亚洲av在线观看| 三上悠亚av全集在线观看| 国产有黄有色有爽视频| 国产精品女同一区二区软件| 99久久人妻综合| 日本av免费视频播放| 亚洲精品日本国产第一区| 深夜精品福利| √禁漫天堂资源中文www| 国产精品久久久久久久久免| 免费观看在线日韩| 美国免费a级毛片| 一级片免费观看大全| 一二三四在线观看免费中文在 | 最近最新中文字幕大全免费视频 | 国产成人免费观看mmmm| 在线免费观看不下载黄p国产| 亚洲经典国产精华液单| 制服人妻中文乱码| 熟妇人妻不卡中文字幕| 亚洲国产精品一区二区三区在线| 欧美另类一区| 内地一区二区视频在线| 天天躁夜夜躁狠狠久久av| 国产1区2区3区精品| av视频免费观看在线观看| 久久精品国产自在天天线| 男人爽女人下面视频在线观看| 母亲3免费完整高清在线观看 | 超碰97精品在线观看| 成人手机av| 亚洲av综合色区一区| 国产免费福利视频在线观看| 免费观看无遮挡的男女| 日韩欧美一区视频在线观看| 免费播放大片免费观看视频在线观看| 午夜福利视频精品| 亚洲性久久影院| 五月玫瑰六月丁香| 18禁观看日本| 免费人妻精品一区二区三区视频| 色婷婷久久久亚洲欧美| 永久网站在线| 丝袜美足系列| 新久久久久国产一级毛片| 国产精品偷伦视频观看了| 亚洲欧美清纯卡通| 亚洲第一av免费看| 久久这里只有精品19| 婷婷色av中文字幕| 国产永久视频网站| 亚洲人与动物交配视频| 最新的欧美精品一区二区| 黄网站色视频无遮挡免费观看| 亚洲成色77777| 日韩成人av中文字幕在线观看| 老女人水多毛片| 91午夜精品亚洲一区二区三区| 免费观看无遮挡的男女| 国产精品不卡视频一区二区| 亚洲国产av新网站| 如何舔出高潮| 国产精品国产三级专区第一集| 亚洲国产精品成人久久小说| 亚洲一码二码三码区别大吗| 国产 精品1| 国产男女超爽视频在线观看| 青春草亚洲视频在线观看| 成年动漫av网址| 亚洲内射少妇av| 999精品在线视频| 日韩av不卡免费在线播放| 中文字幕精品免费在线观看视频 | 9色porny在线观看| 免费女性裸体啪啪无遮挡网站| 少妇熟女欧美另类| 国产一区二区在线观看av| 国产色婷婷99| 最近最新中文字幕免费大全7| 国产精品一二三区在线看| 成人毛片60女人毛片免费| 亚洲精品国产av成人精品| 国产免费现黄频在线看| 国精品久久久久久国模美| 80岁老熟妇乱子伦牲交| 色哟哟·www| 久久精品国产鲁丝片午夜精品| 亚洲精品日韩在线中文字幕| 亚洲国产色片| 亚洲美女黄色视频免费看| 中文字幕免费在线视频6| 观看av在线不卡| 国产69精品久久久久777片| 午夜日本视频在线| 97精品久久久久久久久久精品| 高清黄色对白视频在线免费看| 亚洲精品av麻豆狂野| 黄色一级大片看看| 亚洲,一卡二卡三卡| 亚洲av电影在线进入| 嫩草影院入口| 熟女电影av网| 久久久久精品性色| 亚洲天堂av无毛| av线在线观看网站| 中国三级夫妇交换| 欧美精品高潮呻吟av久久| 亚洲精品色激情综合| 国产成人精品一,二区| 欧美精品国产亚洲| 精品人妻一区二区三区麻豆| 日韩制服骚丝袜av| 国产乱人偷精品视频| 极品人妻少妇av视频| 午夜福利视频在线观看免费| 国产精品人妻久久久影院| 成人综合一区亚洲| 午夜久久久在线观看| 少妇被粗大的猛进出69影院 | 国产成人av激情在线播放| 久久精品国产亚洲av涩爱| 亚洲欧美日韩卡通动漫| 亚洲中文av在线| 夜夜骑夜夜射夜夜干| 交换朋友夫妻互换小说| 欧美精品亚洲一区二区| 日本与韩国留学比较| 精品亚洲成a人片在线观看| 久久ye,这里只有精品| 精品少妇黑人巨大在线播放| 美女内射精品一级片tv| 日本黄色日本黄色录像| 热re99久久国产66热| 一级毛片黄色毛片免费观看视频| 人妻人人澡人人爽人人| 午夜福利在线观看免费完整高清在|