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

    基于微地震數(shù)據(jù)的增強型地?zé)醿訁?shù)及采熱的數(shù)值模擬研究

    2022-11-23 03:12:02馬子涵邢會林靳國棟譚玉陽閆偉超李四海
    水文地質(zhì)工程地質(zhì) 2022年6期
    關(guān)鍵詞:生產(chǎn)井微震均質(zhì)

    馬子涵,邢會林,靳國棟,譚玉陽,閆偉超,李四海

    (1.深海圈層與地球系統(tǒng)前沿科學(xué)中心,海底科學(xué)與探測技術(shù)教育部重點實驗室,中國海洋大學(xué)海洋地球科學(xué)學(xué)院, 山東 青島 266100;2.青島海洋科學(xué)與技術(shù)試點國家實驗室,山東 青島 266237;3.中國海洋大學(xué)海底科學(xué)與工程計算國際中心,山東 青島 266100)

    地?zé)崮茏鳛橐环N清潔、穩(wěn)定、可再生的新型能源,符合綠色環(huán)保可持續(xù)的發(fā)展要求,被譽為21世紀(jì)最具潛力的能源之一[1]。干熱巖型地?zé)豳Y源是指儲存在地下低滲透性高溫巖體中的熱能[2],約占地?zé)豳Y源總量的30%,相當(dāng)于全球石油、天然氣和煤炭能量總和的30 倍[3]。但是,干熱巖儲層具有埋藏深、巖石致密、滲透率極低等特點,使用傳統(tǒng)技術(shù)難以有效開發(fā)其中的熱能[4]。利用增強型地?zé)嵯到y(tǒng)(Enhanced Geothermal Systems, EGS),即使用水力壓裂技術(shù)將原有儲層改造為高滲透性的人工儲層,而后通過水、二氧化碳等工作流體的循環(huán)注采,可實現(xiàn)干熱巖中熱能的高效開采[5]。EGS 技術(shù)作為實現(xiàn)深層干熱巖型地?zé)豳Y源綜合開發(fā)利用的關(guān)鍵技術(shù),已經(jīng)成為各國新能源發(fā)展的重點關(guān)注方向之一[6]。

    利用EGS 技術(shù)提取熱能是一個包含傳熱、滲流等多物理場耦合的復(fù)雜過程[7],因其影響參數(shù)眾多[8],物理性質(zhì)相差懸殊,時空演化不均衡等特點,通常利用數(shù)值模擬技術(shù)分析EGS 采熱過程[9-10]。近年來,國內(nèi)外學(xué)者提出了不同的數(shù)值模型,其中最常用的方法為等效多孔介質(zhì)模型。Jiang 等[6]、Ayling 等[11]建立了用于研究EGS 采熱的熱流耦合模型,將熱儲層等效為多孔介質(zhì),進行模擬并與實際情況進行對比,驗證了模型的有效性。Llanos 等[12]為庫珀盆地Habanero EGS 開發(fā)了一個三維熱流模型,考慮了滲透率各向異性,結(jié)果表明滲透率各向異性是控制儲層流體流動的關(guān)鍵,同時為干熱巖儲層的開發(fā)設(shè)計了不同的井布局,并得出東西向交錯的布局在該模型中為最優(yōu)井布局。Hogarth等[13]在此模型基礎(chǔ)上進行了單井注入和2 口井聯(lián)通、閉環(huán)測試的模擬,將模型計算的生產(chǎn)井溫度曲線與實際記錄溫度進行對比,驗證了模擬結(jié)果的準(zhǔn)確性。Samin 等[14]提出孔隙度均一的三維多孔介質(zhì)THM 模型,儲層分為基巖和裂縫,裂縫區(qū)的滲透率具有各向異性,以此模型優(yōu)化裂縫區(qū)滲透率、流體注入壓力、注采井距等,實現(xiàn)熱發(fā)電性能與成本控制的綜合改善。然而,均勻多孔介質(zhì)模型沒有考慮儲層空間分布上的非均質(zhì)性,過于理想化,不能精確地表征壓裂改造后的儲層滲流特征[15]。此外,另一種比較常用的方法為規(guī)則裂縫或隨機裂縫模型,Sun 等[1]提出二維模型中儲層由巖石基質(zhì)和規(guī)則裂縫組成,并以此模型對地?zé)醿又械牧黧w流動、傳熱等特征進行了綜合模擬,通過敏感性分析研究了EGS 采熱的主控參數(shù),結(jié)果表明THM 耦合效應(yīng)對采熱效率有重要影響。Xu等[2]提出由隨機裂縫和巖石基質(zhì)組成的模型,以此研究EGS 的熱流耦合系統(tǒng),將該模型應(yīng)用于Habanero 熱儲層,對其采熱過程進行模擬,應(yīng)用結(jié)果表明該模型是模擬工業(yè)規(guī)模地?zé)醿觽鳠岬挠行Х椒āang等[4]提出的THM 模型由巖石基質(zhì)塊體和隨機裂縫組成,模擬包含三維隨機裂縫的地?zé)醿拥拈L期地?zé)嵘a(chǎn)過程,并將其應(yīng)用于澳大利亞的Habanero EGS 項目,評估了20 a 期間不同注入和生產(chǎn)壓力下的采熱量,結(jié)果表明,通過適當(dāng)控制注采壓差,可以實現(xiàn)持續(xù)穩(wěn)定的供電。另外有學(xué)者利用單裂縫概念化模型研究采熱過程,即認(rèn)為熱儲層是由基質(zhì)巖體和水平單裂縫組成或?qū)⒘芽p簡化為一維線段嵌入巖石基質(zhì)網(wǎng)絡(luò)中,Yao 等[16]在三維熱流單裂縫模型基礎(chǔ)上,研究區(qū)域流動方向?qū)Σ蔁岬挠绊?,結(jié)果表明,采熱可持續(xù)性隨著區(qū)域流動方向和2 井中心連線的夾角的增加而增加。Li 等[17]在THM 耦合單裂縫模型基礎(chǔ)上進行開發(fā)過程研究,發(fā)現(xiàn)高滲透率或大孔徑會導(dǎo)致過早發(fā)生熱突破,因此要在流動和傳熱之間進行權(quán)衡。然而,在這些模型中單一或者規(guī)則的裂縫分布過于理想,不能體現(xiàn)儲層中裂縫的真實情況;隨機裂縫相對靈活,但自由變量過多、模型數(shù)據(jù)量大、試驗數(shù)據(jù)缺乏。

    盡管已經(jīng)開展了大量干熱巖EGS 多物理場耦合的數(shù)值模擬工作,但由于改造后的滲透率通常極為復(fù)雜,常規(guī)方法難以描述,因此準(zhǔn)確地估計儲層中的熱量提取過程仍然存在大量挑戰(zhàn)。近年來,Xing 等[8-9]利用微震數(shù)據(jù)計算出地?zé)醿拥臐B透率分布,為干熱巖開采設(shè)計了8 井布局,模擬結(jié)果展示了地?zé)醿又辛鲌龊蜏囟葓鲭S時間的演化。李庭樑等[15]結(jié)合地質(zhì)背景和壓裂數(shù)據(jù),計算出滲透率等水動力學(xué)參數(shù)的分布,并對熱儲進行數(shù)值模擬,將其結(jié)果與分塊均質(zhì)熱儲進行比較,結(jié)果表明熱儲的非均質(zhì)性會導(dǎo)致EGS 采熱性能顯著下降。Xu 等[18]利用水力壓裂過程中記錄的微震事件生成了裂縫模型,并將該方法應(yīng)用于Habanero 儲層,證明了 Habanero 井之間的連通性。Fang 等[19]分別利用基于立方定律的等效多孔介質(zhì)方法和離散裂縫網(wǎng)絡(luò)方法估算壓裂改造后的儲層滲透率,結(jié)果表明儲層中滲透率最高的區(qū)域即微震最密集的區(qū)域。以上研究表明水力壓裂誘發(fā)微地震的分布可以較為準(zhǔn)確地表征儲層改造后的滲透率,微地震分布優(yōu)勢的方向也即儲層中流體流動的優(yōu)勢路徑。目前利用微地震數(shù)據(jù)評估滲透率也存在不足,如網(wǎng)格分布和震源精度的限制等??赏ㄟ^優(yōu)化網(wǎng)格劃分、改善定位方法等提高滲透率評估的準(zhǔn)確性。

    基于上述研究,本文以澳大利亞庫珀盆地的Habanero 熱儲為研究區(qū),創(chuàng)新性地構(gòu)建了微地震監(jiān)測數(shù)據(jù)與儲層滲透率之間的等效關(guān)系并建立了不同的熱儲模型,包括分區(qū)均質(zhì)模型和非均質(zhì)模型;利用自主研發(fā)的多物理場耦合軟件PANDAS 對不同模型中滲流場、溫度場和生產(chǎn)井井底溫度的變化等進行了數(shù)值模擬、分析比較,探究采熱性能的差異。

    1 理論與方法

    本文使用的模擬工具為自主研發(fā)的PANDAS,這是一種基于有限元和格子玻爾茲曼方法(Lattice Boltzmann Method,LBM)的數(shù)值模擬程序;其中熱-流動耦合數(shù)值模擬模塊[20]將應(yīng)用于Habanero 地?zé)醿娱_發(fā)中的有限元模擬及評估研究。該模型基本假設(shè)為:

    (1)干熱巖儲層中充滿單相的飽和流體水,且流體在儲層中的流動遵循達西定律;

    (2)流體和巖石間的換熱通過對流和傳導(dǎo)實現(xiàn);

    (3)不考慮流體與礦物的物理化學(xué)反應(yīng)。

    1.1 連續(xù)性方程

    流體在多孔介質(zhì)中的流動可以表示為:

    式中: φ—孔隙度;

    βm—綜合壓縮系數(shù)/Pa-1;

    βf、βr—流體和基巖的壓縮系數(shù)/Pa-1;

    P—流體壓力/Pa;

    ρ—流體密度/(kg·m-3);

    v—流體流速/(m·s-1)。

    遵循達西定律的單相流體的流動控制方程為:

    式中:k—多孔介質(zhì)的滲透率/m2;

    μ—流體的動力黏度/(Pa·s);

    g—重力加速度/(m·s-2);

    D—深度/m。

    1.2 多孔介質(zhì)中的對流-傳導(dǎo)換熱方程

    多孔介質(zhì)中基質(zhì)和流體能量平衡方程分別為:

    式中:C—比熱/(J·kg-1·K-1);

    λ—導(dǎo)熱系數(shù)/(W·m-1·K-1);

    T—溫度/K;

    r—巖石基質(zhì);

    f—流體。

    從宏觀角度出發(fā),假定EGS 各階段中局部處于熱平衡狀態(tài),即Tr=Tf=T。聯(lián)立式(4)和式(5),多孔介質(zhì)中的傳導(dǎo)-對流換熱控制方程為:

    式中:m—基質(zhì)和流體的混合系統(tǒng)。

    2 研究區(qū)概況

    庫珀盆地位于南澳大利亞東北部,如圖1 所示。據(jù)資料顯示,庫珀盆地下方3 500~4 500 m 范圍內(nèi)埋藏有大量的高溫花崗巖,溫度約為250 °C[21],是其他地區(qū)相同深度花崗巖溫度的2 倍?;◢弾r巖體上方覆蓋著約3 500 m 厚的沉積蓋層,其主要由低孔、低滲的碎屑巖和煤系夾層組成[12],是花崗巖基底天然的“絕熱屏障”。

    圖1 庫珀盆地?zé)醿友芯繀^(qū)域(修改自文獻[12,24])Fig.1 Geothermal reservoir region to be studied in Cooper basin (modified from [12, 24])

    Habanero 3 號井(H03)位于1 號井東北方向約560 m處,井深約4 221 m。在2008年對3 號井進行了水力壓裂[22]。并對1 號井和3 號井進行了時長為42 d 的閉環(huán)流動測試。實測數(shù)據(jù)表明,1 號井井口壓力為44.5 MPa,3 號井井口壓力為33.5 MPa,注采流量為15.5 kg/s。

    水力壓裂誘發(fā)的微地震的空間分布表明存在呈近水平結(jié)構(gòu)的Habanero 斷層[23]。測井資料表明,該斷層是構(gòu)造成因,水力壓裂重新激活并擴展了之前存在的裂縫,擴大了儲層高滲區(qū)域[25]。本文將以庫珀盆地現(xiàn)場水力壓裂及相關(guān)聯(lián)通試驗為基礎(chǔ),開展數(shù)值模擬評估研究。

    3 Habanero EGS 數(shù)值模型

    3.1 三維數(shù)值模型

    基于Habanero 地?zé)醿拥牡刭|(zhì)構(gòu)造特征及水力壓裂改造效果,建立了三維數(shù)值模型。該模型在x,y和z方向的模型尺寸分別為1 200,1 500,500 m。模型在水力壓裂區(qū)劃分較密,網(wǎng)格尺寸為2.5 m,其他區(qū)域的網(wǎng)格尺寸由壓裂區(qū)向外逐漸變大,在15~20 m 之間。模型計算深度為地下4 000~4 500 m,如圖2(a)所示。

    圖2 分區(qū)均質(zhì)模型和剖面圖Fig.2 Homogeneous model (a) and homogeneous model profile view (b)

    3.1.1 分區(qū)均質(zhì)模型

    根據(jù)庫珀盆地現(xiàn)場資料及研究成果,將研究區(qū)簡化為從上至下分為3 層的三維模型,如圖2(b)所示。上層的深度為4 000~4 245 m,為未壓裂的致密花崗巖;中間層深度4 245~4 255 m,厚度為10 m,為孔隙度和滲透率較高的改造區(qū);下層的深度為4 255~4 500 m,為未壓裂的致密花崗巖。

    3.1.2 基于微地震數(shù)據(jù)的非均質(zhì)模型

    微震可以反映改造過程中巖石破裂發(fā)生的位置及擴展范圍,故本文利用已知地質(zhì)信息、微震監(jiān)測數(shù)據(jù)構(gòu)建Habanero 非均質(zhì)熱儲模型。本文的數(shù)據(jù)為Geodynamics 公司在Habanero 1 號井水力壓裂過程中記錄的誘發(fā)微震事件。研究表明,微震點的密集程度可以用來表征儲層內(nèi)的滲透率分布,一般認(rèn)為,微震分布密集區(qū)域的滲透率要高于微震分布稀疏區(qū)域[15]。因此,本文基于Habanero 熱儲區(qū)巖石破裂的機理,根據(jù)微震點到模型中各個節(jié)點的距離建立了三維非均質(zhì)儲層滲透率模型。

    基于微地震數(shù)據(jù)的非均質(zhì)模型使用的網(wǎng)格模型與分區(qū)均質(zhì)模型使用的網(wǎng)格相同。在本文的計算方法中,滲透率是網(wǎng)格節(jié)點到微震點距離的指數(shù)函數(shù),各個節(jié)點處的滲透率由2 部分組成:第1 部分為初始滲透率,第2 部分通過計算各個微震事件到該節(jié)點的距離進行評估。具體計算方法引自Xing 等[8]。在有限元計算中對微震事件和滲透率計算節(jié)點的加權(quán)距離評估:

    式中:K—節(jié)點處的滲透率/m2;

    K0—初始滲透率/m2;

    A—震級對滲透率的影響系數(shù)(此處A=1.0×10-15);

    dx、dy、dz—考慮破裂面傾角影響的微震事件與網(wǎng)格節(jié)點的距離;

    B1、B2、B3—微震事件在x、y、z方向上影響范圍的參數(shù)。

    3.2 模型參數(shù)

    根據(jù)測井?dāng)?shù)據(jù),庫珀盆地在3 500~5 000 m 深度范圍內(nèi),初始溫度和壓力隨深度而變化,溫度梯度為31 °C/km,壓力梯度為8.6 MPa/km。將模型上下表面的壓力和溫度設(shè)置為恒定,上表面壓力為71.8 MPa,溫度為240.1 °C,下表面壓力為76.1 MPa,溫度為255.6 °C。根據(jù)Llanos 等[12]和O’Sullivan 等[26]的研究,模型基于所有側(cè)面沒有熱量或質(zhì)量的傳遞建立,因此模型四周設(shè)置為封閉邊界。

    分區(qū)均質(zhì)模型中除滲透率以外,其他的巖石物理性質(zhì)參考Llanos 等[12]的研究,如表1 所示,在基于微地震數(shù)據(jù)的非均質(zhì)模型中,滲透率通過微震監(jiān)測數(shù)據(jù)進行計算,其他巖石物理性質(zhì)參數(shù)與分區(qū)均質(zhì)模型保持一致。

    表1 分區(qū)均質(zhì)模型中巖石物理性質(zhì)參數(shù)Table 1 Parameters of rock in the zonal homogenization model

    在2 個模型中,均使用水作為循環(huán)工作流體。在模型初始狀態(tài)時,假設(shè)儲層中水的溫度與其所處的圍巖溫度相同,注入工作流體溫度為90 ℃。表2 為循環(huán)工作流體的物理性質(zhì)參數(shù)。

    表2 循環(huán)工作流體的物理性質(zhì)參數(shù)Table 2 Parameters of the circulating fluid

    3.3 模型滲透率的確定

    利用2008年8月進行的Habanero 1 號井(H01,注入井)與Habanero 3 號井(H03,生產(chǎn)井)(圖3)的閉環(huán)流動現(xiàn)場測試數(shù)據(jù),對研究區(qū)滲透率進行校正。以現(xiàn)場實測壓力、流量數(shù)據(jù)為基礎(chǔ),確定模型中的井底壓力,調(diào)整模型中改造區(qū)的滲透率,并對比模擬流量與實際流量。

    圖3 模型中H01、H03 位置Fig.3 Locations of H01 and H03 in the model

    3.3.1 分區(qū)均質(zhì)模型滲透率

    Habanero 地?zé)醿邮侵旅芑◢弾r體,地層初始滲透率低,約為1.0×10-20m2。根據(jù)水力壓裂產(chǎn)生的微震點的空間分布,結(jié)合Llanos 等[12]的研究,Habanero 地?zé)醿訅毫押鬂B透率存在各向異性,y方向和x方向的滲透率之比為2∶1,z方向滲透率較低,在此基礎(chǔ)上利用實測流量對滲透率進行調(diào)整。表3 為在不同滲透率下,生產(chǎn)井中的模擬流量與實際流量的對比,由此可得,壓裂后改造區(qū)滲透率的最佳擬合值為x方向1.3×10-13m2,y方向2.6×10-13m2,z方向4.0×10-14m2,10 m 厚的改造區(qū)以外滲透率均為1.0×10-20m2。

    表3 分區(qū)均質(zhì)模型不同滲透率下生產(chǎn)井中的模擬流量與實際流量Table 3 Simulated flow and actual flow of production wells in the homogeneous model with different permeabilities

    3.3.2 非均質(zhì)模型滲透率

    在非均質(zhì)模型中,通過調(diào)整相關(guān)參數(shù),得到不同滲透率下生產(chǎn)井中的模擬流量(表4),并與實際流量進行對比。模型中滲透率的最佳擬合相關(guān)參數(shù)為B1=5.0×10-5、B2=2.0×10-5、B3=2.5(B1、B2、B3是微震事件在x、y、z方向上影響范圍的權(quán)重參數(shù),庫珀盆地水力壓裂微震事件以水平分布為主,因此設(shè)置Habanero地?zé)醿訅毫研Ч麑?jié)點至微地震點x、y方向的距離較為敏感)。該情況下,滲透率分布見圖4 和圖5,圖4 顯示整體模型的滲透率的分布,圖5 為滲透率分布的等值面圖,主要顯示水力壓裂區(qū)的滲透率分布。

    表4 非均質(zhì)模型不同滲透率相關(guān)參數(shù)下生產(chǎn)井中模擬流量與實際流量Table 4 Simulated flow and actual flow of production wells in the heterogeneous model with different permeability related parameters

    圖4 (a)非均質(zhì)模型的滲透率分布圖和(b)滲透率剖面圖Fig.4 Permeability distribution for the heterogeneous model (a)and permeability profile (b)

    圖5 (a)滲透率分布等值面圖和(b)滲透率等值面的剖面圖Fig.5 Contour surface view of permeability distribution (a) and isosurface profile of permeability distribution (b)

    4 結(jié)果與討論

    4.1 流場

    圖6(a)(d)是2 種模型中深度為4 250 m 截面的流速分布,可以看出,2 種滲透率模型中流速較大的是2 口井之間的連線及其附近的區(qū)域,說明儲層內(nèi)流體流動的優(yōu)勢路徑為兩井之間的連通區(qū)域。對比2 種滲透率模型的平面流速,分區(qū)均質(zhì)模型中該深度截面流速較大區(qū)域的范圍要比非均質(zhì)模型更大。原因在于2 種模型滲透率的不同,分區(qū)均質(zhì)模型10 m 厚改造區(qū)內(nèi)是具有各向異性的均一分布滲透率,x方向1.30×10-13m2,y方向2.60×10-13m2,z方向4.0×10-14m2。非均質(zhì)模型中滲透率由微震事件及破裂面決定,如圖7 所示,注入井附近滲透率最大,約為9.4 ×10-13m2,從注入井向生產(chǎn)井方向滲透率逐漸變小,在生產(chǎn)井附近滲透率約為1.0×10-14m2,這是因為微震事件發(fā)生的集中區(qū)域為注入井附近,并由此向外延伸,在微震事件發(fā)生較多的區(qū)域滲透率較高,儲層中2 井連線及其附近的區(qū)域表現(xiàn)出流動優(yōu)勢,遠(yuǎn)離注入井的區(qū)域由于微震事件較少,所以裂隙較少、滲透率較低,流速也小。非均質(zhì)模型中除注入井附近,其他位置滲透率要比分區(qū)均質(zhì)模型中低,因此分區(qū)均質(zhì)模型井間流速大于非均質(zhì)模型。

    圖6 不同模型的流速分布(上為分區(qū)均質(zhì)模型,下為非均質(zhì)模型)Fig.6 Velocity distribution of different models (the upper is the zonal homogeneous model and the bottom is the heterogeneous model)

    圖7 非均質(zhì)模型中z=4 250 m 剖面滲透率分布Fig.7 Permeability distribution at the section z=4 250 m in the heterogeneous model

    由圖6(b)(e)可知,在垂直方向上,分區(qū)均質(zhì)模型的流動區(qū)域限于10 m 厚的改造區(qū)內(nèi),而非均質(zhì)模型在深度方向上流體流動區(qū)域集中于該區(qū)域,但不限于此,這是因為微震事件在垂向上有較大延伸。

    4.2 溫度場

    溫度分布等值面圖,范圍為90~230 °C(圖8),總的來看,在初始運行階段,冷流體自注入井進入壓裂改造區(qū),由于與周圍高溫巖體進行了熱交換,注入井附近巖體溫度接近注入溫度90 °C。隨著注入流體的增多,低溫區(qū)(90~115 °C)逐漸擴大,并沿優(yōu)勢流動通道向生產(chǎn)井?dāng)U散。

    圖8 不同模型溫度隨時間變化的等值面圖Fig.8 Isosurface of temperature variation with time in different models

    對比來看,分區(qū)均質(zhì)模型中的冷流體從注入井向改造區(qū)中流動擴散,逐漸形成一個近似扁橢圓體的低溫區(qū),近y方向上的低溫區(qū)影響范圍大于近x方向的影響范圍。在深度方向上,低溫區(qū)范圍基本限制在深度4 245~4 255 m 的10 m 厚改造區(qū)內(nèi)。

    在非均質(zhì)模型中,注入井兩側(cè)冷流體流動的影響范圍與均質(zhì)模型相比,冷流體向著生產(chǎn)井流動的趨勢更加明顯。近y方向上的低溫區(qū)影響范圍同樣大于近x方向的影響范圍。但深度方向?qū)Ρ葋砜矗蜏貐^(qū)影響范圍不再局限于10 m 厚的范圍。這是因為水力壓裂形成的微震事件在深度上的分布范圍較廣,微震發(fā)生的區(qū)域內(nèi)儲層滲透率也隨之變大。

    4.3 采熱性能

    2 種模型的生產(chǎn)井井底溫度在30 a 間的變化情況見圖9(a)。分區(qū)均質(zhì)模型30 a 間溫度變化趨勢為前期(0~7.9 a)保持不變,之后(7.9~30 a)緩慢下降,非均質(zhì)模型30 a 間溫度變化趨勢為前期(0~8.5 a)保持不變,之后(8.5~30 a)快速下降,且下降速度逐漸變慢。系統(tǒng)運行到30 a,分區(qū)均質(zhì)模型溫度下降到220.9 °C,下降了26.9 °C,非均質(zhì)模型溫度下降到170.0 °C,下降了77.8 °C。

    圖9 分區(qū)均質(zhì)和非均質(zhì)模型生產(chǎn)井井底溫度和采熱量對比Fig.9 Comparison of the bottom hole temperature of production wells and heat extraction comparison in the two models

    對比來看,2 種模型生產(chǎn)井初始溫度保持不變的時間十分接近;后期分區(qū)均質(zhì)模型溫度下降少,非均質(zhì)模型下降多。由此分析,分區(qū)均質(zhì)模型略早下降0.6 a,即熱突破時間稍微早一些,這是分區(qū)均質(zhì)模型生產(chǎn)井附近滲透率比非均質(zhì)模型大導(dǎo)致的。非均質(zhì)模型溫度下降幅度比分區(qū)均質(zhì)模型大,下降速度比非均質(zhì)模型快,這是由于非均質(zhì)模型的高滲區(qū)集中于2 口井連通區(qū)域內(nèi),使得非均質(zhì)模型中優(yōu)勢流動通道顯著,所以溫度下降更快。

    2 種模型單位質(zhì)量流體采熱量的對比關(guān)系,即采出熱量與采出質(zhì)量的比值[27],見圖9(b)。分區(qū)均質(zhì)模型和非均質(zhì)模型在儲層模型運行初期單位質(zhì)量流體采熱量比較接近,約為1.04 MJ/kg,分區(qū)均質(zhì)模型初始時間略高。之后隨著生產(chǎn)井溫度的快速下降,2 種模型的單位質(zhì)量流體采熱量也迅速降低,對比來看,兩種模型單位質(zhì)量流體采熱量下降發(fā)生的時間接近,分區(qū)均質(zhì)模型采熱量下降速度很慢,在模型計算到30 a時,單位質(zhì)量流體采熱量降低至1.003 MJ/kg,降低了3.74%;非均質(zhì)模型采熱量下降速度較快,在模型計算到30 a 時,采熱量降低至0.906 MJ/kg,降低了12.72%。從采熱量對比看,分區(qū)均質(zhì)模型在熱儲層運行期間采熱量變化相對更穩(wěn)定,非均質(zhì)模型采熱量降幅較大。

    對比生產(chǎn)井溫度曲線和采熱量曲線,分區(qū)均質(zhì)模型的模擬結(jié)果相比于非均質(zhì)模型結(jié)果來說,溫度下降幅度小、單位質(zhì)量流體采熱量高,但由微震分布可知實際儲層中的滲透率分布并不均勻,2 井之間優(yōu)勢流動通道效應(yīng)顯著,因此分區(qū)均質(zhì)模型的模擬單位質(zhì)量流體采熱量高于實際單位質(zhì)量流體采熱量。

    原始地?zé)醿訅毫押?,激活并擴展了先存裂縫,在水壓足夠高時產(chǎn)生新的裂縫,提高了儲層中尤其是近井區(qū)域的滲透率,形成優(yōu)勢流動通道,改善了儲層連通性。此外,水力壓裂后儲層滲透率具有顯著的滲透率各向異性,是影響干熱巖開發(fā)動態(tài)的重要因素之一[28]。相比分區(qū)均質(zhì)模型,非均質(zhì)模型中的優(yōu)勢流動通道、滲透率各向異性更加明顯,更加符合壓裂后的儲層實際情況,因此在實際應(yīng)用中,非均質(zhì)模型的模擬結(jié)果對實際工程更具參考意義。

    5 結(jié)論

    (1)水力壓裂誘發(fā)的微震事件集中于近井區(qū)域,導(dǎo)致該區(qū)滲透率增高,因此非均質(zhì)模型中優(yōu)勢流動通道顯著,冷流體從注入井更快流向生產(chǎn)井,溫度下降速度快;分區(qū)均質(zhì)模型中優(yōu)勢流動通道效應(yīng)較弱,溫度下降速度慢。

    (2)深度方向上,分區(qū)均質(zhì)模型僅限于10 m 厚的壓裂改造區(qū)內(nèi)流動;非均質(zhì)模型由于微震事件在深度方向上有較大延伸,流體流動區(qū)域集中于該層但不限于該范圍。

    (3)分區(qū)均質(zhì)模型在熱儲層開采模擬期間采熱量變化相對穩(wěn)定,降幅為3.74%;非均質(zhì)模型單位質(zhì)量流體采熱量降幅較大,為12.72%。

    (4)對比2 種模型的模擬結(jié)果,分區(qū)均質(zhì)模型的生產(chǎn)井溫度下降幅度小、單位質(zhì)量流體采熱量高。但由微震監(jiān)測數(shù)據(jù)可知,實際儲層中的滲透率分布不均,優(yōu)勢流動通道效應(yīng)顯著,分區(qū)均質(zhì)模型的模擬采熱量高于實際采熱量,因此非均質(zhì)模型的模擬結(jié)果在干熱巖資源實際開發(fā)過程中更具指導(dǎo)意義。

    猜你喜歡
    生產(chǎn)井微震均質(zhì)
    淺談KJ768煤礦微震監(jiān)測系統(tǒng)的應(yīng)用
    新疆鋼鐵(2021年1期)2021-10-14 08:45:48
    鄂爾多斯某區(qū)塊致密油產(chǎn)能影響因素分析與優(yōu)化研究
    化工管理(2021年23期)2021-08-25 03:01:00
    長平煤業(yè)5302 綜放工作面頂板巖層移動規(guī)律研究
    魅力中國(2020年46期)2020-02-02 11:50:08
    基于流線的聚驅(qū)雙高油田調(diào)整策略研究
    Orlicz對偶混合均質(zhì)積分
    基于隸屬度分析的回采面沖擊地壓微震能量閾值判定
    中國煤炭(2016年1期)2016-05-17 06:11:34
    EILog生產(chǎn)井技術(shù)平臺的五員“開路先鋒”
    石油知識(2016年2期)2016-02-28 16:20:13
    非均質(zhì)巖心調(diào)堵結(jié)合技術(shù)室內(nèi)實驗
    生產(chǎn)井實時檢測工藝技術(shù)及系統(tǒng)分析
    向斜軸部附近微震活動規(guī)律及對沖擊礦壓的影響
    金屬礦山(2013年5期)2013-03-11 16:53:46
    午夜福利在线观看免费完整高清在| 亚洲美女搞黄在线观看| 久久久久免费精品人妻一区二区| 国产免费一级a男人的天堂| 美女大奶头视频| 国产一区二区亚洲精品在线观看| 麻豆成人午夜福利视频| 亚洲在线观看片| 色网站视频免费| 日本色播在线视频| 国产黄色小视频在线观看| 三级毛片av免费| 亚洲人成网站高清观看| 亚洲电影在线观看av| 婷婷六月久久综合丁香| 亚洲欧美日韩卡通动漫| 久久草成人影院| 亚洲国产日韩欧美精品在线观看| 自拍偷自拍亚洲精品老妇| 日本猛色少妇xxxxx猛交久久| 免费观看的影片在线观看| 色尼玛亚洲综合影院| 少妇熟女aⅴ在线视频| 午夜日本视频在线| 亚洲欧洲日产国产| 人妻少妇偷人精品九色| 日韩电影二区| 天天躁夜夜躁狠狠久久av| 久久人人爽人人爽人人片va| 熟女人妻精品中文字幕| 国产激情偷乱视频一区二区| 亚洲精品视频女| 亚洲内射少妇av| 好男人在线观看高清免费视频| 免费黄频网站在线观看国产| 我的老师免费观看完整版| 亚洲一级一片aⅴ在线观看| 18禁在线无遮挡免费观看视频| 晚上一个人看的免费电影| 亚洲无线观看免费| 国产午夜精品久久久久久一区二区三区| 成人特级av手机在线观看| 成人性生交大片免费视频hd| 国产成人午夜福利电影在线观看| 天天躁夜夜躁狠狠久久av| 少妇的逼好多水| 五月天丁香电影| 久久久久久久亚洲中文字幕| 欧美精品一区二区大全| 亚洲欧洲日产国产| 亚洲精品成人av观看孕妇| 日韩av免费高清视频| 99久久精品国产国产毛片| 又爽又黄a免费视频| 成人亚洲精品一区在线观看 | 九九爱精品视频在线观看| 日本熟妇午夜| 精品一区二区免费观看| 国语对白做爰xxxⅹ性视频网站| 色尼玛亚洲综合影院| 老司机影院毛片| 国产成人91sexporn| 欧美xxxx黑人xx丫x性爽| 丝瓜视频免费看黄片| 国产成人精品福利久久| 亚洲精品成人久久久久久| 搡老妇女老女人老熟妇| 成人亚洲精品av一区二区| 九草在线视频观看| 亚洲国产精品专区欧美| 美女cb高潮喷水在线观看| 在线观看人妻少妇| 国产免费福利视频在线观看| 有码 亚洲区| 乱系列少妇在线播放| 日韩一区二区视频免费看| av网站免费在线观看视频 | 极品少妇高潮喷水抽搐| 精品国产露脸久久av麻豆 | 晚上一个人看的免费电影| 人妻制服诱惑在线中文字幕| www.色视频.com| 成人无遮挡网站| 国产综合懂色| 国产精品一区二区性色av| 高清av免费在线| 国产亚洲av嫩草精品影院| 真实男女啪啪啪动态图| 男的添女的下面高潮视频| 国产一区二区亚洲精品在线观看| 久久草成人影院| 少妇高潮的动态图| 国产成人免费观看mmmm| 好男人在线观看高清免费视频| 日韩一区二区视频免费看| 男人舔奶头视频| 搡老乐熟女国产| 天堂俺去俺来也www色官网 | 中国国产av一级| 99热这里只有是精品在线观看| 国产精品99久久久久久久久| www.色视频.com| 2021少妇久久久久久久久久久| 亚洲人成网站高清观看| 亚洲精品自拍成人| 校园人妻丝袜中文字幕| 99久久人妻综合| 九九在线视频观看精品| 九色成人免费人妻av| 亚洲精品第二区| 国产伦理片在线播放av一区| 日韩欧美精品v在线| 国精品久久久久久国模美| 99热这里只有精品一区| 日本-黄色视频高清免费观看| 国产男人的电影天堂91| 少妇熟女欧美另类| 伦理电影大哥的女人| 久久99热这里只有精品18| 婷婷色综合www| 成人综合一区亚洲| 神马国产精品三级电影在线观看| 国精品久久久久久国模美| 久久久久久久久久久丰满| 夜夜看夜夜爽夜夜摸| 777米奇影视久久| 熟女人妻精品中文字幕| 少妇熟女aⅴ在线视频| 免费av毛片视频| 午夜精品在线福利| 中国美白少妇内射xxxbb| eeuss影院久久| 国产精品久久视频播放| 日韩欧美一区视频在线观看 | 欧美 日韩 精品 国产| 看十八女毛片水多多多| 色吧在线观看| 非洲黑人性xxxx精品又粗又长| 搞女人的毛片| 中文精品一卡2卡3卡4更新| 午夜免费激情av| 成人一区二区视频在线观看| 看非洲黑人一级黄片| 亚洲国产最新在线播放| 国产伦精品一区二区三区四那| 汤姆久久久久久久影院中文字幕 | 国产午夜精品一二区理论片| 91av网一区二区| 大片免费播放器 马上看| 亚洲在线观看片| 久久99精品国语久久久| 欧美成人a在线观看| 十八禁国产超污无遮挡网站| 自拍偷自拍亚洲精品老妇| 亚洲国产精品成人综合色| 亚洲国产精品国产精品| 在线观看av片永久免费下载| 男女视频在线观看网站免费| 欧美潮喷喷水| 国内揄拍国产精品人妻在线| 精品欧美国产一区二区三| 国产成人午夜福利电影在线观看| 日韩不卡一区二区三区视频在线| 国产美女午夜福利| 久久久久国产网址| 亚洲av成人av| 午夜精品国产一区二区电影 | 日产精品乱码卡一卡2卡三| 啦啦啦中文免费视频观看日本| 精品人妻一区二区三区麻豆| 两个人的视频大全免费| 最近中文字幕高清免费大全6| 尤物成人国产欧美一区二区三区| 中文字幕av成人在线电影| 国产av国产精品国产| 又黄又爽又刺激的免费视频.| 男女边摸边吃奶| 国精品久久久久久国模美| 一级爰片在线观看| 亚洲精品日韩av片在线观看| 亚洲欧美中文字幕日韩二区| 国产白丝娇喘喷水9色精品| 五月玫瑰六月丁香| 麻豆乱淫一区二区| 亚洲成人av在线免费| 国产黄a三级三级三级人| 亚洲精品456在线播放app| 九九在线视频观看精品| 成人一区二区视频在线观看| 国产精品久久久久久精品电影| 国产人妻一区二区三区在| 特级一级黄色大片| 五月玫瑰六月丁香| 狠狠精品人妻久久久久久综合| 亚洲第一区二区三区不卡| 搡老乐熟女国产| 伊人久久国产一区二区| 久久99精品国语久久久| 亚洲av免费高清在线观看| 国产精品99久久久久久久久| 国产精品三级大全| 九九爱精品视频在线观看| 久久久久久九九精品二区国产| 久久97久久精品| 亚洲三级黄色毛片| 久久精品久久久久久噜噜老黄| 久久久久网色| 国产成人精品福利久久| 亚洲精品国产成人久久av| 九九爱精品视频在线观看| 天天躁夜夜躁狠狠久久av| 在线观看免费高清a一片| 久久精品熟女亚洲av麻豆精品 | 国产精品嫩草影院av在线观看| 色尼玛亚洲综合影院| 国产精品综合久久久久久久免费| 日韩大片免费观看网站| 老司机影院成人| 伊人久久精品亚洲午夜| 夫妻午夜视频| 精品久久久久久久久久久久久| 亚洲av中文av极速乱| 欧美xxxx黑人xx丫x性爽| 久久久欧美国产精品| 日韩制服骚丝袜av| 美女高潮的动态| 日本黄大片高清| 天美传媒精品一区二区| 美女被艹到高潮喷水动态| 少妇被粗大猛烈的视频| 免费黄色在线免费观看| 精品不卡国产一区二区三区| 日日摸夜夜添夜夜爱| 青春草视频在线免费观看| 精品久久久久久久久亚洲| 老司机影院成人| 欧美不卡视频在线免费观看| 亚洲av一区综合| 免费少妇av软件| 亚洲欧美精品专区久久| 男人爽女人下面视频在线观看| 日本熟妇午夜| 美女xxoo啪啪120秒动态图| 十八禁国产超污无遮挡网站| 亚洲国产av新网站| 成人二区视频| 只有这里有精品99| 成人性生交大片免费视频hd| 老司机影院毛片| 青青草视频在线视频观看| 夜夜爽夜夜爽视频| 高清日韩中文字幕在线| 婷婷色av中文字幕| 91午夜精品亚洲一区二区三区| 亚洲婷婷狠狠爱综合网| 精品少妇黑人巨大在线播放| 一级a做视频免费观看| 欧美成人午夜免费资源| 国产午夜精品久久久久久一区二区三区| 久久久色成人| 欧美区成人在线视频| 校园人妻丝袜中文字幕| 赤兔流量卡办理| 免费电影在线观看免费观看| 欧美高清性xxxxhd video| 亚洲精品视频女| 亚洲四区av| 久久久久久久久久久免费av| 午夜亚洲福利在线播放| 午夜免费男女啪啪视频观看| 中文精品一卡2卡3卡4更新| 九九爱精品视频在线观看| 亚洲真实伦在线观看| 成年人午夜在线观看视频 | 视频中文字幕在线观看| a级一级毛片免费在线观看| 美女黄网站色视频| 亚洲图色成人| 男人舔女人下体高潮全视频| 久久99热6这里只有精品| 两个人的视频大全免费| 亚洲最大成人手机在线| 国产高清不卡午夜福利| 99九九线精品视频在线观看视频| 欧美高清成人免费视频www| 美女黄网站色视频| 亚洲天堂国产精品一区在线| 国产精品久久久久久精品电影| 国产精品久久久久久久久免| 2021天堂中文幕一二区在线观| 在线 av 中文字幕| 国产男女超爽视频在线观看| 日韩av在线大香蕉| www.av在线官网国产| 中国国产av一级| 久久久久久久大尺度免费视频| 波野结衣二区三区在线| 欧美激情在线99| 极品少妇高潮喷水抽搐| 亚洲av国产av综合av卡| 成年人午夜在线观看视频 | 丝袜美腿在线中文| 久久精品国产亚洲av涩爱| 激情五月婷婷亚洲| 美女被艹到高潮喷水动态| 美女黄网站色视频| 国产大屁股一区二区在线视频| 久久精品夜夜夜夜夜久久蜜豆| 天堂俺去俺来也www色官网 | 成人性生交大片免费视频hd| 国产免费一级a男人的天堂| 18禁裸乳无遮挡免费网站照片| 免费大片黄手机在线观看| 青春草亚洲视频在线观看| 亚洲三级黄色毛片| 2018国产大陆天天弄谢| 精品国内亚洲2022精品成人| 777米奇影视久久| 国产大屁股一区二区在线视频| 精品国产三级普通话版| 中文字幕制服av| 亚洲国产精品sss在线观看| 国产乱来视频区| 国产精品综合久久久久久久免费| 午夜爱爱视频在线播放| 1000部很黄的大片| 五月天丁香电影| 你懂的网址亚洲精品在线观看| 亚洲精品影视一区二区三区av| 高清视频免费观看一区二区 | 国产精品.久久久| 视频中文字幕在线观看| 日韩国内少妇激情av| 美女xxoo啪啪120秒动态图| 国产又色又爽无遮挡免| 国产精品爽爽va在线观看网站| 国产欧美日韩精品一区二区| 精品久久久久久久久av| av专区在线播放| 国产国拍精品亚洲av在线观看| 国产综合懂色| 国产精品人妻久久久影院| av黄色大香蕉| 国产精品av视频在线免费观看| 亚洲国产成人一精品久久久| 成人国产麻豆网| 狠狠精品人妻久久久久久综合| 免费观看a级毛片全部| 国产亚洲一区二区精品| 深夜a级毛片| 建设人人有责人人尽责人人享有的 | 国语对白做爰xxxⅹ性视频网站| 老司机影院毛片| 久久久国产一区二区| 国产色爽女视频免费观看| 久久久久久伊人网av| 别揉我奶头 嗯啊视频| 亚洲av成人精品一二三区| 特大巨黑吊av在线直播| 亚洲欧洲国产日韩| 99久久人妻综合| 亚洲国产av新网站| 久久久a久久爽久久v久久| 七月丁香在线播放| 十八禁国产超污无遮挡网站| 国产亚洲91精品色在线| 免费人成在线观看视频色| 久久久国产一区二区| 联通29元200g的流量卡| 精品久久国产蜜桃| 亚洲精华国产精华液的使用体验| 亚洲av电影在线观看一区二区三区 | 国产午夜精品论理片| 欧美性猛交╳xxx乱大交人| 尾随美女入室| 搡老妇女老女人老熟妇| 久久久久久久久久久免费av| 一区二区三区乱码不卡18| 一级毛片黄色毛片免费观看视频| av卡一久久| 夫妻性生交免费视频一级片| 日韩欧美精品免费久久| 大话2 男鬼变身卡| 狂野欧美白嫩少妇大欣赏| 亚洲av国产av综合av卡| 国产单亲对白刺激| 搡女人真爽免费视频火全软件| 精品久久久久久成人av| 久久久久久久久大av| 老师上课跳d突然被开到最大视频| 美女黄网站色视频| 久久精品熟女亚洲av麻豆精品 | 丝袜喷水一区| 久久人人爽人人爽人人片va| 午夜免费男女啪啪视频观看| 视频中文字幕在线观看| 三级经典国产精品| 18禁裸乳无遮挡免费网站照片| 亚洲电影在线观看av| 老司机影院成人| ponron亚洲| 久久人人爽人人片av| 人人妻人人澡欧美一区二区| 欧美丝袜亚洲另类| 91av网一区二区| 欧美极品一区二区三区四区| 国产精品一及| 中文字幕久久专区| 乱人视频在线观看| 国产精品一区二区三区四区免费观看| 一级毛片我不卡| 欧美日韩一区二区视频在线观看视频在线 | av在线播放精品| 国产黄片美女视频| 久久午夜福利片| 亚洲国产高清在线一区二区三| 男人和女人高潮做爰伦理| 边亲边吃奶的免费视频| 国产亚洲一区二区精品| 又大又黄又爽视频免费| 在线观看人妻少妇| 亚洲国产欧美人成| 搞女人的毛片| 欧美成人一区二区免费高清观看| 精品熟女少妇av免费看| 人妻少妇偷人精品九色| 91aial.com中文字幕在线观看| 国产精品蜜桃在线观看| 色吧在线观看| 久久这里有精品视频免费| 嫩草影院精品99| 伊人久久精品亚洲午夜| 精品久久久精品久久久| 永久免费av网站大全| 99热全是精品| 看非洲黑人一级黄片| av黄色大香蕉| 精品久久国产蜜桃| 一个人看视频在线观看www免费| 午夜福利视频精品| 国内精品美女久久久久久| 国产麻豆成人av免费视频| 免费看a级黄色片| 一级av片app| 亚洲欧美成人综合另类久久久| 最近最新中文字幕大全电影3| 日本猛色少妇xxxxx猛交久久| 中文资源天堂在线| 亚洲国产精品成人综合色| 夜夜看夜夜爽夜夜摸| 久99久视频精品免费| 91久久精品国产一区二区成人| 一个人看的www免费观看视频| 一本久久精品| 中文字幕av成人在线电影| 成人鲁丝片一二三区免费| 少妇丰满av| 丝袜喷水一区| 日韩强制内射视频| 少妇高潮的动态图| 男插女下体视频免费在线播放| 国产高潮美女av| 日本黄色片子视频| 亚洲国产精品sss在线观看| 国产乱来视频区| 亚洲av国产av综合av卡| 中文字幕免费在线视频6| 美女主播在线视频| 晚上一个人看的免费电影| 久久国产乱子免费精品| 国产综合懂色| 噜噜噜噜噜久久久久久91| 不卡视频在线观看欧美| 婷婷色综合www| 日韩av在线免费看完整版不卡| a级毛色黄片| 少妇丰满av| 国产成人精品一,二区| a级毛片免费高清观看在线播放| 国产精品99久久久久久久久| 亚洲综合精品二区| 男人舔奶头视频| 18禁动态无遮挡网站| 丰满人妻一区二区三区视频av| 国产在线男女| 国产成人a∨麻豆精品| 身体一侧抽搐| 国产单亲对白刺激| 亚洲精品乱久久久久久| 尾随美女入室| 一级毛片我不卡| 久久久亚洲精品成人影院| 久久草成人影院| 丰满人妻一区二区三区视频av| 亚洲人成网站在线观看播放| 国精品久久久久久国模美| 天堂中文最新版在线下载 | 内地一区二区视频在线| 日韩av免费高清视频| 搞女人的毛片| 亚洲四区av| 久久99蜜桃精品久久| 尤物成人国产欧美一区二区三区| 激情五月婷婷亚洲| 国产白丝娇喘喷水9色精品| 中文乱码字字幕精品一区二区三区 | 国内精品宾馆在线| 黄片无遮挡物在线观看| 精品久久久久久久久久久久久| 久久精品综合一区二区三区| www.av在线官网国产| 老司机影院成人| 深夜a级毛片| av在线老鸭窝| 一级二级三级毛片免费看| 夫妻性生交免费视频一级片| 亚洲精品第二区| 亚洲怡红院男人天堂| 联通29元200g的流量卡| 两个人视频免费观看高清| 日本黄大片高清| 男人狂女人下面高潮的视频| 五月玫瑰六月丁香| 日本熟妇午夜| 精品一区二区三区视频在线| 日本色播在线视频| 大又大粗又爽又黄少妇毛片口| 欧美日韩视频高清一区二区三区二| 久久久久久久久久人人人人人人| 免费在线观看成人毛片| 亚洲伊人久久精品综合| 在线播放无遮挡| 好男人视频免费观看在线| 日韩精品青青久久久久久| 99久久人妻综合| 精品一区在线观看国产| av天堂中文字幕网| 国产精品国产三级国产专区5o| 国产成人91sexporn| 97超碰精品成人国产| 最后的刺客免费高清国语| 亚洲真实伦在线观看| 老师上课跳d突然被开到最大视频| 听说在线观看完整版免费高清| 波野结衣二区三区在线| 99热全是精品| 丰满人妻一区二区三区视频av| 麻豆成人av视频| 2021少妇久久久久久久久久久| 国产 一区精品| 99热这里只有是精品50| 我要看日韩黄色一级片| 美女脱内裤让男人舔精品视频| 久久久久性生活片| 亚洲国产av新网站| 日韩成人伦理影院| 最近手机中文字幕大全| 日本与韩国留学比较| 国产 一区 欧美 日韩| 亚洲av一区综合| 国内精品美女久久久久久| 久久久久久久大尺度免费视频| 精品久久国产蜜桃| 久久久久久久久久人人人人人人| 国产视频内射| 亚洲精品色激情综合| 不卡视频在线观看欧美| 亚洲欧美中文字幕日韩二区| 亚洲美女视频黄频| 中文字幕av在线有码专区| 国产精品久久视频播放| 日本与韩国留学比较| 国产成人精品一,二区| 国产av国产精品国产| 亚洲av福利一区| 伦理电影大哥的女人| 亚洲精品第二区| 美女黄网站色视频| 少妇熟女aⅴ在线视频| 国产美女午夜福利| 亚洲av.av天堂| 久久精品国产亚洲av涩爱| 简卡轻食公司| 亚洲综合色惰| 亚洲在久久综合| 国产精品一区www在线观看| 日韩制服骚丝袜av| 国产成人一区二区在线| 久久久色成人| 婷婷六月久久综合丁香| av在线亚洲专区| 精品一区二区三区视频在线| 亚洲怡红院男人天堂| av在线亚洲专区| 男人爽女人下面视频在线观看| 婷婷六月久久综合丁香| 国产成人一区二区在线| 在线免费十八禁| 99热网站在线观看| 不卡视频在线观看欧美| 狂野欧美白嫩少妇大欣赏| 国产精品国产三级国产av玫瑰| 黄色配什么色好看| 免费在线观看成人毛片| 亚洲av福利一区| 黄色日韩在线| 欧美丝袜亚洲另类| 一级片'在线观看视频| 亚洲经典国产精华液单| 美女高潮的动态| 精品久久国产蜜桃| 在线天堂最新版资源| 午夜精品国产一区二区电影 | 精品久久久久久久末码| 只有这里有精品99|