孫永亮 陳 琦 張 澤 李國(guó)良 房金偉 南雪芬
(①中國(guó)石油渤海鉆探工程有限公司油氣合作開發(fā)分公司;②中國(guó)石油冀東油田公司陸上作業(yè)區(qū))
蘇里格氣田位于鄂爾多斯盆地,是一個(gè)低壓、低孔、低滲透、低豐度,以河流砂體為主體儲(chǔ)層的大面積分布的巖性氣藏,儲(chǔ)層非均質(zhì)性強(qiáng),氣井能量衰竭快,有效波及范圍小、儲(chǔ)量動(dòng)用程度低。蘇X 區(qū)塊經(jīng)過(guò)多年開發(fā),目前已開發(fā)富集區(qū)內(nèi)井網(wǎng)逐漸完善,可供鉆探的目標(biāo)很少。外圍區(qū)地質(zhì)條件差,屬于低品位儲(chǔ)量區(qū),儲(chǔ)量豐度低,儲(chǔ)層橫向相變快,可動(dòng)水飽和度高,氣井產(chǎn)量普遍較低,開發(fā)效果差。國(guó)內(nèi)外開發(fā)實(shí)踐表明,井網(wǎng)加密是致密氣藏提高采收率的有效手段之一。本文圍繞蘇X 區(qū)塊富集區(qū)提高采收率的生產(chǎn)需求,采用數(shù)值模擬技術(shù)手段,論述剩余氣分布、井網(wǎng)加密方案及現(xiàn)場(chǎng)應(yīng)用的研究成果。
蘇X 區(qū)塊位于蘇里格氣田中部,鄂爾多斯盆地伊陜斜坡的西北側(cè),總體為一區(qū)域性西傾大單斜,局部發(fā)育一些鼻狀隆起的小型構(gòu)造。平面上廣覆式生烴、上古生界儲(chǔ)集巖大面積分布,源下儲(chǔ)上配置關(guān)系較好,在二疊系石盒子組和山西組形成了大面積含氣區(qū)[1]。盒8 段、山1 段沉積為緩坡型三角洲沉積體系,具有典型的河流相沉積特征,河道砂體縱向疊置,橫向復(fù)合連片,砂體厚40~60 m,氣層厚平均9 m;儲(chǔ)層非均質(zhì)性較強(qiáng),有效砂體規(guī)模小,含氣性橫向變化快[1]。該區(qū)塊為典型“三低”致密氣藏。低滲:孔隙度5%~14%,滲透率0.01~0.1 mD,裂縫不發(fā)育。低壓:主力層埋深2 800~3 700 m,壓力系數(shù)0.70~0.98。低豐度:(1.1~1.5)×108m3/km2。儲(chǔ)層孔喉結(jié)構(gòu)復(fù)雜,氣體充注程度低,地層水以自由水、滯留水和束縛水的狀態(tài)賦存,含水飽和度相對(duì)較大,平均40%,氣井產(chǎn)量普遍較低。
數(shù)值模擬的研究目標(biāo)是在地質(zhì)模型的基礎(chǔ)上,利用氣藏靜態(tài)及動(dòng)態(tài)數(shù)據(jù),優(yōu)化數(shù)值模擬模型,通過(guò)擬合氣井各個(gè)生產(chǎn)階段的產(chǎn)出特征,模擬氣水運(yùn)動(dòng)規(guī)律,預(yù)測(cè)剩余氣的定量分布,為下一步氣藏開發(fā)調(diào)整方案提供依據(jù)。
研究區(qū)位置:數(shù)值模擬研究區(qū)位于蘇X 工區(qū)中東部,面積76.7 km2,地質(zhì)儲(chǔ)量氣井158口。
網(wǎng)格尺寸及數(shù)量:以盒8 段、山1 段及山2 段為模擬對(duì)象,模型采用角點(diǎn)網(wǎng)格,平面網(wǎng)格大小為50 m×50 m,數(shù)量164×187 個(gè)。縱向上,盒8 段以單砂體作為一個(gè)網(wǎng)格單元,山1 段和山2 段以小層作為一個(gè)網(wǎng)格單元,共劃分14個(gè)單元格,網(wǎng)格總數(shù)為429 352個(gè)。
儲(chǔ)層參數(shù)場(chǎng)模型:對(duì)三維地質(zhì)建模提供的凈毛比、孔隙度、滲透率、飽和度模型進(jìn)行粗化,粗化后的模型作為數(shù)值模擬的儲(chǔ)層參數(shù)場(chǎng)輸入模型。以開發(fā)方案中的原始靜壓數(shù)據(jù)和回歸公式為依據(jù),分別建立盒8段、山1段和山2段的初始地層壓力模型。
流體和高壓物性參數(shù)參考開發(fā)方案研究成果,主要流體參數(shù)包括:地面條件下天然氣相對(duì)密度為0.613,地面條件下水的密度為1.008 g/cm3,地層水的體積系數(shù)為1.029,地層水的壓縮系數(shù)為3.05×10-5MPa-1,原始油藏條件下地層水的粘度為0.283 mPa·s,巖石壓縮系數(shù)為4.4×10-5MPa-1。PVT 曲線由數(shù)值模擬軟件提供的PVT 模塊生成,該模塊能夠根據(jù)天然氣的相對(duì)密度,自動(dòng)生成不同地層壓力下的天然氣膨脹系數(shù)、天然氣粘度等。
相滲曲線:針對(duì)石盒子組和山西組采用不同的相滲曲線,由于不同物性區(qū)域相滲數(shù)據(jù)的不確定性,實(shí)際生產(chǎn)擬合時(shí)會(huì)有所調(diào)整。精細(xì)計(jì)算砂體含氣面積和凈毛比,可使模型儲(chǔ)量接近地質(zhì)儲(chǔ)量,誤差率小于3%。
油藏生產(chǎn)動(dòng)態(tài)數(shù)據(jù)中,氣井產(chǎn)氣量數(shù)據(jù)相對(duì)齊全,在歷史擬合階段,通常用產(chǎn)氣量作為工作制度參與模擬計(jì)算,而將井底壓力等作為擬合指標(biāo)。在分析整理該氣田生產(chǎn)井歷年來(lái)的開發(fā)綜合數(shù)據(jù)、單井生產(chǎn)數(shù)據(jù)、測(cè)壓資料的基礎(chǔ)上,在模擬軟件中建立油藏動(dòng)態(tài)模擬模型。蘇X 區(qū)塊從2006 年12 月投產(chǎn)至今,氣井生產(chǎn)歷史較長(zhǎng),歷史擬合以1個(gè)月為計(jì)算時(shí)間步長(zhǎng),氣藏不存在邊底水且無(wú)產(chǎn)水計(jì)量,主要對(duì)氣井產(chǎn)量和井底壓力進(jìn)行擬合。
模型初始化:常規(guī)初始化方法有2種,平衡初始化和非平衡初始化。平衡初始化受相滲曲線限制,與地質(zhì)模型初始飽和度場(chǎng)數(shù)據(jù)有差異;非平衡初始化直接加載地質(zhì)模型初始飽和度場(chǎng),模型不穩(wěn)定且收斂性差。本次模擬采用了毛管力標(biāo)定地質(zhì)模型飽和度場(chǎng)、平衡啟動(dòng)初始化方法。首先為模型加載地質(zhì)模型提供的飽和度場(chǎng)數(shù)據(jù),然后應(yīng)用Eclipse軟件的端點(diǎn)標(biāo)定功能,計(jì)算每個(gè)網(wǎng)格標(biāo)定的最大毛管力。這樣既保證了每個(gè)網(wǎng)格的初始飽和度為指定數(shù)值,又采用了平衡法啟動(dòng),通過(guò)標(biāo)定保證了重力和毛管力的平衡[2]。
產(chǎn)氣量擬合:氣藏因壓裂等工程因素破壞了氣井射開層段周圍的儲(chǔ)層,增大了儲(chǔ)層的實(shí)際滲透能力。產(chǎn)氣量擬合時(shí)參考?jí)毫呀Y(jié)果,增大氣井射開層段的地層系數(shù),同時(shí)調(diào)整氣井滲流半徑區(qū)域的滲透率和傳導(dǎo)率,通過(guò)模擬壓裂生產(chǎn),使模擬產(chǎn)氣量接近實(shí)際產(chǎn)氣量。采用儲(chǔ)層分類相滲控制的方法,對(duì)不同層段不同類型儲(chǔ)層進(jìn)行分區(qū),分別設(shè)置對(duì)應(yīng)相滲曲線,提高了擬合精度。目前研究區(qū)日產(chǎn)氣59.8×104m3,已累計(jì)產(chǎn)氣36.5×108m3。
壓力擬合:氣藏壓力主要與壓縮系數(shù)、產(chǎn)氣和產(chǎn)水量有關(guān)。地層水壓縮系數(shù)和巖石壓縮系數(shù)可調(diào)范圍不大。由于致密氣藏低滲透特征和壓裂生產(chǎn)方式,壓力擬合過(guò)程中,除了常規(guī)修改儲(chǔ)層相滲曲線、井的完善系數(shù)等參數(shù),還需要結(jié)合動(dòng)態(tài)分析和試井解釋結(jié)果,調(diào)整氣井滲流半徑區(qū)域的傳導(dǎo)率和表皮因子,反復(fù)修改參數(shù),使單井壓力得到較好擬合(圖1)。
圖1 SX-10-11井產(chǎn)量與壓力擬合
偏差分析:研究區(qū)氣井生產(chǎn)歷史擬合符合率為81.6%。擬合偏差的原因包括:節(jié)流器失效或計(jì)量不準(zhǔn)確導(dǎo)致氣井采氣曲線產(chǎn)生陡峭尖峰,數(shù)模軟件由于時(shí)間步長(zhǎng)等參數(shù)設(shè)置因素,無(wú)法完全反映這種非線性流動(dòng);由于氣井產(chǎn)量低、攜液不足等原因,存在不同程度的井筒積液,導(dǎo)致井底壓力波動(dòng)幅度大,不能真實(shí)反映地層能量。
目前國(guó)內(nèi)外研究剩余氣分布的常用方法主要有以下5 種:取心井巖心分析法、剩余氣飽和度測(cè)井法、油藏?cái)?shù)值模擬法、開發(fā)地質(zhì)與動(dòng)態(tài)綜合分析法、水動(dòng)力學(xué)法。本文采用地質(zhì)研究和數(shù)值模擬相結(jié)合、定性研究與定量表征相結(jié)合的方法研究剩余氣分布。
蘇里格氣田剩余氣分布主要受地質(zhì)和開發(fā)兩大因素的影響。地質(zhì)因素主要指沉積微相、儲(chǔ)層特性、微型構(gòu)造、流體性質(zhì)等,表現(xiàn)為對(duì)剩余氣分布的控制;開發(fā)因素主要指采氣速度、井網(wǎng)分布等,影響剩余氣儲(chǔ)量的動(dòng)態(tài)變化[3]。
沉積因素:沉積微相決定儲(chǔ)集砂體的外部形態(tài)、內(nèi)部結(jié)構(gòu)、平面和垂向非均質(zhì)性,控制著氣、水的運(yùn)動(dòng)方向,從而導(dǎo)致剩余氣沿一定的相帶分布。心灘、邊灘微相砂巖厚度較大,泥質(zhì)含量少,剩余氣飽和度較低;河道微相砂巖厚度相對(duì)小,泥質(zhì)含量較少,剩余氣飽和度較高,但是其剩余儲(chǔ)量豐度和剩余可采儲(chǔ)量較小[4]。
儲(chǔ)層因素:儲(chǔ)層平面非均質(zhì)性強(qiáng)、物性差異大,滲透率大的儲(chǔ)層導(dǎo)流能力強(qiáng),氣藏壓降快,剩余儲(chǔ)量變??;縱向上多砂體合采,非主力薄層低滲砂體產(chǎn)量受到抑制,剩余儲(chǔ)量較大;砂體分布范圍小,氣層對(duì)應(yīng)性差,剩余儲(chǔ)量分布比較分散且豐度較高。
開發(fā)因素:剩余氣分布與井網(wǎng)分布、采氣速度有直接關(guān)系。平面上井網(wǎng)不完善的部位,儲(chǔ)量動(dòng)用程度低,是剩余氣的相對(duì)富集區(qū);采氣速度越快,剩余儲(chǔ)量越小。
應(yīng)用剩余地質(zhì)儲(chǔ)量豐度描述地下剩余氣的分布,綜合考慮了有效厚度、孔隙度、含氣飽和度、天然氣體積系數(shù)等多參數(shù)對(duì)計(jì)算地質(zhì)儲(chǔ)量的影響,克服了應(yīng)用剩余氣飽和度描述剩余氣的片面性,可以反映出地下剩余儲(chǔ)量分布狀況??v向上,將每一小層的剩余地質(zhì)儲(chǔ)量豐度相加,可得到整個(gè)氣藏地下剩余地質(zhì)儲(chǔ)量的分布,為今后的開發(fā)調(diào)整提供可靠依據(jù)。本次研究采用數(shù)值模擬和動(dòng)態(tài)分析相結(jié)合的方法,計(jì)算含氣飽和度、地層壓力、儲(chǔ)量豐度的變化情況,定量表征剩余氣分布。
數(shù)值模擬計(jì)算研究區(qū)剩余地質(zhì)儲(chǔ)量為74.83×108m3,盒8 段、山1 段、山2段剩余儲(chǔ)量分別為39.43×108m3、21.81×108m3、13.59×108m3。剩 余儲(chǔ)量富集區(qū)有5個(gè),分別位于北部、中部、南西部、南中部、南東部(圖2),其中:北部富集區(qū)面積3.60 km2,剩余儲(chǔ)量5.10×108m3;中部富集區(qū)面積3.56 km2,剩余儲(chǔ)量4.77×108m3;南西部富集區(qū)面積7.23 km2,剩余儲(chǔ)量13.10×108m3;南中部富集區(qū)面積4.12 km2,剩余儲(chǔ)量6.73×108m3;南東部富集區(qū)面積8.70 km2,剩余儲(chǔ)量11.65×108m3。
圖2 研究區(qū)剩余儲(chǔ)量豐度分布
綜合現(xiàn)有井網(wǎng)分布、剩余氣分布、有效儲(chǔ)層展布,定性劃分剩余氣類型主要為4種:井網(wǎng)未控制型、復(fù)合砂體內(nèi)阻流帶型、水平井漏失型、射孔不完善型[5]。
井網(wǎng)未控制型:有效砂體規(guī)模小,橫向連通性差,縱向發(fā)育頻率低,平面上以孤立分布為主。
復(fù)合砂體內(nèi)阻流帶型:有效砂體規(guī)模較大,橫向氣藏連通性強(qiáng),縱向發(fā)育頻率高,平面連片分布,復(fù)合砂體內(nèi)部不連通,垂直水流方向發(fā)育多個(gè)“阻流帶”,影響砂體滲流能力和儲(chǔ)量動(dòng)用程度,形成一定規(guī)模的剩余氣。
水平井漏失型:水平井利用多段壓裂改造技術(shù)提升了有效砂體的儲(chǔ)量動(dòng)用程度,但多層含氣的地質(zhì)特點(diǎn)使其不可避免地遺漏縱向上部分層段的儲(chǔ)量。
射孔不完善型:受開發(fā)早期直井分層壓裂技術(shù)限制,部分差氣層射孔不完善或壓裂改造不完善形成了剩余氣[6]。
蘇X 區(qū)塊含氣富集區(qū)目前井網(wǎng)井距為600 m×1 000 m,且已基本完善,區(qū)塊中北部和西南部氣水關(guān)系復(fù)雜。隨著蘇X 區(qū)塊進(jìn)一步開發(fā),區(qū)塊井位部署難度逐年增大。為了保持區(qū)塊穩(wěn)產(chǎn),有必要開展蘇X 區(qū)塊富集區(qū)井網(wǎng)加密地質(zhì)技術(shù)研究,為區(qū)塊保持穩(wěn)產(chǎn)提供技術(shù)支撐。
井距、排距:加密井網(wǎng)井距、排距可借鑒蘇里格氣田其他區(qū)塊干擾試驗(yàn)結(jié)果。井距≤400 m,干擾概率為71%,井距>600 m,基本無(wú)干擾,故井距應(yīng)設(shè)置在500 m 左右。排距≤700 m,干擾概率為38.5%,排距≤600 m,干擾概率為44.4%,故排距應(yīng)大于600 m。目前蘇X 區(qū)塊基礎(chǔ)井網(wǎng)600 m×1 000 m,且已較完善,根據(jù)目前現(xiàn)有井網(wǎng)形勢(shì),采用平行四邊形加密方式在排距間加密。
井網(wǎng)未控制型:試驗(yàn)區(qū)位于工區(qū)中東部SX-11-11S井區(qū),地質(zhì)儲(chǔ)量13.57×108m3,采出程度25.8%。加密試驗(yàn)區(qū)砂體分布范圍小,氣層對(duì)應(yīng)性差,剩余儲(chǔ)量分布分散,適合開展井網(wǎng)加密[7]。井區(qū)加密井位部署見(jiàn)圖3,部署加密直井8 口,井網(wǎng)600 m×550 m,井網(wǎng)密度2.80 口/km2。設(shè)置數(shù)值模擬限制條件為年有效生產(chǎn)天數(shù)330 d,井底流壓2.9 MPa,經(jīng)濟(jì)極限產(chǎn)氣量1 000 m3/d,模擬計(jì)算采收率為39.4%,比不加密方案提高了8.2%(圖4)。
圖3 SX-11-11S井區(qū)加密井位部署
圖4 SX-11-11S井區(qū)兩種方案模擬預(yù)測(cè)采收率對(duì)比
復(fù)合砂體內(nèi)阻流帶型:試驗(yàn)區(qū)位于工區(qū)南東部,地質(zhì)儲(chǔ)量12.85×108m3,采出程度26.3%。針對(duì)盒8段復(fù)合連片含氣砂體部署加密直井7 口,井網(wǎng)600 m×(550~650)m,井網(wǎng)密度2.56 口/km2。模擬計(jì)算采收率為37.9%,比不加密方案提高了7%。
水平井漏失型:試驗(yàn)區(qū)位于工區(qū)南東部,針對(duì)水平井開發(fā)井網(wǎng)。預(yù)測(cè)氣層厚度20 m,水平井動(dòng)用8 m,部署加密直叢井4口,采收率預(yù)測(cè)提高了6.0%。
SX-11-11S 井區(qū)位于蘇X 區(qū)塊中東部富集區(qū),目的層盒8段及山西組處于有利沉積相帶且分布基本穩(wěn)定,有效砂體分布范圍有限,有效儲(chǔ)層物性好、厚度大,縱向氣層分散。該井區(qū)已鉆直井15 口,動(dòng)用地質(zhì)儲(chǔ)量8.26×108m3,剩余地質(zhì)儲(chǔ)量5.31×108m3。各氣井投產(chǎn)時(shí)間均為2007 年,單井累產(chǎn)氣量均超過(guò)1 000×104m3,整體生產(chǎn)效果較好。2020 年該井區(qū)實(shí)施了8 口加密井,動(dòng)用地質(zhì)儲(chǔ)量4.36×108m3。單井氣層厚度7.2~32.0 m,平均15.3 m,單井無(wú)阻流量(10.14~37.74)×104m3/d,平均13.5×104m3/d,單井初期產(chǎn)量平均1.0×104m3/d,動(dòng)靜態(tài)綜合評(píng)價(jià)Ⅰ類井比例87.5%。其中7口井效果良好(表1)。
表1 SX-11-11S井區(qū)加密井實(shí)施效果統(tǒng)計(jì)
針對(duì)蘇里格氣田蘇X 區(qū)塊致密儲(chǔ)層低滲透特征和壓裂生產(chǎn)方式,優(yōu)化了數(shù)值模擬方法,采用毛管力標(biāo)定平衡初始化、儲(chǔ)層分段分類相滲控制的方法,調(diào)整滲透率、傳導(dǎo)率、表皮因子等模型參數(shù),精細(xì)擬合氣井產(chǎn)量和井底流壓,進(jìn)而模擬地層壓力和剩余儲(chǔ)量豐度變化,定量表征剩余氣分布,定性劃分剩余儲(chǔ)量類型為井網(wǎng)未控制型、復(fù)合砂體內(nèi)阻流帶型、水平井漏失型、射孔不完善型4種。針對(duì)不同類型剩余儲(chǔ)量,分別制定加密方案,預(yù)測(cè)采收率提高了6.0%~8.2%。優(yōu)選SX-11-11S 試驗(yàn)井區(qū)實(shí)施8 口加密井,綜合評(píng)價(jià)Ⅰ類井占比為87.5%,效果良好。