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

    熱儲(chǔ)上下巖層熱補(bǔ)償作用對(duì)CO2羽流地?zé)嵯到y(tǒng)性能的影響

    2017-12-22 05:36:44李靜巖劉中良周宇李艷霞
    化工學(xué)報(bào) 2017年12期
    關(guān)鍵詞:生產(chǎn)井開(kāi)采量基巖

    李靜巖,劉中良,周宇,李艷霞

    (北京工業(yè)大學(xué)環(huán)境與能源工程學(xué)院,傳熱強(qiáng)化與過(guò)程節(jié)能教育部重點(diǎn)實(shí)驗(yàn)室,北京 100124)

    熱儲(chǔ)上下巖層熱補(bǔ)償作用對(duì)CO2羽流地?zé)嵯到y(tǒng)性能的影響

    李靜巖,劉中良,周宇,李艷霞

    (北京工業(yè)大學(xué)環(huán)境與能源工程學(xué)院,傳熱強(qiáng)化與過(guò)程節(jié)能教育部重點(diǎn)實(shí)驗(yàn)室,北京 100124)

    二氧化碳羽流地?zé)嵯到y(tǒng)(CPGS)能夠在直接開(kāi)采地?zé)岬耐瑫r(shí)實(shí)現(xiàn) CO2的地質(zhì)封存,熱儲(chǔ)上側(cè)和下側(cè)巖層的熱補(bǔ)償作用是影響 CPGS性能的重要因素。建立了三維蓋巖-熱儲(chǔ)-基巖整體模型,研究了熱補(bǔ)償作用對(duì) CPGS采熱性能及熱儲(chǔ)內(nèi)巖石和流體溫度演化過(guò)程的影響,為優(yōu)化CPGS設(shè)計(jì)、提升CPGS的系統(tǒng)經(jīng)濟(jì)性和工程應(yīng)用提供了理論參考。研究結(jié)果表明:熱補(bǔ)償作用減小了生產(chǎn)流體在豎直方向上的溫度差異及系統(tǒng)運(yùn)行后期生產(chǎn)流體溫度下降的速率,有效延長(zhǎng)系統(tǒng)運(yùn)行壽命,使系統(tǒng)獲得更好的采熱性能和更穩(wěn)定的產(chǎn)能輸出;考慮熱補(bǔ)償時(shí)系統(tǒng)總熱開(kāi)采量明顯提高;與蓋巖相比,基巖的熱補(bǔ)償作用更強(qiáng)。

    巖層;熱儲(chǔ);熱補(bǔ)償;二氧化碳;地?zé)嵯到y(tǒng);多孔介質(zhì);數(shù)值模擬

    引 言

    當(dāng)今社會(huì)主要依賴于傳統(tǒng)的化石能源,而化石燃料的大量使用,造成了以二氧化碳(CO2)為主的溫室氣體大量排放,引發(fā)了全球變暖等諸多環(huán)境問(wèn)題[1-2]。二氧化碳地質(zhì)封存是一種有效減少 CO2排放、緩解溫室效應(yīng)的技術(shù)途徑[3-5],然而單一的CO2地質(zhì)封存成本過(guò)高,難以實(shí)現(xiàn)大規(guī)模工程應(yīng)用[6]。針對(duì)這一問(wèn)題,Brown[7]于2000年首次指出CO2具有良好的流動(dòng)特性和傳熱特性,可以替代水作為開(kāi)采地?zé)岬膫鳠峁べ|(zhì),即基于CO2的增強(qiáng)型地?zé)嵯到y(tǒng)(enhanced geothermal systems powered by CO2,CO2-EGS)進(jìn)行地?zé)衢_(kāi)采,系統(tǒng)運(yùn)行過(guò)程中CO2在地下的流量損失可以間接達(dá)到地質(zhì)封存的目的。隨后 Pruess[8-11]進(jìn)一步驗(yàn)證了 CO2用于地?zé)衢_(kāi)采的可行性,并且圍繞CO2-EGS的采熱能力展開(kāi)了研究。但是CO2-EGS存在很大的局限性。首先,人工熱儲(chǔ)層的構(gòu)造需要較高的壓裂費(fèi)用,從而使系統(tǒng)整體的經(jīng)濟(jì)性較低;此外,“人工壓裂”技術(shù)有可能破壞地質(zhì)構(gòu)造,進(jìn)而引起地震活動(dòng)[12]。因此,CO2-EGS的發(fā)展受到了極大的經(jīng)濟(jì)與政治約束[13]。

    針對(duì)CO2-EGS存在的諸多問(wèn)題,Randolph等[14-15]提出了二氧化碳羽流地?zé)嵯到y(tǒng)(CO2plume geothermal system,CPGS)這一全新的概念。與CO2-EGS的目標(biāo)儲(chǔ)層不同,CPGS使用的是自然存在的地質(zhì)儲(chǔ)層,這種天然的儲(chǔ)層具有一定的滲透率和孔隙率,無(wú)須人工壓裂,可以直接用于地?zé)衢_(kāi)采,從而提高了系統(tǒng)的經(jīng)濟(jì)效益;同時(shí)其較大的儲(chǔ)層體積使之具有比CO2-EGS更大的采熱潛力和CO2地質(zhì)封存潛力[16-18]。

    CPGS的目標(biāo)儲(chǔ)層普遍位于地層深處,上、下側(cè)分別被低滲透性的蓋巖與基巖包圍。隨著系統(tǒng)的運(yùn)行,熱儲(chǔ)在與冷流體的換熱過(guò)程中逐漸被冷卻,造成與蓋巖、基巖的溫度差越來(lái)越大,蓋巖和基巖均將對(duì)熱儲(chǔ)產(chǎn)生一定的熱補(bǔ)償作用。考慮到地?zé)嵯到y(tǒng)較長(zhǎng)的運(yùn)行時(shí)間(一般為幾十年),來(lái)自蓋巖與基巖的熱補(bǔ)償作用會(huì)對(duì)熱儲(chǔ)內(nèi)的溫度變化以及系統(tǒng)的采熱性能造成不可忽略的影響[19-20]。

    以往的研究多只關(guān)注于CPGS的熱儲(chǔ)部分,圍繞熱儲(chǔ)滲透率、比熱容、初始鹽度、初始溫度等物性參數(shù)以及系統(tǒng)布井方式等對(duì)系統(tǒng)采熱性能的影響展開(kāi)研究[21-25],而針對(duì)熱補(bǔ)償作用對(duì)采熱性能影響的研究較少。目前,建立模型時(shí),對(duì)熱補(bǔ)償作用的處理主要有兩種方式:一種是忽略熱儲(chǔ)上下側(cè)巖層的熱補(bǔ)償[21,23],即將熱儲(chǔ)上下側(cè)邊界假設(shè)為絕熱邊界,應(yīng)用此類模型計(jì)算得到的系統(tǒng)采熱量中未包含熱補(bǔ)償作用的貢獻(xiàn)量,與系統(tǒng)實(shí)際采熱量相比具有較大偏差;另一種是通過(guò)將熱儲(chǔ)邊界設(shè)置為一種溫度隨時(shí)間變化的傳熱邊界來(lái)考慮上下巖層的熱補(bǔ)償[25-28],這類邊界模型能夠較為準(zhǔn)確地描述熱補(bǔ)償?shù)淖饔眠^(guò)程,與忽略熱補(bǔ)償?shù)哪P拖啾雀系責(zé)衢_(kāi)采的實(shí)際情況。然而,實(shí)際上,蓋巖和基巖的物性盡管可以與熱儲(chǔ)有明顯不同,甚至可以忽略其質(zhì)量傳遞作用,但從熱量傳遞的角度看,它們與熱儲(chǔ)是不可分割的。所以這種通過(guò)設(shè)定溫度邊界的辦法有很大的隨意性,只能模擬考慮熱補(bǔ)償后熱儲(chǔ)內(nèi)的熱量傳遞過(guò)程和巖石的溫度變化,不能得到蓋巖和基巖中的溫度分布,同時(shí)無(wú)法獲知系統(tǒng)總熱開(kāi)采量中熱補(bǔ)償作用的貢獻(xiàn)量。Jiang等[29-30]通過(guò)在熱儲(chǔ)周圍設(shè)置滲透性可忽略不計(jì)的巖石區(qū)域,考慮了周圍巖石的熱補(bǔ)償作用,與上述兩種對(duì)熱補(bǔ)償作用的處理方式相比,Jiang的模型[31-32]更符合地?zé)衢_(kāi)采時(shí)熱補(bǔ)償過(guò)程的實(shí)際情況,但 Jiang的研究主要針對(duì)EGS,且以CO2作為采熱工質(zhì)時(shí)忽略了熱儲(chǔ)內(nèi)初始水的影響。

    基于以上問(wèn)題,本文首先結(jié)合非等溫的 CO2-水兩相滲流原理建立了基于CPGS的,能夠準(zhǔn)確描述熱補(bǔ)償過(guò)程和地?zé)衢_(kāi)采過(guò)程的蓋巖-熱儲(chǔ)-基巖整體模型,隨后通過(guò)對(duì)比考慮熱補(bǔ)償和忽略熱補(bǔ)償兩種情況,針對(duì)熱補(bǔ)償作用對(duì)CPGS采熱性能及熱儲(chǔ)內(nèi)巖石和流體溫度演化過(guò)程的影響展開(kāi)研究,最后對(duì)兩種情況下系統(tǒng)熱開(kāi)采量的組成進(jìn)行了分析,并提供了計(jì)算蓋巖和基巖熱補(bǔ)償作用對(duì)總熱開(kāi)采量的貢獻(xiàn)量的方法。為后期優(yōu)化CPGS設(shè)計(jì)奠定基礎(chǔ),為提升CPGS的系統(tǒng)經(jīng)濟(jì)性提供理論依據(jù),為CPGS的工程應(yīng)用提供參考。

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

    為合理描述 CPGS運(yùn)行過(guò)程中儲(chǔ)層內(nèi)的流動(dòng)與傳熱過(guò)程,基于以下假設(shè)建立了非等溫條件下的兩相滲流模型:①將熱儲(chǔ)視為均勻性質(zhì)的多孔介質(zhì);②忽略水和CO2間的溶解作用,液相中只有水組分,氣相中只有CO2組分;③不考慮流體與巖石之間的化學(xué)反應(yīng);④忽略巖石熱應(yīng)力引起的熱儲(chǔ)孔隙率及滲透率的變化;⑤模型整體滿足局部熱平衡。

    對(duì)于多孔含水層中的水和CO2兩相,其連續(xù)性方程可表示為

    流動(dòng)滿足達(dá)西方程,即

    輔助方程

    巖層內(nèi)水氣巖三相統(tǒng)一的能量守恒方程為

    式中,φ、S、u、K、κr、μ、p、ρ、g、u、h、C、T、λ、q和Qh分別表示巖石孔隙率、飽和度、速度、巖石固有滲透率、流體相對(duì)滲透率、流體動(dòng)力黏度、流體壓力、密度、重力加速度、流體比內(nèi)能、流體比焓、比熱容、溫度、熱導(dǎo)率、質(zhì)量源和熱源。下標(biāo)w、c、s分別表示水、CO2和固體巖石。毛細(xì)壓力及相對(duì)滲透率采用 Brooks-Corey-Burdine(BCB)模型[33]進(jìn)行計(jì)算,具體計(jì)算公式見(jiàn)表1。

    圖1 CO2和水的物性-壓力-溫度關(guān)系Fig.1 Pressure-temperature relationship of CO2 and H2O physical properties

    表1 毛細(xì)壓力及相對(duì)滲透率模型Table 1 Model of capillary force and relative permeability

    數(shù)值求解采用有限元分析軟件COMSOL,通過(guò)自定義偏微分方程(PDE)模塊建立上述方程。模擬過(guò)程中考慮流體的物性變化,由NIST數(shù)據(jù)庫(kù)獲取CO2和水的物性建立物性數(shù)據(jù)差值表,通過(guò)線性內(nèi)插獲得如圖1所示的物性隨溫度和壓力變化的擬合曲面用于計(jì)算。

    2 算例設(shè)置

    參考前人的經(jīng)驗(yàn)[22,34-35],選用五點(diǎn)式布井方案建立研究區(qū)域(圖2),注入井位于儲(chǔ)層中心位置,生產(chǎn)井位于方形區(qū)域的4個(gè)角上,與注入井的井間距均為500 m。由于五點(diǎn)式布井方式的對(duì)稱性,只考慮方形區(qū)域的1/4部分,建立三維幾何模型(圖3)。模型區(qū)域包括注入井、生產(chǎn)井、多孔介質(zhì)結(jié)構(gòu)的熱儲(chǔ)及滲透率可以忽略不計(jì)的蓋巖和基巖,熱儲(chǔ)中心位于地表以下2500 m的深處,熱儲(chǔ)尺寸為353 m × 353 m × 100 m,蓋巖和基巖厚度為500 m,注入井和生產(chǎn)井地下部分只對(duì)熱儲(chǔ)開(kāi)放,井筒半徑均為0.15 m。

    初始時(shí)刻,熱儲(chǔ)孔隙中被水充滿,流體壓力服從靜水壓強(qiáng),巖石的溫度分布隨深度線性增加,地表溫度293.15 K,地溫梯度為4 K/100 m。不考慮井筒中的流動(dòng)與換熱問(wèn)題,采用固定井底壓力的方式進(jìn)行開(kāi)采,注入井井底和生產(chǎn)井井底與儲(chǔ)層之間的壓力差均取2 MPa,注入井井底CO2的注入溫度為313.15 K。模型通過(guò)改變熱儲(chǔ)頂部和底部邊界(熱儲(chǔ)與蓋巖和基巖的交界面)的邊界條件控制考慮熱補(bǔ)償和忽略熱補(bǔ)償兩種情況,考慮熱補(bǔ)償時(shí)設(shè)為開(kāi)邊界,忽略熱補(bǔ)償時(shí)設(shè)為絕熱邊界,兩種情況下模型的其余邊界均設(shè)為絕熱邊界。熱儲(chǔ)、基巖和蓋巖的物性參數(shù)見(jiàn)表2。

    圖2 五點(diǎn)式布井方案Fig.2 Five-spot well pattern

    圖3 儲(chǔ)層模型幾何尺寸Fig.3 Geometry (including geometric dimensions) of a doublet CPGS

    計(jì)算區(qū)域網(wǎng)格劃分采用自由剖分四面體網(wǎng)格,考慮到注入井和生產(chǎn)井的井筒尺寸較小,重點(diǎn)細(xì)化井筒周圍網(wǎng)格。為保證計(jì)算結(jié)果不受網(wǎng)格數(shù)量的影響,對(duì)網(wǎng)格進(jìn)行加密,直至獲得網(wǎng)格無(wú)關(guān)性計(jì)算結(jié)果,最終網(wǎng)格數(shù)為168664,其中蓋巖和基巖的網(wǎng)格數(shù)均為73005,熱儲(chǔ)的網(wǎng)格數(shù)為22654,四面體網(wǎng)格最小單元尺寸為1.45 m,最大單元尺寸為25 m。

    表2 物理性質(zhì)Table 2 Physical properties

    3 結(jié)果與討論

    3.1 巖層溫度分布

    圖4顯示了考慮熱補(bǔ)償和忽略熱補(bǔ)償兩種情況下不同時(shí)刻儲(chǔ)層內(nèi)的溫度分布。注入低溫CO2后,注入井附近巖石被迅速冷卻,形成低溫區(qū)域,隨CPGS運(yùn)行時(shí)間的推移,CO2逐漸向生產(chǎn)井側(cè)擴(kuò)散,低溫區(qū)域逐漸向生產(chǎn)井側(cè)擴(kuò)展。同時(shí),在豎直方向上,重力作用導(dǎo)致底部流體流動(dòng)速度較快,使得儲(chǔ)層底部巖體被冷卻的速度加快。而圖4(b)、(d)、(f)、(h)顯示出,隨時(shí)間的推移,蓋巖與基巖對(duì)熱儲(chǔ)的熱補(bǔ)償作用逐漸明顯,豎直方向上熱補(bǔ)償作用的強(qiáng)度由熱儲(chǔ)頂部與底部向熱儲(chǔ)中部逐漸減弱,從而出現(xiàn)了頂部與底部溫度變化較慢,而中部溫度變化較快的現(xiàn)象。

    考慮熱補(bǔ)償時(shí),隨熱儲(chǔ)溫度的不斷降低,熱儲(chǔ)與蓋巖及基巖的溫差逐漸增大,蓋巖和基巖與熱儲(chǔ)間發(fā)生熱量傳遞,并造成前者的溫度變化。圖5給出了考慮熱補(bǔ)償情況下不同時(shí)刻蓋巖和基巖的溫度分布。由圖中可以看出,由于熱儲(chǔ)內(nèi)流體的對(duì)流傳熱作用,蓋巖和基巖與熱儲(chǔ)交界面附近的巖體首先被冷卻,在注入井附近形成低溫區(qū)域,隨時(shí)間的推移,低溫區(qū)域逐漸向水平及豎直方向擴(kuò)展。由于蓋巖和基巖中無(wú)流體流動(dòng),巖體溫度的變化僅由巖石的導(dǎo)熱作用引起,因此低溫區(qū)域在豎直方向上的擴(kuò)展速度較小。通過(guò)對(duì)比圖 5(a)、(c)、(e)、(g)和圖 5(b)、(d)、(f)、(h)不難發(fā)現(xiàn),系統(tǒng)運(yùn)行至20年時(shí),基巖中低溫區(qū)域在豎直方向上的最大波及范圍位于注入井附近,超過(guò)70 m,而同一時(shí)刻蓋巖中低溫區(qū)域在豎直方向上的最大波及范圍僅為50 m。由此可知,在系統(tǒng)開(kāi)采過(guò)程中基巖的熱量損失更為明顯。

    3.2 生產(chǎn)流體溫度

    圖6展示了考慮熱補(bǔ)償和忽略熱補(bǔ)償兩種情況下儲(chǔ)層出口處生產(chǎn)井內(nèi)流體溫度隨時(shí)間的分布。從圖中可以看出,兩種情況下生產(chǎn)井中生產(chǎn)流體的溫度分布在豎直方向上具有一定差異。在忽略熱補(bǔ)償?shù)那闆r下,由于重力作用儲(chǔ)層底部流體流速較快,導(dǎo)致儲(chǔ)層底部的生產(chǎn)流體首先出現(xiàn)溫度下降,從而產(chǎn)生儲(chǔ)層頂部生產(chǎn)流體溫度較高,底部生產(chǎn)流體溫度較低的現(xiàn)象。在系統(tǒng)運(yùn)行至20年時(shí),生產(chǎn)流體最低溫度位于儲(chǔ)層底部,溫度為340 K左右,最高溫度位于儲(chǔ)層頂部,溫度為380 K左右,頂部與底部溫差較大,溫度波動(dòng)范圍較大。而在考慮熱補(bǔ)償時(shí),熱補(bǔ)償對(duì)儲(chǔ)層底部流體起到了加熱作用,從而減緩了底部生產(chǎn)流體溫度下降的速率,同時(shí)生產(chǎn)流體首先出現(xiàn)溫度下降的位置由儲(chǔ)層底部上移至底部與中部之間。在系統(tǒng)運(yùn)行至20年時(shí),生產(chǎn)流體最低溫度位于儲(chǔ)層中部偏下,最低溫度為360 K左右,最高溫度與最低溫度溫差較小,整體溫度波動(dòng)范圍較小。通過(guò)對(duì)比圖6(a)、(b)可知,熱補(bǔ)償作用有效抑制了儲(chǔ)層底部生產(chǎn)流體溫度的過(guò)早過(guò)快下降,同時(shí)縮小了生產(chǎn)流體溫度在豎直方向上的波動(dòng)范圍,這有助于減小系統(tǒng)運(yùn)行后期生產(chǎn)流體溫度的下降速率,延長(zhǎng)系統(tǒng)運(yùn)行壽命。

    圖4 不同時(shí)刻熱儲(chǔ)溫度分布Fig.4 Temperature distribution of geothermal reservoir in different time

    圖7給出了兩種情況下生產(chǎn)流體平均溫度的變化曲線。由圖中可以看出,在忽略熱補(bǔ)償?shù)那闆r下,生產(chǎn)流體溫度出現(xiàn)明顯下降的時(shí)間較晚,溫度隨時(shí)間的下降速率較快。而考慮熱補(bǔ)償時(shí),生產(chǎn)流體溫度出現(xiàn)明顯下降的時(shí)間較早,溫度隨時(shí)間的下降速率較慢。這是由于考慮熱補(bǔ)償時(shí),熱儲(chǔ)溫度相對(duì)較高,流體流動(dòng)速度較快,因此流體溫度出現(xiàn)明顯下降的時(shí)間較為提前。同時(shí)由于上下巖層對(duì)熱儲(chǔ)的加熱作用,減慢了熱儲(chǔ)溫度的下降速率,使得生產(chǎn)流體溫度隨時(shí)間的下降速率也較為緩慢。

    圖5 不同時(shí)刻蓋巖和基巖溫度分布Fig.5 Temperature distribution of cap and base rocks in different time

    值得指出的是,在實(shí)際工程中,地?zé)嵯到y(tǒng)地表的能量轉(zhuǎn)換設(shè)備一般為汽輪機(jī)或熱泵機(jī)組,較低的生產(chǎn)流體溫度不利于熱能的轉(zhuǎn)換,為提高機(jī)組的能量轉(zhuǎn)換效率,假定以生產(chǎn)流體溫度下降至373 K(100℃)時(shí)作為系統(tǒng)運(yùn)行結(jié)束的時(shí)間,則從圖 7可以看出,考慮熱補(bǔ)償時(shí)系統(tǒng)運(yùn)行壽命可以達(dá)到20年,而忽略熱補(bǔ)償時(shí)只有18年左右,由此可見(jiàn)熱儲(chǔ)上下側(cè)巖層熱補(bǔ)償作用可以有效延長(zhǎng)CPGS的運(yùn)行壽命,從而對(duì)CPGS的采熱性能產(chǎn)生積極影響。

    3.3 系統(tǒng)采熱性能分析

    圖8展示了考慮熱補(bǔ)償和忽略熱補(bǔ)償兩種情況下生產(chǎn)流體質(zhì)量流量隨時(shí)間的變化曲線。運(yùn)行前期儲(chǔ)層內(nèi)生產(chǎn)流體僅有單相水,隨著CO2被逐漸注入儲(chǔ)層并突破生產(chǎn)井,生產(chǎn)井附近壓力梯度逐漸增大,生產(chǎn)流體的速度隨之增大,導(dǎo)致生產(chǎn)流體質(zhì)量流量逐漸上升。CO2突破生產(chǎn)井后,生產(chǎn)流體質(zhì)量流量大幅下降,這是由于CO2密度小于水,在相同條件下的質(zhì)量流量小于水,從而引起生產(chǎn)流體質(zhì)量流量的下降。隨著生產(chǎn)流體中CO2飽和度的增大,CO2的相對(duì)滲透率隨之增大,造成CO2的流速增大,使得生產(chǎn)流體質(zhì)量流量在很長(zhǎng)一段時(shí)間內(nèi)呈上升趨勢(shì)。蓋巖和基巖的熱補(bǔ)償作用對(duì)生產(chǎn)流體質(zhì)量流量的影響主要體現(xiàn)在CO2突破生產(chǎn)井后,從圖8中可以看出,考慮熱補(bǔ)償作用時(shí),生產(chǎn)流體的質(zhì)量流量有一定的提升。這是由于熱補(bǔ)償?shù)淖饔眠^(guò)程比較緩慢,在系統(tǒng)運(yùn)行前期并不明顯,隨時(shí)間的推移熱補(bǔ)償作用導(dǎo)致熱儲(chǔ)溫度較高,流體流動(dòng)性能較強(qiáng),流動(dòng)速度較快,因此質(zhì)量流量較高。

    圖6 儲(chǔ)層出口區(qū)域生產(chǎn)井中流體溫度隨時(shí)間分布Fig.6 Distribution of temperature of production fluid in thermal reservoir

    圖7 生產(chǎn)流體平均溫度曲線Fig.7 Average temperature curves of production fluid

    熱提取率是評(píng)價(jià)系統(tǒng)性能的重要參數(shù)。Pruess[11]的熱提取率計(jì)算公式為

    圖8 生產(chǎn)流體質(zhì)量流量隨時(shí)間變化曲線Fig.8 Mass flow rate time-varying curves of production fluid

    式(7)僅適用于單相流體的計(jì)算,因此采用式(7)的修改公式[22]進(jìn)行熱提取率的計(jì)算,即

    式中,Hpro-CO2為生產(chǎn)井中氣相 CO2的焓值;Hpro-H2O為生產(chǎn)井中水的焓值;Hinj-CO2為注入井中CO2的焓值;Fpro-CO2為生產(chǎn)井中的CO2的流量。CO2和水的焓值由注入和產(chǎn)出時(shí)的溫度和壓力確定,其焓值與壓力及溫度的關(guān)系已在圖1中給出。

    圖9 熱提取率隨時(shí)間變化曲線Fig.9 Time-varying curves of thermal extraction efficiency

    圖9展示的是系統(tǒng)熱提取率隨時(shí)間的變化曲線。從圖中可以看出,熱提取率隨時(shí)間的變化存在兩個(gè)峰,在第2個(gè)峰的峰值點(diǎn)前,熱提取率的變化規(guī)律與質(zhì)量流量的變化一致,這是由于熱提取率在第2次下降之前,生產(chǎn)流體的溫度幾乎不變,熱提取率的變化主要由生產(chǎn)流體質(zhì)量流量決定,當(dāng)生產(chǎn)流體溫度開(kāi)始下降時(shí),生產(chǎn)流體焓值隨之下降,從而引起熱提取率的下降。熱補(bǔ)償對(duì)熱提取率的影響起始于系統(tǒng)運(yùn)行至2年左右,也就是CO2突破生產(chǎn)井時(shí)開(kāi)始,考慮熱補(bǔ)償時(shí)的熱提取率相對(duì)較高,同時(shí)到達(dá)第2個(gè)峰值的時(shí)間較晚。在CO2突破生產(chǎn)井至第 2峰值點(diǎn)之間,考慮熱補(bǔ)償?shù)那闆r下生產(chǎn)流體質(zhì)量流量較高,導(dǎo)致其熱提取率較高。在第2個(gè)峰值點(diǎn)之后,由于考慮熱補(bǔ)償時(shí)生產(chǎn)流體溫度隨時(shí)間變化的速率較慢,因此熱提取率的下降也較慢。

    與 CPGS總的運(yùn)行時(shí)間相比,CO2驅(qū)替水的過(guò)程較短,因此系統(tǒng)正常的工作時(shí)間應(yīng)從CO2突破生產(chǎn)井后開(kāi)始。由圖 9可知,考慮熱補(bǔ)償時(shí)系統(tǒng)的熱提取率較高,同時(shí)熱提取率隨時(shí)間的變化幅度較小,這一結(jié)果表明,熱補(bǔ)償作用能夠有效提高CPGS的熱提取率,同時(shí)減弱熱提取率隨時(shí)間的波動(dòng),使CPGS具有更好的采熱性能和更穩(wěn)定的產(chǎn)能輸出。

    3.4 系統(tǒng)熱開(kāi)采量

    圖10給出了考慮和忽略熱補(bǔ)償兩種情況下各巖層熱開(kāi)采量及系統(tǒng)總熱開(kāi)采量隨時(shí)間的變化曲線??紤]熱補(bǔ)償時(shí),系統(tǒng)總熱開(kāi)采量由3部分組成,其中熱儲(chǔ)開(kāi)采量為系統(tǒng)從熱儲(chǔ)中獲得的熱量,是總熱開(kāi)采量的主要組成部分,蓋巖和基巖開(kāi)采量是熱補(bǔ)償作用對(duì)系統(tǒng)采熱的貢獻(xiàn)量,在總熱開(kāi)采量中所占的比例較小。熱儲(chǔ)、蓋巖、基巖及系統(tǒng)總熱開(kāi)采量通過(guò)式(9)~式(11)進(jìn)行計(jì)算。

    式中,Wi、Wr和W分別為蓋巖或基巖開(kāi)采量、熱儲(chǔ)開(kāi)采量和總熱開(kāi)采量;vi為蓋巖或基巖的巖層體積;Tini和T(t)分別為巖石初始和t時(shí)刻的溫度。

    圖10 熱開(kāi)采量隨時(shí)間變化曲線Fig.10 Time-varying curves of heat production

    從圖 10中可以看出考慮熱補(bǔ)償時(shí)的總熱開(kāi)采量較高,但熱儲(chǔ)開(kāi)采量較少,這是由于蓋巖與基巖的熱補(bǔ)償作用減緩了熱儲(chǔ)溫度下降的速率。蓋巖和基巖的熱補(bǔ)償為總熱開(kāi)采量做出了一定貢獻(xiàn),其中基巖的熱補(bǔ)償作用更為明顯。在系統(tǒng)運(yùn)行至 20年時(shí),考慮熱補(bǔ)償?shù)那闆r下總熱開(kāi)采量達(dá)到了1.93×1015J,比忽略熱補(bǔ)償時(shí)高出1.5×1014J,同時(shí)蓋巖和基巖開(kāi)采量分別達(dá)到了 1.43×1014J和3.18×1014J,分別占總熱開(kāi)采量的7.4%和16.5%,基巖對(duì)總熱開(kāi)采量的貢獻(xiàn)量約為蓋巖的兩倍。造成基巖開(kāi)采量較高的原因主要有兩點(diǎn),一是由于地溫梯度的存在,基巖的溫度高于蓋巖,在相同情況下與熱儲(chǔ)的溫差更大,熱補(bǔ)償作用更強(qiáng);二是由于重力作用,儲(chǔ)層底部流體的流速更快,換熱作用更為劇烈,相同時(shí)間內(nèi)能夠從基巖中獲得更多的熱量。

    4 結(jié) 論

    建立了 CPGS蓋巖-儲(chǔ)層-基巖整體模型,研究了熱儲(chǔ)上下側(cè)巖層熱補(bǔ)償作用對(duì)系統(tǒng)性能的影響。結(jié)果表明,熱補(bǔ)償作用對(duì)CPGS生產(chǎn)流體溫度、質(zhì)量流量、系統(tǒng)熱提取率及熱開(kāi)采量均會(huì)產(chǎn)生一定的影響,在評(píng)估采熱性能時(shí)是一不能簡(jiǎn)單忽略的重要因素。熱補(bǔ)償作用可減小重力作用造成的生產(chǎn)井中生產(chǎn)流體在豎直方向上的溫度差異,同時(shí)減小系統(tǒng)運(yùn)行后期生產(chǎn)流體溫度下降的速率,有效延長(zhǎng)系統(tǒng)運(yùn)行壽命;有助于提高CPGS生產(chǎn)流體質(zhì)量流量和系統(tǒng)熱提取率,并減弱熱提取率隨時(shí)間的波動(dòng),使系統(tǒng)獲得更好的采熱性能和更穩(wěn)定的產(chǎn)能輸出;考慮熱補(bǔ)償時(shí)系統(tǒng)總熱開(kāi)采量明顯提高,但由于熱補(bǔ)償?shù)募訜嶙饔脺p緩了熱儲(chǔ)溫度的下降速率,因此熱儲(chǔ)開(kāi)采量減少。此外,與蓋巖相比,基巖的熱補(bǔ)償作用更為明顯,本文算例中基巖對(duì)總熱開(kāi)采量的貢獻(xiàn)量約為蓋巖的兩倍。

    [1]HOLLOWAY S.Storage of fossil fuel-derived carbon dioxide beneath the surface of the Earth[J].Annual Review Energy Environment,2001,26:145-166.

    [2]WEST J M,PEARCE J,BENTHAM M,et al.Issue profile:environmental issues and the geological storage of CO2[J].European Environment,2005,15:250-259.

    [3]GENTZIS T.Subsurface sequestration of carbon dioxide-an overview from an Alberta (Canada) perspective[J].International Journal of Coal Geology,2000,43:287-305.

    [4]GOUGH C.State of the art in carbon dioxide capture and storage in the UK:an experts’ review[J].International Journal of Greenhouse Gas Control,2008,2:155-168.

    [5]HOLLOWAY S.Underground sequestration of carbon dioxide—a viable greenhouse gas mitigation option[J].Energy,2005,30:2318-2333.

    [6]METZ B,DAVIDSON O R,BOSCH P R,et al.Contribution of working group Ⅲ to the fourth assessment report of the intergovernmental panel on climate change[R].Cambridge:Cambridge University Press,2007.

    [7]BROWN D.A hot dry rock geothermal energy concept utilizing supercritical CO2instead of water[C]// Proceedings of the Twenty-Fifth Workshop on Geothermal Reservoir Engineering,Stanford,2000:233-238.

    [8]PRUESS K.Enhanced geothermal systems (EGS) using CO2as working fluid-a novel approach for generating renewable energy with simultaneous sequestration of carbon[J].Geothermics,2006,35:351-367.

    [9]PRUESS K.Enhanced geothermal systems (EGS) comparing water with CO2as heat transmission fluids[R].Lawrence Berkeley National Laboratory,2007.

    [10]PRUESS K.On the feasibility of using supercritical CO2as heat transmission fluid in an engineered hot dry rock geothermal system thirty-first workshop on geothermal reservoir engineering:SGP-TR-179[R].California:Stanford University,2006.

    [11]PRUESS K.On production behavior of enhanced geothermal systems with CO2as working fluid[J].Energy Conversion and Management,2008,49(6):1446-1454.

    [12]MAJER E L,BARIA R,STARK M,et al.Induced seismicity associated with enhanced geothermal systems[J].Geothermics,2007,36:185-222.

    [13]GLANZ J.Deep in bedrock,clean energy and quake fears[N].New York:New York Times,2009-06-23 (1).

    [14]RANDOLPH J B,SAAR M O.Coupling geothermal energy capture with carbon dioxide sequestration in naturally permeable,porous geologic formations:a comparison with enhanced geothermal systems[J].GRC Trans.,2010,34:433-438.

    [15]RANDOLPH J B,SAAR M O.Combining geothermal energy capture with geologic carbon dioxide sequestration[J].Geophys Res.Lett.,2011,38(10):415-421.

    [16]RANDOLPH J B,SAAR M O.Coupling carbon dioxide sequestration with geothermal energy capture in naturally permeable,porous geologic formations:implications for CO2sequestration[J].Energy Procedia,2011,4:2206-2213.

    [17]ZHANG L,EZEKIEL J,LI D X,et al.Potential assessment of CO2injection for heat mining and geological storage in geothermal reservoirs of China[J].Applied Energy,2014,122:237-246.

    [18]XU T F,FENG G H,SHI Y.On fluid-rock chemical interaction in CO2-based geothermal systems[J].Journal of Geochemical Exploration,2014,144:179-193.

    [19]陳繼良,羅良,蔣方明.熱儲(chǔ)周圍巖石熱補(bǔ)償對(duì)增強(qiáng)型地?zé)嵯到y(tǒng)采熱過(guò)程的影響[J].計(jì)算物理,2013,30(6):862-870.CHEN J L,LUO L,JIANG F M.Thermal compensation of rocks encircling heat reservoir in heat extraction of enhanced geothermal system[J].Chinese Journal of Computational Physics,2013,30(6):862-870.

    [20]陳繼良,蔣方明.增強(qiáng)型地?zé)嵯到y(tǒng)熱開(kāi)采性能的數(shù)值模擬分析[J].可再生能源,2013,31(12):111-117.CHEN J L,JIANG F M.A numerical study on heat extraction performance of enhanced geothermal systems[J].Renewable Energy Resources,2013,31(12):111-117.

    [21]封官宏,李佳琦,許天福,等.二氧化碳羽流地?zé)嵯到y(tǒng)中儲(chǔ)層物性參數(shù)對(duì)熱提取率的影響[J].可再生能源,2013,31(7):85-92.FENG G H,LI J Q,XU T F,et al.Effects of property of reservoir on heat extraction in CO2plume geothermal system[J].Renewable Energy Resources,2013,31(7):85-92.

    [22]魏銘聰,楊冰,許天福,等.二氧化碳羽流地?zé)嵯到y(tǒng)中井間距和儲(chǔ)層滲透率對(duì)熱提取率的影響:以松遼盆地為例[J].地質(zhì)科技情報(bào),2015,34(2):188-193.WEI M C,YANG B,XU T F,et al.Effects of well spacing and reservoir permeability on heat extraction in CO2plume geothermal system:a case study of Songliao basin[J].Geological Science and Technology Information,2015,34(2):188-193.

    [23]任韶然,崔國(guó)棟,李德祥,等.注超臨界CO2開(kāi)采高溫廢棄氣藏地?zé)釞C(jī)制與采熱能力分析[J].中國(guó)石油大學(xué)學(xué)報(bào)(自然科學(xué)版),2016,40(2):91-98.REN S R,CUI G D,LI D X,et al.Development of geothermal energy from depleted high temperature gas reservoirviasupercritical CO2injection[J].Journal of China University of Petroleum(Edition of Natural Science),2016,40(2):91-98.

    [24]GARAPATI N,RANDOLPH J B,JOSE L,et al.CO2-plume geothermal (CPG) heat extraction in multi-layered geologic reservoirs[J].Energy Procedia,2014,63:7631-7643.

    [25]GARAPATI N,RANDOLPH J B,SAAR M O.Brine displacement by CO2,energy extraction rates,and lifespan of a CO2-limited CO2-Plume Geothermal (CPG) system with a horizontal production well[J].Geothermics,2015,55:182-194.

    [26]SHI Y,WANG F G,YANG Y L,et al.Use of a CO2geological storage system to develop geothermal resources:a case study of a sandstone reservoir in the Songliao basin of northeast China[M]//HOU M Z,XIE H P,WERE P.Clean Energy Systems in the Subsurface:Production,Storage and Conversion.Berlin:Springer,2013:89-103.

    [27]ADAMS B M,KUEHN T H,BIELICKI J M,et al.On the importance of the thermosiphon effect in CPG (CO2plume geothermal) power systems[J].Energy,2014,69:409-418.

    [28]ADAMS B M,KUEHN T H,BIELICKI J M,et al.A comparison of electric power output of CO2Plume Geothermal (CPG) and brine geothermal systems for varying reservoir conditions[J].Applied Energy,2015,140:365-377.

    [29]CHEN J L,JIANG F M.A numerical study of EGS heat extraction process based on a thermal non-equilibrium model for heat transfer in subsurface porous heat reservoir[J].Heat and Mass Transfer,2016,52:255-267.

    [30]JIANG F M,CHEN J L,HUANG W B,et al.A three-dimensional transient model for EGS subsurface thermo-hydraulic process[J].Energy,2014,72:300-310.

    [31]CHEN J L,JIANG F M.Designing multi-well layout for enhanced geothermal system to better exploit hot dry rock geothermal energy[J].Renewable Energy,2015,74:37-48.

    [32]CAO W J,HUANG W B,JIANG F M.Numerical study on variable thermophysical properties of heat transfer fluid affecting EGS heat extraction[J].International Journal of Heat and Mass Transfer,2016,92:1205-1217.

    [33]FAGERLUND F F,NIEMI A,ODEN M.Comparison of relative permeability-fluid saturation-capillary pressure relations in the modelling of non-aqueous phase liquid infiltration in variably saturated,layered media[J].Advances in Water Resources,2006,29(11):1705-1730

    [34]楊艷林,靖晶,王福剛,等.CO2增強(qiáng)型地?zé)嵯到y(tǒng)中的井網(wǎng)間距優(yōu)化研究[J].太陽(yáng)能學(xué)報(bào),2014,(7):1130-1137.YANG Y L,JING J,WANG F G,et al.Optimal design of well spacing on CO2enhanced geothermal[J].Acta Energiae Solaris Sinica,2014,(7):1130-1137.

    [35]張俊虎,劉君.煤層氣井網(wǎng)布置優(yōu)化設(shè)計(jì)的探討[J].科技情報(bào)開(kāi)發(fā)與經(jīng)濟(jì),2008,(10):210-212.ZHANG J H,LIU J.Probe into the optimal design of coal-bed methane well network[J].Sci-tech Information Development &Economy,2008,(10):210-212.

    date:2017-04-17.

    Prof.LIU Zhongliang,liuzhl@ bjut.edu.cn

    Influence of thermal compensation of geothermal reservoir rock formation on CO2plume geothermal system performance

    LI Jingyan,LIU Zhongliang,ZHOU Yu,LI Yanxia
    (Key Laboratory of Enhanced Heat Transfer and Energy Conservation,Ministry of Education,College of Environmental and Energy Engineering,Beijing University of Technology,Beijing100124,China)

    Carbon dioxide plume geothermal system (CPGS) can be used to exploit geothermal energy and realize carbon dioxide geological sequestration simultaneously.The thermal compensation from the rock formation around the geothermal reservoir is one of the important factors that affect the performance of CPGS.Based on a three dimensional base and cap rocks enclosed heat reservoir model,the influences of the thermal compensation on the temperature evolutionary process of the rock and fluid in the geothermal reservoir and the heat collection performance of CPGS were studied.The distribution of geothermal reservoir temperature and the temperature of production fluid were compared with that without consideration of the thermal compensation.The results show that the thermal compensation reduces both the production fluid temperature variation along the vertical direction and its temperature decreasing rate in the later period of system operation,therefore extends the lifetime of CPGS and gains better heat collection performance.With consideration of the thermal compensation,the heat production is improved significantly.The results also show that the thermal compensation of the base rocks is stronger than that of the cap rocks.

    rock formation; geothermal reservoir; thermal compensation; carbon dioxide; geothermal system;porous media; numerical simulation

    TK 124

    A

    0438—1157(2017)12—4526—11

    10.11949/j.issn.0438-1157.20170413

    2017-04-17收到初稿,2017-07-15收到修改稿。

    聯(lián)系人:劉中良。

    李靜巖(1992—),男,碩士研究生。

    猜你喜歡
    生產(chǎn)井開(kāi)采量基巖
    鄂爾多斯某區(qū)塊致密油產(chǎn)能影響因素分析與優(yōu)化研究
    化工管理(2021年23期)2021-08-25 03:01:00
    再談河北省灤平縣馬營(yíng)子鄉(xiāng)高鍶天然礦泉水特征與開(kāi)采量估算
    基于流線的聚驅(qū)雙高油田調(diào)整策略研究
    輸水渠防滲墻及基巖滲透系數(shù)敏感性分析
    基于改進(jìn)物元的大壩基巖安全評(píng)價(jià)
    河北省基巖熱儲(chǔ)開(kāi)發(fā)利用前景
    利用統(tǒng)計(jì)分析法預(yù)測(cè)地?zé)崴砷_(kāi)采量應(yīng)注意的問(wèn)題
    中國(guó)新疆石油開(kāi)采量總額增長(zhǎng)
    中亞信息(2016年1期)2016-03-07 01:00:42
    EILog生產(chǎn)井技術(shù)平臺(tái)的五員“開(kāi)路先鋒”
    生產(chǎn)井實(shí)時(shí)檢測(cè)工藝技術(shù)及系統(tǒng)分析
    18禁动态无遮挡网站| 在线观看国产h片| 美女福利国产在线 | 免费av中文字幕在线| 国内揄拍国产精品人妻在线| 黑人高潮一二区| 涩涩av久久男人的天堂| 国产午夜精品一二区理论片| 在现免费观看毛片| 草草在线视频免费看| 少妇精品久久久久久久| 日韩强制内射视频| 国产av精品麻豆| 美女高潮的动态| 激情 狠狠 欧美| www.色视频.com| 最黄视频免费看| 精品一品国产午夜福利视频| 成人无遮挡网站| 久久国产乱子免费精品| 另类亚洲欧美激情| 五月天丁香电影| 啦啦啦在线观看免费高清www| 久久精品人妻少妇| 午夜福利影视在线免费观看| 九草在线视频观看| 亚洲国产av新网站| 嫩草影院入口| 纯流量卡能插随身wifi吗| 国产av精品麻豆| 大片电影免费在线观看免费| 干丝袜人妻中文字幕| 欧美变态另类bdsm刘玥| 成人影院久久| 全区人妻精品视频| 国产爽快片一区二区三区| 免费大片18禁| 国产精品国产三级专区第一集| 亚洲国产精品国产精品| 国产精品欧美亚洲77777| 日本黄色日本黄色录像| 老司机影院成人| 一本—道久久a久久精品蜜桃钙片| 亚洲中文av在线| 特大巨黑吊av在线直播| 男女边吃奶边做爰视频| 国产成人一区二区在线| 麻豆精品久久久久久蜜桃| av一本久久久久| 免费播放大片免费观看视频在线观看| 亚洲电影在线观看av| 日韩视频在线欧美| 青春草国产在线视频| 天堂8中文在线网| 亚洲精品视频女| 国产欧美亚洲国产| 一个人免费看片子| 简卡轻食公司| 日本爱情动作片www.在线观看| 26uuu在线亚洲综合色| 国产熟女欧美一区二区| 日韩视频在线欧美| 麻豆成人av视频| 日韩欧美精品免费久久| 汤姆久久久久久久影院中文字幕| 青春草视频在线免费观看| 亚洲av电影在线观看一区二区三区| 人人妻人人澡人人爽人人夜夜| 欧美xxxx性猛交bbbb| 国产男人的电影天堂91| 国产男女超爽视频在线观看| 亚洲人与动物交配视频| 我的女老师完整版在线观看| 麻豆成人av视频| 亚洲四区av| 亚洲av成人精品一二三区| 97超碰精品成人国产| 精品久久久久久久末码| 欧美精品一区二区大全| 高清在线视频一区二区三区| 天堂8中文在线网| 91久久精品电影网| 久久久久久人妻| 欧美另类一区| 大话2 男鬼变身卡| 亚洲av成人精品一区久久| 亚洲精品成人av观看孕妇| 中文字幕精品免费在线观看视频 | 欧美成人精品欧美一级黄| 久久久国产一区二区| 色综合色国产| 亚洲精品成人av观看孕妇| 国产高清不卡午夜福利| 国产成人免费观看mmmm| 狂野欧美激情性bbbbbb| 另类亚洲欧美激情| 国产精品免费大片| 97在线人人人人妻| 中文字幕av成人在线电影| 激情 狠狠 欧美| 中文字幕精品免费在线观看视频 | 亚洲精品视频女| 人体艺术视频欧美日本| 亚洲,欧美,日韩| 国产精品人妻久久久久久| 午夜视频国产福利| 国产av一区二区精品久久 | 国产又色又爽无遮挡免| 国产日韩欧美在线精品| 国产乱人偷精品视频| 欧美高清性xxxxhd video| 国产视频内射| 女性被躁到高潮视频| 精品久久久精品久久久| av卡一久久| 免费大片18禁| 青青草视频在线视频观看| 国国产精品蜜臀av免费| 亚洲自偷自拍三级| 一本色道久久久久久精品综合| 最新中文字幕久久久久| 国产精品偷伦视频观看了| 日产精品乱码卡一卡2卡三| 亚洲av欧美aⅴ国产| 全区人妻精品视频| 熟女电影av网| 制服丝袜香蕉在线| 最黄视频免费看| 中文精品一卡2卡3卡4更新| 中国国产av一级| 性色avwww在线观看| 国产精品久久久久久精品电影小说 | 我的女老师完整版在线观看| 国产国拍精品亚洲av在线观看| 有码 亚洲区| 男女国产视频网站| 国产精品女同一区二区软件| av卡一久久| 精品一区二区三区视频在线| 黄色视频在线播放观看不卡| 亚洲人与动物交配视频| 妹子高潮喷水视频| 国产乱人偷精品视频| 观看免费一级毛片| 亚洲国产毛片av蜜桃av| 男女无遮挡免费网站观看| 在线播放无遮挡| 免费看光身美女| 国产有黄有色有爽视频| 一级毛片黄色毛片免费观看视频| 午夜免费男女啪啪视频观看| av专区在线播放| 国产毛片在线视频| 全区人妻精品视频| 亚洲精品自拍成人| 色吧在线观看| 99九九线精品视频在线观看视频| 日本欧美国产在线视频| 亚洲国产色片| 国产69精品久久久久777片| av线在线观看网站| 97精品久久久久久久久久精品| 久久久久精品性色| 国产精品一二三区在线看| 欧美成人午夜免费资源| h视频一区二区三区| 日本黄色日本黄色录像| 欧美日韩综合久久久久久| 中文天堂在线官网| 亚洲精品久久久久久婷婷小说| 国产午夜精品久久久久久一区二区三区| av播播在线观看一区| 黄片wwwwww| 亚洲av.av天堂| av一本久久久久| 久久国产亚洲av麻豆专区| 最后的刺客免费高清国语| 欧美日韩国产mv在线观看视频 | 免费看不卡的av| 亚洲国产日韩一区二区| 制服丝袜香蕉在线| 自拍偷自拍亚洲精品老妇| 国产乱人视频| 又大又黄又爽视频免费| 少妇的逼好多水| 少妇裸体淫交视频免费看高清| 亚洲精品久久午夜乱码| 国产精品一及| 国产精品免费大片| 毛片一级片免费看久久久久| 久久久久网色| 国产乱来视频区| 久久精品国产鲁丝片午夜精品| 国产色爽女视频免费观看| 久久久欧美国产精品| 国产精品一区二区三区四区免费观看| 大香蕉久久网| 国产白丝娇喘喷水9色精品| 久久国产精品大桥未久av | 在线观看免费视频网站a站| 男人添女人高潮全过程视频| 亚洲第一av免费看| 中文字幕久久专区| 99re6热这里在线精品视频| 午夜福利在线观看免费完整高清在| 久久99精品国语久久久| 亚洲欧美成人综合另类久久久| 99热网站在线观看| 欧美性感艳星| 18禁在线无遮挡免费观看视频| 女性被躁到高潮视频| 欧美日韩亚洲高清精品| 一级爰片在线观看| 免费看光身美女| 男女免费视频国产| 亚洲真实伦在线观看| 日韩三级伦理在线观看| 最近中文字幕高清免费大全6| 2021少妇久久久久久久久久久| 欧美 日韩 精品 国产| 丝袜喷水一区| 亚洲精品一二三| 亚洲av二区三区四区| 亚洲熟女精品中文字幕| 久久国产乱子免费精品| 国产精品国产三级国产av玫瑰| 亚洲国产欧美人成| 免费大片18禁| 老司机影院成人| 少妇熟女欧美另类| 亚洲国产毛片av蜜桃av| 黄色一级大片看看| 成人毛片a级毛片在线播放| 只有这里有精品99| 亚州av有码| 日韩成人av中文字幕在线观看| 乱系列少妇在线播放| 中文字幕免费在线视频6| h视频一区二区三区| 美女脱内裤让男人舔精品视频| 国产精品偷伦视频观看了| 日韩av免费高清视频| 性色avwww在线观看| 精品99又大又爽又粗少妇毛片| 人人妻人人看人人澡| 身体一侧抽搐| 一级毛片 在线播放| 搡老乐熟女国产| 王馨瑶露胸无遮挡在线观看| 国产精品国产三级国产专区5o| 久久久a久久爽久久v久久| 日韩,欧美,国产一区二区三区| 国产精品国产三级专区第一集| 国产一区亚洲一区在线观看| 91在线精品国自产拍蜜月| 下体分泌物呈黄色| 久久精品国产鲁丝片午夜精品| 秋霞伦理黄片| 日韩欧美 国产精品| 2018国产大陆天天弄谢| 亚洲高清免费不卡视频| 乱系列少妇在线播放| 国产在线一区二区三区精| 看十八女毛片水多多多| av福利片在线观看| 22中文网久久字幕| 99久久精品热视频| av卡一久久| 国产精品麻豆人妻色哟哟久久| 欧美激情极品国产一区二区三区 | 91精品一卡2卡3卡4卡| 国产精品99久久久久久久久| 熟女电影av网| 婷婷色麻豆天堂久久| 18禁裸乳无遮挡动漫免费视频| 日韩伦理黄色片| 国产精品一及| 久久ye,这里只有精品| 国内精品宾馆在线| 久久 成人 亚洲| 日本欧美国产在线视频| 一二三四中文在线观看免费高清| 热99国产精品久久久久久7| 又爽又黄a免费视频| 91精品伊人久久大香线蕉| 丝袜脚勾引网站| 97在线人人人人妻| 欧美+日韩+精品| 欧美精品人与动牲交sv欧美| 日本wwww免费看| 成人漫画全彩无遮挡| 国产 精品1| 亚洲精品一区蜜桃| 国产欧美亚洲国产| 亚洲激情五月婷婷啪啪| 国产精品三级大全| 精品久久久精品久久久| 国产91av在线免费观看| 看免费成人av毛片| av播播在线观看一区| 大码成人一级视频| 国产av国产精品国产| 午夜激情福利司机影院| 特大巨黑吊av在线直播| kizo精华| 国内少妇人妻偷人精品xxx网站| 亚洲国产欧美在线一区| 能在线免费看毛片的网站| 晚上一个人看的免费电影| 91精品一卡2卡3卡4卡| 制服丝袜香蕉在线| 精品久久久久久电影网| 一级毛片久久久久久久久女| 精品一区二区三区视频在线| 内地一区二区视频在线| 97超碰精品成人国产| 男的添女的下面高潮视频| 日韩中文字幕视频在线看片 | 丝袜喷水一区| 网址你懂的国产日韩在线| 婷婷色综合大香蕉| 久久人人爽av亚洲精品天堂 | 欧美区成人在线视频| 成人黄色视频免费在线看| 免费大片18禁| 日韩 亚洲 欧美在线| 秋霞伦理黄片| 1000部很黄的大片| 亚洲精品乱久久久久久| 国模一区二区三区四区视频| 国产一级毛片在线| 亚洲欧美成人精品一区二区| 国产成人精品久久久久久| 97超视频在线观看视频| 卡戴珊不雅视频在线播放| 亚洲精品456在线播放app| 欧美高清成人免费视频www| 久久热精品热| 精品午夜福利在线看| av国产久精品久网站免费入址| a级毛色黄片| 日韩在线高清观看一区二区三区| 国产免费一级a男人的天堂| 色吧在线观看| 亚洲成人手机| 国产久久久一区二区三区| 91精品一卡2卡3卡4卡| 永久免费av网站大全| 国产精品国产三级国产av玫瑰| 久久久久久久精品精品| 少妇人妻精品综合一区二区| 精品人妻视频免费看| 性高湖久久久久久久久免费观看| 久久久亚洲精品成人影院| 精品久久国产蜜桃| 麻豆成人午夜福利视频| 国产免费又黄又爽又色| 欧美xxⅹ黑人| 一本—道久久a久久精品蜜桃钙片| 色视频在线一区二区三区| 人妻系列 视频| av视频免费观看在线观看| 天堂中文最新版在线下载| 97超视频在线观看视频| 黄色视频在线播放观看不卡| 免费看不卡的av| 亚洲精品视频女| 欧美国产精品一级二级三级 | 国产高清国产精品国产三级 | 美女国产视频在线观看| 久久久色成人| 欧美老熟妇乱子伦牲交| 日本色播在线视频| 久久精品久久精品一区二区三区| av在线app专区| 日本黄大片高清| 国产av国产精品国产| 中文字幕久久专区| 亚洲真实伦在线观看| av国产久精品久网站免费入址| av黄色大香蕉| 国产免费福利视频在线观看| 超碰av人人做人人爽久久| 亚洲国产精品国产精品| 一个人看的www免费观看视频| 少妇丰满av| 天堂8中文在线网| 天堂俺去俺来也www色官网| 亚洲丝袜综合中文字幕| 在线看a的网站| 亚州av有码| av卡一久久| 久久久久久人妻| 大香蕉久久网| 深夜a级毛片| 校园人妻丝袜中文字幕| 在现免费观看毛片| 黄色欧美视频在线观看| 三级国产精品欧美在线观看| 亚洲精品乱久久久久久| 国内精品宾馆在线| av专区在线播放| 日本av免费视频播放| 亚洲婷婷狠狠爱综合网| 久久久久久久久久久丰满| 亚洲成色77777| 亚州av有码| 少妇 在线观看| 久久精品熟女亚洲av麻豆精品| 久久99精品国语久久久| 成人漫画全彩无遮挡| 国产精品一区二区在线不卡| 一级毛片 在线播放| 亚洲精品成人av观看孕妇| 精品国产露脸久久av麻豆| 日韩强制内射视频| 亚洲精品国产av蜜桃| 亚洲精品aⅴ在线观看| 亚洲综合精品二区| 成人免费观看视频高清| 一级片'在线观看视频| 九九爱精品视频在线观看| 两个人的视频大全免费| 嫩草影院入口| 国产精品99久久久久久久久| 高清av免费在线| 亚洲国产精品国产精品| 久久99精品国语久久久| 国模一区二区三区四区视频| 小蜜桃在线观看免费完整版高清| 亚洲欧洲国产日韩| 久久久久精品性色| 亚洲欧美日韩东京热| 亚洲国产精品专区欧美| 免费黄色在线免费观看| 久久久久久久久久久免费av| 国产一区亚洲一区在线观看| 在线免费十八禁| 赤兔流量卡办理| 免费高清在线观看视频在线观看| 久久综合国产亚洲精品| 久热久热在线精品观看| 精品熟女少妇av免费看| 日韩av免费高清视频| 欧美 日韩 精品 国产| 国产精品久久久久久久电影| 亚洲国产欧美人成| 国产日韩欧美亚洲二区| 久久热精品热| 日日摸夜夜添夜夜爱| 一级爰片在线观看| 欧美亚洲 丝袜 人妻 在线| 日韩欧美精品免费久久| 国产av一区二区精品久久 | 日韩,欧美,国产一区二区三区| 亚洲综合色惰| 国产 精品1| 欧美激情国产日韩精品一区| 日本vs欧美在线观看视频 | 亚洲,一卡二卡三卡| 直男gayav资源| 欧美三级亚洲精品| 一级毛片电影观看| 国产一级毛片在线| 最近中文字幕2019免费版| 男女边吃奶边做爰视频| 国产乱来视频区| 一个人看的www免费观看视频| 国产一区亚洲一区在线观看| 国产人妻一区二区三区在| 国产成人一区二区在线| 国产精品秋霞免费鲁丝片| 18禁裸乳无遮挡免费网站照片| 亚洲性久久影院| 午夜激情久久久久久久| 乱系列少妇在线播放| 精品国产一区二区三区久久久樱花 | h日本视频在线播放| 一级毛片 在线播放| 亚洲精品中文字幕在线视频 | 又黄又爽又刺激的免费视频.| 久久久午夜欧美精品| 免费观看的影片在线观看| 日日摸夜夜添夜夜添av毛片| 高清日韩中文字幕在线| 久久人人爽人人片av| 久久综合国产亚洲精品| 免费黄频网站在线观看国产| 街头女战士在线观看网站| 夜夜骑夜夜射夜夜干| 亚洲欧美日韩东京热| 国产精品欧美亚洲77777| 熟妇人妻不卡中文字幕| 永久网站在线| 中国国产av一级| 免费黄网站久久成人精品| 欧美日韩综合久久久久久| 免费久久久久久久精品成人欧美视频 | 国产视频内射| 亚洲国产高清在线一区二区三| 国产精品一区二区在线不卡| 国产在线一区二区三区精| 这个男人来自地球电影免费观看 | 精品熟女少妇av免费看| 噜噜噜噜噜久久久久久91| 在现免费观看毛片| 国产色爽女视频免费观看| 大话2 男鬼变身卡| 日韩伦理黄色片| 嫩草影院新地址| 晚上一个人看的免费电影| 国产精品一区二区在线观看99| 一级毛片 在线播放| 国内精品宾馆在线| 在线观看免费视频网站a站| 男人舔奶头视频| 97在线视频观看| 建设人人有责人人尽责人人享有的 | 日本与韩国留学比较| 最近中文字幕高清免费大全6| 熟妇人妻不卡中文字幕| 看免费成人av毛片| 26uuu在线亚洲综合色| 六月丁香七月| 91精品一卡2卡3卡4卡| 多毛熟女@视频| 日本欧美视频一区| 美女国产视频在线观看| 91精品国产九色| 国产亚洲5aaaaa淫片| 能在线免费看毛片的网站| 日本黄色日本黄色录像| 亚州av有码| 成年免费大片在线观看| 国产熟女欧美一区二区| 久久 成人 亚洲| 国产伦精品一区二区三区四那| 亚洲av福利一区| 亚洲av在线观看美女高潮| 免费观看a级毛片全部| 插阴视频在线观看视频| 欧美+日韩+精品| 一本—道久久a久久精品蜜桃钙片| 国产精品国产三级国产专区5o| 看十八女毛片水多多多| 91在线精品国自产拍蜜月| 91狼人影院| 舔av片在线| 亚洲欧美成人精品一区二区| 国产男女内射视频| 亚洲人成网站高清观看| 婷婷色综合大香蕉| 免费看av在线观看网站| av播播在线观看一区| 精品人妻熟女av久视频| 日韩欧美 国产精品| 一本久久精品| 人体艺术视频欧美日本| 国产免费视频播放在线视频| 五月开心婷婷网| 五月玫瑰六月丁香| 蜜桃亚洲精品一区二区三区| 久久久午夜欧美精品| 一区二区三区乱码不卡18| 国产一区二区三区综合在线观看 | 在线观看免费日韩欧美大片 | 天天躁夜夜躁狠狠久久av| 一本一本综合久久| 欧美日韩综合久久久久久| 日本av手机在线免费观看| 欧美日本视频| 欧美另类一区| 午夜免费男女啪啪视频观看| 九色成人免费人妻av| 亚洲人成网站高清观看| 亚洲av男天堂| 精品亚洲乱码少妇综合久久| 精华霜和精华液先用哪个| 国产亚洲一区二区精品| 国产极品天堂在线| 最近2019中文字幕mv第一页| 久久久久久九九精品二区国产| 女性被躁到高潮视频| 日本黄大片高清| 最近2019中文字幕mv第一页| 欧美精品人与动牲交sv欧美| 亚洲欧美中文字幕日韩二区| 一级毛片黄色毛片免费观看视频| 一区二区三区四区激情视频| 中国美白少妇内射xxxbb| 欧美精品国产亚洲| 亚洲图色成人| 一区在线观看完整版| 久久久久国产精品人妻一区二区| 国产男女超爽视频在线观看| 99九九线精品视频在线观看视频| h视频一区二区三区| 一本色道久久久久久精品综合| 久热久热在线精品观看| 日韩欧美 国产精品| 大片免费播放器 马上看| 亚洲国产精品999| 日韩三级伦理在线观看| 香蕉精品网在线| 天天躁夜夜躁狠狠久久av| 最近手机中文字幕大全| 99国产精品免费福利视频| 高清av免费在线| 蜜桃久久精品国产亚洲av| 欧美亚洲 丝袜 人妻 在线| 国产成人91sexporn| 久久婷婷青草| 天美传媒精品一区二区| 少妇高潮的动态图| 婷婷色综合www|