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

    基于巖土失調(diào)溫度限值的土壤源熱泵系統(tǒng)土壤蓄能狀態(tài)評價(jià)

    2015-03-08 06:34:48尹暢昱金逸韜
    關(guān)鍵詞:現(xiàn)值源熱泵能效

    王 勇,尹暢昱,金逸韜

    (1.重慶大學(xué) 三峽庫區(qū)生態(tài)環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室,重慶 400045;2.重慶大學(xué) 低碳綠色建筑國際聯(lián)合研究中心,重慶 400045;3.蘇州工業(yè)園區(qū)設(shè)計(jì)研究院股份有限公司,江蘇 蘇州 215021)

    基于巖土失調(diào)溫度限值的土壤源熱泵系統(tǒng)土壤蓄能狀態(tài)評價(jià)

    王 勇1,2?,尹暢昱1,2,金逸韜3

    (1.重慶大學(xué) 三峽庫區(qū)生態(tài)環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室,重慶 400045;2.重慶大學(xué) 低碳綠色建筑國際聯(lián)合研究中心,重慶 400045;3.蘇州工業(yè)園區(qū)設(shè)計(jì)研究院股份有限公司,江蘇 蘇州 215021)

    基于巖土失調(diào)溫度限值、建筑負(fù)荷特征、傳熱機(jī)理等多參數(shù)耦合的分析方法,提出了系統(tǒng)節(jié)能率、系統(tǒng)總運(yùn)行費(fèi)用現(xiàn)值和系統(tǒng)失效指標(biāo)等評價(jià)指標(biāo),并建立了土壤蓄能狀態(tài)評價(jià)體系.以評價(jià)體系為基礎(chǔ),通過CFD軟件建立三維地下埋管管群數(shù)值計(jì)算模型,構(gòu)建了評價(jià)全壽命周期內(nèi)的土壤源熱泵系統(tǒng)土壤蓄能狀態(tài)的計(jì)算方法.計(jì)算和分析結(jié)果表明,評價(jià)體系和計(jì)算方法不僅可以較客觀地反映土壤的蓄能狀態(tài),同時(shí)能夠預(yù)測基于實(shí)際工程運(yùn)行特性下的巖土失調(diào)溫度限值范圍和全壽命周期內(nèi)節(jié)能性和經(jīng)濟(jì)性的失效時(shí)間點(diǎn).

    土壤源熱泵;全壽命周期;溫度限值;CFD模擬;土壤蓄能狀態(tài)

    目前中國的土壤源熱泵系統(tǒng)發(fā)展越來越快,大規(guī)模的地下埋管換熱器對系統(tǒng)的長期高效運(yùn)行提出了嚴(yán)峻的考驗(yàn),土壤蓄能狀態(tài)問題已經(jīng)成為制約其科學(xué)發(fā)展的關(guān)鍵因素.花莉等[1]利用TRNSYS軟件進(jìn)行熱平衡問題影響分析,研究了冷熱負(fù)荷比較大時(shí)系統(tǒng)運(yùn)行50年間COP的變化,并進(jìn)行單因素分析得出緩解由于冷熱不平衡而引起的土壤溫升的措施,但模擬中未將建筑負(fù)荷比與埋管工況等進(jìn)行多參數(shù)耦合分析,所以其結(jié)論對具體工程指導(dǎo)性不強(qiáng).范蕊等[2]將夏季土壤吸熱負(fù)荷100 kW,冬季土壤放熱負(fù)荷50 kW定義為土壤全年熱不平衡率為50%;冬夏季負(fù)荷相等定義為土壤全年熱不平衡率為0%,其定義沒有考慮土壤的自平衡能力和建筑負(fù)荷與埋管取(釋)熱負(fù)荷的耦合關(guān)系等因素.林東超等[3-4]提出土壤熱失衡問題應(yīng)從冬夏空調(diào)負(fù)荷情況、地埋管換熱器的間距、地埋管換熱器系統(tǒng)構(gòu)成和實(shí)際運(yùn)行情況等各方面進(jìn)行考慮,擬定系統(tǒng)不同負(fù)荷下的運(yùn)行策略,但未能建立以實(shí)際工程為基礎(chǔ)的三維管群模型,也沒有建立建筑負(fù)荷與機(jī)組出力負(fù)荷之間的耦合關(guān)系式.楊衛(wèi)波等[5]建立了部分埋管群的模型,計(jì)算了十年間的土壤溫度變化,得出冷負(fù)荷熱比較大時(shí),土壤溫升較為明顯,且溫升與負(fù)荷比呈正相關(guān),但因?yàn)槟P腿源嬖谝欢ǖ募俣l件,也未以全壽命周期為基礎(chǔ)進(jìn)行分析,所以未能給出土壤溫度對熱泵系統(tǒng)失效的限值.

    由于目前土壤蓄能狀態(tài)評價(jià)模型的缺失,定義的不明確,地下三維管群模型的失真度較高,造成熱平衡理論研究的突破點(diǎn)較少.所以越來越多的研究偏向于從實(shí)際工程出發(fā),基于測試數(shù)據(jù),對土壤蓄能狀態(tài)進(jìn)行分析評判.楊紅輝等[6]實(shí)地測試北京某工程的單U管和雙U管的換熱情況,得出雖然冷熱負(fù)荷差存在一定不匹配,但是大地可以平衡一部分冷熱負(fù)荷差,但是其測試未能取得管群的影響效果,也未給出大地自我調(diào)節(jié)的極限溫度值.范龍華等[7-9]已經(jīng)意識到土壤蓄能失衡的主要原因是冬夏季負(fù)荷不等,但是埋管深度、管間距和建筑負(fù)荷等因素也存在影響,根據(jù)工程實(shí)際提出了分區(qū)運(yùn)行,間歇運(yùn)行,增加冰蓄冷系統(tǒng)等策略,但是仍處于定性分析的階段,未能給出土壤蓄能失調(diào)溫度限值和輔助冷熱源開啟時(shí)機(jī)等參數(shù)的定量結(jié)果.楊昌智等[10]創(chuàng)造性地提出單位井深換熱成本指數(shù),將運(yùn)行費(fèi)用與換熱能力結(jié)合,在一定程度上反映系統(tǒng)的經(jīng)濟(jì)性,但是未將末端和機(jī)組的運(yùn)行費(fèi)用考慮在內(nèi).

    目前研究中用到的土壤源熱泵土壤失效溫度基本取自GB 50366—2005《地源熱泵系統(tǒng)工程技術(shù)規(guī)范》,而規(guī)范中的溫度限值僅僅是保證機(jī)組正常運(yùn)行的低位冷卻水/輔熱水溫度范圍,當(dāng)作為巖土溫度限值時(shí)并不能保證熱泵系統(tǒng)的經(jīng)濟(jì)性和節(jié)能性.缺少對土壤蓄能狀態(tài)的評價(jià)和定義,缺乏多參數(shù)耦合分析,無法建立與實(shí)際工程高度仿真的埋管三維管群模型,這是目前對土壤源熱泵系統(tǒng)熱平衡問題分析的不足之處.對于土壤源熱泵的土壤蓄能狀態(tài)分析僅僅停留在冬夏季負(fù)荷或累計(jì)負(fù)荷比上進(jìn)行估算是不科學(xué)的,必須將建筑負(fù)荷特征,機(jī)組級數(shù)控制,傳熱機(jī)理,埋管工況等多參數(shù)耦合建模,考慮動(dòng)態(tài)能效,進(jìn)行全壽命周期內(nèi)土壤蓄能狀態(tài)的分析[11].

    在已有研究成果的基礎(chǔ)上,針對目前研究的不足,筆者建立了以系統(tǒng)節(jié)能率、系統(tǒng)總運(yùn)行費(fèi)用現(xiàn)值和系統(tǒng)失效指標(biāo)為基礎(chǔ)的土壤蓄能狀態(tài)的評價(jià)模型以及土壤蓄能狀態(tài);引入機(jī)組動(dòng)態(tài)能效,以建筑末端負(fù)荷和非穩(wěn)態(tài)三維管群模型為基礎(chǔ)[12],考慮整個(gè)壽命周期內(nèi)系統(tǒng)運(yùn)行工況,提出了一種分析評價(jià)土壤源熱泵系統(tǒng)的土壤蓄能狀態(tài)的一般計(jì)算方法.

    1 評價(jià)參數(shù)數(shù)學(xué)描述

    1.1 系統(tǒng)節(jié)能率

    為了具體說明土壤源熱泵系統(tǒng)相對常規(guī)空調(diào)系統(tǒng)的節(jié)能特點(diǎn),以常規(guī)空調(diào)系統(tǒng)運(yùn)行的能效比為對比基礎(chǔ),分析土壤源熱泵系統(tǒng)運(yùn)行時(shí)相對于常規(guī)空調(diào)系統(tǒng)的節(jié)能率[13].系統(tǒng)節(jié)能率定義為以使用冷水機(jī)組和燃?xì)忮仩t的常規(guī)空調(diào)系統(tǒng)能效為比較對象,EERs1(t) 為土壤源熱泵系統(tǒng)運(yùn)行的能效比,EERs2(t) 為常規(guī)空調(diào)系統(tǒng)運(yùn)行的能效比,兩者的差值與常規(guī)空調(diào)系統(tǒng)運(yùn)行的能效比相比,得到系統(tǒng)節(jié)能率f.系統(tǒng)節(jié)能率的數(shù)學(xué)描述表達(dá)為:

    (1)

    系統(tǒng)節(jié)能率f值若為正值,則代表t時(shí)刻土壤源熱泵系統(tǒng)的運(yùn)行能效比高于常規(guī)空調(diào)系統(tǒng),反之則相反.

    1.2 系統(tǒng)總運(yùn)行費(fèi)用現(xiàn)值

    系統(tǒng)節(jié)能率說明土壤源熱泵系統(tǒng)相對常規(guī)空調(diào)系統(tǒng)的節(jié)能特點(diǎn),但無法說明節(jié)能量的多少,故需要對兩個(gè)系統(tǒng)全壽命周期內(nèi)的系統(tǒng)總能耗(運(yùn)行費(fèi)用)進(jìn)行計(jì)算.

    1)土壤源熱泵系統(tǒng)的總能耗W的構(gòu)成:熱泵機(jī)組的能耗W1,地埋管側(cè)水泵的能耗W2,用戶側(cè)循環(huán)水泵的能耗W3,水處理儀的能耗W4及末端風(fēng)系統(tǒng)的能耗W5等.

    2)常規(guī)空調(diào)系統(tǒng)的總能耗N的構(gòu)成:冷水機(jī)組與燃?xì)忮仩t的能耗N1,冷卻水泵的能耗N2,用戶側(cè)循環(huán)水泵的能耗N3,水處理儀的能耗N4,末端風(fēng)系統(tǒng)的能耗N5及冷卻塔的能耗N6等.

    3)系統(tǒng)總運(yùn)行費(fèi)用現(xiàn)值的計(jì)算方法:數(shù)值計(jì)算的時(shí)間步長取為1 d,故系統(tǒng)總能耗的計(jì)算結(jié)果也應(yīng)以“d”為最小單位.例如,第n年第td的常規(guī)空調(diào)系統(tǒng)及土壤源熱泵系統(tǒng)的能耗可按式(2),(3)計(jì)算:

    N(t)=[N1(t)+N2(t)+N3(t)+N4(t)+

    N5(t)+N6(t)]×24;

    (2)

    W(t)=[W1(t)+W2(t)+W3(t)+W4(t)+

    W5(t)]×24.

    (3)

    故土壤源熱泵系統(tǒng)全壽命周期內(nèi)的總運(yùn)行費(fèi)用現(xiàn)值(S1)和常規(guī)空調(diào)系統(tǒng)全壽命周期內(nèi)的總運(yùn)行費(fèi)用現(xiàn)值(S2)可分別按式(4),(5)計(jì)算:

    (4)

    (5)

    式中:P為電價(jià);i為折現(xiàn)率,由于系統(tǒng)總運(yùn)行費(fèi)用現(xiàn)值的計(jì)算周期為全壽命周期,不可以忽略資金的時(shí)間價(jià)值(包括通貨膨脹、收回資金的不確定性及機(jī)會(huì)成本)[14],需要考慮折現(xiàn)率,本文采用2%~3%[15].

    1.3 系統(tǒng)失效指標(biāo)

    系統(tǒng)失效指標(biāo)是指土壤溫度T與保證熱泵系統(tǒng)全壽命周期內(nèi)的經(jīng)濟(jì)性和節(jié)能性的土壤溫度范圍M(Tmin,Tmax)的關(guān)系,指標(biāo)失效也即巖土溫度已不適合土壤源熱泵系統(tǒng)正常運(yùn)行,即

    1)當(dāng)T

    2)當(dāng)T∈M時(shí),指標(biāo)即將失效,建議采用相應(yīng)措施來緩解;

    3)當(dāng)T>Tmax時(shí),指標(biāo)已失效,必須采用相應(yīng)措施保證系統(tǒng)的正常運(yùn)行.

    同一建筑對于地埋管換熱器系統(tǒng)中的不同鉆孔,由于幾何位置不同,傳熱條件也不盡相同,本文在計(jì)算地埋管換熱器系統(tǒng)失效問題時(shí)選擇最不利的鉆孔來進(jìn)行判斷.

    GB 50366—2005《地源熱泵系統(tǒng)工程技術(shù)規(guī)范》規(guī)定地埋管地源熱泵系統(tǒng)熱泵機(jī)組正常運(yùn)行的冷熱源溫度范圍:10~40 ℃(制冷),-5~25 ℃(制熱),該溫度限值僅僅是針對保證機(jī)組正常運(yùn)行的低位冷卻水/輔熱水溫度范圍.指標(biāo)失效的主要原因是地下累積的冷(熱)超出了土壤的自平衡能力,失效指標(biāo)是土壤蓄能平衡評價(jià)的關(guān)鍵參數(shù),土壤蓄能是否平衡直接關(guān)聯(lián)指標(biāo)是否失效.本文提出的土壤蓄能狀態(tài)和計(jì)算方法是將全壽命周期內(nèi)的土壤源熱泵系統(tǒng)的節(jié)能性和經(jīng)濟(jì)性進(jìn)行耦合,并考慮全年動(dòng)態(tài)負(fù)荷和大地自平衡能力,可以較準(zhǔn)確地預(yù)測全壽命周期內(nèi)的節(jié)能性和經(jīng)濟(jì)性的失效區(qū)間以及對應(yīng)參數(shù).

    2 土壤源熱泵系統(tǒng)土壤蓄能狀態(tài)

    土壤源熱泵系統(tǒng)土壤蓄能狀態(tài)涉及到多參數(shù)評價(jià),為了定義土壤蓄能狀態(tài),首先要定義巖土失調(diào)溫度限值,而該定義又與各種巖土溫度相關(guān).

    2.1 與巖土溫度失調(diào)限值相關(guān)的各種巖土溫度

    2.1.1 巖土溫度

    實(shí)際情況下,在沿埋管半徑的水平方向上,巖土的溫度會(huì)隨距埋管中心的距離不同而有所變化;在沿埋管深度的豎直方向上,巖土的溫度會(huì)隨埋管深度的不同而有所變化.故巖土溫度的定義對象不能是一個(gè)簡單的點(diǎn)或是一個(gè)面,而是能夠綜合評價(jià)巖土溫度的定義.

    考慮到埋管近管壁處的巖土溫度對地埋管運(yùn)行周期內(nèi)埋管的換熱性能及埋管進(jìn)、出口水溫度的影響較大,從而影響熱泵機(jī)組的能耗及系統(tǒng)的能效比,故將影響系統(tǒng)能耗及能效比的巖土溫度定義為:埋管回填孔壁及其以內(nèi)回填體的巖土溫度[16].巖土溫度t*的數(shù)學(xué)表達(dá)式為:

    (6)

    式中:V為整個(gè)回填區(qū)域體;txyz為V區(qū)域體內(nèi)各單元體積的溫度值.

    2.1.2 系統(tǒng)能效溫度限值的定義

    當(dāng)巖土溫度t*達(dá)到某一值時(shí),此巖土溫度下的埋管出水溫度導(dǎo)致系統(tǒng)節(jié)能率小于零,這個(gè)溫度值就是巖土的能效溫度限值.能效溫度限值表征的是兩個(gè)系統(tǒng)能效比高低的臨界點(diǎn),是瞬態(tài)值,不適用于全壽命周期內(nèi)巖土溫度限值的評價(jià).

    2.1.3 系統(tǒng)運(yùn)行費(fèi)用溫度限值的定義

    當(dāng)巖土溫度t*達(dá)到某一值時(shí),土壤源熱泵系統(tǒng)的累計(jì)總運(yùn)行費(fèi)用現(xiàn)值(S1)超過常規(guī)空調(diào)系統(tǒng)的累計(jì)總運(yùn)行費(fèi)用現(xiàn)值(S2),這個(gè)溫度值就是巖土的運(yùn)行費(fèi)用溫度限值.由于運(yùn)行費(fèi)用溫度限值從能效的累積上對巖土溫度進(jìn)行了分析,故可以根據(jù)它來評價(jià)巖土溫度失調(diào)限值.

    圖1為土壤源熱泵系統(tǒng)與常規(guī)空調(diào)系統(tǒng)運(yùn)行費(fèi)用現(xiàn)值之差(S1-S2)的趨勢圖.從第4年開始,系統(tǒng)的運(yùn)行費(fèi)用就由于熱累積效應(yīng)開始逐漸上升,但是系統(tǒng)仍處于經(jīng)濟(jì)區(qū),第7年開始是震蕩區(qū),這是因?yàn)榇蟮貛r土自身的冷熱調(diào)節(jié),系統(tǒng)運(yùn)行費(fèi)用隨季節(jié)變化的波動(dòng).在第9年左右出現(xiàn)了一個(gè)運(yùn)行費(fèi)用溫度限值點(diǎn),當(dāng)經(jīng)過這個(gè)點(diǎn)以后S1恒大于S2,系統(tǒng)進(jìn)入非經(jīng)濟(jì)區(qū).

    年度

    由于在土壤源熱泵系統(tǒng)全壽命周期的運(yùn)行時(shí)間內(nèi),可能會(huì)出現(xiàn)多個(gè)運(yùn)行費(fèi)用溫度限值點(diǎn),則土壤源熱泵系統(tǒng)的巖土失調(diào)溫度限值范圍應(yīng)為最小值點(diǎn)至最大值點(diǎn).

    2.2 巖土失調(diào)溫度限值

    巖土失調(diào)溫度限值可定義為:在空調(diào)系統(tǒng)正常使用周期內(nèi),以給定的建筑負(fù)荷特征為基礎(chǔ),計(jì)算時(shí)間以年為單位,當(dāng)土壤源熱泵系統(tǒng)與常規(guī)系統(tǒng)的運(yùn)行費(fèi)用之差為零時(shí)求解得到的對應(yīng)溫度,其溫度變化范圍為區(qū)間L(Tmin,Tmax).將區(qū)間L賦值于系統(tǒng)失效指標(biāo)的區(qū)間M,即可得到針對工程實(shí)際的保證系統(tǒng)節(jié)能性和經(jīng)濟(jì)性的系統(tǒng)失效指標(biāo).該參數(shù)可以用系統(tǒng)節(jié)能率f和系統(tǒng)總運(yùn)行費(fèi)用現(xiàn)值S來表征.

    2.3 土壤蓄能狀態(tài)

    在理解巖土溫度失調(diào)限值定義的基礎(chǔ)上,就可以定義土壤蓄能狀態(tài):在空調(diào)系統(tǒng)正常使用周期內(nèi),以給定的建筑負(fù)荷特征為基礎(chǔ),計(jì)算時(shí)間以年為單位,如果系統(tǒng)節(jié)能率恒大于零,且系統(tǒng)失效指標(biāo)不失效,則稱該負(fù)荷特征對應(yīng)下的土壤蓄能處于平衡狀態(tài);如果系統(tǒng)節(jié)能率小于零,或系統(tǒng)失效指標(biāo)失效,則稱該負(fù)荷特征對應(yīng)下的土壤蓄能處于不平衡狀態(tài).

    3 計(jì)算方法、數(shù)學(xué)模型描述

    3.1 建筑逐時(shí)能耗計(jì)算

    本文采用DeST軟件作為建筑能耗模擬軟件.根據(jù)項(xiàng)目的施工圖,按照建筑實(shí)際情況建立模型,室內(nèi)設(shè)計(jì)參數(shù)參考GB 50189-2005《公共建筑節(jié)能設(shè)計(jì)標(biāo)準(zhǔn)》,可以得到示例工程的逐時(shí)負(fù)荷圖.

    3.2 末端負(fù)荷至出力負(fù)荷轉(zhuǎn)換數(shù)學(xué)模型

    根據(jù)熱泵機(jī)組的壓縮機(jī)級數(shù)調(diào)節(jié)和啟??刂撇呗詠硖幚斫ㄖ饡r(shí)負(fù)荷計(jì)算結(jié)果,將末端側(cè)的冷熱負(fù)荷需求量Q1轉(zhuǎn)化為熱泵機(jī)組的實(shí)際冷熱負(fù)荷輸出量Q2.以夏季制冷工況為例,以具體項(xiàng)目的熱泵機(jī)組的壓縮機(jī)實(shí)際調(diào)節(jié)級數(shù)為基礎(chǔ),將建筑逐時(shí)負(fù)荷導(dǎo)入編譯的程序后,得到機(jī)組的出力逐時(shí)負(fù)荷[9].

    為保證土壤源熱泵系統(tǒng)與常規(guī)空調(diào)系統(tǒng)有一個(gè)比較的平臺(tái),假設(shè)冷水機(jī)組在夏季的啟停及運(yùn)行情況的確定方法與土壤源熱泵機(jī)組一致,而冬季燃?xì)忮仩t的啟停及運(yùn)行情況則由建筑負(fù)荷直接決定.根據(jù)工程的實(shí)際冷水機(jī)組及燃?xì)忮仩t選型可以得到常規(guī)空調(diào)系統(tǒng)的逐時(shí)出力負(fù)荷.

    3.3 動(dòng)態(tài)能效比下出力負(fù)荷與進(jìn)出水溫耦合關(guān)系數(shù)學(xué)模型

    3.3.1 土壤源熱泵機(jī)組耦合關(guān)系式

    熱泵機(jī)組的出力負(fù)荷Q2與向大地的放(取)熱量Q3之間的關(guān)系式在文獻(xiàn)[17]中可以查得,該式適用于冬、夏季工況.

    EER=Q2/Php,

    (7)

    (8)

    (9)

    3.3.2 冷水機(jī)組耦合關(guān)系式

    常規(guī)空調(diào)系統(tǒng)夏季制冷依靠冷水機(jī)組和冷卻塔,若需要進(jìn)行冷水機(jī)組的動(dòng)態(tài)能效分析,需要建立冷卻塔模型,以逆流式冷卻塔作為討論對象.冷卻塔的冷卻數(shù)N可以由冷卻塔進(jìn)水溫度t1,出水溫度t2,室外空氣濕球溫度ts三者組成的函數(shù)計(jì)算[18].

    根據(jù)工程實(shí)際數(shù)據(jù)求得冷卻塔的特性數(shù)N′=1.07,再按照圖2所示的流程圖進(jìn)行冷卻塔出水溫度的求解.

    圖2 冷卻塔出水溫度的求解流程圖

    3.4 土壤源熱泵系統(tǒng)及常規(guī)空調(diào)系統(tǒng)的全年能耗數(shù)學(xué)模型

    3.4.1 土壤源熱泵系統(tǒng)全年能耗數(shù)學(xué)模型

    土壤源熱泵系統(tǒng)的能耗模型為:

    (10)

    式中:Q(t)為熱泵機(jī)組的制冷(熱)量;W1(t)為熱泵機(jī)組的運(yùn)行能耗;W2(t)為地埋管側(cè)水泵的運(yùn)行能耗;W3(t)為用戶側(cè)循環(huán)水泵的運(yùn)行能耗;W4(t)為水處理儀能耗;W5(t)為末端風(fēng)系統(tǒng)的運(yùn)行能耗;t為空調(diào)的運(yùn)行時(shí)刻.

    3.4.2 常規(guī)空調(diào)系統(tǒng)全年能耗數(shù)學(xué)模型

    常規(guī)系統(tǒng)的系統(tǒng)能耗數(shù)學(xué)方程式可用式(11)表示:

    EERs2(t)=

    (11)

    式中:Q(t)為夏季制冷工況下為冷水機(jī)組的制冷量,冬季制熱工況下為燃?xì)忮仩t的制熱量;N1(t)為夏季制冷工況下為冷水機(jī)組的運(yùn)行能耗,冬季制熱工況下為燃?xì)忮仩t的運(yùn)行能耗;N2(t)為夏季制冷工況下,冷卻側(cè)水泵的運(yùn)行能耗;N3(t)為用戶側(cè)循環(huán)水泵的運(yùn)行能耗;N4(t)為水處理儀能耗;N5(t)為末端風(fēng)系統(tǒng)的運(yùn)行能耗;N6(t) 為夏季制冷工況下,冷卻塔的運(yùn)行能耗;t為空調(diào)的運(yùn)行時(shí)刻.系統(tǒng)冷源為冷水機(jī)組+冷卻塔,熱源為鍋爐.對于冬季來說,W2(t)=W3(t)=0.

    由式(10)和(11)可以計(jì)算得到兩種系統(tǒng)的能效比,再由式(1)可以計(jì)算土壤源熱泵系統(tǒng)相比于常規(guī)空調(diào)系統(tǒng)的系統(tǒng)節(jié)能率,從而評價(jià)土壤源熱泵系統(tǒng)的節(jié)能性.

    上述能耗模型中各項(xiàng)的詳細(xì)能耗數(shù)學(xué)模型在很多論文中都有提及,模型已比較完善,本文采用課題組建立的已有模型[11-12].

    4 基于三維傳熱管群模型的實(shí)際工程計(jì)算分析

    4.1 工程概況

    該示例工程位于重慶市,空調(diào)系統(tǒng)在夏季為該樓裙房(共4層,第1層主要為食堂;第2層主要為辦公室和會(huì)議室;第3和第4層為專家公寓)供冷,總冷負(fù)荷為231.21kW.冬季為該樓第1層食堂、第2層辦公室和會(huì)議室的低溫輻射地板及副樓羽毛球場供暖,總熱負(fù)荷為230kW.地埋管換熱器布置在建筑旁邊的草坪下面,地下?lián)Q熱器系統(tǒng)由60個(gè)鉆孔組成,鉆孔成6×10矩形布置,鉆孔直徑110mm,鉆孔之間的距離為4m,換熱器采用直徑為25mm,深度為80m的PE管,兩支管間距為50mm,埋管換熱器回水管上設(shè)有60m深的保溫層,管群沒有分區(qū)[12].

    該工程選用2臺(tái)3機(jī)頭熱泵機(jī)組,壓縮機(jī)級數(shù)調(diào)節(jié)分為“0%,-33%,-66%,-100%” 4個(gè)等級,壓縮機(jī)運(yùn)行級數(shù)的選擇由末端負(fù)荷大小決定.熱泵機(jī)組的啟停控制策略見表1,控制策略參數(shù)分為冬季運(yùn)行工況和夏季運(yùn)行工況.

    表1 熱泵機(jī)組啟停控制策略

    4.2 地下巖土三維傳熱管群模型

    數(shù)值計(jì)算的幾何模型按示例工程實(shí)際情況建立,所以該模型是高度仿真的.模型中共60組地埋管,取1/4區(qū)域埋管,其布置示意圖如圖3所示.

    圖3 三維管群模型平面布置示意圖

    為了便于計(jì)算分析,作出以下簡化:

    1)埋管內(nèi)液體的流速在徑向上均勻一致;

    2)假定巖土的導(dǎo)熱系數(shù)、比熱、密度等物性參數(shù)不隨溫度的變化而變化,且是均勻一致的;

    3)地埋管同截面具有相同的溫度和流速;

    4)忽略巖土中水分遷移而引起的熱濕遷移;

    5)忽略巖土表面夜間輻射的得熱、失熱;

    6)忽略回填材料和孔洞的接觸熱阻.

    鉆孔孔深80 m,管間距為4 m,鉆孔直徑130 mm,進(jìn)、出水管為DN25,管內(nèi)流速為0.818 m/s (按所選熱泵機(jī)組流量計(jì)算而得),三維模型遠(yuǎn)邊界延伸至管群外10 m.整個(gè)模型為一個(gè)28 m× 20 m ×80 m(長×寬×高)的長方體,采用Gambit建模,總網(wǎng)格數(shù)約120萬.網(wǎng)格劃分情況如圖4所示.

    將系統(tǒng)能效測試解和數(shù)值解進(jìn)行對比分析,可以得到該模型的系統(tǒng)能效偏差,地埋管溫度偏差都在誤差范圍之內(nèi),故可知所建立的三維管群模型用于地埋管換熱器土壤熱響應(yīng)計(jì)算是可行的[12].

    圖4 三維管群模型網(wǎng)格劃分圖

    5 土壤源熱泵系統(tǒng)土壤蓄能狀態(tài)分析的一般計(jì)算方法流程圖

    土壤源熱泵系統(tǒng)土壤蓄能狀態(tài)分析的一般計(jì)算方法流程圖如圖5所示.

    1)利用DeST軟件計(jì)算得到建筑全年逐時(shí)末端負(fù)荷;

    2)利用末端負(fù)荷至出力負(fù)荷轉(zhuǎn)換數(shù)學(xué)模型求得熱泵機(jī)組及冷水機(jī)組(鍋爐)的全年出力負(fù)荷;

    3)利用出力負(fù)荷與進(jìn)出水溫耦合關(guān)系數(shù)學(xué)模型求得動(dòng)態(tài)能效比下的數(shù)值計(jì)算模型埋管進(jìn)出水溫度;

    4)將埋管進(jìn)、出水溫度的動(dòng)態(tài)耦合關(guān)系式通過UDF作為埋管進(jìn)口的邊界條件,導(dǎo)入高度仿真的數(shù)值計(jì)算模型進(jìn)行全壽命周期的流動(dòng)、換熱計(jì)算,求得全壽命周期內(nèi)地下?lián)Q熱器埋管進(jìn)出水溫度及土壤溫度分布場;

    5)將數(shù)值計(jì)算得到的全壽命周期埋管出水溫度導(dǎo)入土壤源熱泵系統(tǒng)的能耗模型,再代入S1,S2,EERs1,EERs2,f的定義式求得土壤源熱泵系統(tǒng)及常規(guī)空調(diào)系統(tǒng)的全壽命周期運(yùn)行費(fèi)用現(xiàn)值S1和S2和系統(tǒng)節(jié)能率f.

    6 計(jì)算結(jié)果分析

    6.1 土壤蓄能處于平衡狀態(tài)的情況

    對4.1節(jié)中提到的4層建筑所有負(fù)荷進(jìn)行分析計(jì)算.冬夏季累計(jì)負(fù)荷比為1∶1.55.

    按圖5所示流程圖進(jìn)行運(yùn)算,最終系統(tǒng)節(jié)能率和系統(tǒng)運(yùn)行費(fèi)用現(xiàn)值分別如圖6和圖7所示.

    圖5 一般計(jì)算方法流程圖

    年度

    年度

    按照本文所提出的計(jì)算方法可以得知,全壽命周期內(nèi),土壤源熱泵系統(tǒng)的能效比始終比常規(guī)空調(diào)系統(tǒng)高,系統(tǒng)節(jié)能率f>0;同時(shí)S1-S2=0方程無解,土壤源熱泵系統(tǒng)運(yùn)行費(fèi)用現(xiàn)值恒小于常規(guī)空調(diào)系統(tǒng)運(yùn)行費(fèi)用現(xiàn)值,也即不存在運(yùn)行費(fèi)用溫度限值,所以該土壤源熱泵系統(tǒng)在此工況下運(yùn)行不存在巖土溫度失調(diào)限值,區(qū)間為空集?.

    再按照本文所提出的計(jì)算方法可以得知,系統(tǒng)節(jié)能率f>0,且土壤溫度T∈?恒不成立,系統(tǒng)失效指標(biāo)恒不失效,所以該土壤源熱泵系統(tǒng)土壤蓄能處于平衡狀態(tài).綜上可知,該熱泵系統(tǒng)在全壽命周期內(nèi),以4層樓的全部建筑負(fù)荷為基礎(chǔ),以年為單位計(jì)算時(shí)間,系統(tǒng)不存在巖土溫度失調(diào)限值,且大地蓄能是平衡的,不必采用相應(yīng)措施來主動(dòng)調(diào)節(jié)大地蓄能,從4年實(shí)際運(yùn)行效果看,系統(tǒng)運(yùn)行是穩(wěn)定的.

    6.2 土壤蓄能處于不平衡狀態(tài)的情況

    為了尋求土壤蓄能處于不平衡狀態(tài)的情況,現(xiàn)假設(shè)該項(xiàng)目冬季供暖的副樓羽毛球場不使用,即對4.1節(jié)中提到的4層建筑所有負(fù)荷除冬季副樓羽毛球場的供暖負(fù)荷外按照計(jì)算方法對該工程重新進(jìn)行計(jì)算分析.冬夏季累計(jì)負(fù)荷比為1∶4.82.

    經(jīng)流程圖運(yùn)算,最終系統(tǒng)節(jié)能率和系統(tǒng)運(yùn)行費(fèi)用現(xiàn)值分別如圖8和圖9所示.

    按照本文所提出的計(jì)算方法可以得知,全壽命周期內(nèi),可求得S1-S2=0出現(xiàn)了最小,最大解——A和B,由運(yùn)行費(fèi)用溫度限值定義和計(jì)算方法可知,將A,B點(diǎn)對應(yīng)查詢土壤溫度,得到限值區(qū)間L(TA,TB)就是該熱泵系統(tǒng)的巖土失調(diào)溫度限值區(qū)間L(31.22 ℃,34.26 ℃),將區(qū)間L賦值于系統(tǒng)失效指標(biāo)的區(qū)間M,得到系統(tǒng)失效指標(biāo)的土壤溫度范圍為M(31.22 ℃,34.26 ℃).

    按照本文所提出的計(jì)算方法可以得知:

    1)開始運(yùn)行至第1年夏季,土壤源熱泵系統(tǒng)能效比要比常規(guī)冷熱源系統(tǒng)能效比高,且土壤源熱泵系統(tǒng)運(yùn)行費(fèi)用現(xiàn)值要比常規(guī)冷熱源運(yùn)行費(fèi)用現(xiàn)值低,所以該熱泵系統(tǒng)土壤蓄能處于平衡狀態(tài).

    年度

    年度

    2)第1年夏季至第7年秋季,土壤源熱泵系統(tǒng)能效比要比常規(guī)冷熱源系統(tǒng)能效比低,但土壤源熱泵系統(tǒng)運(yùn)行費(fèi)用現(xiàn)值也要比常規(guī)冷熱源運(yùn)行費(fèi)用現(xiàn)值低,建議采用相應(yīng)措施來主動(dòng)調(diào)節(jié)大地蓄能,以緩解土壤蓄能不平衡狀態(tài).

    3)第7年秋季至第8年夏季及以后,土壤源熱泵系統(tǒng)能效比要比常規(guī)冷熱源系統(tǒng)能效比低,且土壤源熱泵系統(tǒng)運(yùn)行費(fèi)用現(xiàn)值要比常規(guī)冷熱源運(yùn)行費(fèi)用現(xiàn)值高,必須采用相應(yīng)措施來主動(dòng)調(diào)節(jié)大地蓄能,以保證系統(tǒng)的正常運(yùn)行.

    7 結(jié) 論

    1) 建立了以系統(tǒng)節(jié)能率、系統(tǒng)總運(yùn)行費(fèi)用現(xiàn)值和系統(tǒng)失效指標(biāo)等參數(shù)為基礎(chǔ)的評價(jià)模型,可以客觀具體地反映土壤蓄能狀態(tài),同時(shí)利用評價(jià)模型對土壤蓄能狀態(tài)進(jìn)行了評價(jià).

    2) 提出了一種能分析評價(jià)土壤源熱泵系統(tǒng)的土壤蓄能狀態(tài)的一般計(jì)算方法,利用該計(jì)算方法能夠?qū)ν寥涝礋岜孟到y(tǒng)的土壤蓄能的平衡狀態(tài)或不平衡狀態(tài)進(jìn)行定量分析.

    3) 利用該計(jì)算方法得出了可以保證系統(tǒng)全壽命周期內(nèi)的節(jié)能性和經(jīng)濟(jì)性的巖土失調(diào)溫度限值范圍,同時(shí)結(jié)合運(yùn)行后的工程實(shí)際,對土壤源熱泵的土壤蓄能狀態(tài)進(jìn)行評價(jià),從而對全壽命周期內(nèi)的節(jié)能性和經(jīng)濟(jì)性的失效時(shí)間點(diǎn)進(jìn)行預(yù)測,以便維護(hù)人員提前做好相應(yīng)措施來滿足系統(tǒng)正常運(yùn)行的需要.

    4) 土壤蓄能平衡研究涉及到的參數(shù)較多,包括地下水流動(dòng)的影響等,沒有在評價(jià)模型中體現(xiàn),在后續(xù)的研究中,將不斷改進(jìn)和完善模型,同時(shí)簡化計(jì)算方法,使其更實(shí)用于工程實(shí)際.

    [1] 花莉,潘毅群,范蕊, 等. 基于TRNSYS的土壤源熱泵熱平衡問題的影響因素分析[J].建筑節(jié)能,2012,40(3):23-29.

    HUA Li,PAN Yi-qun,FAN Rui,etal. Impacting factor on ground heat balance of GSHP system using TRNSYS[J]. Buliding Energy Efficiency,2012,40(3):23-29.(In Chinese)

    [2] 范蕊,高巖,陳旭, 等. 基于熱平衡的土壤源熱泵系統(tǒng)特性分析[J].同濟(jì)大學(xué)學(xué)報(bào):自然科學(xué)版,2011,39(2):282-286.

    FAN Rui,GAO Yan, CHEN Xu,etal. Analysis on influence of yearly heat imbalance on heat transfer characteristic of ground-coupled heat pump system[J].Journal of Tongji University:Natural Science,2011,39(2):282-286.(In Chinese)

    [3] 馬宏權(quán),龍惟定. 地埋管地源熱泵系統(tǒng)的熱平衡[J].暖通空調(diào), 2009,39(1):102-106.

    MA Hong-quan, LONG Wei-ding. Ground heat balance in GSHP[J]. Heating Ventilating and Air Conditioning,2009,39(1):102-106. (In Chinese)

    [4] 林東超,沈劍. 土壤源熱泵熱平衡的驗(yàn)證[J].發(fā)電與空調(diào), 2013, 34(2):44-48.

    LIN Dong-chao, SHEN Jian. Verification of heat balance in ground-source heat pump[J].Refrigeration Air Conditioning and Electric Power Machinery,2013,34(2):44-48. (In Chinese)

    [5] 楊衛(wèi)波,陳振乾,劉光遠(yuǎn). 土壤源熱泵系統(tǒng)地下熱平衡問題分析[C]//吳元煒.中國制冷學(xué)會(huì)2009年學(xué)術(shù)年會(huì)論文集.天津:中國制冷學(xué)會(huì),2009:963-968.

    YANG Wei-bo,CHEN Zhen-qian,LIU Guang-yuan. Analysis on the ground heat balance of ground coupled heat pump system[C]//WU Yuan-wei. Academic Conference Proceedings of Chinese Association of Refrigeration in 2009.Tianjin:Chinese Association of Refrigeration, 2009:963-968. (In Chinese)

    [6] 楊紅輝. 土壤源熱泵系統(tǒng)的地埋管熱平衡分析[J]. 工程建設(shè)與設(shè)計(jì),2012,10:100-104.

    YANG Hong-hui.Heat balance analysis in ground-source heat pump[J]. Construction and Design for Project,2012,10:100-104. (In Chinese)

    [7] 范龍華,王勇. 土壤源熱泵系統(tǒng)中地埋管換熱器熱平衡問題及解決方案[J]. 黑龍江冶金, 2011,31(1) :44-45.

    FAN Long-hua,WANG Yong. Heat balance problem ground heat exchanger and solutions in ground-source heat pump system [J]. Heilongjiang Yejin, 2011,31(1) :44-45.(In Chinese)

    [8] 崔宏楨. 地源熱泵系統(tǒng)土壤熱平衡分析[J]. 工程建設(shè)與設(shè)計(jì), 2011, 11:95-98.

    CUI Hong-zhen. Heat balance analysis of soil for the ground-source heat pump system[J].Construction and Design for Project,2011,11:95-98. (In Chinese)

    [9] 李永安,高亞南,邢滕,等. 土壤源熱泵蓄能系統(tǒng)的研究與應(yīng)用[J]. 湖南大學(xué)學(xué)報(bào):自然科學(xué)版,2009,36(12):40-43.

    LI Yong-an,GAO Ya-nan,XING Teng,etal. Research and application of ground-coupled heat pump and ice storage [J]. Journal of Hunan University:Natural Sciences,2009,36(12):40-43. (In Chinese)

    [10]楊昌智,黃兵.U型管換熱性能影響因素研究[J].湖南大學(xué)學(xué)報(bào):自然科學(xué)版,2009,36(12):44-48.

    YANG Chang-zhi,HUANG Bing.Study on influence factor of heat transfer performance of the U-tube[J].Journal of Hunan University:Natural Sciences,2009,36(12):44-48.(In Chinese)

    [11]金逸韜.夏熱冬冷地區(qū)土壤源熱泵系統(tǒng)的蓄能失調(diào)限值分析[D].重慶:重慶大學(xué)城市建設(shè)與環(huán)境工程學(xué)院,2012:30-33.

    JIN Yi-tao.The energy storage imbalance limitanalysis of ground source heat pumpsystem in hot-summer and cold-winter zone[D].Chongqing: College of Urban Construction and Environmental Engineering,Chongqing University,2012:30-33.(In Chinese)

    [12]賈亞北.典型氣候區(qū)典型建筑土壤源熱泵系統(tǒng)蓄能不平衡率分析[D].重慶:重慶大學(xué)城市建設(shè)與環(huán)境工程學(xué)院,2012:14-19.

    JIA Ya-bei.Analysis on heat storage unbalance rate in ground source heat pump system of the typical buildings of the typical climate area[D].Chongqing: College of Urban Construction and Environmental Engineering,Chongqing University,2012:14-19.(In Chinese)

    [13]王勇,韓傳璞,李文,等. 開式地表水源熱泵系統(tǒng)取水能耗限值確定方法[J]. 重慶大學(xué)學(xué)報(bào),2011,34(3):110-115.

    WANG Yong,HAN Chuan-pu,LI Wen,etal. The determination method of limit energy consumption of water supply in open-loop surface water source heatpump systems[J]. Journal of Chongqing University,2011,34(3): 110-115. (In Chinese)

    [14]王光偉. 貨幣,利率與匯率經(jīng)濟(jì)學(xué)[M].北京:清華大學(xué)出版社,2003:107-111.

    WANG Guang-wei.Currency, interest rate and exchange rate economics[M].Beijing:Tsinghua University Press,2003:107-111.(In Chinese)

    [15]周立,葛耀君. 全壽命經(jīng)濟(jì)分析中折現(xiàn)率的確定[J]. 上海公路,2007(2):51-54.

    ZHOU Li,GE Yao-jun. Discount rates of the life cycle cost analysis[J]. Shanghai Highways,2007(2):51-54. (In Chinese)

    [16] 王勇,金逸韜. 回填空氣間隙對地埋管巖土溫度恢復(fù)性能的影響[J]. 土木建筑與環(huán)境工程,2012,34(4):142-148.

    WANG Yong, JIN Yi-tao. Influence of backfill air gap on soil temperature rescovery in ground-coupled heat pumps[J].Journal of Civil,Architectural & Environmental Engineering,2012, 34(4):142-148.(In Chinese)

    [17]GB 50366-2009 地源熱泵系統(tǒng)工程技術(shù)規(guī)范[S].北京:中國建筑工業(yè)出版社,2009.

    GB 50366-2009 Technical code for ground-sourse heat pump system[S].Beijing:China Architecture & Building Press,2009 .(In Chinese)

    [18]王景昱.重慶市公共建筑節(jié)能75%設(shè)計(jì)技術(shù)路線及冷源系統(tǒng)節(jié)能指標(biāo)研究[D].重慶:重慶大學(xué)城市建設(shè)與環(huán)境工程學(xué)院,2012:32-40.

    WANG Jing-yu.Analysis on technical route of “chongqingdesign for energy efficiency 75% of public building” and energy saving index of the cold source system[D].Chongqing: College of Urban Construction and Environmental Engineering,Chongqing University,2012:32-40.(In Chinese)

    Evaluation of the Ground Energy Storage Condition of Ground Source Heat Pump System Based on Energy Storage Imbalance Temperature Limit

    WANG Yong1,2?,YIN Chang-yu1,2,JIN Yi-tao3

    (1.Key Laboratory of the Three Gorges Reservoir Region's Eco-Environment, Ministry of Education, Chongqing Univ, Chongqing 400045, China; 2.National Centre for International Research of Low-carbon and Green Buildings, Chongqing Univ, Chongqing 400045, China;3. Suzhou Industrial Park Design & Research Institute Co, Ltd, Suzhou, Jiangsu 215021, China)

    Based on the analysis of multi-parameter coupling, which is based on ground source imbalance temperature limit, the building load characteristics and the heat transfer mechanism, this paper put forward an evaluation index, including the energy storage imbalance temperature limit, energy imbalance, and system invalidating indicator, and built the evaluation system on ground energy storage condition. Based on this evaluation system, by using the three-dimensional pipe-group model created by CFD software, a calculation method to accurately analyze the ground energy storage condition of ground source heat pump system was proposed. The result of calculation and analysis show that the evaluation system and calculation method reflect ground energy storage condition objectively and the ground source imbalance temperature limit can be concluded on the basis of the characteristics of practical projects in this method. The failure time of energy saving and efficiency can be predicted in the life cycle.

    geothermal heat pumps; life cycle; temperature limit; CFD simulation; ground energy storage condition

    1674-2974(2015)01-0127-09

    2014-02-13

    國家自然科學(xué)基金資助項(xiàng)目(51178482),National Natural Science Foundation of China(51178482);高等學(xué)校學(xué)科創(chuàng)新引智計(jì)劃資助項(xiàng)目(B13041)

    王 勇(1971-),男,重慶人,重慶大學(xué)教授,博士?通訊聯(lián)系人,E-mail:wyfree1@126.com

    TU381

    A

    猜你喜歡
    現(xiàn)值源熱泵能效
    暖通空調(diào)設(shè)計(jì)中地源熱泵實(shí)踐
    空氣源熱泵用于天然氣加熱的經(jīng)濟(jì)環(huán)保性
    煤氣與熱力(2021年9期)2021-11-06 05:22:54
    空氣源熱泵供暖期耗電量計(jì)算
    煤氣與熱力(2021年6期)2021-07-28 07:21:18
    上海:穩(wěn)中有進(jìn) 能效趨優(yōu)
    關(guān)注能效
    地源熱泵系統(tǒng)的研究與應(yīng)用
    資金時(shí)間價(jià)值中的系數(shù)關(guān)系探析
    資金時(shí)間價(jià)值基礎(chǔ)運(yùn)算解讀
    凈現(xiàn)值法對比煤層氣與常規(guī)天然氣經(jīng)濟(jì)效益
    中國煤層氣(2015年2期)2015-08-22 03:29:15
    利用改進(jìn)凈現(xiàn)值法降低舊工業(yè)建筑投資不確定性研究
    啦啦啦观看免费观看视频高清| 久久久色成人| 韩国av在线不卡| 欧美xxxx性猛交bbbb| 一级a爱片免费观看的视频| 老司机午夜福利在线观看视频| 69av精品久久久久久| 久久婷婷人人爽人人干人人爱| 91狼人影院| 又爽又黄无遮挡网站| 69人妻影院| 赤兔流量卡办理| 国产女主播在线喷水免费视频网站 | 欧美成人免费av一区二区三区| 99久久无色码亚洲精品果冻| 日本一二三区视频观看| 1024手机看黄色片| 欧美xxxx性猛交bbbb| 少妇人妻精品综合一区二区 | 亚洲中文字幕日韩| 1000部很黄的大片| 精品久久久久久久久亚洲 | 亚洲欧美日韩东京热| 男人舔奶头视频| 有码 亚洲区| 久久久色成人| 赤兔流量卡办理| 一区二区三区高清视频在线| 国产黄a三级三级三级人| 国产精华一区二区三区| 可以在线观看毛片的网站| 中文亚洲av片在线观看爽| 国产精品av视频在线免费观看| 天美传媒精品一区二区| 久久精品国产99精品国产亚洲性色| 在线国产一区二区在线| 国产 一区精品| 最近视频中文字幕2019在线8| 乱系列少妇在线播放| 99热这里只有是精品50| 成年女人看的毛片在线观看| 国产精品三级大全| 桃色一区二区三区在线观看| 亚洲最大成人中文| 精品国产三级普通话版| 欧美在线一区亚洲| 亚洲天堂国产精品一区在线| 色综合色国产| 午夜日韩欧美国产| 夜夜爽天天搞| 国产精品女同一区二区软件 | 99热这里只有精品一区| 网址你懂的国产日韩在线| 婷婷色综合大香蕉| 亚洲精品影视一区二区三区av| 99热这里只有精品一区| av福利片在线观看| 五月玫瑰六月丁香| 欧美xxxx性猛交bbbb| 看片在线看免费视频| 国产高潮美女av| 欧美精品国产亚洲| 日韩高清综合在线| 日韩一区二区视频免费看| 91久久精品电影网| 国产高清视频在线观看网站| 在线免费观看不下载黄p国产 | 免费大片18禁| 久久久国产成人免费| 啦啦啦观看免费观看视频高清| 午夜福利在线观看免费完整高清在 | 老司机深夜福利视频在线观看| 精品人妻1区二区| 淫秽高清视频在线观看| 亚洲欧美精品综合久久99| av中文乱码字幕在线| 波多野结衣高清无吗| 午夜福利在线在线| 精品午夜福利在线看| 狂野欧美激情性xxxx在线观看| 搡老妇女老女人老熟妇| 久久午夜亚洲精品久久| 色av中文字幕| 国产欧美日韩精品一区二区| 亚洲欧美日韩东京热| 女同久久另类99精品国产91| 亚洲乱码一区二区免费版| 内地一区二区视频在线| 国产69精品久久久久777片| 欧美极品一区二区三区四区| 淫妇啪啪啪对白视频| 少妇的逼水好多| 日日摸夜夜添夜夜添小说| 欧美日本视频| 在线观看66精品国产| 免费在线观看影片大全网站| 美女黄网站色视频| 国产真实乱freesex| 99久久精品一区二区三区| 国产乱人视频| 一区二区三区四区激情视频 | 国产人妻一区二区三区在| 联通29元200g的流量卡| 欧美一区二区精品小视频在线| 成年女人永久免费观看视频| 男人狂女人下面高潮的视频| 亚洲精品一卡2卡三卡4卡5卡| 丰满的人妻完整版| 国产成人aa在线观看| 午夜免费男女啪啪视频观看 | 黄色视频,在线免费观看| 亚洲人成伊人成综合网2020| 国产毛片a区久久久久| 日本 av在线| 午夜老司机福利剧场| 中文在线观看免费www的网站| av福利片在线观看| 国产精品无大码| 赤兔流量卡办理| 欧美丝袜亚洲另类 | 久久久久久久久久成人| 日韩中文字幕欧美一区二区| 人人妻人人澡欧美一区二区| 一级黄色大片毛片| 男人舔奶头视频| 黄色女人牲交| 国产精品自产拍在线观看55亚洲| 日韩av在线大香蕉| 午夜精品久久久久久毛片777| 亚洲不卡免费看| 中文字幕av在线有码专区| 老熟妇乱子伦视频在线观看| 午夜福利高清视频| 人人妻人人澡欧美一区二区| 亚州av有码| 久久久国产成人免费| 看黄色毛片网站| 深夜a级毛片| 亚洲av第一区精品v没综合| 婷婷丁香在线五月| 国内精品美女久久久久久| 国产一区二区亚洲精品在线观看| 91麻豆av在线| 日本黄大片高清| 久久精品综合一区二区三区| 人妻久久中文字幕网| 变态另类成人亚洲欧美熟女| 亚洲av日韩精品久久久久久密| 色视频www国产| 国产精品国产高清国产av| 国产精品一及| 成人性生交大片免费视频hd| 精品久久久噜噜| 在线天堂最新版资源| 男女下面进入的视频免费午夜| 国产一区二区亚洲精品在线观看| 国产精华一区二区三区| 欧美高清性xxxxhd video| 国产成人一区二区在线| 亚洲av二区三区四区| 大型黄色视频在线免费观看| eeuss影院久久| 麻豆国产97在线/欧美| 欧美成人性av电影在线观看| 18禁黄网站禁片午夜丰满| 国产高清三级在线| 老司机深夜福利视频在线观看| 欧美一区二区精品小视频在线| 久久婷婷人人爽人人干人人爱| 深爱激情五月婷婷| 国产精品嫩草影院av在线观看 | 国产乱人伦免费视频| 少妇人妻一区二区三区视频| a级毛片a级免费在线| 乱系列少妇在线播放| 日韩欧美在线乱码| 91午夜精品亚洲一区二区三区 | 国产精品三级大全| 日日撸夜夜添| 露出奶头的视频| 给我免费播放毛片高清在线观看| 精品午夜福利视频在线观看一区| 国产伦一二天堂av在线观看| 亚洲经典国产精华液单| 久久久午夜欧美精品| 久久99热这里只有精品18| 国产三级中文精品| 亚洲精品影视一区二区三区av| 成人综合一区亚洲| 五月伊人婷婷丁香| 久99久视频精品免费| 亚洲精华国产精华精| 啦啦啦观看免费观看视频高清| 自拍偷自拍亚洲精品老妇| 嫩草影院入口| 搡老熟女国产l中国老女人| 热99在线观看视频| 狂野欧美白嫩少妇大欣赏| 欧美日韩瑟瑟在线播放| 久久久久久久午夜电影| 精品久久久久久成人av| 国产av在哪里看| 天美传媒精品一区二区| av在线老鸭窝| 国产不卡一卡二| 最近在线观看免费完整版| 国产精品亚洲美女久久久| 99热网站在线观看| 久久精品国产99精品国产亚洲性色| 在线天堂最新版资源| 舔av片在线| 色尼玛亚洲综合影院| 成人av一区二区三区在线看| videossex国产| 亚洲中文日韩欧美视频| 国内揄拍国产精品人妻在线| 两人在一起打扑克的视频| 国产老妇女一区| 亚洲av日韩精品久久久久久密| 亚洲男人的天堂狠狠| 91久久精品电影网| 久久久久久久精品吃奶| 欧美日本视频| 长腿黑丝高跟| 成年女人毛片免费观看观看9| 亚洲中文日韩欧美视频| 国产三级在线视频| 精品欧美国产一区二区三| 久久精品国产亚洲av涩爱 | av.在线天堂| av在线天堂中文字幕| 精品久久久久久久久久久久久| 亚洲av中文av极速乱 | 国产爱豆传媒在线观看| 午夜老司机福利剧场| 中文资源天堂在线| 长腿黑丝高跟| 老师上课跳d突然被开到最大视频| 免费看a级黄色片| 久久久久九九精品影院| 色综合色国产| 亚洲av一区综合| 亚洲av成人精品一区久久| 国产一区二区三区av在线 | 俄罗斯特黄特色一大片| 99精品久久久久人妻精品| 久久久久精品国产欧美久久久| 国产欧美日韩精品一区二区| av在线蜜桃| 悠悠久久av| 亚洲国产精品sss在线观看| 国产伦精品一区二区三区视频9| 国产av在哪里看| 亚洲 国产 在线| 婷婷色综合大香蕉| 精品乱码久久久久久99久播| 天堂动漫精品| 国产精品女同一区二区软件 | 又爽又黄a免费视频| 亚洲专区国产一区二区| 精品午夜福利在线看| 好男人在线观看高清免费视频| 午夜精品一区二区三区免费看| 午夜福利在线观看免费完整高清在 | 久久久久性生活片| 国产极品精品免费视频能看的| 99九九线精品视频在线观看视频| 麻豆一二三区av精品| 日日撸夜夜添| 国内久久婷婷六月综合欲色啪| 自拍偷自拍亚洲精品老妇| 99在线视频只有这里精品首页| 网址你懂的国产日韩在线| 国产精品精品国产色婷婷| 天堂网av新在线| 国产精品av视频在线免费观看| 美女高潮的动态| 久久久久久九九精品二区国产| 午夜福利在线观看吧| 欧美+日韩+精品| 亚洲精品粉嫩美女一区| 久久久精品大字幕| 香蕉av资源在线| 国产精品久久久久久久久免| 无遮挡黄片免费观看| 国产人妻一区二区三区在| 熟女电影av网| www.www免费av| 看免费成人av毛片| 观看美女的网站| 成人精品一区二区免费| 精品午夜福利在线看| 国产久久久一区二区三区| 无遮挡黄片免费观看| 午夜激情福利司机影院| 又粗又爽又猛毛片免费看| 禁无遮挡网站| 国产高清不卡午夜福利| 久久午夜福利片| 日本熟妇午夜| 日本a在线网址| 国产亚洲精品综合一区在线观看| 成年女人永久免费观看视频| 免费在线观看日本一区| 国产男人的电影天堂91| 精品久久久噜噜| 久久久久久久久大av| 婷婷丁香在线五月| 精品久久久久久久久亚洲 | 色综合亚洲欧美另类图片| 91精品国产九色| 午夜福利在线在线| 日本撒尿小便嘘嘘汇集6| 22中文网久久字幕| 久9热在线精品视频| 国产视频一区二区在线看| 丝袜美腿在线中文| 99久久九九国产精品国产免费| 日韩欧美免费精品| 成人永久免费在线观看视频| 亚洲性夜色夜夜综合| 又黄又爽又刺激的免费视频.| 村上凉子中文字幕在线| 性欧美人与动物交配| 97碰自拍视频| 国内精品久久久久久久电影| 此物有八面人人有两片| 中文字幕精品亚洲无线码一区| 国产一区二区三区视频了| 欧美性猛交╳xxx乱大交人| 直男gayav资源| 亚洲精品影视一区二区三区av| 国产亚洲欧美98| 天堂网av新在线| 中文字幕av在线有码专区| 久久久久免费精品人妻一区二区| 亚洲av二区三区四区| 99久久久亚洲精品蜜臀av| 国产精品精品国产色婷婷| 国产高清视频在线播放一区| 校园人妻丝袜中文字幕| 精品久久久噜噜| 午夜福利18| 最后的刺客免费高清国语| 午夜精品一区二区三区免费看| av在线观看视频网站免费| 欧美黑人巨大hd| 成人综合一区亚洲| 少妇丰满av| 99久久精品国产国产毛片| 免费av观看视频| 日本 av在线| 最近在线观看免费完整版| av在线亚洲专区| 男女做爰动态图高潮gif福利片| 午夜福利视频1000在线观看| 久久久久久久久久成人| 美女高潮喷水抽搐中文字幕| 欧美色欧美亚洲另类二区| 国产精品1区2区在线观看.| 精品久久久久久久久av| 欧美高清性xxxxhd video| 老熟妇仑乱视频hdxx| 搡老妇女老女人老熟妇| 精品国内亚洲2022精品成人| 国内揄拍国产精品人妻在线| 一进一出抽搐动态| 久久欧美精品欧美久久欧美| 免费在线观看成人毛片| 午夜免费成人在线视频| 麻豆av噜噜一区二区三区| 国产av在哪里看| 欧美色欧美亚洲另类二区| 亚洲精品在线观看二区| 欧美丝袜亚洲另类 | 精品欧美国产一区二区三| 亚洲国产高清在线一区二区三| 午夜老司机福利剧场| 欧美日韩综合久久久久久 | 亚洲无线观看免费| 嫩草影院新地址| 亚洲无线观看免费| 日本a在线网址| 天天一区二区日本电影三级| 99国产精品一区二区蜜桃av| 亚洲专区国产一区二区| 99久久精品国产国产毛片| 成人特级黄色片久久久久久久| 亚洲欧美清纯卡通| 麻豆久久精品国产亚洲av| 日本欧美国产在线视频| 日韩欧美 国产精品| 国产精品嫩草影院av在线观看 | 中亚洲国语对白在线视频| 一个人看的www免费观看视频| 免费人成视频x8x8入口观看| 女同久久另类99精品国产91| 久久午夜亚洲精品久久| 亚洲成人久久性| 成人性生交大片免费视频hd| 午夜免费成人在线视频| 成年女人毛片免费观看观看9| 一级a爱片免费观看的视频| 欧美成人免费av一区二区三区| 欧美bdsm另类| 精品久久久久久久久久久久久| 亚洲av.av天堂| 国产精品一区二区免费欧美| 成年女人永久免费观看视频| 亚洲无线观看免费| 亚洲欧美日韩卡通动漫| 黄色日韩在线| 成年人黄色毛片网站| 三级国产精品欧美在线观看| 欧美xxxx性猛交bbbb| 午夜福利在线在线| 国产精品久久久久久久久免| 久久婷婷人人爽人人干人人爱| 精品人妻一区二区三区麻豆 | 国国产精品蜜臀av免费| 色5月婷婷丁香| 天美传媒精品一区二区| 国语自产精品视频在线第100页| 嫩草影院入口| 永久网站在线| 午夜激情福利司机影院| 91av网一区二区| 亚洲av不卡在线观看| 久久精品国产自在天天线| 欧美高清性xxxxhd video| 香蕉av资源在线| 淫妇啪啪啪对白视频| 国产一区二区激情短视频| 色综合婷婷激情| 超碰av人人做人人爽久久| 少妇熟女aⅴ在线视频| 男女视频在线观看网站免费| 在线a可以看的网站| 毛片女人毛片| 久久欧美精品欧美久久欧美| 极品教师在线视频| 国产精品一及| 免费一级毛片在线播放高清视频| 我要搜黄色片| 亚洲经典国产精华液单| 精品欧美国产一区二区三| 大型黄色视频在线免费观看| 成年女人毛片免费观看观看9| 日本 av在线| 国产色爽女视频免费观看| 变态另类丝袜制服| 国产伦一二天堂av在线观看| 日本 欧美在线| 日本黄色视频三级网站网址| 三级毛片av免费| 精品人妻视频免费看| 亚洲不卡免费看| 88av欧美| 18+在线观看网站| 亚洲av第一区精品v没综合| 亚洲精品一卡2卡三卡4卡5卡| 久久国内精品自在自线图片| 欧美又色又爽又黄视频| 亚洲av电影不卡..在线观看| 免费电影在线观看免费观看| 天堂av国产一区二区熟女人妻| 中文在线观看免费www的网站| 日韩欧美精品v在线| 国内精品久久久久久久电影| 最新中文字幕久久久久| 99热这里只有是精品50| 乱人视频在线观看| 国产精品一区二区三区四区免费观看 | а√天堂www在线а√下载| 免费av不卡在线播放| 午夜精品在线福利| 国产三级在线视频| 久久精品国产鲁丝片午夜精品 | 国产精品亚洲美女久久久| 啪啪无遮挡十八禁网站| 国产老妇女一区| 午夜爱爱视频在线播放| 97人妻精品一区二区三区麻豆| 黄色丝袜av网址大全| 少妇的逼好多水| 韩国av一区二区三区四区| 亚洲av熟女| 国产午夜福利久久久久久| 国产麻豆成人av免费视频| 国产真实伦视频高清在线观看 | 床上黄色一级片| 国产精品99久久久久久久久| 一边摸一边抽搐一进一小说| 成人毛片a级毛片在线播放| 美女高潮的动态| 精品一区二区免费观看| 联通29元200g的流量卡| 国产伦人伦偷精品视频| 成人特级黄色片久久久久久久| 国产麻豆成人av免费视频| 婷婷色综合大香蕉| 亚洲在线观看片| 我要看日韩黄色一级片| 在线观看66精品国产| 久99久视频精品免费| 美女免费视频网站| 久久久精品大字幕| 久久久久久久午夜电影| 午夜日韩欧美国产| 我要看日韩黄色一级片| 综合色av麻豆| 日韩大尺度精品在线看网址| 我的女老师完整版在线观看| 最近最新中文字幕大全电影3| 国产欧美日韩精品一区二区| 精品人妻1区二区| 久久99热6这里只有精品| 国产av在哪里看| 午夜爱爱视频在线播放| 精品不卡国产一区二区三区| 久久久精品欧美日韩精品| 老女人水多毛片| 国产伦精品一区二区三区视频9| 亚洲avbb在线观看| 黄片wwwwww| 国产69精品久久久久777片| 国产免费一级a男人的天堂| 国产欧美日韩一区二区精品| 欧美激情久久久久久爽电影| 在线观看午夜福利视频| 身体一侧抽搐| 99热6这里只有精品| 亚洲自拍偷在线| 国产精品久久电影中文字幕| 国产精品国产三级国产av玫瑰| 国产色爽女视频免费观看| 精品午夜福利在线看| 午夜视频国产福利| 亚洲国产欧美人成| 日日摸夜夜添夜夜添av毛片 | 三级毛片av免费| 久久久久久九九精品二区国产| 好男人在线观看高清免费视频| 国产成人一区二区在线| 观看美女的网站| 亚洲va日本ⅴa欧美va伊人久久| 一夜夜www| 欧美激情久久久久久爽电影| 免费人成视频x8x8入口观看| 在线观看66精品国产| a在线观看视频网站| 亚洲中文字幕日韩| 亚洲欧美日韩高清专用| 久久久久久九九精品二区国产| 非洲黑人性xxxx精品又粗又长| 亚洲真实伦在线观看| 老司机午夜福利在线观看视频| 亚洲精品成人久久久久久| 成人亚洲精品av一区二区| 国产主播在线观看一区二区| 日韩欧美国产在线观看| 亚州av有码| 亚洲av中文字字幕乱码综合| 极品教师在线免费播放| 色综合色国产| 春色校园在线视频观看| 成人毛片a级毛片在线播放| aaaaa片日本免费| 国产蜜桃级精品一区二区三区| 欧美高清成人免费视频www| 最后的刺客免费高清国语| 午夜影院日韩av| 亚洲在线观看片| 亚洲欧美日韩卡通动漫| 在线观看av片永久免费下载| 高清在线国产一区| 舔av片在线| 久久精品91蜜桃| 91麻豆av在线| 男女边吃奶边做爰视频| 伦精品一区二区三区| 丝袜美腿在线中文| 久久精品91蜜桃| 日韩 亚洲 欧美在线| 亚洲乱码一区二区免费版| 免费黄网站久久成人精品| 久久中文看片网| 欧美潮喷喷水| 国产精品女同一区二区软件 | 一个人观看的视频www高清免费观看| 国产高清不卡午夜福利| 国产真实乱freesex| 午夜福利18| 国产av一区在线观看免费| 国产国拍精品亚洲av在线观看| 少妇的逼水好多| 日本免费a在线| 又黄又爽又免费观看的视频| 国产精华一区二区三区| 此物有八面人人有两片| 欧美极品一区二区三区四区| 高清日韩中文字幕在线| 国产av麻豆久久久久久久| 在线观看美女被高潮喷水网站| 婷婷精品国产亚洲av| 亚洲国产色片| 深夜a级毛片| 日本三级黄在线观看| 亚洲熟妇中文字幕五十中出| 尾随美女入室| 午夜福利在线观看吧| 国产爱豆传媒在线观看| 国产精品久久久久久久久免| 丰满人妻一区二区三区视频av| 亚洲中文日韩欧美视频|