• <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.色视频.com| 亚洲国产成人一精品久久久| 亚洲精品国产色婷婷电影| 97超碰精品成人国产| 丁香六月天网| 丝瓜视频免费看黄片| 国产精品99久久99久久久不卡 | 亚洲图色成人| 999精品在线视频| 如何舔出高潮| 九九爱精品视频在线观看| 乱码一卡2卡4卡精品| av在线老鸭窝| 久久婷婷青草| 人妻 亚洲 视频| av视频免费观看在线观看| 欧美+日韩+精品| 久久精品国产亚洲av涩爱| 亚洲国产日韩一区二区| 老司机影院毛片| 久久99蜜桃精品久久| √禁漫天堂资源中文www| 欧美 亚洲 国产 日韩一| 亚洲国产精品国产精品| 大陆偷拍与自拍| 成人黄色视频免费在线看| 国产精品免费大片| 国产片特级美女逼逼视频| 新久久久久国产一级毛片| 90打野战视频偷拍视频| 十分钟在线观看高清视频www| 亚洲成人手机| 男人舔女人的私密视频| 亚洲中文av在线| 69精品国产乱码久久久| 卡戴珊不雅视频在线播放| 国产精品久久久久久精品古装| 午夜免费男女啪啪视频观看| 一级,二级,三级黄色视频| 国产欧美另类精品又又久久亚洲欧美| 中文字幕精品免费在线观看视频 | 亚洲av福利一区| 亚洲精品国产av成人精品| 极品少妇高潮喷水抽搐| 男女免费视频国产| 最近最新中文字幕免费大全7| av卡一久久| 久久久久久人人人人人| 亚洲经典国产精华液单| 成人国语在线视频| 热99国产精品久久久久久7| av在线老鸭窝| 亚洲国产精品一区二区三区在线| 国产国拍精品亚洲av在线观看| av不卡在线播放| 老女人水多毛片| 国产男女超爽视频在线观看| 天天躁夜夜躁狠狠久久av| 色吧在线观看| a级毛片黄视频| av黄色大香蕉| 中文字幕另类日韩欧美亚洲嫩草| 久久人人爽人人爽人人片va| 在线免费观看不下载黄p国产| 久久婷婷青草| 街头女战士在线观看网站| 两个人看的免费小视频| 久久久久久人人人人人| 一级毛片黄色毛片免费观看视频| 2018国产大陆天天弄谢| 麻豆乱淫一区二区| 国产欧美亚洲国产| 亚洲av.av天堂| 久久精品国产综合久久久 | 少妇的逼好多水| 在现免费观看毛片| 久久久欧美国产精品| 最新的欧美精品一区二区| 国产精品欧美亚洲77777| 亚洲欧美一区二区三区国产| 人人妻人人爽人人添夜夜欢视频| 又大又黄又爽视频免费| 9191精品国产免费久久| 精品亚洲成a人片在线观看| 色网站视频免费| 亚洲国产精品专区欧美| 大香蕉97超碰在线| 国产高清三级在线| 高清黄色对白视频在线免费看| 久久久国产一区二区| 9热在线视频观看99| 成年动漫av网址| 精品人妻偷拍中文字幕| 美女xxoo啪啪120秒动态图| 丰满饥渴人妻一区二区三| 国产毛片在线视频| 看非洲黑人一级黄片| 免费观看无遮挡的男女| 精品亚洲成国产av| 22中文网久久字幕| 亚洲av电影在线观看一区二区三区| 超碰97精品在线观看| 国产一区二区三区综合在线观看 | 国产日韩欧美在线精品| 1024视频免费在线观看| 中文字幕av电影在线播放| 18禁观看日本| 大话2 男鬼变身卡| 国产精品女同一区二区软件| 日韩制服丝袜自拍偷拍| 国产成人精品久久久久久| 精品久久国产蜜桃| 999精品在线视频| 亚洲经典国产精华液单| 国产精品三级大全| 香蕉精品网在线| videosex国产| 有码 亚洲区| 一级黄片播放器| 久久久久久久久久成人| 少妇人妻 视频| 人妻人人澡人人爽人人| 午夜免费男女啪啪视频观看| 国产极品粉嫩免费观看在线| 国产片特级美女逼逼视频| 久久毛片免费看一区二区三区| 超色免费av| 久久综合国产亚洲精品| 国产爽快片一区二区三区| 大陆偷拍与自拍| 欧美精品国产亚洲| 麻豆精品久久久久久蜜桃| 久久亚洲国产成人精品v| 国产一区二区三区av在线| 久久国产亚洲av麻豆专区| 伊人久久国产一区二区| 视频在线观看一区二区三区| 啦啦啦视频在线资源免费观看| 亚洲图色成人| 伊人亚洲综合成人网| 欧美精品av麻豆av| 久久精品国产a三级三级三级| 国产有黄有色有爽视频| 制服人妻中文乱码| 激情视频va一区二区三区| 成年女人在线观看亚洲视频| 欧美激情国产日韩精品一区| 1024视频免费在线观看| 久久这里只有精品19| 亚洲,一卡二卡三卡| 91精品伊人久久大香线蕉| 亚洲图色成人| 成人亚洲精品一区在线观看| 观看av在线不卡| 黄色怎么调成土黄色| 亚洲欧美清纯卡通| 精品一区二区三区视频在线| 老女人水多毛片| 在线精品无人区一区二区三| 人体艺术视频欧美日本| 国产一区二区三区av在线| 国产成人av激情在线播放| 韩国av在线不卡| 成人综合一区亚洲| 亚洲成色77777| 欧美日韩精品成人综合77777| 亚洲精品乱久久久久久| 青春草亚洲视频在线观看| 国产片特级美女逼逼视频| 乱人伦中国视频| 精品人妻熟女毛片av久久网站| 在线看a的网站| av片东京热男人的天堂| 亚洲 欧美一区二区三区| 国产精品一二三区在线看| 天天影视国产精品| 99精国产麻豆久久婷婷| 国产一区二区在线观看日韩| 亚洲精品国产色婷婷电影| 在线 av 中文字幕| 成人综合一区亚洲| 欧美日韩一区二区视频在线观看视频在线| 黑人欧美特级aaaaaa片| 日本欧美视频一区| 26uuu在线亚洲综合色| 色婷婷久久久亚洲欧美| av有码第一页| 有码 亚洲区| 中文字幕亚洲精品专区| 久久免费观看电影| 一区二区三区乱码不卡18| 国产又色又爽无遮挡免| 色5月婷婷丁香| av免费观看日本| 亚洲av免费高清在线观看| 18禁国产床啪视频网站| 国产一区二区在线观看av| 欧美xxxx性猛交bbbb| 成人免费观看视频高清| 亚洲国产精品一区三区| 国产又爽黄色视频| 日韩制服骚丝袜av| 啦啦啦啦在线视频资源| 美女大奶头黄色视频| 啦啦啦视频在线资源免费观看| 亚洲三级黄色毛片| 国产一区二区在线观看av| 2021少妇久久久久久久久久久| a级片在线免费高清观看视频| 精品久久久久久电影网| 国产 精品1| 五月开心婷婷网| 国产激情久久老熟女| 中文字幕制服av| 纵有疾风起免费观看全集完整版| 精品卡一卡二卡四卡免费| 日韩一区二区视频免费看| 美女主播在线视频| 欧美精品亚洲一区二区| 全区人妻精品视频| 少妇猛男粗大的猛烈进出视频| 精品视频人人做人人爽| 我的女老师完整版在线观看| 天天躁夜夜躁狠狠躁躁| 亚洲精品久久久久久婷婷小说| av国产久精品久网站免费入址| 午夜免费男女啪啪视频观看| 啦啦啦在线观看免费高清www| 亚洲精品456在线播放app| 男女免费视频国产| 一二三四中文在线观看免费高清| 亚洲精品一二三| 我的女老师完整版在线观看| 男人舔女人的私密视频| 亚洲熟女精品中文字幕| 一级黄片播放器| 熟女电影av网| av卡一久久| 日韩在线高清观看一区二区三区| 亚洲内射少妇av| 亚洲伊人色综图| 久久99蜜桃精品久久| 国产精品熟女久久久久浪| 岛国毛片在线播放| 97在线视频观看| 99久久人妻综合| 日产精品乱码卡一卡2卡三| 欧美人与性动交α欧美软件 | 国产麻豆69| 成人18禁高潮啪啪吃奶动态图| 一区在线观看完整版| 午夜久久久在线观看| 国产精品久久久久久av不卡| 欧美成人精品欧美一级黄| 人人妻人人澡人人爽人人夜夜| 美女国产视频在线观看| 好男人视频免费观看在线| 国产色婷婷99| 国产精品 国内视频| 一级a做视频免费观看| 黄色视频在线播放观看不卡| 天天影视国产精品| 国产成人精品在线电影| 极品人妻少妇av视频| 男男h啪啪无遮挡| 亚洲少妇的诱惑av| 精品少妇黑人巨大在线播放| 69精品国产乱码久久久| 嫩草影院入口| 精品国产一区二区三区四区第35| 亚洲少妇的诱惑av| 亚洲国产欧美日韩在线播放| 久久精品熟女亚洲av麻豆精品| 日本免费在线观看一区| 国产精品免费大片| www日本在线高清视频| 久久久a久久爽久久v久久| 亚洲成人av在线免费| 视频区图区小说| 美国免费a级毛片| 亚洲欧美色中文字幕在线| 在线看a的网站| 男女高潮啪啪啪动态图| 激情视频va一区二区三区| 成人免费观看视频高清| 国产欧美另类精品又又久久亚洲欧美| 国产一区有黄有色的免费视频| 午夜日本视频在线| 欧美精品av麻豆av| 秋霞在线观看毛片| 久久久久久久久久人人人人人人| xxx大片免费视频| 亚洲天堂av无毛| 91精品伊人久久大香线蕉| 午夜福利乱码中文字幕| 久久99蜜桃精品久久| 国产精品熟女久久久久浪| 中国美白少妇内射xxxbb| 草草在线视频免费看| 在现免费观看毛片| 不卡视频在线观看欧美| 亚洲综合色惰| 99热6这里只有精品| 成人亚洲精品一区在线观看| 亚洲美女视频黄频| 精品一区二区三区视频在线| 成年人免费黄色播放视频| 免费播放大片免费观看视频在线观看| 又黄又粗又硬又大视频| 国产不卡av网站在线观看| 最新中文字幕久久久久| 国产成人aa在线观看| 性高湖久久久久久久久免费观看| 免费久久久久久久精品成人欧美视频 | 91午夜精品亚洲一区二区三区| 免费观看性生交大片5| 极品少妇高潮喷水抽搐| 日韩av在线免费看完整版不卡| 久久青草综合色| 久久国产精品男人的天堂亚洲 | 久久久a久久爽久久v久久| 人成视频在线观看免费观看| 久久午夜综合久久蜜桃| 精品人妻一区二区三区麻豆| 交换朋友夫妻互换小说| 黄片播放在线免费| 熟妇人妻不卡中文字幕| 亚洲一码二码三码区别大吗| 黑人猛操日本美女一级片| 在线天堂中文资源库| 亚洲一区二区三区欧美精品| 久久久久视频综合| av不卡在线播放| 99国产精品免费福利视频| 成人国产av品久久久| a级毛色黄片| 十八禁网站网址无遮挡| 精品午夜福利在线看| 少妇 在线观看| 欧美日本中文国产一区发布| 中文字幕人妻熟女乱码| 2021少妇久久久久久久久久久| 国产成人精品福利久久| 爱豆传媒免费全集在线观看| 国产无遮挡羞羞视频在线观看| 最近中文字幕2019免费版| 亚洲av.av天堂| 国产精品99久久99久久久不卡 | 午夜激情av网站| 99re6热这里在线精品视频| 最近最新中文字幕免费大全7| a级毛片黄视频| 久久久久久久久久久免费av| 97在线视频观看| 伊人久久国产一区二区| 国产精品 国内视频| 亚洲精品av麻豆狂野| 亚洲伊人色综图| 汤姆久久久久久久影院中文字幕| 老司机亚洲免费影院| 成年动漫av网址| 最黄视频免费看| 精品卡一卡二卡四卡免费| 久久婷婷青草| 亚洲美女视频黄频| 久久久久精品性色| 欧美老熟妇乱子伦牲交| 日韩制服骚丝袜av| 少妇 在线观看| 成人亚洲精品一区在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 国产欧美日韩一区二区三区在线| 久久久亚洲精品成人影院| 热re99久久国产66热| 97超碰精品成人国产| 日本vs欧美在线观看视频| 国产精品一区www在线观看| 欧美+日韩+精品| 日韩三级伦理在线观看| 狠狠精品人妻久久久久久综合| 精品福利永久在线观看| 毛片一级片免费看久久久久| 亚洲欧美成人精品一区二区| 日韩三级伦理在线观看| 精品亚洲乱码少妇综合久久| 亚洲天堂av无毛| 熟女人妻精品中文字幕| 黄色配什么色好看| 男女国产视频网站| 成年av动漫网址| 在线观看免费日韩欧美大片| 亚洲国产精品成人久久小说| 欧美xxⅹ黑人| 精品久久久精品久久久| 亚洲经典国产精华液单| 免费av中文字幕在线| 18禁观看日本| 成人毛片a级毛片在线播放| 热re99久久国产66热| 午夜福利,免费看| 国产亚洲欧美精品永久| 国产欧美日韩一区二区三区在线| 日韩人妻精品一区2区三区| 一级毛片我不卡| 国产精品国产av在线观看| 国产成人91sexporn| 久久久久久人妻| 婷婷色综合大香蕉| 97人妻天天添夜夜摸| 久久国产精品大桥未久av| 久久久久人妻精品一区果冻| 国产精品久久久av美女十八| 亚洲色图 男人天堂 中文字幕 | 久久国产精品男人的天堂亚洲 | 久久久久久久大尺度免费视频| 90打野战视频偷拍视频| 亚洲欧洲精品一区二区精品久久久 | 乱码一卡2卡4卡精品| 亚洲av中文av极速乱| 亚洲精品国产av成人精品| 久久久国产精品麻豆| 成人影院久久| 亚洲伊人色综图| 亚洲精品乱码久久久久久按摩| 啦啦啦视频在线资源免费观看| 精品一区二区三卡| 人人妻人人澡人人看| 黄网站色视频无遮挡免费观看| 久久99一区二区三区| 成人无遮挡网站| 黄色一级大片看看| 美女xxoo啪啪120秒动态图| 秋霞伦理黄片| 国产精品人妻久久久影院| 在线亚洲精品国产二区图片欧美| 天天躁夜夜躁狠狠久久av| 嫩草影院入口| 91精品三级在线观看| 久久久精品94久久精品| 天天影视国产精品| 色吧在线观看| 丝瓜视频免费看黄片| 汤姆久久久久久久影院中文字幕| 亚洲四区av| 久久久久久久久久人人人人人人| 中文字幕av电影在线播放| 尾随美女入室| 欧美国产精品一级二级三级| 久久久国产欧美日韩av| 少妇的逼水好多| 黄片无遮挡物在线观看| 热99国产精品久久久久久7| 国产成人精品在线电影| 两个人看的免费小视频| 久久这里有精品视频免费| 爱豆传媒免费全集在线观看| 自线自在国产av| 国产精品无大码| 精品一区二区三区视频在线| 国产av码专区亚洲av| av国产久精品久网站免费入址| 丁香六月天网| av在线观看视频网站免费| 国产 精品1| 午夜福利网站1000一区二区三区| 日本欧美国产在线视频| 哪个播放器可以免费观看大片| 丝袜脚勾引网站| 香蕉丝袜av| videossex国产| 国产在线视频一区二区| 亚洲欧洲精品一区二区精品久久久 | kizo精华| 国产视频首页在线观看| 老女人水多毛片| 国产欧美日韩综合在线一区二区| 免费看av在线观看网站| 蜜臀久久99精品久久宅男| 精品一区二区三区视频在线| 日本-黄色视频高清免费观看| 高清欧美精品videossex| 欧美xxⅹ黑人| 9热在线视频观看99| 亚洲一码二码三码区别大吗| 一级毛片我不卡| 久久久国产精品麻豆| 国精品久久久久久国模美| 亚洲成av片中文字幕在线观看 | 成人国语在线视频| 亚洲伊人色综图| 少妇人妻久久综合中文| 女的被弄到高潮叫床怎么办| 男女啪啪激烈高潮av片| 国产黄色视频一区二区在线观看| 亚洲av日韩在线播放| 热99国产精品久久久久久7| 熟女电影av网| 热99国产精品久久久久久7| 人成视频在线观看免费观看| 婷婷色综合www| 成年美女黄网站色视频大全免费| 亚洲综合精品二区| 天堂中文最新版在线下载| 2021少妇久久久久久久久久久| 一级毛片黄色毛片免费观看视频| 精品国产国语对白av| 免费高清在线观看日韩| 午夜免费鲁丝| 熟女av电影| 国产男女内射视频| 午夜影院在线不卡| 色吧在线观看| av在线app专区| 国产麻豆69| 精品人妻一区二区三区麻豆| www.色视频.com| 日日啪夜夜爽| av又黄又爽大尺度在线免费看| 在线观看人妻少妇| 美女国产视频在线观看| 乱码一卡2卡4卡精品| 各种免费的搞黄视频| 亚洲欧美色中文字幕在线| 一本—道久久a久久精品蜜桃钙片| 精品亚洲成国产av| 七月丁香在线播放| 亚洲国产毛片av蜜桃av| 五月玫瑰六月丁香| 欧美97在线视频| 日韩三级伦理在线观看| 国产色爽女视频免费观看| 欧美日韩视频精品一区| 少妇人妻精品综合一区二区| 免费人妻精品一区二区三区视频| 妹子高潮喷水视频| av免费观看日本| 狠狠精品人妻久久久久久综合| 少妇人妻久久综合中文| 夜夜爽夜夜爽视频| 一级毛片电影观看| 午夜福利乱码中文字幕| 欧美精品国产亚洲| 国产成人午夜福利电影在线观看| 久久人人爽av亚洲精品天堂| 亚洲在久久综合| 26uuu在线亚洲综合色| 卡戴珊不雅视频在线播放| 全区人妻精品视频| 亚洲欧美一区二区三区国产| 欧美亚洲 丝袜 人妻 在线| 国产无遮挡羞羞视频在线观看| 欧美日韩成人在线一区二区| 欧美日韩av久久| 人体艺术视频欧美日本| 日韩免费高清中文字幕av| 天天操日日干夜夜撸| 免费黄色在线免费观看| 超色免费av| 日韩制服丝袜自拍偷拍| 女人精品久久久久毛片| 亚洲第一av免费看| 多毛熟女@视频| 日韩一区二区视频免费看| 国产日韩一区二区三区精品不卡| 蜜桃在线观看..| 在线观看美女被高潮喷水网站| 午夜福利,免费看| 天美传媒精品一区二区| 久久国产精品大桥未久av| 国产精品熟女久久久久浪| 国产黄色免费在线视频| 一边摸一边做爽爽视频免费| 亚洲欧美精品自产自拍| 又黄又粗又硬又大视频| 精品福利永久在线观看| 天堂俺去俺来也www色官网| 中文字幕人妻熟女乱码| 亚洲少妇的诱惑av| 各种免费的搞黄视频| 久久精品aⅴ一区二区三区四区 | 国产精品一国产av| 99久久人妻综合| 18禁国产床啪视频网站| 成人影院久久| 午夜福利,免费看| 99久久综合免费| 天堂中文最新版在线下载| 欧美成人午夜精品| av福利片在线| 日本91视频免费播放| 夫妻午夜视频| 免费大片黄手机在线观看| 好男人视频免费观看在线| 大话2 男鬼变身卡| 久久99热6这里只有精品| 亚洲美女视频黄频| 巨乳人妻的诱惑在线观看| 中文字幕最新亚洲高清| 亚洲中文av在线| 一级毛片我不卡| 成人毛片60女人毛片免费| 亚洲五月色婷婷综合| 最黄视频免费看| 午夜福利视频在线观看免费| 欧美精品av麻豆av| 色94色欧美一区二区| 在现免费观看毛片| 在线观看免费日韩欧美大片| 国产日韩欧美亚洲二区| 一级a做视频免费观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲精品,欧美精品| 成年av动漫网址| 日韩中文字幕视频在线看片| 两个人免费观看高清视频|