華文靜, 李治平
(1.中國地質(zhì)大學(xué)(北京)能源學(xué)院,北京 100083; 2.中國地質(zhì)大學(xué)(北京)非常規(guī)天然氣地質(zhì)評價與開發(fā)北京市重點實驗室,北京 100083)
中國致密儲層分布廣泛,儲量豐富,其勘探開發(fā)對中國能源行業(yè)有著重要意義[1-2],但也廣泛存在著注水開發(fā)注不進、采不出等問題。延長油田作為中國最大的低滲透油田,歷經(jīng)多年勘探開發(fā),天然能量衰竭,大部分地區(qū)進入注水開發(fā)中后期階段,油水滲流能力下降,產(chǎn)能低,需要采取有效措施補充地層能量。CO2驅(qū)油最先在美國、加拿大等地取得成功[3],中國也在吉林油田、大慶油田等地取得多次礦場先導(dǎo)性試驗的成功[4]。廣泛的研究表明,CO2驅(qū)油技術(shù)對于低滲透和稠油油藏有良好的增產(chǎn)作用,為水驅(qū)效果不明顯或水驅(qū)施行難的致密油藏提供了解決方案。
隨著對CO2驅(qū)油技術(shù)的研究,其局限性也逐漸呈現(xiàn)。不同的驅(qū)替方式,如氣水交替驅(qū)、連續(xù)氣驅(qū)與CO2吞吐,其驅(qū)油結(jié)果不盡相同,即使在同一驅(qū)替方式下,工作制度的不同,也極大地影響著CO2驅(qū)替的效果[5-8]。注入量、注入時機、燜井時間以及吞吐輪次是CO2吞吐開采的重要參數(shù),其取值將直接影響油藏開發(fā)以及最終的CO2吞吐效果。因此,確定合理的CO2吞吐工作參數(shù),對于提高油田采收率有著重要的意義。
文獻[9]利用Elm Coulee Field油田的數(shù)據(jù)進行了不同注氣方式的對比研究,發(fā)現(xiàn)混相驅(qū)替能極大程度提高致密油的采收率;文獻[10]采用數(shù)值模擬的方法對巴肯油田進行注氣開發(fā)研究,認為井網(wǎng)設(shè)計是提高采收率的關(guān)鍵;文獻[11]同樣采用數(shù)值模擬方法,對致密油的不同的注氣方式進行研究,并研究了水平井中CO2吞吐參數(shù)對開采效果的影響。文獻[12]通過某區(qū)塊某口井CO2吞吐實驗與數(shù)值模擬的方法,得到了5個敏感因素對吞吐效果的影響與現(xiàn)場實施應(yīng)選取的參數(shù)范圍;文獻[13]結(jié)合模糊綜合評判方法與神經(jīng)網(wǎng)絡(luò)技術(shù),進行CO2吞吐候選井的優(yōu)選與綜合評價,并從CO2吞吐實例中優(yōu)選油藏參數(shù)及最優(yōu)參數(shù)組合。
這些研究已經(jīng)證實CO2驅(qū)油工作制度以及地層因素是影響CO2驅(qū)油效果的重要因素,將直接影響開發(fā)效益以及最終的采收率。目前主要通過室內(nèi)實驗和數(shù)值模擬兩種方式進行參數(shù)優(yōu)化。室內(nèi)實驗多是利用巖心驅(qū)替實驗,既可以使用人造巖心,也可利用實際巖心,對不同條件下CO2驅(qū)油效果進行觀察,對比驅(qū)油結(jié)果,從而進行驅(qū)替方式以及工作參數(shù)等的優(yōu)選[14-16]。數(shù)值模擬方法是根據(jù)油田實際數(shù)據(jù)建立數(shù)值模型,控制變量對開發(fā)指標進行計算對比,從而得到各因素對CO2驅(qū)油效果的影響[17-20],這兩種方法適用于不同類型的油藏。許多學(xué)者結(jié)合二者進行相互驗證,并應(yīng)用到礦場實際,提出現(xiàn)場實施應(yīng)選取的參數(shù)范圍。近年也有學(xué)者利用試井分析等方法進行研究[21-23]。
以延長致密油藏吳定地區(qū)SHH區(qū)塊長8儲層為研究對象,該區(qū)儲層厚度分布穩(wěn)定,構(gòu)造為由東向西傾斜的平緩單斜,地層走向北北東(NNE),北部區(qū)域發(fā)育小型低幅度鼻狀隆起,與區(qū)域整體地勢相吻合。該區(qū)長8儲層油源巖為三疊系湖侵形成的厚層深湖暗色泥巖,儲層受水下分流河道控制,屬于扇三角洲前緣沉積,發(fā)育水下分流河道、席狀砂、分流間灣等沉積微相,扇三角洲前緣亞相是主要儲集層。研究區(qū)長8面孔率一般為3%~7%,樣品中最大面孔率為7%,平均面孔率為3.7%,孔隙直徑大部分分布在0.03~0.15 mm,平均喉道半徑分布在0.1~0.54 μm,孔隙大,分選不好,分選系數(shù)分布在1.41~5.22。長8裂縫產(chǎn)狀為垂直裂縫與高角度裂縫,裂縫寬度小于1 mm,裂縫縫長大多分布在10~30 cm,高度在10 m左右,多有充填。
根據(jù)儲層資料,建立該區(qū)SHH區(qū)塊地質(zhì)模型,通過組分模型模擬CO2吞吐,采用正交試驗設(shè)計方法設(shè)計模擬方案,以累積產(chǎn)油量為評價指標,對工作參數(shù)進行優(yōu)化設(shè)計。
對4個主要注采因素進行優(yōu)化,分別如下。
注入速度越快,指進現(xiàn)象越嚴重。較快的注入速度有利于深入油藏內(nèi)部,與更多原油接觸,有利于進行反應(yīng)。但當注入速度過快時,CO2與原油接觸時間過短,原油膨脹、萃取等作用機理降低,同時對注入壓力的要求增大,注入壓力也受到地層破裂壓力的限制。設(shè)計注入速度為10 000、20 000、30 000 m3/d 3個水平。
隨著CO2注入時間增長,注入量增多,驅(qū)油效果先增強,注入時間過長時,CO2將油推離井底較遠,增產(chǎn)量將降低,并且CO2注入量受經(jīng)濟效益約束。
注氣時長共設(shè)計了15 d、1個月、2個月、3個月、4個月5水平,初始優(yōu)化所選水平為2個月、3個月、4個月。
燜井時間過短,CO2不能與原油充分接觸,不能充分溶于原油之中,原油膨脹與酸化解堵效果不明顯,輕質(zhì)組分萃取效果不佳,會降低CO2吞吐效果,采收率提升不大;燜井時間過長,CO2會擴散到更遠的地層中去,能量散失,產(chǎn)能降低,從而導(dǎo)致產(chǎn)量降低甚至無法產(chǎn)油。
涉及燜井時間的水平有5 d、7 d、15 d、1個月、2個月、3個月、4個月,初始方案設(shè)計所選3水平為2個月、3個月、4個月。
隨著周期數(shù)的增加,累積產(chǎn)油量會相應(yīng)增加,但是每一周期累積產(chǎn)油量的增速會放緩,因此并不是周期數(shù)越多經(jīng)濟效益就會越高,應(yīng)計算各輪次的累積產(chǎn)油量,找出采收率較高同時每一輪次的累積產(chǎn)油量較大的周期數(shù)。吞吐輪次設(shè)計為2周期、3周期、4周期。
以上參數(shù)范圍是根據(jù)長8儲層實際參數(shù)水平確定的。各實驗因素與設(shè)計水平值如表1所示。
表1 CO2吞吐參數(shù)與水平
采用正交試驗設(shè)計方法設(shè)計方案,正交試驗以概率論、數(shù)理統(tǒng)計等為理論基礎(chǔ),從綜合測試中選擇一些具有均勻分散性和均質(zhì)性特點的點組合成方案進行測試,這些方案具有代表性、科學(xué)性,能正確分析試驗結(jié)果,定性定量地確定參數(shù)對指標的影響趨勢、主次順序及顯著程度。根據(jù)四因素三水平正交試驗設(shè)計表格,得到工作參數(shù)的實驗方案,共9個實驗方案,如表2所示。
表2 工作參數(shù)正交試驗方案
CO2吞吐涉及油、氣、水三相,采用Eclipse組分模型進行模擬。采用延長SHH區(qū)塊粗化所得地層模型,如圖1所示,網(wǎng)格維數(shù)為129×155×18,平面網(wǎng)格尺寸為18 m×7 m。設(shè)定油藏壓力為100 MPa,達不到該區(qū)混相壓力190 MPa,采用吞吐開發(fā)與定井底壓力的方式進行開采。根據(jù)物性測試結(jié)果,將長8的油藏流體擬合并為8個組分:CO2、C1、C2、C3、C4-6、C7+1、C7+2、C7+3,各組分初始含量如表3所示。并在實驗室條件下測得各組分的臨界溫度、臨界壓力、壓縮因子、偏差因子、質(zhì)量分數(shù)、組分比容等物性參數(shù),采用三參數(shù)彭-羅賓遜(Peng Robinson,PR)方程進行相平衡計算。根據(jù)長8儲層原油組成和流體性質(zhì),通過室內(nèi)模擬實驗獲取油水相對滲透率、油氣相對滲透率、傳導(dǎo)率等參數(shù)。進行參數(shù)優(yōu)化的模型不考慮天然裂縫與人工壓裂的情形。
圖1 模擬所用地質(zhì)模型Fig.1 The geological model used in the simulation
組分含量/%組分含量/%CO2C1C2C30.017.2110.6719.7C4-6C7+1C7+2C7+320.5914.7516.0111.06
根據(jù)以上條件進行CO2吞吐組分模型建立,通過E300模塊模擬計算。
將以上設(shè)計好的實驗方案輸入Eclipse的生產(chǎn)制度模塊,在Excel中統(tǒng)計各方案運行得到的累積產(chǎn)油量,各方案結(jié)果如表4所示。采用直觀評價法進行參數(shù)的初步優(yōu)化,即繪制方案與產(chǎn)油量關(guān)系曲線,優(yōu)選出產(chǎn)量最高的一組方案。在該方案的基礎(chǔ)上,繼續(xù)運用Eclipse進行研究,控制除所要研究的變量以外參數(shù)為最優(yōu)方案參數(shù)值,對研究參數(shù)進行修改,統(tǒng)計運行結(jié)果,繪制各參數(shù)與累積產(chǎn)油量關(guān)系曲線,找出累積產(chǎn)油量變化明顯的參數(shù)值,并考慮時間成本,優(yōu)選單一工作參數(shù),作為參數(shù)最優(yōu)值,獲得最佳的生產(chǎn)制度。
表4 累積產(chǎn)油量
根據(jù)以上各方案的累積產(chǎn)油量結(jié)果,首先優(yōu)選出累積產(chǎn)油量最高的方案4作為進一步優(yōu)化的基礎(chǔ)工作制度。
控制注氣時長2個月、燜井時間3個月、吞吐輪次2輪,依次設(shè)計注入速度為10 000、20 000、30 000 m3/d。根據(jù)累積產(chǎn)油量結(jié)果(圖2),注入速度從10 000 m3/d變動到20 000 m3/d時,累積產(chǎn)油量上升很快,從388 842.742 m3增長到61 353.156 m3,當注入速度增大到30 000 m3/d時,累積產(chǎn)油量不再增加,甚至略微下降,累積產(chǎn)油量為60 919.195 m3,因此最優(yōu)的注氣速度是20 000 m3/d。
圖2 注入速度與累積產(chǎn)油量關(guān)系曲線Fig.2 Relationship curve between injection rate and cumulative oil production
研究注氣時長時,控制注入速度為20 000 m3/d,燜井時間為3個月,吞吐輪次為2周期,注氣時長依次設(shè)計為2個月、3個月、4個月,累積產(chǎn)油量數(shù)值在60 000 m3浮動,數(shù)值較為穩(wěn)定,沒有突變的趨勢。增加注氣1個月與注氣15 d兩個工作制度,繪制累積產(chǎn)油量關(guān)系曲線(圖3)。當注氣15 d變化到注氣1個月時,累積產(chǎn)油量的數(shù)值有較大幅度的上升,約10 000 m3。因此最佳注氣時長為1個月。由圖3可知,注氣時長小于1個月,地層能量補充不足,當注氣時長大于3個月時,由于注氣量比較大,油被推離井底較遠,生產(chǎn)效果并不會隨之變好。
圖3 注氣時長與累積產(chǎn)油量關(guān)系曲線Fig.3 Relationship curve between gas injection time and cumulative oil production
控制注入速度為20 000 m3/d,注氣時長為2個月,吞吐輪次2周期不改變,燜井時間為2、3個月的累積產(chǎn)油量結(jié)果差別不大,燜井時間4個月時產(chǎn)油量反而下降很快,即燜井時間過長,地層能量散失。初始三水平無法得到最佳燜井時間,增加1個月、15 d、7 d、5 d 4個水平,統(tǒng)計數(shù)值模擬結(jié)果(圖4)。燜井時間從5 d增加到7 d時累積產(chǎn)油量數(shù)值增加明顯,達到62 271.004 m3/d,此后隨著時間的延長,累積產(chǎn)油量圍繞該水平上下浮動,沒有較大幅度改善。因此燜井時間選擇7 d是一個合理的值。
圖4 燜井時間與累積產(chǎn)油量關(guān)系曲線Fig.4 Relationship curve between well time and cumulative oil production
將生產(chǎn)制度設(shè)定為在注入速度為20 000 m3/d、注氣2個月、燜井3個月的條件生產(chǎn)4個周期,導(dǎo)出Eclipse中累積產(chǎn)油量的數(shù)據(jù)與曲線(圖5),將各個輪次下截至各年的累積產(chǎn)油量統(tǒng)計成表格(表5),分析數(shù)據(jù)可知,前兩輪的產(chǎn)油量增加明顯,而進入第3輪時,因為產(chǎn)能逸失嚴重,累積產(chǎn)油量開始停滯不前,而第4輪雖然產(chǎn)油量增加,但是增幅過小,3、4兩周期明顯無法達到盈利目的,因此吞吐周期數(shù)在兩輪比較合適。
圖5 吞吐輪次與累積產(chǎn)油量關(guān)系曲線Fig.5 Relationship curve between throughput and cumulative oil production
周期累積產(chǎn)油量/m3周期累積產(chǎn)油量/m3134 639.350 218.9360 830.760 830.7255 985.460 823.2462 798.965 112.2
(1)從累積產(chǎn)油量與工作參數(shù)關(guān)系可知,累積產(chǎn)油量的值隨CO2吞吐工作參數(shù)變化有明顯的變化,工作參數(shù)對CO2吞吐效果具有重要影響作用。
(2)基于Eclipse數(shù)值模擬結(jié)果,得出長8儲層注入速度最優(yōu)值為20 000 m3/d,注氣時長最優(yōu)值為1個月,燜井時間最優(yōu)值為7 d,吞吐輪次以2周期最為合適。