鄂爾多斯盆地鐵邊城地區(qū)長(zhǎng)2油層組剩余油模擬

張曉明, 吳勤博, 鄒焰, 王皓立, 韓明, 郭龍飛, 狄曉磊, 屈樂

(1.中國(guó)石油長(zhǎng)慶油田分公司第八采油廠, 西安 710069; 2.西安市致密油(頁(yè)巖油)重點(diǎn)實(shí)驗(yàn)室(西安石油大學(xué)), 西安 710065)

低滲透和特低滲透致密砂巖油氣藏分布廣泛,但開發(fā)難度很大[1-2]。近年來,常規(guī)油氣儲(chǔ)量急劇下降,低滲透油氣藏日益受到重視[3-4]。鄂爾多斯盆地是一個(gè)含多套油層的克拉通盆地,其中三疊系延長(zhǎng)組含油性能最好。延長(zhǎng)組發(fā)育大型陸相碎屑巖沉積體系,致密砂巖儲(chǔ)層非均質(zhì)性強(qiáng)[5-6]。在低滲透和特低滲透儲(chǔ)層背景下,強(qiáng)非均質(zhì)儲(chǔ)層剩余油分布規(guī)律和影響因素是一個(gè)研究熱點(diǎn)和難點(diǎn)[7-8]。此外,隨著油氣勘探開發(fā)的滾動(dòng)發(fā)展,獲取剩余油分布特征動(dòng)態(tài)變化規(guī)律,可以為油田可持續(xù)發(fā)展提供重要依據(jù)[9-10]。對(duì)于老油田來說,如何延長(zhǎng)其生命周期,獲得更高的采出程度是一個(gè)科學(xué)問題。要解決該科學(xué)問題,需要結(jié)合靜態(tài)和動(dòng)態(tài)數(shù)據(jù),采用數(shù)值模擬方法對(duì)致密油藏地質(zhì)結(jié)構(gòu)及剩余油分布進(jìn)行精細(xì)定量表征。

鄂爾多斯盆地延長(zhǎng)組頻繁互層的砂泥巖地質(zhì)結(jié)構(gòu)造成致密砂巖儲(chǔ)層具有強(qiáng)非均質(zhì)性,油氣分布規(guī)律復(fù)雜[1-3]。延長(zhǎng)組長(zhǎng)1段到長(zhǎng)10段均有油氣資源分布,但隨著長(zhǎng)時(shí)間的油氣開發(fā),很多地區(qū)的主力產(chǎn)層開始出現(xiàn)高含水,原油以剩余油的形式分布在特低-超低滲儲(chǔ)層中[4-5]。如何有效預(yù)測(cè)和開發(fā)剩余油資源是油田增儲(chǔ)上產(chǎn)及老井挖潛的重要工作[6-7]。鄂爾多斯盆地西北部姬塬油田鐵邊城地區(qū)學(xué)3區(qū)塊的主力產(chǎn)層為長(zhǎng)2油層組,采用注水開發(fā)生產(chǎn)方式進(jìn)行油氣開發(fā)。但是,今年來,長(zhǎng)2油層組逐漸顯現(xiàn)出一些開發(fā)矛盾。例如,受頻繁欠注影響,油藏壓力保持水平下降明顯,導(dǎo)致供液能力下降[8-9]。目前,學(xué)3區(qū)塊長(zhǎng)2油層組平均單井注入壓力13.7 MPa,吸水能力變差,水井頻繁堵塞。此外,局部油層開始出現(xiàn)暴性水淹,治理難度大。近年來,油藏局部地層堵塞頻繁,低產(chǎn)井比例不斷增加,且井間產(chǎn)能差異大。

因此,現(xiàn)利用大量地質(zhì)、測(cè)井資料系統(tǒng)開展長(zhǎng)2油層組剩余油分布數(shù)值模擬研究,以期獲得剩余油分布規(guī)律新認(rèn)識(shí),進(jìn)而為下一步油藏開發(fā)措施調(diào)整提供指導(dǎo)。

1 地質(zhì)背景

圖1 研究區(qū)區(qū)域構(gòu)造位置示意圖

2 長(zhǎng)2油層組儲(chǔ)層特征

2.1 巖性及沉積結(jié)構(gòu)

研究區(qū)長(zhǎng)2段粒度較細(xì),結(jié)構(gòu)成熟度較高,代表遠(yuǎn)源沉積[10-11]。鏡下顯微組分鑒定結(jié)果顯示,所發(fā)育的巖性包括巖屑長(zhǎng)石砂巖及長(zhǎng)石巖屑砂巖(圖2)。學(xué)3區(qū)長(zhǎng)2碎屑成分以長(zhǎng)石為主,其含量為18.0%～56.4%,平均為37.8%;其次為石英,含量為15.5%～48.3%,平均為26.12%。粒徑一般為0.1～0.3 mm,最大粒徑在0.25～0.45 mm。

圖2 學(xué)3區(qū)長(zhǎng)2段巖石類型三角圖

對(duì)長(zhǎng)2油層組巖心的沉積結(jié)構(gòu)進(jìn)行了系統(tǒng)的觀察,部分巖心樣品的照片結(jié)果見圖3。

圖3 長(zhǎng)2油層組砂巖儲(chǔ)層沉積構(gòu)造特征

水平層理代表穩(wěn)定的沉積環(huán)境[圖3(a)],可見明顯具有不同粒度的沉積所形成的層理結(jié)構(gòu)[12-13]。該類砂體通常形成于安靜的水體中,表明同沉積期構(gòu)造活動(dòng)不活躍;板狀交錯(cuò)層理代表相對(duì)淺層砂體發(fā)生不斷遷移而形成[圖3(b)]。圖3(b)中的板狀結(jié)構(gòu)層理的角度約為35°,砂體遷移速率較高。

此外,部分巖性中可見蟲孔構(gòu)造[圖3(c)],蟲孔構(gòu)造主要垂直于層理方向擴(kuò)展,代表該層處于一定氧化區(qū),生物活動(dòng)較為活躍;目的層中還可見沙紋交錯(cuò)層理,如圖3(d)所示。沙紋結(jié)構(gòu)呈羽狀透鏡狀擴(kuò)展,組成砂泥頻繁互層的復(fù)雜結(jié)構(gòu)[14-15]。沙紋結(jié)構(gòu)代表該區(qū)域水動(dòng)力條件較強(qiáng),其通常不利于有機(jī)質(zhì)的保存,但該區(qū)域砂體的沖刷強(qiáng)度較高,砂體的物性可能會(huì)得到一定程度的改善。復(fù)雜的古沉積環(huán)境造成了目的層砂體豐富的沉積構(gòu)造特征。

2.2 物性及孔隙類型

圖4 長(zhǎng)2段孔隙類型分布圖

原生粒間孔主要發(fā)育與石英顆粒之間,且外圍多包括有綠泥石[圖5(a)]。目的層中的溶蝕作用也較為強(qiáng)烈,由圖5(b)和圖5(c)可見溶蝕作用在長(zhǎng)石顆粒中極為強(qiáng)烈,溶蝕作用破壞了顆粒內(nèi)部結(jié)構(gòu),形成連通性較強(qiáng)的孔隙團(tuán)[16-17]。因而,這類孔隙的連通性通常較好。

圖5 長(zhǎng)2段孔隙組合類型

另外,儲(chǔ)層內(nèi)部填隙物中也可觀察到大量微孔隙[圖5(d)],這類孔隙尺寸較小,且聯(lián)通性差;只有當(dāng)發(fā)育一些微裂縫時(shí),孔隙間的聯(lián)通性才有一定程度提高[17-18]。

3 三維地質(zhì)建模

3.1 建模工作流程

3.2 三維地質(zhì)結(jié)構(gòu)模型的建立

確定學(xué)3地區(qū)地質(zhì)模型的5個(gè)模擬單元。建模是從深度域做起,建模邊界條件:精細(xì)模型工區(qū)面積6.42 km2,工區(qū)內(nèi)井?dāng)?shù)128口。

3.3 儲(chǔ)層屬性模型

采用相控參數(shù)場(chǎng)隨機(jī)模擬方法,即在建立巖相邊界控制下進(jìn)行儲(chǔ)層屬性隨機(jī)模擬(圖6)。砂體展布采用序貫指示方法實(shí)現(xiàn),儲(chǔ)層物性參數(shù)則采用序貫高斯方法來實(shí)現(xiàn)。

圖6 目的層的巖相模型

3.3.1 巖相模型

利用測(cè)井曲線中劃分的巖相,劃分的標(biāo)準(zhǔn)是對(duì)于有儲(chǔ)集能力和滲透能力的砂巖劃分為砂層,無(wú)儲(chǔ)集能力和滲透能力的劃分為泥巖[22-25]。通過這樣劃分可以很好地反映儲(chǔ)層的空間分布,有利于控制儲(chǔ)層物性模型的計(jì)算。

依據(jù)學(xué)3地區(qū)塊測(cè)井綜合處理后的小層成果可以較好地劃分開砂巖、泥巖(圖6)。對(duì)網(wǎng)格數(shù)據(jù)作粗化處理,即為網(wǎng)格單元產(chǎn)生一個(gè)有代表性的數(shù)據(jù)值。

孔隙度、含油飽和度為標(biāo)量,并且在三維空間中的變化不大,因而粗化方法采用算術(shù)平均值;而滲透率為空間矢量,并且其數(shù)據(jù)在三維空間下變化劇烈,故使用幾何平均進(jìn)行數(shù)據(jù)粗化[25-27]。

變差函數(shù)分析被用來確定對(duì)砂體縱橫向展布進(jìn)行約束,并確定砂巖層變程方向的物源方向,變差類型則選擇球型,而對(duì)于塊金值人為的給定為0。這樣可以使計(jì)算結(jié)果更符合地質(zhì)規(guī)律。

3.3.2 儲(chǔ)層物性模型

目的層物性參數(shù)的垂向分布特征如圖7所示。協(xié)克里金函數(shù)被用來進(jìn)行物性參數(shù)差值處理。為了提高所建立模型的準(zhǔn)確度,利用測(cè)井解釋的孔、滲及含水飽和度結(jié)果對(duì)模型進(jìn)行約束。

圖7 學(xué)76-75井～學(xué)76-86井物性參數(shù)模擬剖面

3.4 模型驗(yàn)證

根據(jù)地層沉積特征、砂體展布特征、孔、滲等儲(chǔ)層參數(shù)分布特征等地質(zhì)認(rèn)識(shí)在縱橫向的分布規(guī)律,通過不斷調(diào)節(jié)模型使得模擬結(jié)果與實(shí)際地質(zhì)認(rèn)識(shí)最為匹配。

對(duì)比了模擬的孔、滲、含水飽和度結(jié)果與模擬結(jié)果之間的誤差,結(jié)果見圖8,其中包括輸入的測(cè)井曲線值、離散后的測(cè)井解釋結(jié)果及模擬值,這3個(gè)值之間存在一定誤差。

圖8 目的層各物性參數(shù)的模擬結(jié)果與測(cè)井解釋結(jié)果的對(duì)比

測(cè)井解釋結(jié)果是經(jīng)過實(shí)測(cè)值校正之后的結(jié)果,具有高精度特征,是單井縱向儲(chǔ)層識(shí)別的主要依據(jù)。測(cè)井離散值代表對(duì)模型網(wǎng)格賦值后各單元的實(shí)際值,其通常略高于或低于實(shí)際測(cè)井值。而模擬值代表了數(shù)值模擬后的預(yù)測(cè)結(jié)果。

將模擬結(jié)果與原始測(cè)井值進(jìn)行了對(duì)比,結(jié)果顯示,兩者的誤差小于5%(圖8)。因此,兩個(gè)結(jié)果極為一致,表明模型預(yù)測(cè)結(jié)果是可靠的。

4 剩余油分布規(guī)律研究

4.1 網(wǎng)格劃分

采用ECLIPSE軟件開展油藏?cái)?shù)值模擬。數(shù)值模擬網(wǎng)格為笛卡爾坐標(biāo)系下的等距網(wǎng)格。所劃分的網(wǎng)格密度為20 m×20 m,縱向網(wǎng)格大小為1.5 m,劃分為32個(gè)。學(xué)3區(qū)模型形成了一個(gè)165×137×32的網(wǎng)格系統(tǒng),網(wǎng)格節(jié)點(diǎn)總數(shù)為714 523個(gè),且所選的網(wǎng)格類型為角點(diǎn)網(wǎng)格。

4.2 生產(chǎn)動(dòng)態(tài)指標(biāo)擬合

油藏?cái)?shù)值模擬需要進(jìn)行產(chǎn)能的生產(chǎn)歷史擬合。主全區(qū)的產(chǎn)油及產(chǎn)水?dāng)M合結(jié)果見圖9,可以看出,擬合結(jié)果與實(shí)測(cè)結(jié)果基本上是一致的。擬合過程中,充分考慮儲(chǔ)層的特性,調(diào)整砂體、物性等參數(shù)到合理的范圍內(nèi)。

圖9 研究區(qū)整體日產(chǎn)油及含水率擬合曲線

4.3 剩余油分布規(guī)律

圖10 長(zhǎng)層及長(zhǎng)層剩余油分布圖

目前, 隨著中國(guó)大多數(shù)油田進(jìn)入開發(fā)中后期,如何有效開發(fā)剩余油成為下一步增儲(chǔ)上產(chǎn)的重要研究方向,其也是當(dāng)前油田實(shí)現(xiàn)穩(wěn)產(chǎn)的熱點(diǎn)問題。在利用大量地質(zhì)及測(cè)井解釋結(jié)果基礎(chǔ)上,利用三維模擬方法實(shí)現(xiàn)了相控建模及剩余油的精細(xì)刻畫。進(jìn)而得出剩余油的分布規(guī)律:微相控制剩余油分布方向、物性+注水+井網(wǎng)控制剩余油富集差異的認(rèn)識(shí)。

基于上述認(rèn)識(shí),建議在剩余油分布區(qū)部署加密井,提高剩余油的動(dòng)用程度。對(duì)于初期產(chǎn)量相對(duì)較高,目前多為低產(chǎn)井的剩余油分布區(qū),宜對(duì)油井實(shí)施小規(guī)模壓裂,提高單井產(chǎn)量,進(jìn)而使得縱向剩余油得到充分動(dòng)用。因此,研究對(duì)后期致密油藏開發(fā)方案的科學(xué)制定具有極為重要的意義。

5 結(jié)論

(2) 根據(jù)相控原則,建立了主力產(chǎn)層的物理屬性模型及流體模型;進(jìn)而,開展了剩余油分布數(shù)值模擬研究。物性模擬結(jié)果與測(cè)井結(jié)果誤差低于5%,模擬結(jié)果可靠。通過基于實(shí)測(cè)產(chǎn)水及產(chǎn)油量約束的歷史生產(chǎn)曲線擬合,獲得了剩余油的縱橫向分布結(jié)果。

(3)研究顯示,剩余油分布受砂體邊界、物性差異及注采關(guān)系的影響。剩余油整體上成片分布,但是受油藏物性差異的影響,注水多呈不均勻受效,因而造成局部剩余油富集。剩余油主要分布在砂體邊界處,這些區(qū)域無(wú)井網(wǎng)控制,因而水驅(qū)程度低,剩余油分布較多。受長(zhǎng)期注采的影響,井間會(huì)形成高滲區(qū),這些其余類似高速通道,會(huì)發(fā)生快速水淹,其周邊則形成剩余油滯留區(qū)。