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

    跨越飽和/非飽和土壤條件下豎直地埋管換熱器傳熱特性研究

    2020-04-15 06:41:30蔡學(xué)龍王松慶
    可再生能源 2020年4期
    關(guān)鍵詞:土壤層型管非飽和

    蔡學(xué)龍,王松慶

    (東北林業(yè)大學(xué) 土木工程學(xué)院,黑龍江 哈爾濱 150040)

    0 引言

    受到占地面積和換熱穩(wěn)定性等因素的限制,土壤源熱泵系統(tǒng)的室外地下?lián)Q熱器多采用豎直地埋管換熱器(以下簡稱為地埋管)。地埋管通常會跨越不同的地質(zhì)層,即一部分位于無滲流區(qū),另一部分位于地下水位以下。 對于前者,地埋管與周圍土壤之間的傳熱過程是一個溫度梯度下的熱量傳遞和土、水勢梯度下的水分遷移相互耦合的復(fù)雜傳熱傳質(zhì)過程;對于后者,地埋管與周圍土壤之間的傳熱過程是一個溫度梯度下的熱傳導(dǎo)和水力梯度下的對流換熱相互耦合的傳熱過程。 對于實(shí)際工程中的地埋管,如果將其與土壤之間的傳熱過程等效為純導(dǎo)熱問題,并單純地選用線熱源模型或柱熱源模型等進(jìn)行分析和設(shè)計,會造成設(shè)計結(jié)果與實(shí)際結(jié)果偏差較大,影響土壤源熱泵系統(tǒng)的運(yùn)行性能及其初期投 資[1],[2]。

    地埋管的換熱性能不僅與其周圍土壤的物性參數(shù)密切相關(guān),還會受到地下水滲流和分層作用的影響[3]。 Chiasson[4]通過有限差分法,分析了地下水滲流對地埋管換熱性能的影響,分析結(jié)果表明,若忽略地下水在土壤中的遷移,會導(dǎo)致地埋管幾何尺寸設(shè)計得偏大,使得初期投資隨之增加。 Angelotti[5]研究了地下水滲流速度對土壤傳熱性能以及熱干擾的影響,分析結(jié)果表明,地下水滲流速度越大,地埋管周圍土壤的等溫線沿滲流方向的偏移度越大。 Capozza[6]研究了地下水滲流對于土壤中熱量遷移情況的影響,分析結(jié)果表明,若考慮地下水滲流,則會導(dǎo)致地埋管總長度的預(yù)測值減小。

    實(shí)際情況下,土壤的分層情況以及物性的不均勻情況均會對地埋管換熱器的傳熱情況產(chǎn)生影響,但影響規(guī)律尚不明確,值得深入探討[7]。 景政[8]通過數(shù)值方法和實(shí)驗(yàn)方法,研究了滲流作用對地埋管傳熱情況的影響,分析結(jié)果表明,滲流作用會大幅度地增加土壤的導(dǎo)熱系數(shù)以及地埋管與周圍土壤之間的換熱量。張琳琳[9]建立了同時考慮分層和滲流作用的地埋管傳熱解析模型,模擬結(jié)果表明,部分土壤層存在滲流時的熱作用距離比全部土壤層存在滲流時的熱作用距離減小了43%。 可見,地下水滲流和分層作用對地埋管傳熱性能的影響是不容忽視的。

    上述研究雖然考慮了滲流作用和分層作用,但沒有考慮滲流作用情況下,非飽和土壤和飽和土壤相鄰間隔分層布置時,地埋管的傳熱特性。本文同時考慮分層作用和滲流作用,并采用數(shù)值仿真方法研究了地埋管跨越非飽和土壤和飽和土壤時的傳熱規(guī)律,為土壤源熱泵系統(tǒng)的設(shè)計提供理論依據(jù)。

    1 控制方程

    1.1 地埋管換熱器傳熱模型

    由于地埋管換熱器內(nèi)循環(huán)液(工質(zhì))的流線會在U 型管的彎管處發(fā)生較大的變化,因此本文采用Realizable k-ε 模型描述地埋管換熱器內(nèi)循環(huán)液的流動和傳熱情況。

    地埋管換熱器內(nèi)循環(huán)液的連續(xù)性方程為

    式中:ρ 為地埋管換熱器內(nèi)循環(huán)液的密度,kg/m3;分別為地埋管換熱器內(nèi)循環(huán)液在 x 軸、y 軸方向上的流速,m/s;xj為y 軸方向上地埋管換熱器內(nèi)循環(huán)液的流動距離,m。

    地埋管換熱器內(nèi)循環(huán)液的動量方程為[10]

    式中:μe為地埋管內(nèi)循環(huán)液分子擴(kuò)散造成的動力粘性;xk為z 軸方向上地埋管換熱器內(nèi)循環(huán)液的流動距離,m;為地埋管換熱器內(nèi)循環(huán)液在z 軸方向上的流速,m/s。

    地埋管換熱器內(nèi)循環(huán)液的能量方程為

    式中:xi為x 軸方向上地埋管換熱器內(nèi)循環(huán)液流動的距離,m;k 為地埋管內(nèi)循環(huán)液的湍流脈動動能,J;μt為地埋管換熱器內(nèi)循環(huán)液的湍流粘性系數(shù);μ 為地埋管換熱器內(nèi)循環(huán)液的動力粘度,(N·s)/m2;ε 為地埋管換熱器的效能;ν 為地埋管換熱器內(nèi)循環(huán)液分子的粘性;Cε1,Cε2均為系數(shù),Cε1=1.2,Cε2=1.44;σk,σε均為常量,σk=1.0,σε=1.2[11]。

    1.2 土壤傳熱模型

    1.2.1 飽和土壤傳熱模型

    對于有、無滲流的飽和土壤,一個物質(zhì)系統(tǒng)或空間體積內(nèi)含有固體和液體2 個部分。因此,對于選取的任意一個控制體積而言,令其孔隙率為φ,那么單位體積內(nèi)液體和固體骨架占據(jù)的體積分?jǐn)?shù)分別為 φ 和1-φ ,于是本文基于體積平均法,構(gòu)建出多孔介質(zhì)中各部分的能量方程[12]。

    多孔介質(zhì)中液體相的能量方程為

    式中:ρf為多孔介質(zhì)中液體相的密度,Kg/m3;cp為多孔介質(zhì)中液體相的定壓比熱,J/(Kg﹒℃);Tf為多孔介質(zhì)中液體相的溫度,℃;λf為多孔介質(zhì)中液體相的導(dǎo)熱系數(shù),W·(m·℃)-1;V 為土壤的體積,m3;qf為多孔介質(zhì)中單位體積液體相的內(nèi)熱源產(chǎn)生的熱量,J/m3。

    多孔介質(zhì)中固體相的能量方程為

    式中:ρ 為土壤的干密度,kg·m-3;c 為土壤的比熱,J/(kg·℃);Ts為土壤的溫度,℃;λs為多孔介質(zhì)中固體骨架的導(dǎo)熱系數(shù),W/(m·℃);qs為多孔介質(zhì)中單位體積固體骨架相的內(nèi)熱源產(chǎn)生的熱量,J/m3。

    通常可以假設(shè)流體和固體可以在瞬間達(dá)到局部熱平衡狀態(tài),此外,本文將流體、固體的熱容和導(dǎo)熱系數(shù)均視為常數(shù),則基于式(5),(6)可得到有滲流的飽和土壤的能量方程為

    式中:λm為土壤的表觀導(dǎo)熱系數(shù),W/(m·℃);T 為土壤的溫度,℃;qm為土壤的表觀內(nèi)熱源產(chǎn)生率;cm為土壤的表觀容積比熱,J/(kg·℃)。

    式(7)中 ρcm,λm,qm的計算式分別為

    1.2.2 非飽和土壤傳熱模型

    在土壤非飽和區(qū),非等溫條件下,土壤水分的遷移可以由擴(kuò)散達(dá)西定律進(jìn)行描述[13]。其中:液相質(zhì)流由對流質(zhì)流和表面擴(kuò)散質(zhì)流組成; 蒸汽質(zhì)流由蒸汽對流質(zhì)流和蒸汽在多孔介質(zhì)中產(chǎn)生的擴(kuò)散質(zhì)流組成;空氣質(zhì)流由對流質(zhì)流和擴(kuò)散質(zhì)流組成。

    液相、蒸汽和空氣的質(zhì)流方程分別為

    式中:JL,Jv,Ja分別為液態(tài)水、蒸汽、空氣的質(zhì)流,kg/s;DTL為水的熱質(zhì)擴(kuò)散系數(shù);DTD為吸附傳遞系數(shù);D 為分子擴(kuò)散系數(shù);DθL為等溫質(zhì)擴(kuò)散系數(shù);DTv為蒸汽熱質(zhì)擴(kuò)散系數(shù);Dv,a為 Knudsen 擴(kuò)散系數(shù);θL為液相的百分比;K 為水力傳導(dǎo)系數(shù);g 為重力加速度,N/kg;ρL為液態(tài)水的密度,kg/m3; ρv為蒸汽的密度,kg/m3;ρa(bǔ)為空氣的密度,kg/m3。

    土壤非飽和區(qū)的質(zhì)量守恒方程為

    土壤非飽和區(qū)的能量守恒方程為

    式中:ε 為土壤的孔隙率;hL,hv,ha分別為液體、蒸汽、空氣的焓值;qv為土壤內(nèi)熱源;λeff當(dāng)量導(dǎo)熱系數(shù),W/(m·℃)。

    λeff的計算式為

    式中:θs,θL,θg分別為固相、液相和氣相的百分比。

    θs,θL,θg的關(guān)系式為

    2 數(shù)值仿真研究

    2.1 物理模型建立與網(wǎng)格劃分

    本文采用相似理論建立幾何尺寸小于實(shí)際情況的物理模型。 在考慮分層和滲流同時作用的情況下,若按照實(shí)際幾何尺寸設(shè)計模型,則會導(dǎo)致在仿真過程中難以控制水在分層土壤中滲流的均勻性,仿真計算也很難收斂。利用幾何尺寸小于實(shí)際情況的物理模型能規(guī)避上述問題。

    本文的模型為長方體,幾何尺寸為800 mm×800 mm×900 mm。 利用3 個等分點(diǎn)將該模型上表面的一條對角線等分成4 份,以該對角線中點(diǎn)兩側(cè)的等分點(diǎn)為中心,向下布置了兩根U 型管,每根U 型管的內(nèi)徑均為8 mm,外徑均為10 mm,長均為830 mm,U 型管兩支管的管心距均為 40 mm。 由于U 型管與土壤之間的接觸面存在復(fù)雜性,因此,該模型采用非結(jié)構(gòu)化網(wǎng)格進(jìn)行劃分,并將U 型管管壁設(shè)置為邊界層,此外,對于U 型管由內(nèi)到外適當(dāng)?shù)卦龃缶W(wǎng)格的尺寸,以減少網(wǎng)格數(shù)量,提高計算效率。

    為了確保計算結(jié)果的準(zhǔn)確性,須定量地確定模型的最優(yōu)網(wǎng)格數(shù),因此,本文進(jìn)行了網(wǎng)格的無關(guān)性分析??紤]模擬過程中的各種約束條件,分別計算了4 種不同的網(wǎng)格劃分模型,網(wǎng)格數(shù)量分別為700 000,900 000,1 100 000 和 1 400 000,網(wǎng)格質(zhì)量均為 0.5~1。 監(jiān)測的溫度:U 型管壁面的最高溫度、平均溫度,以及U 型管出口工質(zhì)的平均溫度。

    每種工況各處監(jiān)測溫度以及各工況之間相對誤差的具體情況如表1 所示。 表1 中的相對誤差為相鄰工況測量結(jié)果之間的偏差。

    表1 網(wǎng)格獨(dú)立性驗(yàn)證Table 1 Verification of grid independence

    由表1 可以看出,當(dāng)網(wǎng)格數(shù)量不低于110 萬時,各項(xiàng)測量結(jié)果大體上保持不變,相對誤差均小于0.1%。因此,本文采用的網(wǎng)格數(shù)量為110 萬,這樣既能夠保證計算精度,又能夠保證計算效率。

    2.2 仿真參數(shù)與測點(diǎn)設(shè)置

    2.2.1 砂土熱物性參數(shù)

    本文對砂土進(jìn)行水平分層布置,自下而上分別為飽和礫砂層、非飽和中砂層和飽和粗砂層,厚度均為300 mm。由于中砂的保水性優(yōu)于礫砂和粗砂,故將中砂設(shè)置成非飽和砂土,這樣會使得非飽和土壤層的含水率更加均衡。

    各層砂土的熱物性參數(shù)如表2 所示。

    表2 各層砂土的熱物性參數(shù)Table 2 Thermophysical parameters of each layer of sand

    2.2.2 邊界條件設(shè)置

    模型的各個壁面均設(shè)置為絕熱壁面,U 型管入口、出口和土壤滲流入口、出口均設(shè)置為速度入口、出口。U 型管壁面及其與土壤層之間的交界面均設(shè)置為耦合傳熱。此外,設(shè)置模型中的工質(zhì)為多孔介質(zhì)混合流,由于土壤中滲流速度較低,因此,該混合流也被視為線性流。U 型管的材質(zhì)為銅,其內(nèi)部的工質(zhì)為水。本文研究了2 種工況下,土壤溫度的變化特性。這2 種工況均為夏季工況。各工況下,模型的仿真參數(shù)如表3 所示。

    表3 2 種工況下,模型的仿真參數(shù)Table 3 Simulation parameters of the model under 2 working conditions

    2.2.3 測點(diǎn)設(shè)置

    圖1 為土壤中各測點(diǎn)的布置情況。

    圖1 土壤中各測點(diǎn)的布置情況Fig.1 Arrangement of measuring points in soil

    由圖1 可知,在每個土壤層的x 軸、z 軸方向上各布置了 5 個測點(diǎn)。 對于下土壤層,測點(diǎn) 1,2,3以及測點(diǎn) 1,6,7 的間距均為 100 mm;測點(diǎn) 3,4,5以及測點(diǎn)7,8,9 的間距均為150 mm;其他土壤層中各測點(diǎn)的布置方式與下土壤層相同。

    2.3 案例分析

    2.3.1 無滲流條件下土壤溫度變化特性研究

    地埋管中的熱量會向土壤的四周傳遞,傳遞情況和傳遞規(guī)律將對該地埋管的設(shè)計有重要影響。 當(dāng)土壤層中無滲流(工況1)時,中土壤層、下土壤層中各測點(diǎn)處的溫度隨時間的變化情況如圖2 所示。

    圖2 當(dāng)土壤層中無滲流時,中土壤層、下土壤層中各測點(diǎn)處的溫度隨時間的變化情況Fig.2 When there is no seepage in the soil layer,the temperature changes with time at each measuring point in the middle soil layer and the lower soil layer

    由圖2 可知,在下土壤層x 軸、z 軸兩個方向上,均是靠近U 型管測點(diǎn)處的溫度升高速率較快,遠(yuǎn)離U 型管測點(diǎn)處的溫度升高速率較慢,且兩個方向上離U 型管等距離測點(diǎn)處溫度的升高幅度大致相等。其中,測點(diǎn)12 與測點(diǎn)16 的溫度變化曲線幾乎重合,x 軸、z 軸方向上其他對應(yīng)測點(diǎn)溫度的變化情況也大致相同,這說明在該土壤層中x 軸、z 軸方向上的傳熱情況相同,因此,在無滲流情況下,分析土壤層x 軸方向上的傳熱情況可以確定土壤層z 軸方向上的傳熱情況。此外,通過研究發(fā)現(xiàn),各土壤層中各測點(diǎn)處溫度的變化趨勢相同。

    為了更全面地明確各測點(diǎn)處土壤溫度的變化規(guī)律,對工況1 下各土壤層中各測點(diǎn)的溫度變化特性進(jìn)行深入分析。 由圖2 可知,下土壤層各測點(diǎn)處的溫度均呈現(xiàn)出逐漸升高的變化趨勢,這說明熱源傳遞至各測點(diǎn)處的熱量始終大于各測點(diǎn)處向外傳遞的熱量 (上土壤層也有相似的傳熱特性)。 中土壤層的傳熱情況不同于下土壤層,雖然該土壤層各測點(diǎn)處的溫度也呈現(xiàn)出逐漸升高的變化趨勢,但在溫度升高的過程中,各測點(diǎn)處的溫度出現(xiàn)波動現(xiàn)象,溫度上升趨勢不穩(wěn)定,這說明熱源傳遞到各測點(diǎn)的熱量有時大于各測點(diǎn)處向外傳遞的熱量,有時小于各測點(diǎn)處向外傳遞的熱量。 這是由于非飽和砂土中的水分在地埋管運(yùn)行過程中會向溫度高的區(qū)域移動,水分的遷移會造成砂土中各區(qū)域的含水率隨時間發(fā)生變化,此外,砂土的導(dǎo)熱系數(shù)與砂土含水率有關(guān),當(dāng)含水率發(fā)生變化時,會引起砂土的導(dǎo)熱系數(shù)隨之發(fā)生變化[14]。綜上可知,飽和土壤中的熱流量比非飽和土壤中的熱流量更穩(wěn)定。

    當(dāng)土壤層中無滲流,運(yùn)行時間為10 h 時,各土壤層x 軸方向上各測點(diǎn)處的溫度相對于初始溫度的升高幅度,如表4 所示。

    表4 當(dāng)土壤層中無滲流,運(yùn)行時間為10 h 時,各土壤層x 軸方向上各測點(diǎn)處的溫度相對于初始溫度的升高幅度Table 4 When there is no seepage in the soil layer and the operation time is 10 hours,the temperature of each measuring point in the x direction of each soil layer increases with respect to the initial temperature

    由表4 可知,對于不同的土壤層,距離U 型管中心越遠(yuǎn),其溫度相對于初始溫度的升高幅度越小,這是由于距離U 型管越遠(yuǎn),相鄰砂土之間的溫度差越小,傳熱量也越小。 因此,土壤中的熱流量與其和U 型管中心之間的徑向距離成反比關(guān)系。

    2.3.2 滲流條件下土壤溫度變化特性研究

    當(dāng)土壤中有滲流時,中土壤層、下土壤層中各測點(diǎn)處的溫度隨時間的變化情況如圖3 所示。

    圖3 當(dāng)土壤中有滲流時,中土壤層、下土壤層中各測點(diǎn)處的溫度隨時間的變化情況Fig.3 When there is seepage in the soil,the temperature changes with time at each measuring point in the middle soil layer and the lower soil layer

    由圖3 可知,當(dāng)土壤中有滲流時,下土壤層x軸方向上測點(diǎn)2~5 和z 軸方向上測點(diǎn)6~9 的溫度變化情況存在較大差異,這表明對于存在滲流的下土壤層,其內(nèi)部x 軸、z 軸方向上有著不同傳熱的特性。 此外,經(jīng)過研究發(fā)現(xiàn),上土壤層也有類似傳熱特性。 由圖3 還可看出,當(dāng)土壤中有滲流時,中土壤層x 軸方向上測點(diǎn)12~15 和z 軸方向上測點(diǎn)16~19 的溫度變化情況相似,這表明對于存在滲流的中土壤層,其x 軸、z 軸方向上有著相同的傳熱特性,綜上可知,當(dāng)飽和土壤層中有滲流時,其x 軸、z 軸方向上有不同的傳熱特性,對于非飽和土壤層以及沒有滲流的飽和土壤層,x 軸、z 軸方向上的傳熱特性相同。

    工況 1,2 下,上土壤層 x 軸、z 軸方向上各測點(diǎn)處的溫度隨時間的變化情況分別如圖4,5 所示。

    圖4 工況1,2 下,上土壤層x 方向上各測點(diǎn)處的溫度隨時間的變化情況Fig.4 Under working condition 1 and 2,the change of temperature with time at each measuring point in x direction of upper soil layer

    圖5 工況1,2 下,上土壤層z 方向上各測點(diǎn)處溫度隨時間的變化情況Fig.5 Under working condition 1 and 2,the change of temperature with time at each measuring point in z direction of upper soil layer

    由圖4,5 可知,兩種工況下,測點(diǎn) 22,23 處的溫度與初始溫度之間的差值會隨著運(yùn)行時間的增加而增大。當(dāng)運(yùn)行時間為10 h 時,工況1,2 下,測點(diǎn)22 處的溫度比初始溫度分別升高了28.8%,40.3%,兩者相差11.5%;由于上土壤層各測點(diǎn)處溫度的上升曲線無交點(diǎn),且變化趨勢穩(wěn)定,故本文在運(yùn)行時間為10 h 的條件下,分析上土壤層的徑向傳熱特性。

    當(dāng)運(yùn)行時間為10 h 時,兩種工況下,上土壤層各測點(diǎn)處的溫度相對于初始溫度的升高幅度,以及不同工況下升高幅度的差值如表5 所示。 表中:ΔT1,ΔT2分別為工況 1,2 下上土壤層各測點(diǎn)處的溫度相對于初始溫度的升高幅度;ΔT 為ΔT1與ΔT2的差值,表征了土壤層的傳熱能力。

    表5 當(dāng)運(yùn)行時間為10 h 時,兩種工況下,上土壤層各測點(diǎn)處的溫度相對于初始溫度的升高幅度,以及不同工況下升高幅度的差值Table 5 When the operation time is 10 hours,the increase range of the temperature at each measuring point of the upper soil layer relative to the initial temperature and the difference of the increase range under different working conditions

    由表5 可知,工況2 下,上土壤層x 軸方向上各測點(diǎn)處的溫度始終高于工況1 下各測點(diǎn)處的溫度; 上土壤層z 軸方向上大部分測點(diǎn)處的溫度始終高于工況1 下各測點(diǎn)處的溫度(測點(diǎn)21 除外),這說明當(dāng)土壤中存在滲流作用時,滲流方向上土壤的傳熱能力被強(qiáng)化,垂直滲流方向上土壤的傳熱能力被弱化。 通過計算得到,滲流方向上ΔT 的平均值約為7%,垂直滲流方向上ΔT 的平均值約為-3.8%,這是由于滲流作用能夠加速滲流方向上的工質(zhì)的熱量傳導(dǎo)速率,從而使更多的熱量向滲流方向擴(kuò)散,導(dǎo)致垂直滲流方向上熱量傳遞得較少。

    3 結(jié)論

    ①地埋管換熱器與飽和土壤之間的熱傳導(dǎo)比其與非飽和土壤之間的熱傳導(dǎo)更穩(wěn)定。 土壤中的熱流量與其和U 型管中心的徑向距離成反比關(guān)系。

    ②滲流對土壤的徑向傳熱有很大影響,當(dāng)土壤層中無滲流時,熱量在U 型管四周均勻地向外傳遞,而當(dāng)土壤層中存在滲流時,滲流方向上ΔT約為 7%,垂直滲流方向上 ΔT 的平均值約為-3.8%。

    猜你喜歡
    土壤層型管非飽和
    雄安新區(qū)上游油松林土壤層物理性質(zhì)研究
    東江中下游流域森林土壤有機(jī)碳空間分布特征研究*
    b型管板與筒體溫差應(yīng)力的分析計算和評定
    化工管理(2021年7期)2021-05-13 00:46:04
    內(nèi)置加勁環(huán)T型管節(jié)點(diǎn)抗沖擊承載力計算
    Study on the Letter Words and Letter-word Phrases under Language Contact Theory
    Montgomery T型管置入術(shù)治療聲門下氣管狹窄失敗三例報道并文獻(xiàn)復(fù)習(xí)
    非飽和原狀黃土結(jié)構(gòu)強(qiáng)度的試驗(yàn)研究
    非飽和多孔介質(zhì)應(yīng)力滲流耦合分析研究
    非飽和土基坑剛性擋墻抗傾覆設(shè)計與參數(shù)分析
    非飽和地基土蠕變特性試驗(yàn)研究
    国产免费一区二区三区四区乱码| 国产精品.久久久| 午夜福利高清视频| 啦啦啦在线观看免费高清www| 黄色怎么调成土黄色| 亚洲av免费高清在线观看| 丝瓜视频免费看黄片| 亚洲性久久影院| 一级毛片久久久久久久久女| 精品一区二区三区视频在线| 老女人水多毛片| 欧美激情极品国产一区二区三区 | 黄色视频在线播放观看不卡| 国产成人a区在线观看| 成人亚洲精品一区在线观看 | 亚洲精品日本国产第一区| 亚洲不卡免费看| 国产伦在线观看视频一区| 边亲边吃奶的免费视频| 爱豆传媒免费全集在线观看| 欧美日韩国产mv在线观看视频 | 国产精品国产三级国产av玫瑰| 久久97久久精品| 国产视频内射| 亚洲性久久影院| 亚洲av国产av综合av卡| av一本久久久久| 国产精品久久久久成人av| 亚洲不卡免费看| 亚洲国产日韩一区二区| 五月伊人婷婷丁香| 国产av国产精品国产| 国产美女午夜福利| 成人综合一区亚洲| 亚洲av欧美aⅴ国产| 亚洲自偷自拍三级| 一级毛片我不卡| 国产爽快片一区二区三区| 欧美高清成人免费视频www| 久久久久久九九精品二区国产| 国产精品人妻久久久影院| av国产久精品久网站免费入址| 嘟嘟电影网在线观看| 中国美白少妇内射xxxbb| 日韩一本色道免费dvd| 亚洲av福利一区| 婷婷色综合www| 尤物成人国产欧美一区二区三区| 熟女电影av网| 国产精品人妻久久久影院| 免费观看的影片在线观看| 尤物成人国产欧美一区二区三区| 欧美xxⅹ黑人| 亚洲国产欧美人成| 51国产日韩欧美| h视频一区二区三区| 国产精品一区二区性色av| 日韩中文字幕视频在线看片 | freevideosex欧美| 国产精品女同一区二区软件| 久久久久性生活片| 日本黄色日本黄色录像| 九九久久精品国产亚洲av麻豆| 久久精品久久精品一区二区三区| 久久久久久人妻| a级一级毛片免费在线观看| 九草在线视频观看| 亚洲综合色惰| 国产亚洲精品久久久com| 老师上课跳d突然被开到最大视频| 免费av不卡在线播放| www.色视频.com| 插阴视频在线观看视频| 国产 一区 欧美 日韩| 人体艺术视频欧美日本| 国产精品三级大全| 妹子高潮喷水视频| 黄色日韩在线| 大码成人一级视频| 免费在线观看成人毛片| 亚洲,欧美,日韩| 日本av免费视频播放| 综合色丁香网| 又爽又黄a免费视频| 人体艺术视频欧美日本| 各种免费的搞黄视频| 免费观看的影片在线观看| 日本欧美视频一区| 一区二区三区精品91| 欧美一区二区亚洲| 内地一区二区视频在线| 日本黄大片高清| 少妇的逼好多水| 国产精品久久久久久精品电影小说 | 国产精品国产av在线观看| 国产乱人视频| 午夜免费鲁丝| 三级国产精品片| 国产伦理片在线播放av一区| 亚洲精品,欧美精品| 在线观看美女被高潮喷水网站| 欧美3d第一页| av天堂中文字幕网| 国产亚洲午夜精品一区二区久久| 韩国高清视频一区二区三区| 十八禁网站网址无遮挡 | 偷拍熟女少妇极品色| 国产69精品久久久久777片| 国产伦理片在线播放av一区| 国产伦精品一区二区三区视频9| 啦啦啦视频在线资源免费观看| 免费黄网站久久成人精品| 波野结衣二区三区在线| 国产亚洲欧美精品永久| 国产亚洲一区二区精品| 最近中文字幕2019免费版| 成人国产av品久久久| 国产精品偷伦视频观看了| xxx大片免费视频| 欧美成人一区二区免费高清观看| 内地一区二区视频在线| 岛国毛片在线播放| 欧美亚洲 丝袜 人妻 在线| 高清av免费在线| 欧美xxxx黑人xx丫x性爽| 国产黄色视频一区二区在线观看| 青春草国产在线视频| 直男gayav资源| 国产伦理片在线播放av一区| 国产精品一及| 超碰av人人做人人爽久久| 男人和女人高潮做爰伦理| 亚洲精品中文字幕在线视频 | 下体分泌物呈黄色| 老熟女久久久| 日本黄色日本黄色录像| 日日啪夜夜撸| av国产免费在线观看| 自拍欧美九色日韩亚洲蝌蚪91 | 精品久久久精品久久久| 国产精品一二三区在线看| 日韩国内少妇激情av| 国产精品秋霞免费鲁丝片| 能在线免费看毛片的网站| 国产成人精品一,二区| 中文字幕av成人在线电影| 激情 狠狠 欧美| 亚洲成人一二三区av| 色综合色国产| 一边亲一边摸免费视频| 午夜激情福利司机影院| 网址你懂的国产日韩在线| 国产欧美另类精品又又久久亚洲欧美| 三级国产精品片| 国产欧美亚洲国产| 777米奇影视久久| 水蜜桃什么品种好| 日韩国内少妇激情av| 精品久久久久久电影网| 简卡轻食公司| 性高湖久久久久久久久免费观看| 日韩av免费高清视频| 成年免费大片在线观看| 国产精品一区二区三区四区免费观看| 日韩大片免费观看网站| 99热这里只有是精品在线观看| 国产成人freesex在线| www.色视频.com| 五月天丁香电影| 国产又色又爽无遮挡免| 亚洲av成人精品一二三区| 久久国内精品自在自线图片| 不卡视频在线观看欧美| 一级片'在线观看视频| 三级国产精品片| 久久久久久久久久久免费av| 日韩强制内射视频| 中文在线观看免费www的网站| 国内少妇人妻偷人精品xxx网站| 青青草视频在线视频观看| 熟妇人妻不卡中文字幕| 男女免费视频国产| 欧美zozozo另类| 国产一区有黄有色的免费视频| 男人舔奶头视频| 嫩草影院新地址| 涩涩av久久男人的天堂| 日本av免费视频播放| 欧美xxxx黑人xx丫x性爽| av专区在线播放| 最黄视频免费看| 99久久精品国产国产毛片| kizo精华| 一区二区av电影网| 久久人人爽人人片av| 赤兔流量卡办理| 国产精品一区二区性色av| 免费大片18禁| 精品一区二区三卡| 国产亚洲5aaaaa淫片| 大片电影免费在线观看免费| 哪个播放器可以免费观看大片| 中文欧美无线码| 成人影院久久| 亚洲人与动物交配视频| 欧美97在线视频| 国产精品久久久久久久电影| 又大又黄又爽视频免费| 色视频在线一区二区三区| 韩国av在线不卡| 小蜜桃在线观看免费完整版高清| 亚洲av中文av极速乱| 精品久久久久久久久av| 我要看黄色一级片免费的| 老师上课跳d突然被开到最大视频| av天堂中文字幕网| 亚洲国产精品专区欧美| 91精品伊人久久大香线蕉| 高清毛片免费看| 人人妻人人澡人人爽人人夜夜| 又大又黄又爽视频免费| 啦啦啦视频在线资源免费观看| 欧美成人精品欧美一级黄| 伦精品一区二区三区| 亚洲av男天堂| 午夜精品国产一区二区电影| 男女免费视频国产| 天天躁日日操中文字幕| 人妻一区二区av| 亚洲av男天堂| 狠狠精品人妻久久久久久综合| 男的添女的下面高潮视频| 久久国产亚洲av麻豆专区| 久热这里只有精品99| 日韩中字成人| 人妻系列 视频| 99视频精品全部免费 在线| 欧美国产精品一级二级三级 | 中文字幕久久专区| 菩萨蛮人人尽说江南好唐韦庄| 99热6这里只有精品| 插逼视频在线观看| 亚洲,一卡二卡三卡| a级一级毛片免费在线观看| 日本wwww免费看| 久久久久性生活片| 我的老师免费观看完整版| 一二三四中文在线观看免费高清| 亚洲无线观看免费| 国产精品人妻久久久影院| 中文字幕制服av| 久久久亚洲精品成人影院| 国产免费一区二区三区四区乱码| 久久亚洲国产成人精品v| 爱豆传媒免费全集在线观看| 国产久久久一区二区三区| 亚洲一区二区三区欧美精品| 中文乱码字字幕精品一区二区三区| 精品酒店卫生间| 亚洲精品国产色婷婷电影| 最近最新中文字幕大全电影3| 国产精品无大码| 久久久久精品性色| 国产精品一区二区三区四区免费观看| 免费大片18禁| 成人午夜精彩视频在线观看| 精品一区二区免费观看| 如何舔出高潮| 少妇人妻一区二区三区视频| 日韩一区二区三区影片| 成人毛片a级毛片在线播放| 国产精品女同一区二区软件| 人人妻人人看人人澡| 各种免费的搞黄视频| 亚洲欧美精品专区久久| 亚洲欧美中文字幕日韩二区| 三级国产精品欧美在线观看| 最后的刺客免费高清国语| 国产成人freesex在线| 欧美国产精品一级二级三级 | 在线天堂最新版资源| 久久99热这里只频精品6学生| 免费人成在线观看视频色| 亚洲av成人精品一二三区| 99热这里只有是精品在线观看| 亚洲久久久国产精品| 亚洲av成人精品一二三区| 麻豆国产97在线/欧美| 18禁在线播放成人免费| 亚洲精品中文字幕在线视频 | 日韩一区二区三区影片| 国产精品一区二区三区四区免费观看| 亚洲自偷自拍三级| 少妇人妻精品综合一区二区| 日韩在线高清观看一区二区三区| 国产精品无大码| 国产免费视频播放在线视频| 国产深夜福利视频在线观看| 亚洲经典国产精华液单| 免费不卡的大黄色大毛片视频在线观看| 男人爽女人下面视频在线观看| 国产精品国产av在线观看| 欧美日韩精品成人综合77777| 欧美xxxx黑人xx丫x性爽| 黄色欧美视频在线观看| 国产永久视频网站| 在线观看美女被高潮喷水网站| 亚洲美女搞黄在线观看| 制服丝袜香蕉在线| 久久久久人妻精品一区果冻| 亚洲精品国产av蜜桃| 天堂中文最新版在线下载| 一级毛片久久久久久久久女| 亚洲第一av免费看| 欧美bdsm另类| 韩国av在线不卡| 插逼视频在线观看| 中国国产av一级| 伊人久久国产一区二区| 男人舔奶头视频| 91精品国产国语对白视频| 男女无遮挡免费网站观看| 日韩精品有码人妻一区| 久久人人爽人人片av| 国产精品久久久久久久电影| 国精品久久久久久国模美| av专区在线播放| 日日撸夜夜添| 欧美日韩亚洲高清精品| 三级国产精品欧美在线观看| 久久ye,这里只有精品| 中国国产av一级| 日韩成人av中文字幕在线观看| 国产精品久久久久久精品古装| 女性生殖器流出的白浆| 日韩三级伦理在线观看| 少妇人妻精品综合一区二区| 高清视频免费观看一区二区| 色婷婷av一区二区三区视频| 一本一本综合久久| 尾随美女入室| 欧美亚洲 丝袜 人妻 在线| 一级毛片电影观看| 久久精品久久久久久久性| 久久这里有精品视频免费| 日韩av不卡免费在线播放| 亚洲,一卡二卡三卡| 国产成人a区在线观看| 久久精品人妻少妇| 在线观看免费高清a一片| 国产伦精品一区二区三区四那| 五月伊人婷婷丁香| 99热这里只有是精品在线观看| 22中文网久久字幕| 免费在线观看成人毛片| 欧美亚洲 丝袜 人妻 在线| 午夜激情福利司机影院| av在线观看视频网站免费| 久久人人爽人人片av| 亚洲精品乱码久久久久久按摩| 欧美精品亚洲一区二区| 国产永久视频网站| 99热全是精品| 亚洲美女黄色视频免费看| 天堂俺去俺来也www色官网| 亚洲精品久久久久久婷婷小说| 51国产日韩欧美| 亚洲久久久国产精品| 少妇人妻久久综合中文| 秋霞在线观看毛片| 夜夜看夜夜爽夜夜摸| 黄色怎么调成土黄色| 男人添女人高潮全过程视频| 26uuu在线亚洲综合色| 99热6这里只有精品| 黄色配什么色好看| 午夜激情久久久久久久| 精品99又大又爽又粗少妇毛片| 国内少妇人妻偷人精品xxx网站| 久久韩国三级中文字幕| a级一级毛片免费在线观看| 国产成人91sexporn| 麻豆精品久久久久久蜜桃| 高清午夜精品一区二区三区| 一级毛片久久久久久久久女| 国产免费又黄又爽又色| 一区二区三区精品91| 国产精品熟女久久久久浪| 国产成人freesex在线| 夫妻午夜视频| 国产深夜福利视频在线观看| 国产精品福利在线免费观看| 日本一二三区视频观看| 久久鲁丝午夜福利片| 国产精品三级大全| 日韩在线高清观看一区二区三区| 久久久久人妻精品一区果冻| 国产探花极品一区二区| 大片免费播放器 马上看| 欧美日韩综合久久久久久| 国语对白做爰xxxⅹ性视频网站| 欧美日韩一区二区视频在线观看视频在线| 美女cb高潮喷水在线观看| 久久人妻熟女aⅴ| 国产片特级美女逼逼视频| 边亲边吃奶的免费视频| 99热这里只有是精品在线观看| 国产精品偷伦视频观看了| 亚洲精品456在线播放app| 在线观看一区二区三区| 中文在线观看免费www的网站| 黑人高潮一二区| 久久婷婷青草| 18禁在线播放成人免费| 女性生殖器流出的白浆| 99久久精品国产国产毛片| 1000部很黄的大片| 久久精品国产亚洲网站| 精品人妻视频免费看| 如何舔出高潮| 久久久久久久久大av| 国产大屁股一区二区在线视频| 亚洲精品乱久久久久久| 成人午夜精彩视频在线观看| 视频中文字幕在线观看| av天堂中文字幕网| 国产片特级美女逼逼视频| 亚洲激情五月婷婷啪啪| 韩国高清视频一区二区三区| 少妇丰满av| 精品亚洲成a人片在线观看 | 精品99又大又爽又粗少妇毛片| 久久99热这里只有精品18| 亚洲图色成人| 又粗又硬又长又爽又黄的视频| 国产av一区二区精品久久 | 午夜激情福利司机影院| 丝袜脚勾引网站| 亚洲不卡免费看| 色综合色国产| 大片电影免费在线观看免费| 51国产日韩欧美| h日本视频在线播放| 少妇裸体淫交视频免费看高清| 国产爽快片一区二区三区| 国产亚洲5aaaaa淫片| 国产片特级美女逼逼视频| 成人亚洲精品一区在线观看 | 国产老妇伦熟女老妇高清| 伦理电影大哥的女人| 欧美一区二区亚洲| 另类亚洲欧美激情| 国产成人精品婷婷| 精品亚洲成国产av| 亚洲怡红院男人天堂| 国产精品久久久久久久久免| 国产精品偷伦视频观看了| 建设人人有责人人尽责人人享有的 | 日本vs欧美在线观看视频 | 精品国产露脸久久av麻豆| 精品久久久精品久久久| 日韩一本色道免费dvd| 精品一区二区三区视频在线| 亚洲欧美精品专区久久| 国产乱人偷精品视频| av又黄又爽大尺度在线免费看| 国产av一区二区精品久久 | 亚洲欧美中文字幕日韩二区| 汤姆久久久久久久影院中文字幕| 内地一区二区视频在线| 一边亲一边摸免费视频| 亚洲国产精品999| 久久 成人 亚洲| av.在线天堂| 亚洲精品第二区| 亚洲天堂av无毛| 乱码一卡2卡4卡精品| 免费av中文字幕在线| av专区在线播放| 欧美另类一区| 亚洲av中文字字幕乱码综合| 观看av在线不卡| 99re6热这里在线精品视频| 18禁裸乳无遮挡动漫免费视频| 亚洲av不卡在线观看| 亚洲,欧美,日韩| 一边亲一边摸免费视频| 高清在线视频一区二区三区| 26uuu在线亚洲综合色| 丝袜脚勾引网站| 黄色日韩在线| 男女边吃奶边做爰视频| 成年av动漫网址| 日韩精品有码人妻一区| 晚上一个人看的免费电影| 观看av在线不卡| 国产一级毛片在线| 中文欧美无线码| 亚洲精品日本国产第一区| 在线精品无人区一区二区三 | 国产精品久久久久久av不卡| 国产欧美日韩精品一区二区| 久久国产乱子免费精品| 在线免费观看不下载黄p国产| 精品酒店卫生间| 国产亚洲91精品色在线| 国产 一区 欧美 日韩| 免费高清在线观看视频在线观看| 国产高清有码在线观看视频| 久久久久久久大尺度免费视频| 亚洲国产高清在线一区二区三| 综合色丁香网| 免费看不卡的av| 国产一区二区在线观看日韩| 亚洲成人一二三区av| 色婷婷av一区二区三区视频| 尤物成人国产欧美一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91 | 另类亚洲欧美激情| 777米奇影视久久| 国产毛片在线视频| 超碰av人人做人人爽久久| 久久国产精品大桥未久av | 婷婷色综合www| 欧美精品一区二区大全| 黄色视频在线播放观看不卡| 欧美日韩视频精品一区| 亚洲成人一二三区av| 日日啪夜夜撸| 精品亚洲乱码少妇综合久久| 亚洲av.av天堂| 国产淫语在线视频| 久久久久视频综合| 国产精品人妻久久久久久| 成人18禁高潮啪啪吃奶动态图 | 99热网站在线观看| 在线观看国产h片| 深爱激情五月婷婷| av在线老鸭窝| 日本av免费视频播放| 久久久久精品性色| 久久鲁丝午夜福利片| 亚洲欧洲国产日韩| 国产视频内射| 国内少妇人妻偷人精品xxx网站| 国产 精品1| 成人高潮视频无遮挡免费网站| 免费观看a级毛片全部| 99久国产av精品国产电影| 亚洲色图av天堂| av在线观看视频网站免费| 日韩制服骚丝袜av| 少妇 在线观看| 黄色一级大片看看| 成人18禁高潮啪啪吃奶动态图 | 国产精品久久久久久精品电影小说 | 大话2 男鬼变身卡| 在线观看国产h片| 国产亚洲av片在线观看秒播厂| 免费高清在线观看视频在线观看| 最近的中文字幕免费完整| 国产成人精品一,二区| 日韩一本色道免费dvd| 亚洲va在线va天堂va国产| 午夜福利网站1000一区二区三区| 能在线免费看毛片的网站| 你懂的网址亚洲精品在线观看| 91精品伊人久久大香线蕉| 91精品一卡2卡3卡4卡| 丰满乱子伦码专区| 伊人久久国产一区二区| av免费在线看不卡| 干丝袜人妻中文字幕| 亚洲av日韩在线播放| 一区二区三区四区激情视频| 青春草视频在线免费观看| 日韩亚洲欧美综合| 国产精品偷伦视频观看了| 久久精品国产鲁丝片午夜精品| 午夜激情久久久久久久| 国产精品不卡视频一区二区| 亚洲美女视频黄频| 亚洲精品456在线播放app| 免费观看在线日韩| 色婷婷久久久亚洲欧美| 国产伦理片在线播放av一区| 各种免费的搞黄视频| 色视频www国产| 街头女战士在线观看网站| 亚洲国产最新在线播放| 国产探花极品一区二区| 国精品久久久久久国模美| 在线播放无遮挡| 精品久久久久久久久亚洲| 高清不卡的av网站| 人人妻人人看人人澡| 亚洲一区二区三区欧美精品| 97热精品久久久久久| 看十八女毛片水多多多| 男女国产视频网站| 国产伦理片在线播放av一区| 亚洲色图综合在线观看| 久久亚洲国产成人精品v| 亚洲欧美日韩卡通动漫| 国产在线男女| 日韩强制内射视频| 国产高清三级在线| av一本久久久久| 亚洲欧美日韩东京热| 亚洲伊人久久精品综合| 亚洲国产成人一精品久久久| 国产男女超爽视频在线观看| 欧美另类一区| 日日摸夜夜添夜夜添av毛片| 最黄视频免费看|