• <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
    久久中文字幕人妻熟女| 国产一卡二卡三卡精品| 国产成人av激情在线播放| 男女之事视频高清在线观看| 少妇 在线观看| 国产熟女午夜一区二区三区| 日韩精品青青久久久久久| 欧美日韩亚洲国产一区二区在线观看| 亚洲av美国av| netflix在线观看网站| 久久精品91无色码中文字幕| 涩涩av久久男人的天堂| 我的亚洲天堂| 丝袜美腿诱惑在线| 久久香蕉精品热| 亚洲成国产人片在线观看| 岛国视频午夜一区免费看| 亚洲国产精品合色在线| 在线观看日韩欧美| 男人舔女人下体高潮全视频| 最新美女视频免费是黄的| 亚洲av成人av| 美女国产高潮福利片在线看| 女人被狂操c到高潮| 淫秽高清视频在线观看| 午夜两性在线视频| 女警被强在线播放| 一区二区日韩欧美中文字幕| 99国产精品99久久久久| 后天国语完整版免费观看| 黄片大片在线免费观看| 19禁男女啪啪无遮挡网站| 国产精品久久久av美女十八| 亚洲av成人不卡在线观看播放网| 亚洲熟女毛片儿| 曰老女人黄片| 丝袜美足系列| 日日摸夜夜添夜夜添小说| 亚洲欧洲精品一区二区精品久久久| 一个人免费在线观看的高清视频| 成人国产一区最新在线观看| av电影中文网址| 精品久久久久久久久久免费视频 | 日韩欧美一区视频在线观看| 亚洲自偷自拍图片 自拍| 中出人妻视频一区二区| 久久亚洲真实| 一进一出好大好爽视频| 麻豆久久精品国产亚洲av | 久久狼人影院| 国产成人欧美| 国产精品乱码一区二三区的特点 | a级毛片在线看网站| 婷婷六月久久综合丁香| 韩国av一区二区三区四区| 精品久久蜜臀av无| 精品久久久久久久毛片微露脸| 91成人精品电影| 水蜜桃什么品种好| 欧美久久黑人一区二区| 在线观看日韩欧美| 满18在线观看网站| 宅男免费午夜| 欧美在线黄色| 亚洲成人久久性| 高清欧美精品videossex| 麻豆一二三区av精品| 成熟少妇高潮喷水视频| 日本vs欧美在线观看视频| 欧美日韩乱码在线| 亚洲,欧美精品.| 天天躁狠狠躁夜夜躁狠狠躁| 999久久久精品免费观看国产| 欧美成人午夜精品| 侵犯人妻中文字幕一二三四区| 久久久久国产精品人妻aⅴ院| 91国产中文字幕| 成人三级黄色视频| 看免费av毛片| 多毛熟女@视频| 久久精品国产99精品国产亚洲性色 | 他把我摸到了高潮在线观看| 国产精品综合久久久久久久免费 | 可以免费在线观看a视频的电影网站| aaaaa片日本免费| 老司机午夜十八禁免费视频| 国产精品秋霞免费鲁丝片| 69精品国产乱码久久久| 大香蕉久久成人网| 嫁个100分男人电影在线观看| 新久久久久国产一级毛片| 午夜91福利影院| 高清在线国产一区| 国产精品久久久久久人妻精品电影| 亚洲美女黄片视频| 久久精品亚洲熟妇少妇任你| 国产日韩一区二区三区精品不卡| 丝袜人妻中文字幕| 国产欧美日韩精品亚洲av| 精品一区二区三卡| 不卡av一区二区三区| 国产97色在线日韩免费| 亚洲精品国产精品久久久不卡| 99re在线观看精品视频| 国产免费男女视频| 欧美中文日本在线观看视频| 午夜福利在线免费观看网站| 丁香欧美五月| 手机成人av网站| 久久伊人香网站| 麻豆av在线久日| 亚洲中文av在线| 国产深夜福利视频在线观看| 亚洲色图综合在线观看| 人人妻人人澡人人看| 最近最新中文字幕大全免费视频| 国产主播在线观看一区二区| 激情在线观看视频在线高清| 国产国语露脸激情在线看| 天堂中文最新版在线下载| 少妇被粗大的猛进出69影院| 欧美日韩av久久| 欧美中文综合在线视频| 久久久久九九精品影院| 午夜福利在线免费观看网站| 97碰自拍视频| 午夜福利影视在线免费观看| 国产色视频综合| 国产欧美日韩精品亚洲av| 精品国产超薄肉色丝袜足j| 久久国产精品影院| av中文乱码字幕在线| 在线观看www视频免费| 脱女人内裤的视频| av网站在线播放免费| cao死你这个sao货| 久久天堂一区二区三区四区| 老汉色∧v一级毛片| 欧美精品亚洲一区二区| 成人影院久久| 成人手机av| 男女之事视频高清在线观看| 国产精品久久久人人做人人爽| 亚洲精品在线美女| 国产精品野战在线观看 | 91在线观看av| 麻豆av在线久日| 国产一区二区三区在线臀色熟女 | 99在线人妻在线中文字幕| 精品久久久精品久久久| 99在线视频只有这里精品首页| 国产三级黄色录像| 在线观看一区二区三区| 少妇 在线观看| 亚洲精品久久午夜乱码| 中出人妻视频一区二区| 丰满人妻熟妇乱又伦精品不卡| 日韩人妻精品一区2区三区| 亚洲精品美女久久久久99蜜臀| 精品久久久久久,| 国产精品1区2区在线观看.| 午夜激情av网站| 国产精品乱码一区二三区的特点 | 多毛熟女@视频| 亚洲欧洲精品一区二区精品久久久| 亚洲精品国产色婷婷电影| 国产片内射在线| 50天的宝宝边吃奶边哭怎么回事| 日日爽夜夜爽网站| 美女国产高潮福利片在线看| 欧美不卡视频在线免费观看 | 精品欧美一区二区三区在线| 亚洲免费av在线视频| 亚洲一卡2卡3卡4卡5卡精品中文| 在线观看免费午夜福利视频| 亚洲人成电影免费在线| 无人区码免费观看不卡| 十分钟在线观看高清视频www| 人人妻人人澡人人看| 亚洲五月婷婷丁香| 成人亚洲精品一区在线观看| 欧美黑人欧美精品刺激| 亚洲国产毛片av蜜桃av| 男人舔女人下体高潮全视频| 亚洲第一av免费看| av欧美777| 国产成年人精品一区二区 | 欧美日韩视频精品一区| 成年人黄色毛片网站| 免费看十八禁软件| 麻豆av在线久日| 视频在线观看一区二区三区| 日韩人妻精品一区2区三区| 亚洲国产欧美一区二区综合| 狠狠狠狠99中文字幕| 国产成人精品无人区| 久久久久久久久久久久大奶| 中文欧美无线码| 久久香蕉激情| 亚洲人成电影免费在线| 女人被狂操c到高潮| 精品久久蜜臀av无| 91大片在线观看| 久久人人爽av亚洲精品天堂| 制服人妻中文乱码| 黑人欧美特级aaaaaa片| 最近最新中文字幕大全电影3 | 欧美日韩精品网址| 身体一侧抽搐| 日韩成人在线观看一区二区三区| 国产黄a三级三级三级人| 国产av精品麻豆| 男女之事视频高清在线观看| av在线天堂中文字幕 | 久久精品91无色码中文字幕| 视频在线观看一区二区三区| 后天国语完整版免费观看| 啪啪无遮挡十八禁网站| 在线观看免费午夜福利视频| 熟女少妇亚洲综合色aaa.| 色综合婷婷激情| 9热在线视频观看99| 18禁观看日本| 午夜福利一区二区在线看| 国产精品日韩av在线免费观看 | 91精品三级在线观看| 俄罗斯特黄特色一大片| 超色免费av| 桃色一区二区三区在线观看| 欧美大码av| 午夜激情av网站| 搡老熟女国产l中国老女人| 夜夜看夜夜爽夜夜摸 | 国产高清激情床上av| 波多野结衣高清无吗| www.自偷自拍.com| 身体一侧抽搐| 亚洲专区国产一区二区| 99久久人妻综合| 999久久久精品免费观看国产| 大型黄色视频在线免费观看| 久久精品91蜜桃| 麻豆成人av在线观看| 老熟妇仑乱视频hdxx| 在线观看www视频免费| 一级a爱片免费观看的视频| 欧美+亚洲+日韩+国产| 色哟哟哟哟哟哟| 丝袜在线中文字幕| 欧美激情 高清一区二区三区| 黑人巨大精品欧美一区二区mp4| 中文字幕精品免费在线观看视频| 欧美精品啪啪一区二区三区| 咕卡用的链子| 久久久久久久久中文| 一本综合久久免费| 丁香六月欧美| 日日爽夜夜爽网站| 怎么达到女性高潮| 国产高清videossex| 亚洲在线自拍视频| 亚洲国产中文字幕在线视频| 久久久国产成人精品二区 | 精品福利观看| 亚洲中文字幕日韩| 国产精品1区2区在线观看.| 又紧又爽又黄一区二区| 亚洲精品一卡2卡三卡4卡5卡| 不卡av一区二区三区| 亚洲欧美日韩另类电影网站| 一边摸一边抽搐一进一出视频| 国产亚洲精品一区二区www| 欧美精品啪啪一区二区三区| 人人澡人人妻人| 男人操女人黄网站| 国产成年人精品一区二区 | 国产男靠女视频免费网站| 国产精品野战在线观看 | 在线播放国产精品三级| 色精品久久人妻99蜜桃| 黄色女人牲交| 国产亚洲av高清不卡| 久久国产亚洲av麻豆专区| 一进一出抽搐动态| 黑丝袜美女国产一区| av在线播放免费不卡| 国产乱人伦免费视频| 免费观看精品视频网站| 久久 成人 亚洲| 国产亚洲精品一区二区www| 女性被躁到高潮视频| 欧美日韩av久久| 999久久久国产精品视频| 神马国产精品三级电影在线观看 | 757午夜福利合集在线观看| 精品福利观看| 激情在线观看视频在线高清| av免费在线观看网站| 91字幕亚洲| 国产一区二区三区视频了| www.自偷自拍.com| 大陆偷拍与自拍| 丝袜人妻中文字幕| 欧美人与性动交α欧美精品济南到| 少妇的丰满在线观看| 男男h啪啪无遮挡| 亚洲精品中文字幕在线视频| 国产精品香港三级国产av潘金莲| av在线天堂中文字幕 | 高清在线国产一区| 国产精品99久久99久久久不卡| 亚洲成人精品中文字幕电影 | 两性夫妻黄色片| 成人18禁高潮啪啪吃奶动态图| 变态另类成人亚洲欧美熟女 | 欧美日韩精品网址| 亚洲 国产 在线| 国产免费av片在线观看野外av| 日日干狠狠操夜夜爽| av片东京热男人的天堂| xxx96com| 亚洲精品成人av观看孕妇| av天堂在线播放| 88av欧美| 久久亚洲真实| 狂野欧美激情性xxxx| 国产伦一二天堂av在线观看| 成人国产一区最新在线观看| 国产高清激情床上av| 亚洲狠狠婷婷综合久久图片| www.www免费av| 中文字幕人妻熟女乱码| 成人三级做爰电影| 午夜91福利影院| 不卡一级毛片| 国产精品国产av在线观看| 最近最新中文字幕大全电影3 | 国产一区二区三区视频了| 国产xxxxx性猛交| 日韩视频一区二区在线观看| 极品人妻少妇av视频| 91老司机精品| av免费在线观看网站| 久久国产乱子伦精品免费另类| 女生性感内裤真人,穿戴方法视频| 一边摸一边抽搐一进一小说| 亚洲男人的天堂狠狠| 亚洲午夜精品一区,二区,三区| 久热这里只有精品99| 在线观看午夜福利视频| 成人18禁高潮啪啪吃奶动态图| 精品一区二区三区视频在线观看免费 | 亚洲色图av天堂| 亚洲激情在线av| 免费女性裸体啪啪无遮挡网站| 别揉我奶头~嗯~啊~动态视频| 成熟少妇高潮喷水视频| 午夜91福利影院| 亚洲伊人色综图| 亚洲狠狠婷婷综合久久图片| 免费在线观看亚洲国产| 亚洲av片天天在线观看| 自线自在国产av| 欧美激情高清一区二区三区| 人妻久久中文字幕网| 亚洲专区国产一区二区| 宅男免费午夜| 欧美精品一区二区免费开放| 欧美最黄视频在线播放免费 | 亚洲成人精品中文字幕电影 | 性少妇av在线| 热99国产精品久久久久久7| av在线播放免费不卡| 18美女黄网站色大片免费观看| x7x7x7水蜜桃| 成人18禁高潮啪啪吃奶动态图| 国产精品久久久久成人av| 免费在线观看影片大全网站| 日韩 欧美 亚洲 中文字幕| 午夜福利在线免费观看网站| 黄片大片在线免费观看| 一级片'在线观看视频| 国产熟女xx| 欧美日韩亚洲高清精品| 1024香蕉在线观看| 中文字幕另类日韩欧美亚洲嫩草| 久久久精品国产亚洲av高清涩受| 国产精品一区二区在线不卡| 99国产精品一区二区三区| 亚洲专区中文字幕在线| 深夜精品福利| 亚洲精华国产精华精| 国产亚洲欧美精品永久| 美女福利国产在线| 国产高清视频在线播放一区| 80岁老熟妇乱子伦牲交| 最近最新中文字幕大全免费视频| 国产极品粉嫩免费观看在线| 亚洲在线自拍视频| av免费在线观看网站| 亚洲国产看品久久| 俄罗斯特黄特色一大片| 黄色片一级片一级黄色片| 级片在线观看| 国产精品1区2区在线观看.| 精品久久久久久,| 午夜福利在线免费观看网站| 国产欧美日韩一区二区精品| 亚洲情色 制服丝袜| 亚洲av五月六月丁香网| 搡老岳熟女国产| 大陆偷拍与自拍| 亚洲欧美激情综合另类| 国产99久久九九免费精品| 999久久久精品免费观看国产| 国产免费av片在线观看野外av| 国产精品成人在线| 色婷婷久久久亚洲欧美| www.999成人在线观看| 午夜福利在线观看吧| 欧美黄色淫秽网站| 久久久久久久久中文| 一进一出抽搐gif免费好疼 | 国产精品九九99| 女生性感内裤真人,穿戴方法视频| 亚洲精品美女久久av网站| 亚洲精品粉嫩美女一区| 色综合站精品国产| 99精国产麻豆久久婷婷| 国产亚洲精品久久久久5区| 日韩av在线大香蕉| 电影成人av| 一区在线观看完整版| 中文字幕人妻熟女乱码| 性欧美人与动物交配| 国产成人啪精品午夜网站| 黄色怎么调成土黄色| 丝袜美足系列| 国产99白浆流出| 久久久久亚洲av毛片大全| av欧美777| 亚洲国产欧美一区二区综合| 欧洲精品卡2卡3卡4卡5卡区| 97人妻天天添夜夜摸| 一区福利在线观看| 大型av网站在线播放| 精品日产1卡2卡| 黄色成人免费大全| 亚洲片人在线观看| 亚洲人成电影观看| 精品久久久久久久久久免费视频 | 亚洲第一av免费看| 中文字幕av电影在线播放| 国产成人欧美| 久久人妻熟女aⅴ| 欧美一级毛片孕妇| 一级黄色大片毛片| 亚洲第一青青草原| 后天国语完整版免费观看| 亚洲色图综合在线观看| 久久久久久免费高清国产稀缺| 午夜福利,免费看| 操出白浆在线播放| 女性被躁到高潮视频| 长腿黑丝高跟| 免费在线观看亚洲国产| 色哟哟哟哟哟哟| 亚洲一码二码三码区别大吗| 国产99久久九九免费精品| 村上凉子中文字幕在线| 国产色视频综合| 757午夜福利合集在线观看| 国产精品九九99| www.熟女人妻精品国产| www.自偷自拍.com| 国产成人影院久久av| 老司机深夜福利视频在线观看| 国产精品一区二区三区四区久久 | 日韩大尺度精品在线看网址 | 日韩免费av在线播放| 黄频高清免费视频| 午夜免费鲁丝| 天天影视国产精品| 男人操女人黄网站| 水蜜桃什么品种好| 一个人免费在线观看的高清视频| 亚洲欧洲精品一区二区精品久久久| 999久久久精品免费观看国产| 国产亚洲精品久久久久5区| 国产在线精品亚洲第一网站| 国产亚洲精品第一综合不卡| 国产黄色免费在线视频| 国产亚洲精品综合一区在线观看 | 色在线成人网| 18禁裸乳无遮挡免费网站照片 | 成人av一区二区三区在线看| 男女之事视频高清在线观看| 精品一区二区三区av网在线观看| 日韩成人在线观看一区二区三区| 久9热在线精品视频| 老熟妇乱子伦视频在线观看| 久久久国产精品麻豆| 久久人人爽av亚洲精品天堂| 日韩欧美一区视频在线观看| 精品欧美一区二区三区在线| 麻豆久久精品国产亚洲av | 午夜91福利影院| 日本精品一区二区三区蜜桃| 久久伊人香网站| 村上凉子中文字幕在线| 在线观看免费午夜福利视频| 成人特级黄色片久久久久久久| 国产伦人伦偷精品视频| 91在线观看av| 欧美午夜高清在线| 国产精品秋霞免费鲁丝片| 国产精品亚洲一级av第二区| 满18在线观看网站| 女人爽到高潮嗷嗷叫在线视频| 啦啦啦 在线观看视频| 国产高清videossex| 美女午夜性视频免费| 很黄的视频免费| 日韩中文字幕欧美一区二区| а√天堂www在线а√下载| 在线国产一区二区在线| 成年人免费黄色播放视频| 淫妇啪啪啪对白视频| 精品一品国产午夜福利视频| 最新在线观看一区二区三区| 黄色 视频免费看| 夜夜爽天天搞| 老司机亚洲免费影院| 久久精品人人爽人人爽视色| 久久午夜综合久久蜜桃| 亚洲美女黄片视频| 一级作爱视频免费观看| 变态另类成人亚洲欧美熟女 | 国产一区二区在线av高清观看| 久久人妻av系列| 久久精品aⅴ一区二区三区四区| 99久久国产精品久久久| 亚洲av五月六月丁香网| 精品乱码久久久久久99久播| 一级片'在线观看视频| 美女午夜性视频免费| 精品无人区乱码1区二区| 色老头精品视频在线观看| 亚洲精品美女久久av网站| 午夜免费鲁丝| 欧美激情 高清一区二区三区| 欧美另类亚洲清纯唯美| 韩国精品一区二区三区| 麻豆国产av国片精品| 午夜两性在线视频| 色精品久久人妻99蜜桃| 手机成人av网站| 一夜夜www| 久久精品亚洲熟妇少妇任你| 午夜亚洲福利在线播放| 中文字幕精品免费在线观看视频| av在线天堂中文字幕 | 欧美成人午夜精品| 夜夜爽天天搞| 亚洲午夜精品一区,二区,三区| 中文字幕人妻丝袜制服| 欧美乱码精品一区二区三区| 美女高潮到喷水免费观看| 国产精品电影一区二区三区| 国产欧美日韩精品亚洲av| 人人妻人人添人人爽欧美一区卜| 亚洲第一青青草原| 国产在线精品亚洲第一网站| 黄片小视频在线播放| 黄色丝袜av网址大全| 国产免费av片在线观看野外av| 一级片'在线观看视频| 国产一区在线观看成人免费| 成年女人毛片免费观看观看9| 亚洲自拍偷在线| 首页视频小说图片口味搜索| 久久婷婷成人综合色麻豆| 高清毛片免费观看视频网站 | 人人妻,人人澡人人爽秒播| 精品无人区乱码1区二区| 精品人妻1区二区| 最新在线观看一区二区三区| 最好的美女福利视频网| 少妇的丰满在线观看| 日本欧美视频一区| 国产亚洲欧美精品永久| av在线播放免费不卡| 国产精品 欧美亚洲| 少妇被粗大的猛进出69影院| 欧美成人免费av一区二区三区| 91成人精品电影| 国产亚洲欧美98| 国产男靠女视频免费网站| 美女 人体艺术 gogo| 自线自在国产av| 人人妻人人澡人人看| 纯流量卡能插随身wifi吗| 69精品国产乱码久久久| 12—13女人毛片做爰片一| 一区二区三区精品91| 久久香蕉精品热| 真人一进一出gif抽搐免费| 精品国产超薄肉色丝袜足j| 嫁个100分男人电影在线观看| 国产成人欧美| 久久香蕉激情| 国产高清视频在线播放一区| 欧美人与性动交α欧美精品济南到| 99久久久亚洲精品蜜臀av| 久久久国产成人免费|