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

    水源熱泵制冷季間歇運(yùn)行的地溫場變化研究

    2020-04-06 04:10劉愛華李娟王哲郭艷春
    城市地質(zhì) 2020年4期
    關(guān)鍵詞:含水層熱泵源熱泵

    劉愛華 李娟 王哲 郭艷春

    摘 要:地下水換熱系統(tǒng)實(shí)驗(yàn)場地位于北京市平原區(qū),地層巖性屬河流相沉積。建立水源熱泵抽灌區(qū)地下水滲流與熱量運(yùn)移的三維耦合數(shù)值模型,采用COMSOL Multiphysics有限元軟件模擬了熱泵單一制冷季運(yùn)行條件下,地下水與巖土體的熱平衡過程。結(jié)果發(fā)現(xiàn):地下水滲流速度隨水井抽灌狀態(tài)改變呈間歇變化,不同含水地層水井熱影響半徑及熱運(yùn)移速度不同,水源熱泵換熱區(qū)地溫在一個(gè)制冷季和一個(gè)間歇期后未恢復(fù)至初始狀態(tài),出現(xiàn)地層熱堆積和抽灌井區(qū)熱貫通現(xiàn)象。

    關(guān)鍵詞:地下水源熱泵;地溫場;數(shù)值模擬;熱堆積

    Abstract: Taking the groundwater source heat pump system experiment field conditions, a coupled 3D numerical model of groundwater seepage and thermal migration has been developed. The finite element software COMSOL Multiphysics was used to simulate thermal equilibrium process of groundwater and rock-soil body during single refrigerating season. The results show that, the seepage field changes intermittently with the changes of groundwater source heat pump daily running state, and thermal diffusion radius as well as migration rate are different in different aquifers. The ground temperatures can not be recovered to the initial state after one year of cooling season and intermission, and produce strata heat stack as well as heat transfixion phenomenon in the pumping and injecting wells place .

    Keywords: groundwater source heat pump;ground thermal field;numerical modeling;heat stack

    地下水地源熱泵系統(tǒng)又稱水源熱泵系統(tǒng),它通過熱泵空調(diào)系統(tǒng)實(shí)現(xiàn)建筑物與地下水之間的熱量交換達(dá)到制冷和供暖的目的。水源熱泵空調(diào)系統(tǒng)的運(yùn)行能效(COP值)一般可達(dá)到3.5~4.0,效能好(劉立才等,2008;趙輝,2008)。然而抽灌井之間的熱貫通效應(yīng)以及井周圍的熱堆積效應(yīng)致使地下水采能區(qū)地溫場發(fā)生熱失衡,影響熱泵系統(tǒng)的運(yùn)行效率和使用壽命(王慧玲等,2011;胡繼華,2009;叢曉春等,2008)。預(yù)測地下水滲流場與溫度場的變化,避免冷熱貫通現(xiàn)象與熱堆積問題,是合理進(jìn)行抽、灌井的布局設(shè)計(jì)的先行工作(周世玲等,2014)。利用數(shù)值模擬手段開展地溫研究,了解地下水流動與地溫場演化之間的關(guān)系是科學(xué)設(shè)計(jì)地下水源熱泵系統(tǒng)的手段之一,具有重要的應(yīng)用價(jià)值(張遠(yuǎn)東等,2006;王慧玲,2010)。當(dāng)前有多井、雙井以及同井回灌數(shù)值模擬研究(王家樂等,2019;余期沖等,2015;李鳳昱等,2020),模擬抽灌方式包括場地內(nèi)雙季節(jié)工況和季節(jié)交替工況。雙工況條件下地溫場基本均衡,季節(jié)交替工況下地溫變化優(yōu)于正常抽灌模式,然而場地條件下夏季單工況間歇運(yùn)行時(shí)地溫變化研究尚未開展。本文通過有限元數(shù)值模擬方法對實(shí)驗(yàn)場地承壓抽灌對井溫度場開展計(jì)算分析,探究水源熱泵制冷季間歇運(yùn)行時(shí)含水層熱量運(yùn)移特點(diǎn)。

    1 實(shí)驗(yàn)場地概況

    地下水換熱系統(tǒng)實(shí)驗(yàn)場地位于北京平原區(qū),場地占地面積18000 m2。第四系為溫榆河沖積物,厚約200 m。地層巖性主要為砂、砂礫石與黏性土互層,以河流相沉積為主(圖1)。

    第四系孔隙水為弱承壓水,含水層厚度較大,補(bǔ)給條件好,富水性較好。根據(jù)地層測井結(jié)果,地層可概略為6層。19~36 m和40~58 m為砂卵石層和粗砂層、58~75 m為細(xì)砂層,這3個(gè)層位均為含水層,其余層位巖性為黏土—粉土層,為弱透水層。潛水層主要接受大氣降水入滲和地下水側(cè)向徑流補(bǔ)給,潛水水位8.5 m。承壓水的補(bǔ)給方式主要為地下水側(cè)向徑流和潛水含水層越流補(bǔ)給。

    2 地下水地源熱泵系統(tǒng)運(yùn)行概況

    在場地內(nèi)鉆鑿兩眼水井,井深均為100 m,兩口井相距70 m,通過開展試驗(yàn)場地抽水試驗(yàn)發(fā)現(xiàn),自然狀態(tài)下,抽水井處水位低,回灌井處水位高,地下水自回灌井流入抽水井,地下水流向與抽水井—回灌井連線方向呈45°斜交(圖2)。試驗(yàn)場地地下水背景溫度在15.0 ℃左右,整體比較穩(wěn)定。地下水水質(zhì)達(dá)到飲用水標(biāo)準(zhǔn),滿足地下水回灌的要求。為了進(jìn)一步觀察熱泵系統(tǒng)換熱對地層溫度的影響,在場地內(nèi)距離抽水井上游方向50 m處設(shè)置溫度采集點(diǎn),觀察地溫變化。

    實(shí)驗(yàn)場地內(nèi)兩地下水井,即是兩眼抽水試驗(yàn)井,又是地下水地源熱泵系統(tǒng)夏季制冷的備用井。由于場地內(nèi)冬季采用地?zé)崴┡?,夏季室?nèi)制冷采用地下水源熱泵系統(tǒng),水源熱泵系統(tǒng)每個(gè)夏季制冷時(shí)間是100 d,非全天連續(xù)運(yùn)行,為節(jié)能,每運(yùn)行12 h,間歇12 h,每年運(yùn)行單一個(gè)制冷季,共100 d,間歇周期為265 d。地下水抽取層位為砂卵礫石層和細(xì)砂層。采用重力回灌方式進(jìn)行回灌,回灌層位與抽水層位相同。

    3 數(shù)值模擬模型建立

    3.1 基本假設(shè)

    抽灌水井的模擬中采用以下假設(shè)條件:含水層非均質(zhì)、各向同性,產(chǎn)狀水平,厚度不變;地下水水流為水平徑向流,和回灌井井流連續(xù),水流通過含水層徑向流服從達(dá)西定律;抽取交換完熱量的水進(jìn)入回灌井中,不考慮井損;回水溫度在換熱時(shí)間與熱泵間歇期保持一恒定值;忽略上部弱透水層的貯水性和彈性釋水量。

    3.2 數(shù)學(xué)模型

    地層多孔介質(zhì)為飽和介質(zhì),均質(zhì)且各向同性,局部瞬態(tài)換熱平衡,能量平衡方程為:

    式中:ρ為流體密度,kg/m3;CP為定壓比熱容,J/(kg·K);q為熱通量,W/m2 ;u為地下水流速,m/s;keff為有效導(dǎo)熱系數(shù),W/(m·K);keff=(1-φ)ks+φkl,ks和kl分別為固體和液體的導(dǎo)熱系數(shù),φ為孔隙率,%;Q為熱源(熱匯),W/m3。

    地下水流動遵循達(dá)西定律,計(jì)算公式為:

    式中:?p為孔隙體積,m3;Q為質(zhì)量源(匯);κ為地層滲透率,m2;μ為動力黏度,Pa·s;p為水頭壓力,Pa。其它符號意義同前。

    3.3 幾何模型和分層參數(shù)

    兩眼井相距70 m,根據(jù)本地地下水補(bǔ)徑排系統(tǒng)劃分,井115 m處的側(cè)向邊界設(shè)為定水頭邊界,因此,模型幾何體長、寬均為300 m,幾何體高與井深相同,設(shè)為100 m。地下水流向垂直于幾何體一邊,與回灌井—抽水井連線方向斜交,如圖3所示。

    各層位對應(yīng)的熱物性、密度、孔隙率如表1所示。不同深度地溫設(shè)置為隨深度變化的實(shí)測值。頂面、底面溫度設(shè)定為恒溫邊界。地下水運(yùn)移中,加入水頭邊界。水頭邊界在模型中為重要的邊界,它影響抽水時(shí)的地層熱運(yùn)移速度,不同含水層水頭根據(jù)場地內(nèi)抽水試驗(yàn)計(jì)算得到。

    3.4 網(wǎng)格剖分

    采用COMSOL Multiphysics有限元軟件模擬水井在抽水、回灌中地下水與巖土體的換熱過程,研究熱-滲耦合作用下的滲流、地溫變化特征。模型網(wǎng)格采用超細(xì)化尺寸的自由四面體網(wǎng)格剖分,網(wǎng)格單元共274297個(gè)。計(jì)算時(shí)間為1個(gè)完整年度(365 d),計(jì)算步長為1 d,其中0~100 d為抽水-回灌時(shí)間,模擬中加入熱泵啟停時(shí)間,白天啟動、夜晚關(guān)閉。100~365 d為停止抽灌、自然恢復(fù)時(shí)間。抽水量、回灌量均為60 m3/h時(shí),回灌時(shí)間與抽水時(shí)間相同。實(shí)驗(yàn)場地夏季工況下多次測量回灌水溫均值約為20 ℃,因而模型回灌水溫度設(shè)置為20 ℃。

    4 模擬結(jié)果分析討論

    4.1 滲流變化特征

    隨著抽水回灌時(shí)間的變化,地層壓力呈現(xiàn)一定的規(guī)律性。由于抽水和回灌時(shí)間在前100 d內(nèi)存在關(guān)停,即壓力有間歇性恢復(fù),故地層壓力在抽水回灌時(shí)壓力變換較大,壓力在抽水井處降低,而回灌井處壓力升高,如圖4左圖所示。抽水回灌暫停時(shí),壓力略微恢復(fù),如圖4右圖所示,由于殘余壓力的存在,壓力未恢復(fù)至初始值。隨著抽水回灌的結(jié)束,地層壓力恢復(fù)至自然狀態(tài)的初始值。

    在滲透率一定時(shí),根據(jù)達(dá)西定律,水頭壓力與流速成正比。抽水回灌層位達(dá)西速度與壓力變化規(guī)律一致。抽水、回灌同時(shí)進(jìn)行時(shí),地下水流場形成“一源一匯”,達(dá)西速度場變化較快。抽水、回灌停止后,地下水流動速度減緩,達(dá)西速度場逐漸恢復(fù)至自然流動分布。在抽—灌中,隨著地層壓力的升高和減小,流速呈現(xiàn)交替增大—變小的規(guī)律性變化,與熱泵間歇運(yùn)行時(shí)間一致,如圖5的第0 d、第20 d、第25 d所示。在100 d停止抽水后,流速恢復(fù)至自然狀態(tài),流線變平直。

    4.2 地溫變化特征

    場地內(nèi)地溫采集點(diǎn)距離抽水井50 m,距離回灌井近90 m。通過采集抽水試驗(yàn)期間場地內(nèi)的地層溫度,獲得了地層30 m和75 m深度處的地溫?cái)?shù)據(jù)(圖6)。從圖中可見,地層實(shí)際采集深度30 m處由于受到氣溫影響,溫度約有0.4 ℃下降,70 m深度處地層溫度幾乎沒有波動。從模擬計(jì)算數(shù)據(jù)來看,未考慮氣溫影響下,30 m深度處溫度變化較小,在熱泵間歇期恢復(fù)時(shí),地溫略微下降,與實(shí)測溫度相差0.5 ℃。70 m處溫度幾乎沒有變化,與地層實(shí)際溫度一致。基于以上溫度對比結(jié)果,認(rèn)為本次模擬計(jì)算結(jié)果有效。

    分別截取位于30 m、70 m處3個(gè)深度地溫切片,研究在第0 d、第100 d、第200 d、第365 d的6 h間點(diǎn)的地溫變化,見圖7和圖8。

    在地層30 m處,第0 d時(shí),抽水、回灌初始階段,熱輸入開始,尚未形成熱團(tuán)。隨著抽灌時(shí)間增長,回灌井熱擴(kuò)散范圍增大,向抽水井發(fā)生偏移。在回灌至第100 d時(shí),回灌井周圍地溫升高至20 ℃,熱云圖匯向抽水井?;毓嗑?8~20 ℃的熱場圖呈橢圓形,熱半徑約為30 m。受回灌水傳熱的影響,位于抽水井上游,靠近回灌井一側(cè)的地溫約升高2 ℃。

    第200 d時(shí),熱場向下游隨地下水流動而偏移,回灌井下游熱半徑大于上游方向,地溫略微恢復(fù)。由于自然流場的原因,熱云團(tuán)持續(xù)向下游偏移,恢復(fù)程度升高。第365 d,溫度持續(xù)恢復(fù),但并未恢復(fù)至初始值,產(chǎn)生熱堆積現(xiàn)象。由于回灌井附近的熱擴(kuò)散路徑隨地下水徑流方向延展,與抽水井方位有一定角度偏差,所以抽水井處地溫在形成熱貫通后,未呈現(xiàn)較大的溫度波動。

    70 m深度處,地層巖性為細(xì)砂,回灌水量較小,注入后熱半徑較小,由于地下水流速相比粗砂層慢,所以熱偏移也較小,在100 d內(nèi)的抽灌中,未形成熱突破?;毓嘟Y(jié)束后,熱半徑向下游偏移,且地溫逐漸恢復(fù)。由于熱擴(kuò)散較慢、偏移速度也較慢,在第365 d時(shí),回灌井處地溫未恢復(fù)至初始值,產(chǎn)生熱堆積。

    5 結(jié)論

    (1)地下水地源熱泵熱-滲耦合數(shù)值計(jì)算結(jié)果表明,存在地下水天然流動時(shí),含水層熱運(yùn)移方向主要與地下水徑流方向一致;地層壓力和地下水滲流隨著間歇運(yùn)行出現(xiàn)強(qiáng)—弱相間的變化規(guī)律。

    (2)由于抽灌井距離較近,在換熱一段時(shí)間之后發(fā)生熱貫通現(xiàn)象,抽水井處溫度略微升高,回灌井下游方向在經(jīng)過一個(gè)夏季制冷和一個(gè)間歇期之后,地溫場并未恢復(fù),發(fā)生熱堆積。

    (3)在地下水地源熱泵應(yīng)用中,應(yīng)考慮雙季節(jié)換熱,以利于地層熱負(fù)荷均衡,同時(shí)應(yīng)加大井距,避免熱貫通效應(yīng),提高地下水地源熱泵熱系統(tǒng)換熱效率。

    參考文獻(xiàn):

    叢曉春,楊文斐,2008. 地下承壓含水層水-熱運(yùn)移特性的模擬研究[J]. 太陽能學(xué)報(bào),29(11):1390-1394.

    胡繼華,2009. 地下水源熱泵水力學(xué)機(jī)理及其對地下溫度場影響研究[D]. 長春:吉林大學(xué).

    劉立才,張遠(yuǎn)東,李世君,等,2008. 冷熱負(fù)荷失衡條件下采能區(qū)地溫場的模擬研究[J]. 太陽能學(xué)報(bào),29(9):1072-1077.

    李鳳昱,許天福,封官宏,等,2020. T2Well單井地下水源熱泵水-熱耦合數(shù)值模擬研究[J]. 太陽能學(xué)報(bào),41(4):278-286.

    王慧玲,2010. 區(qū)域流場作用下地下水地源熱泵系統(tǒng)抽灌模式對地溫場的影響[J]. 地下水,32(4):42-44+50.

    王慧玲,王文峰,王峰,等,2011. 地下水地源熱泵系統(tǒng)應(yīng)用對地溫場的影響[J]. 水文地質(zhì)工程地質(zhì),38(3):134-138.

    王家樂,賈寶杰,王志剛,2019. 巖溶含水層地下水源熱泵運(yùn)行對地下水溫度場的影響[J]. 水電能源科學(xué),37(2):55-58.

    余期沖,祝曉彬,吳吉春,等,2015. 場地地下水源熱泵系統(tǒng)抽灌實(shí)驗(yàn)及數(shù)值模擬研究[J]. 地下水,37(3):24-27+62.

    趙輝,2008. 地下水源熱泵工程對淺部地溫場的影響研究[D]. 阜新:遼寧工程技術(shù)大學(xué).

    張遠(yuǎn)東,魏加華,王光謙,2006. 區(qū)域流場對含水層采能區(qū)地溫場的影響[J]. 清華大學(xué)學(xué)報(bào)(自然科學(xué)版),46(9):1518-1521.

    周世玲,駱祖江,于丹丹,等,2014. 地下水源熱泵系統(tǒng)熱平衡預(yù)測三維數(shù)值模擬[J]. 中國煤炭地質(zhì),26(9):34-39.

    猜你喜歡
    含水層熱泵源熱泵
    地源熱泵新型技術(shù)的應(yīng)用
    應(yīng)用空氣能熱泵干燥設(shè)備提高農(nóng)產(chǎn)品加工的效率和品質(zhì)
    熱泵——迸發(fā)生活新動力
    地源熱泵在建筑空調(diào)系統(tǒng)設(shè)計(jì)中的應(yīng)用分析
    基于地層及水化學(xué)特征分析采煤對地下水環(huán)境的影響
    掘進(jìn)巷道遇含水層施工方法研究
    國家質(zhì)檢總局抽查:6批次熱泵熱水機(jī)(器)不合格
    寧夏某礦地下水補(bǔ)給來源分析
    中高溫?zé)岜脵C(jī)組研究應(yīng)用
    一本大道久久a久久精品| 99国产精品99久久久久| 国产精品一区二区精品视频观看| 黄色 视频免费看| 久久久国产欧美日韩av| 成人亚洲精品一区在线观看| 久久精品人人爽人人爽视色| 在线看a的网站| 国产成人免费观看mmmm| 国产精品偷伦视频观看了| 日本vs欧美在线观看视频| 在线播放国产精品三级| 亚洲一区中文字幕在线| 亚洲伊人色综图| 欧美黑人欧美精品刺激| 色老头精品视频在线观看| 我要看黄色一级片免费的| 美女视频免费永久观看网站| 国产一卡二卡三卡精品| 免费人妻精品一区二区三区视频| 国产aⅴ精品一区二区三区波| 一本大道久久a久久精品| 丝袜在线中文字幕| 99在线人妻在线中文字幕 | www.熟女人妻精品国产| 国产精品1区2区在线观看. | 啦啦啦 在线观看视频| 国产精品自产拍在线观看55亚洲 | 一级毛片电影观看| 国产精品久久久久久精品电影小说| 丝瓜视频免费看黄片| 亚洲成av片中文字幕在线观看| 一区二区三区精品91| 国产激情久久老熟女| 丰满少妇做爰视频| 少妇粗大呻吟视频| 99在线人妻在线中文字幕 | 丰满迷人的少妇在线观看| 亚洲一区中文字幕在线| 午夜福利一区二区在线看| 免费观看人在逋| 午夜老司机福利片| 欧美日韩国产mv在线观看视频| 欧美亚洲日本最大视频资源| 亚洲精品国产精品久久久不卡| 免费人妻精品一区二区三区视频| 黄色怎么调成土黄色| 一区二区三区激情视频| 80岁老熟妇乱子伦牲交| 欧美亚洲 丝袜 人妻 在线| 男女高潮啪啪啪动态图| 亚洲欧美激情在线| 国产精品一区二区免费欧美| 9热在线视频观看99| 国产亚洲精品第一综合不卡| 精品国产乱码久久久久久男人| 操出白浆在线播放| netflix在线观看网站| 在线亚洲精品国产二区图片欧美| 久久精品人人爽人人爽视色| 一区二区av电影网| 99九九在线精品视频| 日韩欧美国产一区二区入口| 男女之事视频高清在线观看| 桃花免费在线播放| 久久久久网色| 热99久久久久精品小说推荐| 欧美乱妇无乱码| 人人澡人人妻人| 精品一区二区三区视频在线观看免费 | 久久国产精品大桥未久av| 美女主播在线视频| 黄色毛片三级朝国网站| 久久精品国产a三级三级三级| 大香蕉久久网| 亚洲国产av影院在线观看| 91精品国产国语对白视频| 亚洲精品在线美女| 国产精品一区二区在线观看99| 99精品久久久久人妻精品| av天堂久久9| 少妇 在线观看| 999久久久国产精品视频| 好男人电影高清在线观看| 日韩成人在线观看一区二区三区| 五月开心婷婷网| 国产精品成人在线| 久久精品国产99精品国产亚洲性色 | 久久人人97超碰香蕉20202| 国产精品.久久久| 午夜老司机福利片| av不卡在线播放| 国产男靠女视频免费网站| 亚洲第一av免费看| 久久久水蜜桃国产精品网| 久久中文字幕一级| 最黄视频免费看| 国产精品一区二区在线不卡| 欧美人与性动交α欧美软件| 一区二区三区激情视频| 国产野战对白在线观看| 大型av网站在线播放| www.精华液| 国产福利在线免费观看视频| 国产成人啪精品午夜网站| 九色亚洲精品在线播放| 欧美老熟妇乱子伦牲交| 亚洲熟女精品中文字幕| 精品卡一卡二卡四卡免费| 乱人伦中国视频| 亚洲中文av在线| 无人区码免费观看不卡 | 精品福利永久在线观看| av天堂久久9| 飞空精品影院首页| 老司机影院毛片| tocl精华| 亚洲欧洲精品一区二区精品久久久| 香蕉国产在线看| 亚洲精品av麻豆狂野| 欧美+亚洲+日韩+国产| 天天添夜夜摸| 午夜成年电影在线免费观看| 亚洲五月色婷婷综合| 国产精品免费一区二区三区在线 | 久久久久国内视频| 精品高清国产在线一区| 成年动漫av网址| 99riav亚洲国产免费| 天堂8中文在线网| 动漫黄色视频在线观看| 国产日韩欧美在线精品| 午夜福利视频精品| 在线播放国产精品三级| 久久久久网色| 精品亚洲成a人片在线观看| 精品乱码久久久久久99久播| 精品免费久久久久久久清纯 | av超薄肉色丝袜交足视频| 麻豆国产av国片精品| 午夜久久久在线观看| 精品久久久精品久久久| 超碰97精品在线观看| 亚洲第一av免费看| 国产精品自产拍在线观看55亚洲 | 十分钟在线观看高清视频www| 真人做人爱边吃奶动态| 亚洲av日韩在线播放| 精品一区二区三卡| 欧美日韩成人在线一区二区| 高清在线国产一区| 老熟妇乱子伦视频在线观看| 日本av手机在线免费观看| 成年人黄色毛片网站| 国产精品免费大片| 午夜福利影视在线免费观看| 丁香六月天网| 亚洲av日韩精品久久久久久密| 亚洲熟女精品中文字幕| 人人妻人人添人人爽欧美一区卜| 精品一区二区三区视频在线观看免费 | 国产精品香港三级国产av潘金莲| 五月开心婷婷网| 一级,二级,三级黄色视频| 一级毛片电影观看| 欧美日韩精品网址| 国产xxxxx性猛交| 99国产综合亚洲精品| 日本wwww免费看| 成人影院久久| 一区二区三区激情视频| www.999成人在线观看| 十八禁人妻一区二区| 极品人妻少妇av视频| 亚洲欧美日韩高清在线视频 | 久久久久久久久久久久大奶| 久久久久国产一级毛片高清牌| 国产成人精品久久二区二区免费| 国产一区二区三区视频了| 国产一卡二卡三卡精品| 色94色欧美一区二区| 啦啦啦 在线观看视频| 午夜福利视频在线观看免费| 精品一品国产午夜福利视频| 欧美黑人欧美精品刺激| 99精品久久久久人妻精品| 亚洲情色 制服丝袜| 欧美成人午夜精品| 一二三四社区在线视频社区8| 国产福利在线免费观看视频| 国产高清videossex| 黄网站色视频无遮挡免费观看| 国产成人系列免费观看| 国产精品久久久久久人妻精品电影 | 亚洲av电影在线进入| av片东京热男人的天堂| 老司机在亚洲福利影院| 18禁国产床啪视频网站| 两性夫妻黄色片| 午夜福利视频精品| 国产日韩欧美在线精品| 视频在线观看一区二区三区| 亚洲性夜色夜夜综合| 怎么达到女性高潮| 少妇的丰满在线观看| 亚洲精品久久成人aⅴ小说| 丝袜美足系列| 欧美日韩福利视频一区二区| 一进一出抽搐动态| 国产av国产精品国产| 91成人精品电影| 18禁美女被吸乳视频| 婷婷丁香在线五月| 欧美精品啪啪一区二区三区| 亚洲视频免费观看视频| 国产国语露脸激情在线看| 精品卡一卡二卡四卡免费| 19禁男女啪啪无遮挡网站| 国产三级黄色录像| 欧美一级毛片孕妇| 久久久水蜜桃国产精品网| 欧美日韩亚洲高清精品| 精品欧美一区二区三区在线| 成人手机av| 久久 成人 亚洲| 色在线成人网| 一区二区三区精品91| 国产1区2区3区精品| 国产精品麻豆人妻色哟哟久久| 纵有疾风起免费观看全集完整版| 最近最新中文字幕大全电影3 | 不卡一级毛片| 黄频高清免费视频| 亚洲欧洲精品一区二区精品久久久| 午夜视频精品福利| 人成视频在线观看免费观看| 亚洲色图综合在线观看| 国产淫语在线视频| 18禁黄网站禁片午夜丰满| 免费日韩欧美在线观看| 国产在线视频一区二区| av网站在线播放免费| 动漫黄色视频在线观看| 在线看a的网站| 久久青草综合色| 国产一区有黄有色的免费视频| 高清毛片免费观看视频网站 | 精品国产乱码久久久久久小说| 亚洲专区字幕在线| 十八禁网站免费在线| 日韩中文字幕视频在线看片| 黄色 视频免费看| 国产av精品麻豆| 丰满人妻熟妇乱又伦精品不卡| 婷婷成人精品国产| 久久久久久久精品吃奶| 啦啦啦视频在线资源免费观看| 国产日韩欧美在线精品| 一进一出抽搐动态| 黑人巨大精品欧美一区二区mp4| 成人免费观看视频高清| 正在播放国产对白刺激| 日韩成人在线观看一区二区三区| 在线观看免费日韩欧美大片| 国产精品熟女久久久久浪| 高清视频免费观看一区二区| 美女扒开内裤让男人捅视频| 丝瓜视频免费看黄片| av超薄肉色丝袜交足视频| 丁香六月欧美| 1024香蕉在线观看| 中文字幕色久视频| 中文字幕制服av| 国产精品 国内视频| 国产成人免费无遮挡视频| 日韩 欧美 亚洲 中文字幕| 久久久国产成人免费| 色婷婷久久久亚洲欧美| 日韩熟女老妇一区二区性免费视频| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久久亚洲精品国产蜜桃av| 五月天丁香电影| 叶爱在线成人免费视频播放| 变态另类成人亚洲欧美熟女 | 在线十欧美十亚洲十日本专区| 久久av网站| 国产成人欧美在线观看 | 美女高潮喷水抽搐中文字幕| 美女扒开内裤让男人捅视频| 午夜免费成人在线视频| 亚洲色图 男人天堂 中文字幕| 国产精品国产高清国产av | 欧美日韩福利视频一区二区| 日韩一区二区三区影片| 亚洲欧美一区二区三区黑人| 成人精品一区二区免费| 人妻一区二区av| 老司机午夜福利在线观看视频 | 亚洲av国产av综合av卡| 免费在线观看视频国产中文字幕亚洲| videosex国产| 黑人欧美特级aaaaaa片| avwww免费| 亚洲国产av新网站| 超碰97精品在线观看| 久久99一区二区三区| 男男h啪啪无遮挡| kizo精华| 精品午夜福利视频在线观看一区 | 国产黄色免费在线视频| 成人永久免费在线观看视频 | 咕卡用的链子| 久久久水蜜桃国产精品网| 高清在线国产一区| 日本欧美视频一区| 亚洲伊人久久精品综合| 悠悠久久av| 一级,二级,三级黄色视频| 欧美乱码精品一区二区三区| 国产欧美亚洲国产| 久久国产精品人妻蜜桃| 亚洲精品久久成人aⅴ小说| 桃花免费在线播放| 国产欧美日韩一区二区三区在线| 免费人妻精品一区二区三区视频| 国产精品1区2区在线观看. | 久久久水蜜桃国产精品网| 色综合婷婷激情| 日本一区二区免费在线视频| 久久精品亚洲av国产电影网| 又大又爽又粗| 视频区欧美日本亚洲| 国产单亲对白刺激| 好男人电影高清在线观看| 高清毛片免费观看视频网站 | 桃花免费在线播放| 日韩制服丝袜自拍偷拍| 久久久国产成人免费| 亚洲三区欧美一区| 亚洲 国产 在线| 国产av国产精品国产| 中文字幕制服av| 亚洲午夜精品一区,二区,三区| 国产在线免费精品| 色视频在线一区二区三区| 久久亚洲精品不卡| 国产亚洲欧美精品永久| 日韩成人在线观看一区二区三区| 一进一出抽搐动态| 国产一区二区三区视频了| 一级片免费观看大全| 一区二区三区乱码不卡18| 啦啦啦 在线观看视频| 亚洲国产欧美在线一区| 久久婷婷成人综合色麻豆| 午夜福利在线观看吧| 黄色毛片三级朝国网站| 他把我摸到了高潮在线观看 | 新久久久久国产一级毛片| 欧美日韩视频精品一区| 免费日韩欧美在线观看| 女性生殖器流出的白浆| 精品少妇久久久久久888优播| 亚洲视频免费观看视频| 国产高清国产精品国产三级| 日韩欧美一区二区三区在线观看 | svipshipincom国产片| 国产亚洲精品第一综合不卡| 日韩一卡2卡3卡4卡2021年| 黄色 视频免费看| 欧美亚洲 丝袜 人妻 在线| 女性被躁到高潮视频| 91成年电影在线观看| 亚洲色图综合在线观看| 国产成人精品在线电影| 亚洲av片天天在线观看| 国产免费现黄频在线看| 亚洲中文日韩欧美视频| 成在线人永久免费视频| 一本久久精品| 少妇裸体淫交视频免费看高清 | 欧美日韩视频精品一区| 日韩免费av在线播放| 一进一出抽搐动态| 国产野战对白在线观看| 少妇的丰满在线观看| 国产xxxxx性猛交| 欧美日韩亚洲国产一区二区在线观看 | 久久国产精品影院| 嫁个100分男人电影在线观看| 亚洲欧美日韩另类电影网站| 国产精品偷伦视频观看了| 亚洲中文字幕日韩| av福利片在线| 99久久人妻综合| 国产成人免费观看mmmm| 精品一区二区三卡| av天堂久久9| 欧美激情 高清一区二区三区| 我的亚洲天堂| 国产精品.久久久| 男男h啪啪无遮挡| 国产av精品麻豆| 超碰成人久久| 午夜福利在线观看吧| 亚洲成人免费电影在线观看| 欧美日韩视频精品一区| 日韩欧美国产一区二区入口| 日韩大片免费观看网站| 国产真人三级小视频在线观看| 亚洲avbb在线观看| 国产伦理片在线播放av一区| 免费在线观看日本一区| 国产免费现黄频在线看| 不卡av一区二区三区| 满18在线观看网站| 免费日韩欧美在线观看| 五月开心婷婷网| 欧美日韩亚洲国产一区二区在线观看 | 在线观看免费视频日本深夜| 两个人免费观看高清视频| 久久香蕉激情| 香蕉久久夜色| tocl精华| 国产精品熟女久久久久浪| 亚洲,欧美精品.| 最黄视频免费看| 国产成+人综合+亚洲专区| 欧美乱码精品一区二区三区| 精品少妇内射三级| 99国产极品粉嫩在线观看| 亚洲熟女精品中文字幕| av天堂久久9| 国产精品98久久久久久宅男小说| 自线自在国产av| 在线观看一区二区三区激情| 男女之事视频高清在线观看| 午夜成年电影在线免费观看| 亚洲专区中文字幕在线| 亚洲av国产av综合av卡| 一区二区三区激情视频| 国产成人av教育| 妹子高潮喷水视频| 亚洲av国产av综合av卡| 国产一区二区三区视频了| 99久久国产精品久久久| 国产一区有黄有色的免费视频| 法律面前人人平等表现在哪些方面| 在线永久观看黄色视频| 亚洲精品av麻豆狂野| 国产精品一区二区在线观看99| 欧美 亚洲 国产 日韩一| 亚洲欧美一区二区三区黑人| 国产精品免费视频内射| 精品熟女少妇八av免费久了| 国产亚洲精品第一综合不卡| 丰满人妻熟妇乱又伦精品不卡| 日本精品一区二区三区蜜桃| 中亚洲国语对白在线视频| 香蕉久久夜色| 亚洲一区二区三区欧美精品| 欧美亚洲日本最大视频资源| bbb黄色大片| 久久久国产欧美日韩av| 国产人伦9x9x在线观看| 19禁男女啪啪无遮挡网站| 一本久久精品| 又大又爽又粗| 黑人巨大精品欧美一区二区蜜桃| 高清av免费在线| 国产精品香港三级国产av潘金莲| 欧美久久黑人一区二区| 久久国产精品大桥未久av| 欧美精品啪啪一区二区三区| 久久九九热精品免费| 韩国精品一区二区三区| 老司机福利观看| 亚洲视频免费观看视频| 色综合婷婷激情| 纵有疾风起免费观看全集完整版| 蜜桃国产av成人99| 夜夜夜夜夜久久久久| 国产精品香港三级国产av潘金莲| 欧美久久黑人一区二区| 99国产综合亚洲精品| 欧美人与性动交α欧美软件| 亚洲国产欧美一区二区综合| 亚洲午夜精品一区,二区,三区| 满18在线观看网站| 老司机午夜十八禁免费视频| 国产成人精品久久二区二区免费| 欧美黄色片欧美黄色片| av一本久久久久| 欧美在线黄色| 一本—道久久a久久精品蜜桃钙片| 欧美 日韩 精品 国产| 一级,二级,三级黄色视频| 激情在线观看视频在线高清 | 成年人免费黄色播放视频| 免费在线观看日本一区| 日韩欧美免费精品| 成人特级黄色片久久久久久久 | 国产精品一区二区精品视频观看| 亚洲一卡2卡3卡4卡5卡精品中文| 国产又爽黄色视频| 精品少妇久久久久久888优播| 中文字幕另类日韩欧美亚洲嫩草| 建设人人有责人人尽责人人享有的| 精品久久久久久电影网| 在线观看免费视频网站a站| 欧美久久黑人一区二区| 久久人妻福利社区极品人妻图片| 女同久久另类99精品国产91| 视频区图区小说| 午夜福利乱码中文字幕| 少妇猛男粗大的猛烈进出视频| 69av精品久久久久久 | 日韩中文字幕视频在线看片| 久久中文看片网| 王馨瑶露胸无遮挡在线观看| 国产成人精品久久二区二区91| 国产视频一区二区在线看| 精品久久久久久电影网| 亚洲视频免费观看视频| 国产成人免费无遮挡视频| 欧美在线一区亚洲| 久久久精品区二区三区| 日韩一卡2卡3卡4卡2021年| 国产精品久久久人人做人人爽| 亚洲中文日韩欧美视频| 中国美女看黄片| av天堂在线播放| 美女扒开内裤让男人捅视频| 日韩欧美三级三区| 美国免费a级毛片| 久久精品亚洲av国产电影网| 丁香六月天网| 最近最新中文字幕大全电影3 | 肉色欧美久久久久久久蜜桃| 妹子高潮喷水视频| 亚洲国产欧美在线一区| 欧美老熟妇乱子伦牲交| 啦啦啦视频在线资源免费观看| 少妇裸体淫交视频免费看高清 | 多毛熟女@视频| 久久久国产精品麻豆| 手机成人av网站| 中文字幕人妻熟女乱码| 啦啦啦免费观看视频1| 亚洲色图综合在线观看| 国产精品 国内视频| 老熟女久久久| 久久国产精品人妻蜜桃| 丝袜美腿诱惑在线| 日本黄色视频三级网站网址 | 757午夜福利合集在线观看| 成年人黄色毛片网站| tube8黄色片| 一级毛片女人18水好多| 精品免费久久久久久久清纯 | 国产欧美日韩精品亚洲av| 亚洲自偷自拍图片 自拍| 欧美成狂野欧美在线观看| 国产精品麻豆人妻色哟哟久久| 一级毛片女人18水好多| 天天操日日干夜夜撸| 午夜精品国产一区二区电影| 亚洲国产中文字幕在线视频| 亚洲av电影在线进入| 亚洲中文字幕日韩| 婷婷丁香在线五月| 欧美日韩中文字幕国产精品一区二区三区 | 日韩欧美国产一区二区入口| 制服人妻中文乱码| 国内毛片毛片毛片毛片毛片| 热re99久久精品国产66热6| 国产精品亚洲av一区麻豆| 久久久久精品人妻al黑| 国产精品98久久久久久宅男小说| 亚洲成a人片在线一区二区| 欧美黑人欧美精品刺激| 国产精品久久久久久精品电影小说| 精品国产一区二区久久| 国产一区有黄有色的免费视频| 美女扒开内裤让男人捅视频| 亚洲国产成人一精品久久久| 午夜福利一区二区在线看| 久久久久视频综合| avwww免费| 最新美女视频免费是黄的| 天堂动漫精品| 中文字幕av电影在线播放| 99精国产麻豆久久婷婷| 久久这里只有精品19| 欧美日韩黄片免| 国产淫语在线视频| 色综合欧美亚洲国产小说| 欧美精品高潮呻吟av久久| 大片电影免费在线观看免费| 精品人妻在线不人妻| 怎么达到女性高潮| 老司机影院毛片| 手机成人av网站| 欧美日韩视频精品一区| 亚洲少妇的诱惑av| 91麻豆av在线| kizo精华| 国产成人精品久久二区二区免费| 亚洲av第一区精品v没综合| 黄色视频在线播放观看不卡| 一进一出抽搐动态| 亚洲av美国av| 国产精品98久久久久久宅男小说| 首页视频小说图片口味搜索| 国产一区二区三区综合在线观看|