• <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)研究
    91成人精品电影| 欧美+日韩+精品| 久久精品国产亚洲av天美| 综合色丁香网| 国语对白做爰xxxⅹ性视频网站| av女优亚洲男人天堂| 午夜福利视频在线观看免费| 日韩,欧美,国产一区二区三区| 国产女主播在线喷水免费视频网站| 亚洲怡红院男人天堂| 黑人欧美特级aaaaaa片| 亚洲精品成人av观看孕妇| 少妇 在线观看| 午夜91福利影院| 久久久久久久亚洲中文字幕| 国产精品一二三区在线看| 九色亚洲精品在线播放| 91在线精品国自产拍蜜月| 色婷婷久久久亚洲欧美| 我要看黄色一级片免费的| 美女cb高潮喷水在线观看| 亚洲精品乱久久久久久| 精品国产一区二区久久| 亚洲熟女精品中文字幕| 26uuu在线亚洲综合色| 2022亚洲国产成人精品| 成人黄色视频免费在线看| 观看av在线不卡| 日韩制服骚丝袜av| 国产欧美另类精品又又久久亚洲欧美| 一级a做视频免费观看| 91久久精品国产一区二区三区| 亚洲av电影在线观看一区二区三区| 97在线视频观看| 免费大片18禁| 性色av一级| a级片在线免费高清观看视频| 日韩av不卡免费在线播放| 少妇的逼好多水| 一本一本综合久久| 少妇的逼水好多| 欧美精品人与动牲交sv欧美| 亚洲国产精品一区三区| 18禁观看日本| 精品亚洲成国产av| 99久久综合免费| 国产不卡av网站在线观看| 十八禁高潮呻吟视频| 久久精品国产亚洲av天美| 丝瓜视频免费看黄片| av黄色大香蕉| 日韩电影二区| 国产成人freesex在线| 久久久久久久久久成人| 国产精品嫩草影院av在线观看| 国产老妇伦熟女老妇高清| 亚洲av免费高清在线观看| 亚洲成人一二三区av| 国产亚洲av片在线观看秒播厂| 18禁观看日本| 九九爱精品视频在线观看| 国产黄色视频一区二区在线观看| 国产精品熟女久久久久浪| 国产精品一区www在线观看| 人妻系列 视频| 高清毛片免费看| 免费少妇av软件| 男女无遮挡免费网站观看| 夜夜爽夜夜爽视频| 人妻一区二区av| 久久国产精品大桥未久av| 成年美女黄网站色视频大全免费 | 久久久久视频综合| 免费大片18禁| 美女cb高潮喷水在线观看| 少妇 在线观看| 午夜福利视频在线观看免费| 伊人久久精品亚洲午夜| 大话2 男鬼变身卡| 如日韩欧美国产精品一区二区三区 | 99re6热这里在线精品视频| 久久亚洲国产成人精品v| 亚洲av二区三区四区| 午夜91福利影院| 99热全是精品| 欧美精品高潮呻吟av久久| 久久精品人人爽人人爽视色| 久久久亚洲精品成人影院| 日日爽夜夜爽网站| 一区二区三区乱码不卡18| 蜜桃国产av成人99| 午夜免费观看性视频| 国产白丝娇喘喷水9色精品| 日本与韩国留学比较| 亚洲内射少妇av| 精品久久久精品久久久| 999精品在线视频| 少妇的逼水好多| 高清午夜精品一区二区三区| 欧美日韩视频精品一区| 国产视频首页在线观看| 中文乱码字字幕精品一区二区三区| 午夜精品国产一区二区电影| 亚洲欧美成人精品一区二区| 成人毛片a级毛片在线播放| 美女国产视频在线观看| 国产亚洲午夜精品一区二区久久| 2021少妇久久久久久久久久久| 水蜜桃什么品种好| 三上悠亚av全集在线观看| 18+在线观看网站| 三级国产精品欧美在线观看| 精品一区二区免费观看| 青青草视频在线视频观看| 国产精品久久久久成人av| 亚洲国产精品国产精品| 亚洲在久久综合| 日韩一本色道免费dvd| 七月丁香在线播放| av在线播放精品| 久久久久精品性色| 亚洲少妇的诱惑av| 国产欧美亚洲国产| 超色免费av| 国产精品欧美亚洲77777| 中文字幕精品免费在线观看视频 | 五月天丁香电影| 国产视频内射| 亚洲三级黄色毛片| 成年美女黄网站色视频大全免费 | 乱码一卡2卡4卡精品| 精品99又大又爽又粗少妇毛片| 色婷婷久久久亚洲欧美| 熟女电影av网| 大话2 男鬼变身卡| 18在线观看网站| 日韩伦理黄色片| 自拍欧美九色日韩亚洲蝌蚪91| 一级毛片 在线播放| 久久久a久久爽久久v久久| 国产淫语在线视频| 亚洲成人av在线免费| 国产精品无大码| 国产精品三级大全| 国产精品三级大全| 天堂8中文在线网| 综合色丁香网| 一区二区三区四区激情视频| 国产永久视频网站| 亚洲欧美成人综合另类久久久| 国产亚洲欧美精品永久| 精品少妇内射三级| 精品人妻熟女av久视频| 亚洲国产成人一精品久久久| 亚洲人成网站在线观看播放| 国产成人精品在线电影| 欧美日韩精品成人综合77777| 亚洲情色 制服丝袜| 丰满饥渴人妻一区二区三| 91久久精品电影网| 久久亚洲国产成人精品v| 熟女人妻精品中文字幕| 亚洲图色成人| 久久精品人人爽人人爽视色| 国产成人免费观看mmmm| av国产精品久久久久影院| 99热这里只有是精品在线观看| 国产av一区二区精品久久| 91久久精品电影网| 五月玫瑰六月丁香| 国产 精品1| 成年av动漫网址| 少妇的逼水好多| 男女啪啪激烈高潮av片| 国产精品国产三级国产av玫瑰| 久久久久人妻精品一区果冻| 亚洲人与动物交配视频| 国产av一区二区精品久久| 18在线观看网站| 亚洲综合色网址| 老司机亚洲免费影院| 内地一区二区视频在线| 欧美变态另类bdsm刘玥| 高清午夜精品一区二区三区| 亚洲国产av影院在线观看| 免费大片18禁| 国产精品一区www在线观看| 国产一区有黄有色的免费视频| 亚洲中文av在线| 91国产中文字幕| 免费久久久久久久精品成人欧美视频 | 精品久久久噜噜| 伦理电影大哥的女人| 亚洲av.av天堂| 美女内射精品一级片tv| 亚洲国产最新在线播放| 亚洲欧美成人综合另类久久久| 国产一区亚洲一区在线观看| 啦啦啦中文免费视频观看日本| 中文字幕制服av| 国产欧美亚洲国产| 国产国语露脸激情在线看| 男男h啪啪无遮挡| 自线自在国产av| 国产精品一区二区在线观看99| 精品久久久久久电影网| 一区二区三区乱码不卡18| 欧美精品高潮呻吟av久久| 夫妻性生交免费视频一级片| 日本色播在线视频| 精品少妇久久久久久888优播| 久热久热在线精品观看| 欧美97在线视频| 美女脱内裤让男人舔精品视频| 最近最新中文字幕免费大全7| 大片电影免费在线观看免费| 亚洲四区av| 精品亚洲乱码少妇综合久久| 秋霞伦理黄片| 亚洲欧美精品自产自拍| 亚洲国产av影院在线观看| 久久 成人 亚洲| 国产视频内射| 欧美老熟妇乱子伦牲交| 久久久亚洲精品成人影院| 免费久久久久久久精品成人欧美视频 | 日韩视频在线欧美| 老司机影院毛片| 少妇人妻 视频| 美女福利国产在线| 精品99又大又爽又粗少妇毛片| 亚洲国产日韩一区二区| 国产一区二区三区av在线| 日日摸夜夜添夜夜爱| 又黄又爽又刺激的免费视频.| 伊人亚洲综合成人网| 中国三级夫妇交换| 简卡轻食公司| 两个人免费观看高清视频| 国产精品人妻久久久久久| 国产伦理片在线播放av一区| h视频一区二区三区| 亚洲色图 男人天堂 中文字幕 | 日本-黄色视频高清免费观看| 国产精品免费大片| 久久精品久久久久久久性| 亚洲精品美女久久av网站| 亚洲综合色惰| 国产精品一区www在线观看| 国产欧美日韩一区二区三区在线 | 人妻一区二区av| 精品久久久久久电影网| 亚洲精品国产av成人精品| 新久久久久国产一级毛片| 欧美日韩综合久久久久久| 国产一区二区在线观看av| av电影中文网址| 亚洲av成人精品一区久久| 在线天堂最新版资源| 我的女老师完整版在线观看| 一区二区三区四区激情视频| 久久久久久久大尺度免费视频| 久久精品国产亚洲av涩爱| 这个男人来自地球电影免费观看 | 最近的中文字幕免费完整| 成人午夜精彩视频在线观看| 亚洲精品日韩在线中文字幕| 日本猛色少妇xxxxx猛交久久| a级毛片黄视频| 丝袜喷水一区| 熟女av电影| 国产亚洲av片在线观看秒播厂| 国产成人精品在线电影| 欧美人与性动交α欧美精品济南到 | 亚洲国产欧美在线一区| 热99国产精品久久久久久7| 精品一区二区三卡| 伊人亚洲综合成人网| av天堂久久9| a级毛色黄片| 欧美+日韩+精品| 一区二区av电影网| 国模一区二区三区四区视频| 久久久久久人妻| 国产熟女欧美一区二区| 黑人高潮一二区| 在现免费观看毛片| 成人午夜精彩视频在线观看| 欧美日韩视频精品一区| 久久99热这里只频精品6学生| 国产色爽女视频免费观看| 亚洲欧美成人综合另类久久久| 韩国高清视频一区二区三区| 久久99精品国语久久久| 人妻系列 视频| 22中文网久久字幕| 亚洲国产色片| 亚洲av日韩在线播放| 亚洲国产最新在线播放| 亚洲成人av在线免费| 亚洲精品av麻豆狂野| 久久av网站| 亚洲精品国产色婷婷电影| 亚洲三级黄色毛片| 国产欧美另类精品又又久久亚洲欧美| av在线老鸭窝| 欧美人与性动交α欧美精品济南到 | 国产又色又爽无遮挡免| 极品人妻少妇av视频| 三上悠亚av全集在线观看| 欧美精品一区二区免费开放| 青春草国产在线视频| 日韩一本色道免费dvd| 女性被躁到高潮视频| 夜夜看夜夜爽夜夜摸| 91精品伊人久久大香线蕉| 十八禁网站网址无遮挡| 国产高清有码在线观看视频| 亚洲在久久综合| 国产日韩一区二区三区精品不卡 | 国产一区亚洲一区在线观看| 午夜影院在线不卡| 啦啦啦在线观看免费高清www| 老女人水多毛片| 亚洲国产精品一区三区| 国产色婷婷99| 大片电影免费在线观看免费| 亚洲第一av免费看| av免费在线看不卡| 午夜福利视频在线观看免费| 精品久久蜜臀av无| 精品一区二区三卡| 亚洲欧洲日产国产| 久热这里只有精品99| 亚洲国产欧美在线一区| 久久精品熟女亚洲av麻豆精品| 爱豆传媒免费全集在线观看| 免费少妇av软件| 爱豆传媒免费全集在线观看| 精品一区二区三卡| 精品久久国产蜜桃| 久久久欧美国产精品| 欧美另类一区| 日韩成人伦理影院| 久久综合国产亚洲精品| 国产探花极品一区二区| 免费黄色在线免费观看| 中文字幕久久专区| 色网站视频免费| 伊人久久精品亚洲午夜| 亚洲国产av新网站| 一边摸一边做爽爽视频免费| 人人妻人人澡人人爽人人夜夜| 啦啦啦在线观看免费高清www| 亚洲国产欧美日韩在线播放| 22中文网久久字幕| 一区二区三区乱码不卡18| 日韩视频在线欧美| tube8黄色片| 欧美日本中文国产一区发布| 欧美精品国产亚洲| 纯流量卡能插随身wifi吗| 两个人的视频大全免费| 人人妻人人添人人爽欧美一区卜| 国产高清三级在线| 3wmmmm亚洲av在线观看| 最近的中文字幕免费完整| 肉色欧美久久久久久久蜜桃| 黄片播放在线免费| av黄色大香蕉| 这个男人来自地球电影免费观看 | 男男h啪啪无遮挡| 国产一区有黄有色的免费视频| 91久久精品电影网| 国产免费一区二区三区四区乱码| 尾随美女入室| 一区在线观看完整版| 嫩草影院入口| 如何舔出高潮| 美女脱内裤让男人舔精品视频| 亚洲欧美成人综合另类久久久| 精品卡一卡二卡四卡免费| 国产av一区二区精品久久| 日韩大片免费观看网站| 国产 精品1| 精品酒店卫生间| 一级毛片电影观看| 伊人久久国产一区二区| 亚洲不卡免费看| 国产免费现黄频在线看| 欧美另类一区| 国产一区有黄有色的免费视频| 精品久久久久久久久亚洲| 日韩 亚洲 欧美在线| 欧美日韩综合久久久久久| 青春草亚洲视频在线观看| 国产精品99久久99久久久不卡 | 日本av手机在线免费观看| 天天影视国产精品| 边亲边吃奶的免费视频| 久久国产亚洲av麻豆专区| 国产又色又爽无遮挡免| 中国三级夫妇交换| 天堂俺去俺来也www色官网| 人体艺术视频欧美日本| 国产免费视频播放在线视频| 18+在线观看网站| 色吧在线观看| 欧美丝袜亚洲另类| 边亲边吃奶的免费视频| 91精品伊人久久大香线蕉| 亚洲欧美成人精品一区二区| 日韩中文字幕视频在线看片| 三上悠亚av全集在线观看| 亚洲av中文av极速乱| 亚洲四区av| 成人亚洲精品一区在线观看| 9色porny在线观看| 黄色欧美视频在线观看| 亚洲怡红院男人天堂| 久久久精品区二区三区| 久久免费观看电影| 2021少妇久久久久久久久久久| 中文字幕人妻丝袜制服| 欧美 亚洲 国产 日韩一| 日韩,欧美,国产一区二区三区| 丰满乱子伦码专区| 麻豆成人av视频| 国产精品99久久99久久久不卡 | 高清视频免费观看一区二区| 777米奇影视久久| 成年av动漫网址| 在线观看免费日韩欧美大片 | 三级国产精品欧美在线观看| 大片免费播放器 马上看| 久久精品国产亚洲av涩爱| 丰满乱子伦码专区| 少妇精品久久久久久久| 看非洲黑人一级黄片| 中文乱码字字幕精品一区二区三区| 丰满乱子伦码专区| 又黄又爽又刺激的免费视频.| 久久精品久久久久久噜噜老黄| 性色av一级| 日韩一本色道免费dvd| 欧美 亚洲 国产 日韩一| 久久精品人人爽人人爽视色| 国产综合精华液| 国产不卡av网站在线观看| 久久久精品94久久精品| 日韩欧美一区视频在线观看| 亚洲av.av天堂| 国产欧美日韩一区二区三区在线 | 国产乱来视频区| 只有这里有精品99| 九九爱精品视频在线观看| 国产欧美另类精品又又久久亚洲欧美| 精品视频人人做人人爽| 观看美女的网站| 久久国产精品大桥未久av| 热99国产精品久久久久久7| 亚洲激情五月婷婷啪啪| 亚洲精品国产av蜜桃| 亚洲一级一片aⅴ在线观看| 久久久精品免费免费高清| 天美传媒精品一区二区| 99视频精品全部免费 在线| 黄色视频在线播放观看不卡| 欧美日韩视频精品一区| 久久久国产欧美日韩av| 午夜福利网站1000一区二区三区| 国产精品熟女久久久久浪| 欧美精品亚洲一区二区| 午夜福利在线观看免费完整高清在| 一区二区三区四区激情视频| 欧美+日韩+精品| 黄片播放在线免费| a级片在线免费高清观看视频| 777米奇影视久久| 久久久a久久爽久久v久久| 九九久久精品国产亚洲av麻豆| 国产欧美日韩一区二区三区在线 | 一区二区三区乱码不卡18| 色94色欧美一区二区| 少妇猛男粗大的猛烈进出视频| 欧美激情极品国产一区二区三区 | 精品少妇内射三级| 国产精品一区二区在线观看99| av卡一久久| 国产男女超爽视频在线观看| av卡一久久| 久久国产精品大桥未久av| 国语对白做爰xxxⅹ性视频网站| 波野结衣二区三区在线| 男男h啪啪无遮挡| 九草在线视频观看| 黄色毛片三级朝国网站| 精品酒店卫生间| 国产精品久久久久久av不卡| .国产精品久久| 好男人视频免费观看在线| 国产精品一区二区在线观看99| 美女内射精品一级片tv| 老女人水多毛片| 精品国产露脸久久av麻豆| 嫩草影院入口| 成人亚洲精品一区在线观看| 亚洲,一卡二卡三卡| 国产一区二区三区综合在线观看 | 99九九线精品视频在线观看视频| 熟女av电影| 一级二级三级毛片免费看| 亚洲国产精品999| 人人澡人人妻人| 欧美激情 高清一区二区三区| av一本久久久久| a级片在线免费高清观看视频| 国产精品女同一区二区软件| 欧美人与善性xxx| 观看美女的网站| 岛国毛片在线播放| 美女福利国产在线| 最近中文字幕2019免费版| 日韩欧美精品免费久久| 日本-黄色视频高清免费观看| 国产片特级美女逼逼视频| 亚洲欧美一区二区三区黑人 | 少妇人妻 视频| 99久久人妻综合| 久久久久久伊人网av| av在线老鸭窝| 国产精品免费大片| a级片在线免费高清观看视频| 人人妻人人添人人爽欧美一区卜| 黄色一级大片看看| 国产免费一区二区三区四区乱码| 国产av精品麻豆| 麻豆精品久久久久久蜜桃| 色吧在线观看| 亚洲在久久综合| 色婷婷久久久亚洲欧美| 26uuu在线亚洲综合色| 日韩欧美一区视频在线观看| 99久久综合免费| 最新的欧美精品一区二区| 国产亚洲精品第一综合不卡 | 飞空精品影院首页| 九色成人免费人妻av| 看免费成人av毛片| 亚洲欧美一区二区三区国产| 少妇人妻久久综合中文| 蜜桃在线观看..| 少妇猛男粗大的猛烈进出视频| 夜夜看夜夜爽夜夜摸| 曰老女人黄片| 国产精品国产av在线观看| 亚洲av国产av综合av卡| 久久精品国产亚洲网站| 日本与韩国留学比较| 亚洲不卡免费看| 最近最新中文字幕免费大全7| 午夜免费观看性视频| 人人澡人人妻人| 久久国产精品大桥未久av| av女优亚洲男人天堂| av播播在线观看一区| 亚洲图色成人| 少妇的逼好多水| 男女无遮挡免费网站观看| 日韩av免费高清视频| 久久av网站| 午夜91福利影院| av网站免费在线观看视频| 国产精品久久久久久精品古装| 亚洲av成人精品一二三区| 日韩视频在线欧美| 国产精品欧美亚洲77777| 看免费成人av毛片| 亚洲精品乱久久久久久| 欧美 日韩 精品 国产| 水蜜桃什么品种好| 哪个播放器可以免费观看大片| 亚洲综合精品二区| 亚洲精品美女久久av网站| 亚洲欧美日韩卡通动漫| 国产白丝娇喘喷水9色精品| 一本久久精品| 国产成人精品久久久久久| 久久精品夜色国产| 国产国语露脸激情在线看| 免费黄网站久久成人精品| 国产极品粉嫩免费观看在线 | 国产高清三级在线| 一级a做视频免费观看| 久久av网站| 精品少妇内射三级| 高清黄色对白视频在线免费看| 99热这里只有是精品在线观看| 看非洲黑人一级黄片| av在线老鸭窝| 国产视频首页在线观看| 少妇高潮的动态图| 亚洲四区av| 国产精品国产三级专区第一集| 午夜精品国产一区二区电影| 97超视频在线观看视频| 少妇被粗大的猛进出69影院 | 久久久国产一区二区| 尾随美女入室| 亚洲婷婷狠狠爱综合网| 午夜福利网站1000一区二区三区| 久久99热6这里只有精品| 少妇丰满av| 久久狼人影院| 日韩成人av中文字幕在线观看|