郭肖，劉瑞璇*，高振東，王玥，王鵬鯤

(1.西南石油大學(xué)，油氣藏地質(zhì)及開發(fā)工程國家重點(diǎn)實(shí)驗(yàn)室，成都 610500；2.延長油田股份有限公司，延安 716000)

中國致密油藏資源豐富，資源儲量約為44×108t，居世界第三位。目前，中國致密油的開發(fā)正處于初級階段[1-2]，一些油田針對致密油低孔低滲的特點(diǎn)對水平井進(jìn)行壓裂開采，產(chǎn)能均得到了提升，但是效果并不理想。因此需要加快對致密油藏的研究與探索，保障國家能源安全。致密油藏物性差，儲層絕對滲透率一般小于2×10-3μm2，在這樣的儲層條件下不進(jìn)行壓裂幾乎無工業(yè)性油流。國內(nèi)外開發(fā)實(shí)踐證明，水平井和大型水力壓裂技術(shù)是成功開發(fā)致密油藏的關(guān)鍵[3-8]。壓裂后的水平井增加了裂縫滲透率及泄油面積，從而提高了油藏產(chǎn)能[9]。

致密油水平井彈性開采存在初期能量遞減快、采收率低、能量不足等問題，針對這些問題[10-13]，有大量學(xué)者進(jìn)行了相關(guān)研究。朱杰等[14]針對新疆油田環(huán)瑪湖區(qū)塊，開展室內(nèi)研究實(shí)驗(yàn)，優(yōu)化吞吐參數(shù)，對該區(qū)塊致密油增產(chǎn)提供了指導(dǎo)。候廣[15]采用實(shí)驗(yàn)與數(shù)值模擬相結(jié)合的方法，優(yōu)化吞吐參數(shù)，解決了致密油產(chǎn)能低的問題。史曉東等[16]針對松遼盆地扶余油層高產(chǎn)期短、產(chǎn)量遞減快等特點(diǎn)，將室內(nèi)模擬實(shí)驗(yàn)與數(shù)值模擬相結(jié)合，確定了最佳CO2吞吐參數(shù)。周拓等[17]通過物理模擬實(shí)驗(yàn)對分段壓裂水平井CO2吞吐機(jī)理進(jìn)行了研究，得到了對CO2吞吐的影響因素。曾佳[18]基于非線性滲流力學(xué)理論,針對長7致密油藏建立了致密油藏復(fù)合三線性滲流數(shù)學(xué)模型，得到了最優(yōu)的開發(fā)方案參數(shù)。郎慧慧[19]針對馬嶺油田M區(qū)進(jìn)行地質(zhì)儲層分析，研究得到體積壓裂水平井CO2吞吐方案。上述研究對提高致密油采收率取得了一定的進(jìn)展，但針對致密油吞吐方式、井網(wǎng)及生產(chǎn)參數(shù)整體優(yōu)化和選擇最優(yōu)的吞吐方案一直沒有很好的解答。因此，基于致密油開發(fā)現(xiàn)狀，通過建立鄂爾多斯長7儲層數(shù)值模擬模型，對吞吐方式、井網(wǎng)參數(shù)、生產(chǎn)參數(shù)等進(jìn)行整體優(yōu)化，并給出合適的吞吐方案，為該儲層提高采收率提供方法支持和借鑒。

1 區(qū)域地質(zhì)概況

鄂爾多斯盆地主要劃分為伊盟隆起、渭北隆起、陜北斜坡、晉西撓褶帶、天環(huán)坳陷和西緣沖斷構(gòu)造帶六個一級構(gòu)造單元[20]，總面積大約8.7×104km2，如圖1(a)所示。延長組長7段是致密油的主要區(qū)塊。長7儲層自上而下共劃分為長71、72、73三個亞段[21]，長71、72段主要以細(xì)砂巖、泥質(zhì)粉砂巖為主，73段主要以頁巖為主，如圖1(b)所示。

圖1 鄂爾多斯盆地長7儲層構(gòu)造及地層示意圖Fig.1 Structure and stratigraphic diagram of Chang 7 reservoir in Ordos Basin

2 數(shù)值模型建立

根據(jù)鄂爾多斯盆地長7儲層的原油常規(guī)物性分析，其原油密度為0.802～0.863 g/cm3，平均為0.832 g/cm3；動力黏度(50 ℃)為2.93～12.59 mPa·s，平均為7.03 mPa·s；凝固點(diǎn)為2～31 ℃，平均為17.5 ℃。長7段致密油儲層壓力系數(shù)較低，分布范圍只有0.64～0.87。鑒于此，現(xiàn)以鄂爾多斯盆地長7儲層為研究對象建立儲層基礎(chǔ)模型，模型構(gòu)建流程如圖2所示。

圖2 致密油水平井補(bǔ)充能量模型流程圖Fig.2 Flowchart of supplementary energy model for tight oil horizontal wells

運(yùn)用CMG(computer modelling group)組分模塊建立了水平井分段壓裂數(shù)值模擬模型。網(wǎng)格建立采用笛卡爾坐標(biāo)系，網(wǎng)格步長均為20 m×20 m×3 m，網(wǎng)格數(shù)總計72×50×8，水平井分段壓裂，基本參數(shù)如表1所示。

表1 模型基本參數(shù)Table 1 Model basic parameters

3 不同注入方式

致密油在彈性開采過程中能量衰竭較快，在保持地層能量的條件下以較高的采油速度開采則需要及時的補(bǔ)充能量。利用CMG的組分模型對不同注入方式CO2、N2、20 ℃水和80 ℃水吞吐進(jìn)行模擬研究，將采收率及累計產(chǎn)油作為分析評定依據(jù)，數(shù)據(jù)模擬結(jié)果如圖3、圖4所示，累計產(chǎn)油量變化趨勢如圖中的折線所示。結(jié)果表明CO2吞吐的開發(fā)效果最好，20年生產(chǎn)末期采收4.692 78%，N2吞吐次之，達(dá)到了3.997 13%，20 ℃水吞吐的效果最差，最終采收率僅為3.316 2%。

圖3 不同注入方式下的采收率Fig.3 Recovery efficiency under different injection modes

圖4 不同注入方式下的累計產(chǎn)油Fig.4 Cumulative oil production under different injection modes

對比CO2、N2、20 ℃水、80 ℃水吞吐采收率，可以發(fā)現(xiàn)CO2在生產(chǎn)前三年的效果并不突出，與其他3種注入方式下的采收效果基本一樣，但是隨著生產(chǎn)時間的增加，CO2吞吐的優(yōu)勢逐漸體現(xiàn)出來。生產(chǎn)模擬時間為20年，CO2吞吐的累計產(chǎn)油量達(dá)到了93 033.7 m3，對比20 ℃水的產(chǎn)油量提高了34.29%。結(jié)果表明在致密油水平井分段壓裂模式下，氣體相較于液體具有更好的流動性，更能充分發(fā)揮裂縫優(yōu)勢，與地層中的原油互溶反應(yīng)，從而起到補(bǔ)充能量的作用以達(dá)到較高的產(chǎn)能。

不同注入方式下，壓力變化如圖5所示，含油飽和度變化如圖6所示，CO2黏度變化如圖7所示。對比變化趨勢可以看出，CO2吞吐過后裂縫和靠近水平井的區(qū)域能量補(bǔ)充及原油排驅(qū)效果較好。說明水平井周圍壓力隨著CO2不斷注入地層而逐漸增加。CO2在較高的地層壓力下驅(qū)替近井端的的可流動原油，同時發(fā)生指進(jìn)而進(jìn)入遠(yuǎn)處地層。在注入階段，地層能量得到補(bǔ)充，生產(chǎn)井控制范圍內(nèi)的壓力也不斷升高。注入之后開始燜井，一般燜井時間與注入時間、流體性質(zhì)相關(guān)。在這一階段中，注入的CO2與原油充分反應(yīng)發(fā)生萃取作用，主要提取原油中的輕質(zhì)組分，達(dá)到降黏、原油膨脹、降低界面張力的作用。燜井結(jié)束之后，進(jìn)行吞吐的最后一步采油，原油會隨注入的CO2流向井筒。整個吞吐的過程達(dá)到補(bǔ)充地層能量、降低含油飽和度的目的。因此，優(yōu)化選取CO2吞吐的注入方式。

圖5 不同注入方式在開采兩年之后的壓力變化Fig.5 Pressure changes of different injection methods after two years of mining

圖6 不同注入方式下的含油飽和度Fig.6 Oil saturation under different injection methods

圖7 不同開采年份CO2黏度變化Fig.7 CO2 viscosity change in different mining years

4 模擬結(jié)果

4.1 井網(wǎng)形式

采用CO2吞吐的開發(fā)方式，利用數(shù)值模擬對不同的井網(wǎng)形式進(jìn)行優(yōu)化分析，生產(chǎn)時間為20年。分別模擬了五點(diǎn)交錯井網(wǎng)、七點(diǎn)正對井網(wǎng)、七點(diǎn)交錯井網(wǎng)、九點(diǎn)正對井網(wǎng)、九點(diǎn)交錯井網(wǎng)，如圖8所示。比較在相同參數(shù)下，不同井網(wǎng)形式的采收率、累積產(chǎn)油量及換油率。開發(fā)效果如圖9所示，結(jié)合累計產(chǎn)油、換油率作為評價依據(jù)。換油率為產(chǎn)出油質(zhì)量與注入氣質(zhì)量之比。累計產(chǎn)油量隨著水平井的增加而增加，總體分析得到九點(diǎn)交錯井網(wǎng)獲得最優(yōu)的生產(chǎn)效果。

圖8 不同井網(wǎng)類型Fig.8 Different well pattern types

圖9 不同井網(wǎng)的累產(chǎn)油及換油率Fig.9 Accumulated oil production and oil exchange rate of different well patterns

4.2 井網(wǎng)參數(shù)優(yōu)化

4.2.1 水平井長度

在其他參數(shù)一定的條件下，選擇水平段長度分別為210、270、330、390、450、510 m進(jìn)行數(shù)值模擬預(yù)測，以采收率和累計產(chǎn)油量作為評定分析標(biāo)準(zhǔn)。數(shù)值計算結(jié)果如圖10所示，采收率隨水平井長度的變化趨勢如圖中的折線所示。在生產(chǎn)初期，水平井長度與采收率、累計產(chǎn)油量成正比，但是當(dāng)水平井長度達(dá)到一定值后，產(chǎn)能增加趨勢逐漸變緩。這是因?yàn)榫矁?nèi)流體所受到的摩擦阻力會隨著水平井長的增加而增加。因此，水平井長度存在最優(yōu)值，從圖10中可以看出：當(dāng)水平段長度達(dá)到330 m后，產(chǎn)量增加幅度明顯減??；優(yōu)化結(jié)果取水平井長度為330 m為宜。

圖10 水平井長度與采收率關(guān)系Fig.10 Relationship between horizontal well length and recovery factor

4.2.2 裂縫半長

模擬不同裂縫半長下水平井分段壓裂的產(chǎn)能變化，對裂縫半長進(jìn)行優(yōu)化。裂縫半長分別取20、40、60、80 m，模擬開采時間為20年，對不同裂縫半長情況下的水平井生產(chǎn)情況進(jìn)行模擬。不同裂縫半長隨累計產(chǎn)油量的變化趨勢如圖中的折線所示。從圖11、圖12看出，在其他參數(shù)一定的情況下，不同的裂縫半長對分段壓裂水平井的采收率及累積產(chǎn)量影響較大。但壓裂后的累積產(chǎn)量并不是隨著裂縫半長增加而線性増長，產(chǎn)能增加幅度隨著裂縫半長的增加逐漸減小。這是因?yàn)殡S著水力壓裂裂縫半長的增大，由于邊界作用使得裂縫波及范圍減小，導(dǎo)致水平井累積產(chǎn)量增幅越來越小，考慮投入產(chǎn)出比，則存在一個相對最佳的裂縫半長值，優(yōu)化結(jié)果表明取裂縫半長為60 m為宜。

圖11 裂縫半長與采收率關(guān)系Fig.11 Relationship between fracture half-length and recovery factor

圖12 裂縫半長與累計產(chǎn)油量關(guān)系Fig.12 Relationship between fracture half-length and cumulative oil production

4.2.3 井排距

排距指的是兩排水平井之間的距離如圖13所示。以井距為180 m進(jìn)行建模，分別取排距的值為160、180、200、220、240 m進(jìn)行數(shù)值模擬預(yù)測，以采收率和累計產(chǎn)油量作為評定分析標(biāo)準(zhǔn)，數(shù)據(jù)結(jié)果如圖14所示。當(dāng)排距小于220 m時，產(chǎn)能隨排距增大而增大，但當(dāng)排距大于220 m時，產(chǎn)能與排距之間呈現(xiàn)負(fù)相關(guān)。排距過小時，井排之間容易產(chǎn)生干擾，而排距過大，會存在大片未波及區(qū)域從而降低采出程度。優(yōu)化結(jié)果表明取井排距為220 m為宜。

圖13 井距及排距示意圖Fig.13 Well spacing and row spacing diagram

圖14 井排距與累計產(chǎn)油量、采收率關(guān)系Fig.14 Relationship between well spacing and cumulative oil production and recovery factor

4.3 吞吐參數(shù)優(yōu)化

4.3.1 注氣速度

保持其他參數(shù)不變，每口水平井選取注氣速度分別為2 000、2 500、3000、3 500、4 000、4 500、5 000 m3/d，開展CO2吞吐的開采模擬研究，以采收率和累計產(chǎn)油量作為評定分析標(biāo)準(zhǔn)，模擬結(jié)果如圖15所示。

圖15 注氣速度與累計產(chǎn)油量、采收率關(guān)系Fig.15 Relationship between gas injection rate and cumulative oil production and recovery factor

當(dāng)日注氣量從每口井2 500 m3/d增大到5 000 m3/d時，采收率與累計產(chǎn)油量成線性增長，注入的CO2越多，則地層能量補(bǔ)充越充足，使得地層壓力在注入階段之后恢復(fù)得更好，以獲得更高的產(chǎn)能。從提高采收率的角度來看，CO2注入越多則開采越好。

但是在實(shí)際過程種，注氣量受到注入壓力的限制不能過大，并且生產(chǎn)過程中隨原油產(chǎn)出的CO2會隨著注入CO2量的增加而增加，因此注入地層的CO2并不能全部發(fā)揮作用。實(shí)際工程中注氣速度根據(jù)現(xiàn)場情況而定，優(yōu)化結(jié)果取最佳注氣速度為5 000 m3/d。

4.3.2 燜井時間

1)短周期燜井

燜井時間是指從注氣階段到開井生產(chǎn)這一時間段過程。CO2注入地層后主要集中水平井附近區(qū)域的裂縫中，在燜井階段中，CO2在分子擴(kuò)散的作用下擴(kuò)散至裂縫更深處，與更多原油充分反應(yīng)，補(bǔ)充地層能量。但是較長的燜井時間會導(dǎo)致CO2從原油中分離出來，降低作用效果?？梢姞F井時間對產(chǎn)能有一定的影響。為了研究燜井時間對分段壓裂水平井CO2吞吐開發(fā)效果的影響，在注入速度均為每口井5 000 m3、注入時間10 d的情況下設(shè)置0、3、5、7、9、11、13 d的關(guān)井時間，對不同燜井時間進(jìn)行數(shù)值模擬研究。以采收率和累計產(chǎn)油量作為評定分析標(biāo)準(zhǔn)，模擬結(jié)果如圖16所示，累計產(chǎn)油量隨不同燜井時間的變化趨勢如圖中的折線所示。

依據(jù)數(shù)據(jù)分析圖對比可以看出燜井時間對生產(chǎn)具有一定的影響。在燜井過程中，注入地層近井端的CO2隨著燜井時間的增加，由裂縫向基質(zhì)更遠(yuǎn)處擴(kuò)散，與原油充分接觸反應(yīng)。但較長的燜井時間會降低產(chǎn)能，如圖16所示。當(dāng)燜井時間超過3 d時，累計產(chǎn)油量及采收率與燜井時間呈現(xiàn)負(fù)相關(guān)。

圖16 不同燜井時間與累計產(chǎn)油量及采收率的關(guān)系Fig.16 The relationship between different soaking time and cumulative oil production and recovery factor

2)長周期燜井

為了對比長、短燜井周期對產(chǎn)能的影響，增加燜井時間，同時將注入時間由短周期的10 d增加至30 d，其他參數(shù)保持不變，設(shè)置5、10、15、20、25、30 d關(guān)井時間，對不同燜井時間進(jìn)行數(shù)值模擬研究。如圖17、圖18所示，長周期吞吐有更高的累計產(chǎn)油、采收率，提高了0.567%，且日產(chǎn)油整體比短周期吞吐穩(wěn)定。

圖17 燜井時間與累計產(chǎn)油量、采收率的關(guān)系Fig.17 Relationship between soak time and cumulative oil production and recovery factor

圖18 日產(chǎn)油與燜井周期的關(guān)系Fig.18 Relationship between daily oil production and soak period

生產(chǎn)開采過程中，CO2首先進(jìn)入裂縫，然后沿裂縫進(jìn)入裂縫周圍區(qū)域，待裂縫周圍區(qū)域壓力升高以后，逐步向壓力未波及區(qū)域擴(kuò)展。如圖19為開采兩年之后長周期燜井時間與地層壓力變化，可以看出在較短時間的關(guān)井時間下，地層壓力還未得到及時的恢復(fù)就要進(jìn)行下一輪的吞吐，由此可見頻繁多次的吞吐方式不能有助于地層能量補(bǔ)充。在燜井時間較短的情況下，地層壓力下降較快且能量不能得到補(bǔ)充。因此，優(yōu)化結(jié)果選擇長周期吞吐方式，燜井時間為10 d。

圖19 不同燜井時間下的壓力變化Fig.19 Pressure change under different soaking time

4.3.3 生產(chǎn)時間

生產(chǎn)時間對于CO2吞吐效果也有一定的影響。設(shè)定單周期生產(chǎn)時間分別為100、150、200、250 d，保持注入時間、燜井時間等其他參數(shù)不變，開展CO2吞吐的開采模擬研究，模擬結(jié)果如圖20所示。隨著生產(chǎn)時間的增加，累計產(chǎn)油及采收率在不斷降低。結(jié)果顯示生產(chǎn)時間為280 d時采收率最高。

圖20 不同生產(chǎn)時間與累計產(chǎn)油、采收率的關(guān)系Fig.20 Relationship between different production time and cumulative oil production and recovery factor

圖21為生產(chǎn)前4年的壓力變化，說明生產(chǎn)時間對產(chǎn)能及地層壓力變化具有一定的影響。在注入及燜井時間不變的情況下，生產(chǎn)時間越長則產(chǎn)能越低，是因?yàn)榈貙訅毫﹄S著生產(chǎn)時間的增加不斷降低，并且長時間的開采是在過快的消耗地層能量，從而影響后期的開采。

圖21 不同生產(chǎn)時間與壓力變化Fig.21 Different production time and pressure change

4.3.4 吞吐時機(jī)

油藏在轉(zhuǎn)吞吐之前衰竭開采的時間長短決定了轉(zhuǎn)吞吐時地層壓力的大小，分別設(shè)計衰竭開采0、1、2、3年后轉(zhuǎn)吞吐開采，數(shù)值模擬結(jié)果如圖22所示，模擬結(jié)果表明在衰竭開采一年之后進(jìn)行吞吐開采的開發(fā)效果較好，產(chǎn)能較高。

圖22 不同吞吐時機(jī)與累產(chǎn)油的關(guān)系Fig.22 Relationship between different huff and puff time and cumulative oil production

根據(jù)模擬結(jié)果可以看出，吞吐時機(jī)對于累計產(chǎn)油量具有一定的影響。原油初始組分中輕質(zhì)組分占比較大，較早的進(jìn)行吞吐，可以發(fā)揮CO2的萃取作用，相同的，轉(zhuǎn)注時機(jī)較晚，原油中的輕質(zhì)組分含量下降，就會削弱CO2增產(chǎn)降黏的作用。

5 結(jié)論

(1)鄂爾多斯長7儲層屬于致密區(qū)塊，優(yōu)選開發(fā)方式為CO2吞吐，20年生產(chǎn)末期采收4.692 78%。

(2)通過數(shù)值模擬研究證明，九點(diǎn)交錯井網(wǎng)具有更好的開發(fā)效果，優(yōu)化確定了水平井長度330 m，裂縫半長60 m，井排距220 m。

(3)通過吞吐參數(shù)敏感性分析優(yōu)化，最優(yōu)的生產(chǎn)制度為：注氣速度為5 000 m3/d，燜井時間10 d，生產(chǎn)時間280 d，吞吐時機(jī)為衰竭開采1年之后。

(4)通過數(shù)值模擬優(yōu)化分析，水平井分段壓裂儲層改造及合理的CO2吞吐參數(shù)設(shè)計有助于致密油提高產(chǎn)能。