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

    基于遺傳算法的巖土及填料熱物性參數(shù)評(píng)估*

    2017-05-15 03:30:22張泉馬小威張林鋒
    關(guān)鍵詞:模型

    張泉,馬小威,張林鋒

    (湖南大學(xué) 土木工程學(xué)院,湖南 長(zhǎng)沙 410082)

    ?

    基于遺傳算法的巖土及填料熱物性參數(shù)評(píng)估*

    張泉?,馬小威,張林鋒

    (湖南大學(xué) 土木工程學(xué)院,湖南 長(zhǎng)沙 410082)

    巖土及填料的熱物性參數(shù)可直接影響地埋管換熱器的性能設(shè)計(jì).針對(duì)現(xiàn)場(chǎng)熱響應(yīng)測(cè)試中填料熱物性參數(shù)需要預(yù)先測(cè)量,而無(wú)法通過(guò)計(jì)算直接評(píng)估的難題,在無(wú)限長(zhǎng)線熱源模型及管壁熱阻修正復(fù)合介質(zhì)線熱源模型的基礎(chǔ)上,應(yīng)用遺傳算法對(duì)巖土及填料的熱物性參數(shù)同時(shí)進(jìn)行評(píng)估,并結(jié)合熱響應(yīng)測(cè)試實(shí)驗(yàn),驗(yàn)證了該方法的準(zhǔn)確性.鉆孔熱阻、巖土導(dǎo)熱系數(shù)和填料導(dǎo)熱系數(shù)的評(píng)估值與實(shí)測(cè)值的相對(duì)誤差分別為3.47%,1.42%和4.93%.2種模型計(jì)算所得流體平均溫度與實(shí)測(cè)值的均方根誤差分別為0.050 5 ℃和0.172 ℃.研究結(jié)果對(duì)地埋管換熱器的設(shè)計(jì)具有重要參考價(jià)值.

    填料熱物性;遺傳算法;復(fù)合介質(zhì)線熱源;熱響應(yīng)

    巖土及填料的熱物性參數(shù)是地埋管換熱器設(shè)計(jì)的關(guān)鍵參數(shù),準(zhǔn)確地評(píng)估巖土及填料的熱物性參數(shù)可以解決地埋管換熱器設(shè)計(jì)不當(dāng)?shù)膯?wèn)題.目前,巖土熱物性參數(shù)主要采用斜率法和參數(shù)估計(jì)法進(jìn)行評(píng)估計(jì)算.相比于斜率法,參數(shù)估計(jì)法可有效降低測(cè)試中熱流不恒定對(duì)評(píng)估參數(shù)值的影響,其評(píng)估精度較高[1].因此,在ASHRAE手冊(cè)中,推薦采用參數(shù)估計(jì)法進(jìn)行巖土熱物性參數(shù)的評(píng)估[2].

    參數(shù)估計(jì)法評(píng)估巖土熱物性參數(shù)是一個(gè)典型的逆?zhèn)鳠釂?wèn)題.該方法以巖土熱物性參數(shù)為決策變量,以各時(shí)刻流體平均溫度的計(jì)算值和實(shí)測(cè)值的平方和誤差(SSE)或均方根誤差(RMSE)為目標(biāo)函數(shù),通過(guò)最小化尋優(yōu)算法求解巖土熱物性參數(shù).但目前所采用的參數(shù)估計(jì)法存在以下2個(gè)方面的問(wèn)題:1) 主要運(yùn)用單純形法、模式搜索法等算法尋優(yōu)求解[3],其收斂速度較慢,并受迭代初值的影響,導(dǎo)致評(píng)估結(jié)果可能存在較大誤差.2) 參數(shù)估計(jì)法評(píng)估巖土熱物性參數(shù)時(shí),常以鉆孔壁溫度為計(jì)算耦合點(diǎn),忽略了鉆孔內(nèi)填料的熱物性,無(wú)法采用參數(shù)估計(jì)模型直接評(píng)估填料熱物性參數(shù).

    為克服參數(shù)估計(jì)法的上述缺點(diǎn),本文結(jié)合地埋管無(wú)限長(zhǎng)線源模型和管壁熱阻修正的復(fù)合介質(zhì)線熱源模型,利用遺傳算法同時(shí)評(píng)估鉆孔熱阻、巖土導(dǎo)熱系數(shù)和熱擴(kuò)散系數(shù)以及填料導(dǎo)熱系數(shù)和熱擴(kuò)散系數(shù)等5個(gè)參數(shù),并結(jié)合熱響應(yīng)測(cè)試實(shí)驗(yàn),驗(yàn)證了該方法的準(zhǔn)確性.

    1 地埋管傳熱模型

    1.1 無(wú)線長(zhǎng)線熱源模型

    Kelvin[4]提出了豎直地埋管鉆孔外一維無(wú)限長(zhǎng)線熱源傳熱模型,該模型通過(guò)鉆孔壁溫度耦合鉆孔內(nèi)穩(wěn)態(tài)導(dǎo)熱和鉆孔外非穩(wěn)態(tài)導(dǎo)熱.當(dāng)傳熱時(shí)間大于10 h時(shí),通過(guò)該簡(jiǎn)化模型可求得流體平均溫度為[5-6]:

    (1)

    式中:t為時(shí)間,s;T0為未擾動(dòng)土壤溫度,℃;ql為線熱源的強(qiáng)度,W/m;Rb為鉆孔熱阻,(m·K)/W;λs為巖土的導(dǎo)熱系數(shù),W/(m·K);rb為鉆孔半徑,m;as為巖土的熱擴(kuò)散系數(shù),m2/s.

    1.2 復(fù)合介質(zhì)線熱源模型

    在無(wú)限長(zhǎng)線熱源模型的基礎(chǔ)上,Li等[7]提出了一種全時(shí)間范圍內(nèi)的地埋管溫度響應(yīng)G函數(shù)模型.該模型將U型管管壁作為溫度耦合點(diǎn),以U型管外壁面A和B兩點(diǎn)的平均溫度代表U型管管壁的平均溫度,如圖1所示.

    圖1 單U型地埋管管壁Fig.1 Pipe wall for single U-shaped tube

    對(duì)于單U型地埋管,單位線熱源強(qiáng)度的管壁溫度響應(yīng)G函數(shù)為[7]:

    (2)

    式中:λb為填料的導(dǎo)熱系數(shù),W/(m·K);ab為填料的熱擴(kuò)散系數(shù),m2/s;rA和rB分別為點(diǎn)A和B的徑向坐標(biāo),m;r′為線熱源的徑向坐標(biāo),m;J2m為階數(shù)2m的第一類貝塞爾函數(shù).無(wú)量綱變量φ,ψ,f,g的計(jì)算式為:

    (3)

    根據(jù)式(3),結(jié)合流體至鉆孔外壁的熱阻Rp,可得流體平均溫度為[7]:

    Tf,CMLS(t)=T0+ql(GCMLS(t)+Rp).

    (4)

    由于傳熱初期鉆孔內(nèi)為非穩(wěn)態(tài)傳熱,若管壁熱阻Rp采用穩(wěn)態(tài)熱阻進(jìn)行計(jì)算,所得的流體平均溫度將產(chǎn)生較大誤差.為此,利用鉆孔熱阻Rb對(duì)管壁熱阻Rp進(jìn)行修正,管壁熱阻Rp的修正系數(shù)為[9]:

    η=2Rb-2Rbp-Rp.

    (5)

    結(jié)合式(4)和(5)即可得出修正管壁熱阻后流體平均溫度的計(jì)算式為:

    Tf,CMLS(t)=T0+ql(GCMLS(t)+Rp+η)=

    T0+ql(GCMLS(t)+Rp+2Rb-2Rbp-Rp).

    (6)

    復(fù)合介質(zhì)線熱源模型由于考慮了填料熱容和鉆孔埋管幾何尺度對(duì)換熱性能的影響,在一定時(shí)間范圍內(nèi),可較準(zhǔn)確地計(jì)算流體平均溫度.

    2 基于遺傳算法的熱物性參數(shù)評(píng)估

    遺傳算法(GeneticAlgorithm)是一種模仿自然界生物進(jìn)化,并帶有隨機(jī)性的全局搜索方法.該算法由可能潛在解集的種群開始計(jì)算,通過(guò)優(yōu)勝劣汰的準(zhǔn)則產(chǎn)生適應(yīng)度較高的解集(種群),在每一代進(jìn)化過(guò)程中,選擇種群中適應(yīng)度水平較高的個(gè)體進(jìn)行交叉重組以及變異,并產(chǎn)生子種群,由此通過(guò)一定代數(shù)的進(jìn)化即可得到最優(yōu)解的個(gè)體[10].與傳統(tǒng)優(yōu)化算法相比,遺傳算法具有較好的全局收斂性、計(jì)算速度快、不受目標(biāo)函數(shù)的約束、可并行計(jì)算等諸多優(yōu)點(diǎn),在工程計(jì)算中得到了廣泛應(yīng)用.

    本文遺傳算法采用謝菲爾德大學(xué)所開發(fā)的Matlab遺傳算法工具箱(gatbx)[11].求解中,設(shè)定鉆孔熱阻Rb,巖土導(dǎo)熱系數(shù)λs,巖土熱擴(kuò)散系數(shù)as,填料導(dǎo)熱系數(shù)λb和填料熱擴(kuò)散系數(shù)ab等5個(gè)參數(shù)為決策變量.以式(1)和式(6)計(jì)算所得的流體平均溫度與實(shí)測(cè)值的綜合均方根誤差為目標(biāo)函數(shù)[12],如式(7)所示.該目標(biāo)函數(shù)最小時(shí)所對(duì)應(yīng)的決策變量即為最終評(píng)估結(jié)果.

    (7)

    式中:Tf,ILS為無(wú)限長(zhǎng)線熱源模型計(jì)算的各時(shí)刻流體平均溫度,℃;Tf,CMLS為復(fù)合介質(zhì)線熱源模型計(jì)算的各時(shí)刻流體平均溫度,℃;Tf,EX為實(shí)測(cè)各時(shí)刻流體平均溫度,℃,計(jì)算中取各時(shí)刻U型管流體進(jìn)出口溫度的算術(shù)平均值;n為用于無(wú)限長(zhǎng)線熱源模型計(jì)算的測(cè)量數(shù)據(jù)組數(shù);m為用于復(fù)合介質(zhì)線熱源模型計(jì)算的測(cè)量數(shù)據(jù)組數(shù).

    在算法優(yōu)化過(guò)程中,幾個(gè)關(guān)鍵遺傳操作設(shè)置如下:

    1)編碼:為提高算法效率和求解精度,本文采用實(shí)數(shù)型編碼的方式對(duì)種群進(jìn)行編碼.

    2)適應(yīng)度函數(shù):目標(biāo)函數(shù)至適應(yīng)度函數(shù)的轉(zhuǎn)換,采用非線性排序算法且選擇壓差為2[11].

    3)交叉:采用隨機(jī)遍歷抽樣(SUS)選擇父體進(jìn)行交叉運(yùn)算.由于采用實(shí)數(shù)型編碼,交叉概率的取值不影響離散重組產(chǎn)生子代個(gè)體的表征,故將其設(shè)為1[11].

    4)變異:為防止算法局部收斂,變異運(yùn)算采用自適應(yīng)算法[13].

    5)遷移:在子種群遷移中,將子種群定義成網(wǎng)狀結(jié)構(gòu),并根據(jù)其適應(yīng)度值進(jìn)行各子種群間的個(gè)體遷移[14].

    算法中各控制參數(shù)設(shè)置及優(yōu)化流程分別如表1和圖2所示.

    表1 遺傳算法初始參數(shù)設(shè)置Tab.1 Control parameters definition of genetic algorithm

    圖2 遺傳算法流程圖Fig.2 Flow chart of genetic algorithm

    3 算法應(yīng)用及分析

    3.1 沙箱熱響應(yīng)測(cè)試

    Beier等[15]在實(shí)驗(yàn)室中搭建了一個(gè)較大的沙箱(sandbox)進(jìn)行熱響應(yīng)測(cè)試實(shí)驗(yàn).沙箱周圍采用木板結(jié)構(gòu)固定,里面填充潮濕沙土,其長(zhǎng)和寬分別為18.32 m和1.83 m.中心處設(shè)置內(nèi)徑為12.6 cm的鋁管構(gòu)建鉆孔壁,其內(nèi)水平安裝單U型地埋管換熱器,并用質(zhì)量分?jǐn)?shù)為20%的泥漿膨潤(rùn)土進(jìn)行回填.

    熱響應(yīng)測(cè)試實(shí)驗(yàn)中利用2臺(tái)電加熱器對(duì)流體進(jìn)行加熱,加熱器輸入功率不確定性為±1%;熱敏電阻測(cè)量U型管進(jìn)出口流體溫度,其測(cè)量不確定性為±0.03 ℃;流量計(jì)測(cè)量流體體積流量,其測(cè)量不確定性為±5%.總測(cè)試時(shí)長(zhǎng)為48 h,其主要測(cè)試參數(shù)見(jiàn)表2.

    表2 熱響應(yīng)測(cè)試主要參數(shù)[15]Tab.2 Major parameters of the TRT

    3.2 目標(biāo)函數(shù)進(jìn)化

    圖3給出了綜合均方根誤差(RSME)最小值及平均值隨遺傳代數(shù)的變化.可以看出,目標(biāo)函數(shù)(綜合均方根誤差)在40代以后變化較緩,此時(shí)計(jì)算收斂,最小值和平均值分別為0.111和0.118 ℃.

    遺傳代數(shù)圖3 綜合均方根誤差進(jìn)化曲線Fig.3 The evolution curve of composite RMSE

    當(dāng)目標(biāo)函數(shù)最小時(shí),2個(gè)模型分別計(jì)算所得的均方根誤差遺傳代數(shù)的變化如圖4所示.當(dāng)算法逐漸收斂之后,復(fù)合介質(zhì)線熱源模型計(jì)算所得均方根誤差比無(wú)限長(zhǎng)線熱源模型大,收斂之后其差異平均值為0.124 ℃.這是因?yàn)椴捎脧?fù)合介質(zhì)模型計(jì)算流體平均溫度時(shí),未舍去熱響應(yīng)測(cè)試前10 h的溫度數(shù)據(jù),其計(jì)算時(shí)長(zhǎng)即為熱響應(yīng)測(cè)試總時(shí)長(zhǎng).由于傳熱前期鉆孔內(nèi)非穩(wěn)態(tài)傳熱對(duì)流體溫度的影響較大,根據(jù)復(fù)合介質(zhì)模型計(jì)算的流體平均溫度和實(shí)際值的溫度差也相應(yīng)較大.

    遺傳代數(shù)圖4 流體平均溫度均方根誤差進(jìn)化曲線Fig.4 The RMSE evolution curve for fluid average temperature

    3.3 參數(shù)評(píng)估結(jié)果及分析

    基于地埋管無(wú)限長(zhǎng)線源模型和管壁熱阻修正的復(fù)合介質(zhì)線熱源模型,分別采用模式搜索法(Hooke-Jeeves算法)、單純形法和遺傳算法進(jìn)行巖土及填料熱物性參數(shù)評(píng)估,結(jié)果見(jiàn)表3.表中鉆孔熱阻的實(shí)測(cè)值是通過(guò)單獨(dú)測(cè)量鉆孔壁平均溫度計(jì)算所得.模式搜索法由初始值開始交替實(shí)施軸向搜索和模式搜索,直至找到最優(yōu)值.單純形算法則是基于線性規(guī)劃問(wèn)題的可行域多面凸集原理,通過(guò)查找該凸集的某頂點(diǎn)得出最優(yōu)值.可以發(fā)現(xiàn),遺傳算法計(jì)算所得的綜合均方根誤差最小,為0.111 ℃,且鉆孔熱阻和導(dǎo)熱系數(shù)的評(píng)估相對(duì)誤差均小于5%.

    表3 熱物性參數(shù)評(píng)估結(jié)果及對(duì)比Tab.3 The evaluation results and comparison of thermal parameters

    圖5所示為綜合均方根誤差最小時(shí)所對(duì)應(yīng)的流體平均溫度隨時(shí)間的變化曲線.由圖可知,模型計(jì)算流體平均溫度與實(shí)測(cè)值均吻合較好,所得流體平均溫度與實(shí)測(cè)值的均方根誤差分別為0.0505和0.172 ℃.同時(shí)可見(jiàn),復(fù)合介質(zhì)模型計(jì)算流體平均溫度與實(shí)測(cè)值的差異主要集中在前10 h以內(nèi).

    時(shí)間/h圖5 流體平均溫度對(duì)比Fig.5 Comparison of fluid average temperature

    4 結(jié) 論

    本文采用管壁熱阻修正的復(fù)合介質(zhì)線熱源模型與無(wú)限長(zhǎng)線熱源模型并聯(lián)求解流體平均溫度,并以流體平均溫度的綜合均方根誤差為目標(biāo)函數(shù),采用自適應(yīng)變異型多種群遺傳算法尋優(yōu)求解了鉆孔熱阻、巖土導(dǎo)熱系數(shù)及熱擴(kuò)散系數(shù)、填料導(dǎo)熱系數(shù)及熱擴(kuò)散系數(shù)等5個(gè)熱物性參數(shù).與傳統(tǒng)參數(shù)估計(jì)法相比,通過(guò)管壁熱阻修正復(fù)合介質(zhì)線熱源模型可較準(zhǔn)確地計(jì)算流體平均溫度,同時(shí)可評(píng)估填料熱物性參數(shù).在沙箱熱響應(yīng)測(cè)試中應(yīng)用遺傳算法進(jìn)行巖土及填料熱物性參數(shù)評(píng)估,通過(guò)計(jì)算可以得出以下結(jié)論:

    1)尋優(yōu)過(guò)程中,遺傳進(jìn)化到40代左右時(shí)計(jì)算收斂,所得最終均方根誤差的最小值為0.111 ℃.

    2)與其他尋優(yōu)算法相比,應(yīng)用遺傳算法評(píng)估巖土及填料的熱物性參數(shù),精度較高,鉆孔熱阻、巖土導(dǎo)熱系數(shù)、填料導(dǎo)熱系數(shù)的相對(duì)誤差分別為3.47%,1.42%和4.93%.

    3)通過(guò)評(píng)估所得巖土及填料熱物性參數(shù)計(jì)算流體平均溫度,與實(shí)測(cè)值相比,根據(jù)線熱源模型和復(fù)合模型計(jì)算所得流體平均溫度的均方根誤差分別為0.050 5 ℃和0.172 ℃.

    [1] SPITLER J D,GEHLIN S E A. Thermal response testing for ground source geat pump systems-an historical review[J]. Renewable and Sustainable Energy Reviews,2015, 50:1125-1137.

    [2] ASHRAE. ASHRAE handbook: HVAC applications [M]. Atlanta,GA: American Society of Heating,Refrigerating and Air-Conditioning Engineers,Inc, 2011:34.13-34.17.

    [3] ZHANG C,GUO Z, LIU Y,etal. A review on thermal response test of ground-coupled heat pump systems[J]. Renewable and Sustainable Energy Reviews,2014,40(8):851-867.

    [4] KELVIN T W. Mathematical and physical papers II[M]. London:Cambridge University Press,1884:41-60.

    [5] 刁乃仁,方肇洪.地埋管地源熱泵技術(shù)[M].北京:高等教育出版社,2006:47-50.

    DIAO Nairen,F(xiàn)ANG Zhaohong.Ground-coupled heat pump technology [M]. Beijing:Higher Education Press,2006:47-50.(In Chinese)

    [6] ZHANG L F,ZHANG Q,HUANG G S,etal. AP(t)-linear average method to estimate the thermal parameters of the borehole heat exchangers for in situ thermal response test[J]. Applied Energy,2014,131(7):211-221.

    [7] LI M,LI P,CHAN V,etal. Full-scale temperature response function (G-function) for heat transfer by borehole ground heat exchangers (GHEs) from sub-hour to decades [J]. Applied Energy,2014,136:197-205.

    [8] LI M, LAI A C K. Analytical model for short-time responses of ground heat exchangers with U-shaped tubes: Model development and validation [J]. Applied Energy,2013,104:510-516.

    [9] ZHANG L F,ZHANG Q,HUANG G S. A transient quasi-3D entire time scale line source model for the fluid and ground temperature prediction of vertical ground heat exchangers[J]. Applied Energy,2016,170:65-75.

    [10]張泉,杜亞星,張林峰,等. 基于遺傳算法的單U地源熱泵鉆孔內(nèi)熱阻研究[J]. 湖南大學(xué)學(xué)報(bào):自然科學(xué)版,2014,41(9):100-105.

    ZHANG Quan,DU Yaxing,ZHANG Linfeng,etal. A study on the heat resistance of single U pipe in the ground source heat pump system based on genetic algorithm[J].Journal of Hunan University: Natural Sciences, 2014,41(9):100-105. (In Chinese)

    [11]ANDREW C,PETER F,HARTMUT P,etal. Genetic algorithm TOOLBOX for use with MATLAB,user’s guide v1.2[R]. Sheffield:University of Sheffield,2003:9-36.

    [12]FUJII H,OKUBO H,NISHI K,etal. An improved thermal response test for U-tube ground heat exchanger based on optical fiber thermometers [J]. Geothermics,2009,38(4):399-406.

    [13]SRINIVAS M,PATNALK L M. Adaptive probabilities of crossover and mutation in genetic algorithms[J]. IEEE Transactions on Systems,Man and Cybernetics,1994,24(4):656-667.

    [14]雷英杰.MATLAB遺傳算法工具箱及應(yīng)用[M].西安:西安電子科技大學(xué)出版社,2005:68-69.

    LEI Yingjie. Genetic algorithm toolbox and application[M].Xi’an:Xidian University Press,2005:68-69. (In Chinese)

    [15]BEIER R A,SMITH M D,SPITLER J D. Reference data sets for vertical borehole ground heat exchanger models and thermal response test analysis [J]. Geothermics,2011,40(1):79-85.

    Properties Based on Genetic Algorithm

    ZHANG Quan?,MA Xiaowei,ZHANG Linfeng

    (College of Civil Engineering,Hunan University,Changsha 410082,China)

    The parameters for thermal properties of ground and grout have directly effect on the performance of ground heat exchanger (GHE). For the problem that grout thermal properties of in-situ thermal response test (TRT) need to be measured in advance and cannot be evaluated directly by calculation, this paper applied genetic algorithm to evaluate the parameters of ground and grout thermal properties simultaneously on the basis of infinite line source (ILS) model and composite media line source (CMLS) model with wall thermal resistance correction. The method was validated by in situ TRT. Compared with the measured values, the calculation relative error of borehole thermal resistance, ground thermal conductivity and diffusivity was 3.47%, 1.42% and 4.93%, respectively. The root mean square error (RMSE) of the calculated average temperature by these two models and the measured value is 0.0505 and 0.172 ℃, respectively. The results provide important reference for GHE design.

    grout thermal properties;genetic algorithms;composite media line source;thermal response test

    2016-03-29

    廣東省教育部產(chǎn)學(xué)研項(xiàng)目(2010B090400301),Project on the Integration of Industry,Education and Research of Guangdong Province and Ministry of Education(2010B090400301);國(guó)家國(guó)際科技合作專項(xiàng)項(xiàng)目(2015DFA61170),International Science and Technology Cooperation Program of China(2015DFA61170)

    張泉(1970-),男,江蘇淮安人,湖南大學(xué)教授,博士生導(dǎo)師 ?通訊聯(lián)系人,E-mail:zhangquan_2002@sina.com

    1674-2974(2017)03-0151-06

    10.16339/j.cnki.hdxbzkb.2017.03.019

    TU831.6

    AParameter Estimation of Soil and Grout Thermal

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務(wù)本地化模型
    適用于BDS-3 PPP的隨機(jī)模型
    提煉模型 突破難點(diǎn)
    函數(shù)模型及應(yīng)用
    p150Glued在帕金森病模型中的表達(dá)及分布
    函數(shù)模型及應(yīng)用
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    3D打印中的模型分割與打包
    婷婷丁香在线五月| 狠狠婷婷综合久久久久久88av| 亚洲人成电影观看| 免费久久久久久久精品成人欧美视频| 天天躁日日躁夜夜躁夜夜| 国产一区二区三区视频了| 久久精品aⅴ一区二区三区四区| 叶爱在线成人免费视频播放| 国产精品美女特级片免费视频播放器 | 99热国产这里只有精品6| 久久午夜综合久久蜜桃| 亚洲精品国产区一区二| 亚洲熟妇中文字幕五十中出 | 精品午夜福利视频在线观看一区| 无限看片的www在线观看| 精品国产超薄肉色丝袜足j| 热99久久久久精品小说推荐| 国产成人影院久久av| 久久久国产一区二区| 看黄色毛片网站| 亚洲国产欧美日韩在线播放| 成年人免费黄色播放视频| 一区二区三区精品91| 亚洲avbb在线观看| 国产精品一区二区精品视频观看| 久久天堂一区二区三区四区| 国产精品偷伦视频观看了| 午夜成年电影在线免费观看| 黑人巨大精品欧美一区二区mp4| av欧美777| 亚洲国产精品一区二区三区在线| 婷婷丁香在线五月| 午夜激情av网站| 啦啦啦免费观看视频1| 欧美日韩乱码在线| 老司机福利观看| 高清欧美精品videossex| 午夜免费成人在线视频| 免费女性裸体啪啪无遮挡网站| 免费不卡黄色视频| 又黄又粗又硬又大视频| 国产xxxxx性猛交| 久久人妻熟女aⅴ| 国产成人免费观看mmmm| 国产免费男女视频| 搡老岳熟女国产| 久久人妻福利社区极品人妻图片| 亚洲第一欧美日韩一区二区三区| 久久久久国产精品人妻aⅴ院 | 中文字幕最新亚洲高清| 麻豆成人av在线观看| 一区福利在线观看| 又大又爽又粗| 国产精品久久电影中文字幕 | 91av网站免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 又黄又粗又硬又大视频| 国产精品香港三级国产av潘金莲| 视频在线观看一区二区三区| 亚洲情色 制服丝袜| 亚洲欧美激情在线| 宅男免费午夜| 后天国语完整版免费观看| 久久久久久久久免费视频了| 久久婷婷成人综合色麻豆| 免费在线观看完整版高清| 欧美亚洲 丝袜 人妻 在线| 欧美精品av麻豆av| 欧美成人午夜精品| 成年人免费黄色播放视频| 久99久视频精品免费| 啦啦啦在线免费观看视频4| 免费在线观看影片大全网站| 久久香蕉国产精品| 成人国语在线视频| 女性被躁到高潮视频| 国产不卡一卡二| 男女下面插进去视频免费观看| 校园春色视频在线观看| 侵犯人妻中文字幕一二三四区| 18禁国产床啪视频网站| 成人黄色视频免费在线看| 国产99久久九九免费精品| 国产精品久久久久久精品古装| 91精品三级在线观看| 操出白浆在线播放| 亚洲一区二区三区不卡视频| 一区二区日韩欧美中文字幕| 国产精品成人在线| 色播在线永久视频| 国产真人三级小视频在线观看| 男女下面插进去视频免费观看| 欧洲精品卡2卡3卡4卡5卡区| 亚洲成人免费电影在线观看| 一个人免费在线观看的高清视频| 人妻久久中文字幕网| 激情在线观看视频在线高清 | 日日爽夜夜爽网站| 大型黄色视频在线免费观看| 俄罗斯特黄特色一大片| 国产精品永久免费网站| 黄色 视频免费看| 国产免费现黄频在线看| 婷婷丁香在线五月| 国产无遮挡羞羞视频在线观看| 一二三四社区在线视频社区8| 久9热在线精品视频| 中文字幕av电影在线播放| 午夜免费成人在线视频| av超薄肉色丝袜交足视频| 欧美大码av| 精品一区二区三区视频在线观看免费 | 国产xxxxx性猛交| 欧美精品人与动牲交sv欧美| 一级毛片女人18水好多| av中文乱码字幕在线| 亚洲av第一区精品v没综合| 熟女少妇亚洲综合色aaa.| 亚洲欧美精品综合一区二区三区| 男女床上黄色一级片免费看| 女人被狂操c到高潮| 老司机在亚洲福利影院| 十八禁网站免费在线| 国产欧美日韩综合在线一区二区| 中文字幕精品免费在线观看视频| 亚洲成人国产一区在线观看| 免费在线观看亚洲国产| 国产精品一区二区精品视频观看| 国产一区二区激情短视频| 99热网站在线观看| 精品福利观看| 一二三四社区在线视频社区8| 国产精品秋霞免费鲁丝片| 久久精品亚洲熟妇少妇任你| 99精品欧美一区二区三区四区| 别揉我奶头~嗯~啊~动态视频| 王馨瑶露胸无遮挡在线观看| 国产又爽黄色视频| 欧美精品一区二区免费开放| 欧美乱码精品一区二区三区| 亚洲成人国产一区在线观看| 在线观看免费日韩欧美大片| 99国产精品99久久久久| 日韩欧美一区视频在线观看| 别揉我奶头~嗯~啊~动态视频| 欧美乱码精品一区二区三区| 国产av一区二区精品久久| tocl精华| 在线观看舔阴道视频| 国产又爽黄色视频| 亚洲黑人精品在线| 亚洲 欧美一区二区三区| 新久久久久国产一级毛片| 国产精品久久久久久人妻精品电影| 国产高清videossex| 国产成人欧美| 国产精品自产拍在线观看55亚洲 | a级片在线免费高清观看视频| 别揉我奶头~嗯~啊~动态视频| 美国免费a级毛片| 欧美黑人精品巨大| 香蕉国产在线看| 久久精品国产综合久久久| 在线观看免费视频网站a站| 国产色视频综合| 99国产精品99久久久久| 亚洲视频免费观看视频| 操美女的视频在线观看| 一区二区三区激情视频| 人人妻人人澡人人爽人人夜夜| 自拍欧美九色日韩亚洲蝌蚪91| 天堂俺去俺来也www色官网| 女性生殖器流出的白浆| 99久久综合精品五月天人人| 不卡一级毛片| 国内毛片毛片毛片毛片毛片| 精品第一国产精品| 在线视频色国产色| 在线观看日韩欧美| 下体分泌物呈黄色| 欧美日韩中文字幕国产精品一区二区三区 | 国产不卡一卡二| 久久午夜亚洲精品久久| 欧美性长视频在线观看| 欧美激情极品国产一区二区三区| 国产男女内射视频| 国产精品综合久久久久久久免费 | 在线观看免费视频网站a站| 男人操女人黄网站| 手机成人av网站| 十八禁网站免费在线| 国产成人欧美在线观看 | 国产麻豆69| 亚洲午夜精品一区,二区,三区| 91精品三级在线观看| 国产单亲对白刺激| 久久久久久久国产电影| 久久久久国内视频| 老司机深夜福利视频在线观看| 一本一本久久a久久精品综合妖精| 两人在一起打扑克的视频| 欧美精品亚洲一区二区| 大陆偷拍与自拍| 男女高潮啪啪啪动态图| videos熟女内射| 一边摸一边抽搐一进一小说 | 最近最新中文字幕大全免费视频| 亚洲色图综合在线观看| 国产av一区二区精品久久| 午夜老司机福利片| 亚洲色图av天堂| 欧美黑人精品巨大| 黄色视频,在线免费观看| 国产精品久久久人人做人人爽| 99国产综合亚洲精品| 国产激情久久老熟女| 黑人巨大精品欧美一区二区mp4| 国产一卡二卡三卡精品| 99riav亚洲国产免费| 脱女人内裤的视频| 十八禁高潮呻吟视频| 最近最新免费中文字幕在线| 99精品在免费线老司机午夜| 婷婷精品国产亚洲av在线 | 法律面前人人平等表现在哪些方面| 黑人操中国人逼视频| 自线自在国产av| 亚洲性夜色夜夜综合| 99久久国产精品久久久| 淫妇啪啪啪对白视频| 国产野战对白在线观看| 国产精品偷伦视频观看了| 亚洲专区字幕在线| 成人免费观看视频高清| 国产aⅴ精品一区二区三区波| 无遮挡黄片免费观看| 在线观看午夜福利视频| 欧美大码av| 无限看片的www在线观看| 国产精品一区二区在线不卡| 久久久久国产精品人妻aⅴ院 | 18禁裸乳无遮挡免费网站照片 | 老司机午夜福利在线观看视频| 免费在线观看日本一区| 亚洲欧美精品综合一区二区三区| 无遮挡黄片免费观看| 老司机影院毛片| 亚洲一卡2卡3卡4卡5卡精品中文| 久久人人爽av亚洲精品天堂| 黑人巨大精品欧美一区二区mp4| 国产免费av片在线观看野外av| 欧美在线一区亚洲| 国产亚洲欧美98| 亚洲欧美一区二区三区久久| 日韩制服丝袜自拍偷拍| 黑丝袜美女国产一区| 夜夜夜夜夜久久久久| 757午夜福利合集在线观看| 成年人免费黄色播放视频| 丰满人妻熟妇乱又伦精品不卡| 涩涩av久久男人的天堂| 人人妻人人澡人人看| 亚洲欧美日韩另类电影网站| 亚洲av熟女| 99久久99久久久精品蜜桃| 人人妻人人爽人人添夜夜欢视频| 男女高潮啪啪啪动态图| av欧美777| 国产成人免费观看mmmm| 亚洲伊人色综图| 制服诱惑二区| 麻豆乱淫一区二区| 老鸭窝网址在线观看| 9191精品国产免费久久| 天堂√8在线中文| 欧美中文综合在线视频| 一级a爱片免费观看的视频| 男人操女人黄网站| 国产99白浆流出| 欧美在线一区亚洲| 成年人午夜在线观看视频| videosex国产| 99国产精品99久久久久| 国产精品影院久久| av福利片在线| 久久久久久人人人人人| 一区在线观看完整版| 九色亚洲精品在线播放| 免费女性裸体啪啪无遮挡网站| 99国产精品99久久久久| 国产亚洲欧美在线一区二区| 国产精品电影一区二区三区 | 国产成人免费无遮挡视频| 久久精品国产清高在天天线| 欧洲精品卡2卡3卡4卡5卡区| 国产精品久久久av美女十八| 国产成人免费观看mmmm| 国产日韩欧美亚洲二区| 一级作爱视频免费观看| 欧美 日韩 精品 国产| 久久精品成人免费网站| 侵犯人妻中文字幕一二三四区| 久久精品国产亚洲av高清一级| 午夜91福利影院| 亚洲精品美女久久av网站| 男女高潮啪啪啪动态图| 少妇被粗大的猛进出69影院| 国产欧美日韩综合在线一区二区| av不卡在线播放| 国产主播在线观看一区二区| 久久天躁狠狠躁夜夜2o2o| 怎么达到女性高潮| 亚洲国产欧美日韩在线播放| 精品人妻在线不人妻| 在线观看日韩欧美| 午夜精品国产一区二区电影| 亚洲人成电影免费在线| 丝袜在线中文字幕| 久久九九热精品免费| 九色亚洲精品在线播放| 午夜福利在线观看吧| 亚洲一区中文字幕在线| 欧美在线一区亚洲| 777久久人妻少妇嫩草av网站| 亚洲欧美一区二区三区黑人| 女人爽到高潮嗷嗷叫在线视频| www.999成人在线观看| 制服人妻中文乱码| 亚洲视频免费观看视频| 国产日韩一区二区三区精品不卡| 久久精品熟女亚洲av麻豆精品| 在线av久久热| 久久久国产欧美日韩av| 欧美黄色淫秽网站| 亚洲美女黄片视频| 中文字幕色久视频| 日本五十路高清| 女警被强在线播放| 乱人伦中国视频| 一级作爱视频免费观看| 天天躁夜夜躁狠狠躁躁| 日本wwww免费看| 女人被躁到高潮嗷嗷叫费观| 日本一区二区免费在线视频| 黑人巨大精品欧美一区二区蜜桃| aaaaa片日本免费| 国产深夜福利视频在线观看| x7x7x7水蜜桃| 亚洲成人免费av在线播放| 日本vs欧美在线观看视频| www.自偷自拍.com| x7x7x7水蜜桃| 亚洲美女黄片视频| 99久久综合精品五月天人人| 国产精品香港三级国产av潘金莲| 多毛熟女@视频| 午夜福利乱码中文字幕| 国产片内射在线| 精品人妻在线不人妻| 欧美中文综合在线视频| 国产不卡av网站在线观看| 国产片内射在线| 国产免费av片在线观看野外av| 法律面前人人平等表现在哪些方面| 日韩精品免费视频一区二区三区| 男女高潮啪啪啪动态图| 午夜福利,免费看| 黑人欧美特级aaaaaa片| 亚洲成av片中文字幕在线观看| 精品久久久久久久久久免费视频 | www.999成人在线观看| 91老司机精品| 99在线人妻在线中文字幕 | 日本a在线网址| 亚洲专区字幕在线| 久久久久久久精品吃奶| 国产99久久九九免费精品| 久久狼人影院| 免费在线观看影片大全网站| 在线看a的网站| 男女高潮啪啪啪动态图| 午夜福利视频在线观看免费| 久久午夜亚洲精品久久| 在线十欧美十亚洲十日本专区| 久久国产精品男人的天堂亚洲| 亚洲人成77777在线视频| 熟女少妇亚洲综合色aaa.| 色在线成人网| 在线观看免费日韩欧美大片| 黄色a级毛片大全视频| 国产在线观看jvid| 亚洲色图 男人天堂 中文字幕| 午夜激情av网站| 香蕉丝袜av| 日日夜夜操网爽| 99精品欧美一区二区三区四区| 大片电影免费在线观看免费| 国产在线一区二区三区精| 成人手机av| 九色亚洲精品在线播放| 黄色视频不卡| 高清av免费在线| 纯流量卡能插随身wifi吗| 亚洲一区二区三区欧美精品| 91九色精品人成在线观看| 女性生殖器流出的白浆| 成人永久免费在线观看视频| 在线看a的网站| 777久久人妻少妇嫩草av网站| 免费看十八禁软件| 老司机靠b影院| 中文字幕人妻熟女乱码| videos熟女内射| 自拍欧美九色日韩亚洲蝌蚪91| 国产高清激情床上av| 国产精品永久免费网站| 99国产精品免费福利视频| 高潮久久久久久久久久久不卡| 欧美激情极品国产一区二区三区| 好看av亚洲va欧美ⅴa在| 香蕉丝袜av| 亚洲国产精品合色在线| a在线观看视频网站| 欧美亚洲日本最大视频资源| 黄色片一级片一级黄色片| 国产在线一区二区三区精| 18禁裸乳无遮挡动漫免费视频| 国产又色又爽无遮挡免费看| 91大片在线观看| 日本黄色日本黄色录像| 久久精品国产综合久久久| 精品免费久久久久久久清纯 | 精品国产一区二区三区久久久樱花| 两性夫妻黄色片| 亚洲色图av天堂| 欧美亚洲 丝袜 人妻 在线| 欧美亚洲日本最大视频资源| 国产国语露脸激情在线看| 91成人精品电影| 欧美另类亚洲清纯唯美| 人人妻人人爽人人添夜夜欢视频| 亚洲成人手机| 亚洲精品中文字幕一二三四区| 老熟妇乱子伦视频在线观看| 欧美日韩成人在线一区二区| 纯流量卡能插随身wifi吗| 99热网站在线观看| 两个人免费观看高清视频| 久久人妻av系列| 亚洲av成人av| 好看av亚洲va欧美ⅴa在| 亚洲黑人精品在线| 婷婷成人精品国产| 99久久精品国产亚洲精品| 久久狼人影院| 人人妻人人添人人爽欧美一区卜| 精品久久久久久久毛片微露脸| 一本大道久久a久久精品| 一级片免费观看大全| 午夜免费鲁丝| 久久中文字幕人妻熟女| 亚洲专区国产一区二区| 久99久视频精品免费| 久久久国产欧美日韩av| 亚洲精品国产色婷婷电影| 99国产精品免费福利视频| 9191精品国产免费久久| 老鸭窝网址在线观看| 国产精品久久久久久精品古装| x7x7x7水蜜桃| 久久精品aⅴ一区二区三区四区| 精品久久久久久久久久免费视频 | cao死你这个sao货| 老汉色∧v一级毛片| 在线观看日韩欧美| 免费少妇av软件| 91大片在线观看| 日韩熟女老妇一区二区性免费视频| 一级毛片女人18水好多| 国产欧美亚洲国产| 日日爽夜夜爽网站| 亚洲人成77777在线视频| 国产成人精品久久二区二区免费| 成人免费观看视频高清| 丁香六月欧美| 91成年电影在线观看| 亚洲午夜理论影院| 黄色视频,在线免费观看| 一本综合久久免费| 中文字幕人妻丝袜制服| 香蕉久久夜色| 欧美日韩乱码在线| 日韩免费高清中文字幕av| 欧美日韩亚洲综合一区二区三区_| 在线观看免费高清a一片| 国产男女超爽视频在线观看| 亚洲中文日韩欧美视频| 国产精品.久久久| 精品国产乱子伦一区二区三区| 另类亚洲欧美激情| 亚洲欧洲精品一区二区精品久久久| 日本欧美视频一区| 操美女的视频在线观看| 韩国精品一区二区三区| 精品国产一区二区久久| 纯流量卡能插随身wifi吗| 一级作爱视频免费观看| 久久久久久人人人人人| 亚洲人成电影免费在线| 纯流量卡能插随身wifi吗| 亚洲人成电影免费在线| 三上悠亚av全集在线观看| 国产精品美女特级片免费视频播放器 | 国产无遮挡羞羞视频在线观看| 香蕉国产在线看| 国产精品久久久久久精品古装| 国产区一区二久久| 国产91精品成人一区二区三区| 王馨瑶露胸无遮挡在线观看| 久久精品成人免费网站| 天堂√8在线中文| 窝窝影院91人妻| 国产成人精品久久二区二区91| 激情视频va一区二区三区| 精品久久久久久电影网| 丰满的人妻完整版| 黑人操中国人逼视频| 免费看十八禁软件| 精品人妻1区二区| 中文字幕最新亚洲高清| 国产精品av久久久久免费| 国产精品免费视频内射| 亚洲欧美激情综合另类| 国产精品国产av在线观看| 亚洲av片天天在线观看| 99riav亚洲国产免费| 新久久久久国产一级毛片| www.熟女人妻精品国产| 黑丝袜美女国产一区| 麻豆乱淫一区二区| 女人被狂操c到高潮| 久久中文看片网| 久久精品国产综合久久久| 日韩大码丰满熟妇| 亚洲美女黄片视频| 夜夜躁狠狠躁天天躁| 男女之事视频高清在线观看| 久久久久久免费高清国产稀缺| 久久久久视频综合| 日本黄色日本黄色录像| 下体分泌物呈黄色| 日本wwww免费看| 在线永久观看黄色视频| 亚洲中文日韩欧美视频| 18禁裸乳无遮挡免费网站照片 | 国产欧美日韩综合在线一区二区| 久久亚洲精品不卡| 99久久99久久久精品蜜桃| 国产精品久久视频播放| videos熟女内射| 精品一区二区三区av网在线观看| 最近最新中文字幕大全免费视频| 久久久久久久国产电影| 国产片内射在线| 国产有黄有色有爽视频| 丝袜在线中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 午夜福利在线免费观看网站| 亚洲成人免费电影在线观看| 老司机午夜福利在线观看视频| 亚洲熟妇熟女久久| 亚洲欧美精品综合一区二区三区| 性少妇av在线| 亚洲久久久国产精品| 不卡av一区二区三区| 午夜精品久久久久久毛片777| 亚洲精品久久成人aⅴ小说| 亚洲国产欧美网| 久久性视频一级片| 国产又色又爽无遮挡免费看| 亚洲成人手机| 久久性视频一级片| 亚洲国产看品久久| 少妇猛男粗大的猛烈进出视频| 18在线观看网站| 亚洲欧美色中文字幕在线| 免费在线观看日本一区| 亚洲五月天丁香| 亚洲美女黄片视频| 99精品在免费线老司机午夜| 色精品久久人妻99蜜桃| 久久久国产欧美日韩av| 啦啦啦在线免费观看视频4| 一级片免费观看大全| 国产野战对白在线观看| 亚洲 欧美一区二区三区| 亚洲精品自拍成人| 亚洲国产精品sss在线观看 | 久久精品亚洲熟妇少妇任你| bbb黄色大片| 国产不卡av网站在线观看| 亚洲va日本ⅴa欧美va伊人久久| 日韩有码中文字幕| a级片在线免费高清观看视频| www日本在线高清视频| 免费在线观看黄色视频的| 国产一卡二卡三卡精品| 欧美日韩瑟瑟在线播放| 亚洲人成伊人成综合网2020| a级片在线免费高清观看视频| 岛国在线观看网站| 色婷婷av一区二区三区视频| 一二三四社区在线视频社区8| 黄色成人免费大全| 成年女人毛片免费观看观看9 |