喻 鵬,盧宗盛
(中國地質大學,湖北 武漢 430074)
?
高凝油油藏剩余油分布規(guī)律研究
喻 鵬,盧宗盛
(中國地質大學,湖北 武漢 430074)
針對常規(guī)地質統(tǒng)計技術及網(wǎng)格粗化近似技術在油藏剩余油飽和度求取過程中存在較大誤差的問題,優(yōu)選性能穩(wěn)、算速快以及無網(wǎng)格定向效應的流線數(shù)值模擬方法,利用水驅開發(fā)油藏三維流線模型研究工區(qū)流線分布與剩余油配置關系。由研究結果可知,流線推進特征時間同層系調(diào)整時間吻合,小層剩余油分布與水線推進規(guī)律匹配,井組剩余油分布同沉積、斷層以及溫度等多因素協(xié)同作用相關。為此,建議進一步完善工區(qū)注采井網(wǎng)系統(tǒng),并考慮在局部井組采用針對性的熱水、蒸汽驅替,克服原油溫敏效應,以提高采收率。通過研究,明確出油及來水方向,找到當前注采流線未能波及的剩余油富集潛力區(qū),為下一步調(diào)整井網(wǎng)挖掘剩余油指明了方向。
流線模擬;配置關系;數(shù)值仿真;剩余油分布;流線分布;高凝油;靜觀2塊
高凝油油藏作為諸多油藏類型中較為獨特的一類,具有高凝固點、高含蠟、高析蠟溫度的“三高”特性,同時原油物性對溫度因子極其敏感,故研究其水驅過程的同時,還應考慮原油自身的性質[1-5]。目前,大部分常規(guī)水驅開發(fā)高凝油油藏均存有溫度敏感導致的滲流場紊亂以及剩余油規(guī)律認識不清的情況,而此次的研究區(qū)塊——靜觀2塊也是靜安堡復式油氣聚集區(qū)帶中段的一個主力高凝油斷塊,原油含蠟量為30.0%~45.0%,平均為35.0%,凝固點為40~50 ℃,含硫量小于0.1%。根據(jù)旋回特點及巖、電特性,將油藏縱向劃為12個油層組,40個砂巖組,93個小層。
油藏自投產(chǎn)以來,共經(jīng)歷5次重大調(diào)整,開發(fā)層系及井網(wǎng)的調(diào)整使得油水關系變得更加復雜,剩余油分布規(guī)律差,大分層量導致粗化模型的網(wǎng)格節(jié)點增多,給仿真計算過程帶來了較大阻力。由于常規(guī)地質統(tǒng)計技術及網(wǎng)格粗化近似技術在油藏剩余油飽和度求取過程中可能存在較大誤差,尤其針對該類多節(jié)點、大計算量的數(shù)據(jù)體,流線模擬器較傳統(tǒng)模擬器具有性能穩(wěn)、算速快以及無網(wǎng)格定向效應優(yōu)勢[6-13],故優(yōu)選流線模擬手段進行工區(qū)注采流線分布與剩余油的配置關系研究,旨在明確出油及來水方向,明確區(qū)塊及井網(wǎng)剩余油分布情況,為油田下一步開發(fā)調(diào)整提供一定的指導。
充分考慮高凝油油藏獨特的開發(fā)特征,結合實驗測試結果可知,常規(guī)水驅進程中,蠟晶、瀝青的析出及毛管作用使得滲流體系阻力劇增,故耦合流線-黑油模型時,將流體滲流機理一并考慮在內(nèi),同時給出流線模型數(shù)學模型的假設:油藏流體滲流遵循達西規(guī)律;滲流進程溫度恒定;巖石及流體皆不可壓縮;氣相對驅替進程干擾權重偏小,對其忽略。
引入毛管及重力項的油水兩相流線法滲流數(shù)學模型[14-16]:
▽{[ρoKKro▽(po-ρogD)]/μo}+qo=φ[?(ρoso)/?t]
(1)
▽{[ρwKKrw▽(pw-ρwgD)]/μw}+qw=φ[?(ρwsw)/?t]
(2)
式中:ρo、ρw分別為油、水密度,kg/m3;K為儲層絕對滲透率,10-3μm2;Kro、Krw分別為油相、水相相對滲透率;μo、μw分別為油相、水相黏度,mPa·s;D為與重力加速度方向一致并按某一基準面的測深,m;qo、qw分別為油、水相單位體積、時間內(nèi)巖石質量流量,kg/(m3·s)。
2.1 物性及動態(tài)參數(shù)設置
結合擬合精度及實際生產(chǎn)需求,在精細地質建?;A上,選取52×27×93=130 572的模型格架(網(wǎng)格步長為50 m),數(shù)據(jù)準備如下。
(1) 模型基本參數(shù)及其初始化。根據(jù)物性表征模型及工區(qū)地質信息,運用端點標定、初始含水調(diào)用結合法平衡初始場,配置壓力參數(shù),孔滲參數(shù)直接對接調(diào)用,壓縮系數(shù)為6.872×10-4MPa-1。
(2) 實際地質信息匯總。清繪工區(qū)儲層及滲流特征的圖件(如構造、砂厚等),同時建立以井點為源的模型信息庫。
(3) 流體特征參數(shù)。結合高壓物性(含原油、地層水)資料成果,架構起特征參數(shù)模擬模型。
(4) 滲流特征參數(shù)。包含毛管及相滲數(shù)據(jù),進而結合生產(chǎn)動態(tài),在計算模擬進程中微調(diào)擬合相滲。
(5) 動態(tài)參數(shù)預處理。前處理生產(chǎn)數(shù)據(jù)時,結合報錯調(diào)整非法步長,完善壓力、射孔。
2.2 流線模擬參數(shù)擬合調(diào)節(jié)
合理的參數(shù)設置可將模擬進程中的不收斂及平衡誤差故障降低,故結合經(jīng)驗及報錯等對若干關鍵字進行調(diào)整。設置制約輸出結果,干預飽和度和壓力調(diào)節(jié),交互式控制誤差及壓頻,反復調(diào)整增加默認平衡誤差,設置最大牛頓迭代次數(shù),保證了模型計算進程的高效性。
2.3 流線模型歷史擬合
擬合階段設置定油生產(chǎn),完成產(chǎn)油量擬合,同時以定壓邊界和定液注入井的方法確保流線計算的穩(wěn)定性。
擬合步驟:①區(qū)塊擬合。計算過程結合非平衡法穩(wěn)定含水的同時,給定壓—深關系模擬壓力系統(tǒng),儲量自然擬合,模擬石油地質儲量為724.58×104t,同工區(qū)實際誤差為1.60%,擬合效果理想。同時,進行含水及產(chǎn)液的擬合,綜合含水擬合狀況較好,見水時間與產(chǎn)液的擬合效果也較好;②結合地質資料反復微調(diào)相滲、NTG及壓力等不確定參數(shù),使單井最終擬合率增至83.10%,達到擬合精度。
3.1 油藏流線分布與剩余油配置關系
靜觀2塊投入開發(fā)至今,共歷經(jīng)5次重大調(diào)整,開發(fā)層系以及井網(wǎng)密度等變化使得模擬驅替進程當中的流線特征變化明顯。如1990至1991年第2次部署期間,為緩解油藏層間矛盾,以原二層系為基礎,將下層系進一步細分成2個開發(fā)層系,加密井網(wǎng)。從部署后期流線特征上看(圖1),下層系水線推進明顯增強,主力片區(qū)掃油效率提高。
圖1 靜觀2塊不同開發(fā)部署期水驅注采流線分布
3.2 小層流線分布與剩余油配置關系
結合水淹等資料,隨機選取流線模型小層模擬結果,以檢驗剩余油與流線的配置關系。高含水層中,若干區(qū)域未貫穿流線或密集程度偏低,預測含水飽和度與水淹情況匹配度較高,印證了流線模型的合理性。結合井史資料分析,主力井區(qū)內(nèi),由于井網(wǎng)射孔不完善,若干層非正常剩余油富集,使得目的層進液不足。工區(qū)東北高部原油富集,為斷塊構造使水驅流線抬升范圍變窄造成,南部和西北片的富集油則由于斷層遮擋、干擾了掃油率造成。
3.3 井組流線分布與剩余油配置關系
工區(qū)各井組和剩余油有不同的配置關系和分布規(guī)律,疊加所有關系即構成整個注采系統(tǒng)和剩余油的配置關系。以上述模擬層為例,觀察各井組同剩余油的分布關系。
0017井注水流線呈傾向性聚簇狀偏南西(圖2)。從構造上看,距離該井北東向不遠存在斷塊的3號斷層,這種特殊的聚簇特征主要受控于斷層封隔作用,也正是由于該作用,使得此井組附近屋檐處有一定量剩余油賦存。
圖2 0017井組水驅流線及對應剩余油分布
由圖3可見,0002—0031井均分布于水下河道,且井距較小,注水驅替構成滲流優(yōu)勢通道,水驅過程中,優(yōu)勢流線簇構成橢圓狀水線,故井區(qū)北東部剩余油富集程度明顯高于南西部。
圖3 0002井組水驅流線及其對應剩余油分布
0022井是工區(qū)儲層傷害實驗選井,其注水域內(nèi)儲層冷傷害嚴重,冷水驅替使原油敏感,溫度降低,析出蠟晶,將孔道堵塞,干擾水驅進程。從模擬結果看(圖4),0022井水驅流線前緣以0.07 m/d的速度推進,但仍未在0012井突破,極高的溫度敏感性致使驅油速度受阻,井區(qū)周邊剩余油零星富集,尤其在未波及井周圍,剩余油存在連片富集態(tài)。
通過以上分析,將工區(qū)油藏仿真流線與剩余油的配置情況歸納為五大類:第1類為流線密集度較高的低剩余油飽和度分布;第2類為井網(wǎng)不完善造成若干小層的剩余油富集;第3類為構造條件差流線難以波及的區(qū)域剩余油富集;第4類為相控條件下層中矛盾突出造成的局部剩余油富集;第5類為高凝油冷采過程中溫敏效應導致流線波及受阻的剩余油富集。結合工區(qū)流線分布及與剩余油的配置關系,可以后續(xù)挖掘潛剩余油潛力區(qū),以此可以提高油藏的最終采收率。
圖4 0022井組水驅流線及對應剩余油分布
(1) 結合流線數(shù)值模擬技術手段,研究分析高凝油油藏注采流線分布與剩余油的配置關系,明確出油及來水方向,找到當前注采流線未能波及的剩余油富集潛力區(qū),為下一步調(diào)整現(xiàn)有注采井網(wǎng)、挖掘剩余油指明方向。
(2) 模型流線推進變化特征時間同層系調(diào)整時間吻合,隨機小層剩余油分布及水線推進規(guī)律匹配,總體擬合度高,最終流線模型合理。通過井組流線分布與剩余油配置關系分析,明確了剩余油分布與沉積、斷層以及溫度等多因素協(xié)同作用的關系,可進一步完善注采井網(wǎng)。
[1] 姜彬,丘凌,劉向東,等.固相沉積模型在高凝油藏注水開發(fā)中的應用[J].石油學報,2015,36(1):101-105.
[2] 陽曉燕,楊勝來,吳向紅,等.高凝油油藏注水開采溫度場變化規(guī)律[J].特種油氣藏,2011,18(4):87-89.
[3] 謝文彥,李曉光,陳振巖,等.遼河油區(qū)稠油及高凝油勘探開發(fā)技術綜述[J].石油學報,2007,28(4):145-150.
[4] 鄧惠,楊勝來,康銘輝,等.高凝油常規(guī)冷采時井筒溫度分布分析[J].特種油氣藏,2008,15(5):91-93.
[5] 聶向榮,楊勝來.高凝油油藏冷傷害特征數(shù)值模擬[J].石油鉆探技術,2014,42(1):100-104.
[6] 鄭強,戴雄軍,張有印,等.中高含水期邊底水砂巖油藏開采政策研究——以哈國Aryskum油田M-II層北區(qū)為例[J].新疆石油天然氣,2013,9(3):35-38.
[7] 侯玉培.整裝油藏流場重整提高采收率的數(shù)值模擬研究[D].青島:中國石油大學,2011.
[8] 朱麗紅.特高含水期厚油層剩余油綜合描述及分布特征[J] .大慶石油地質與開發(fā),2013,32(6):58-62.
[9] 黃有泉,周志軍,劉志軍,等.葡南三斷塊特高含水期油藏數(shù)值模擬精度提高方法[J] .油氣地質與采收率,2014,21(5):65-68.
[10] 陳曉彬,石為為.高臺子油層河口壩砂體精細識別及剩余油分析[J].大慶石油地質與開發(fā),2014,33(4):66-69.
[11] 張世明,吳曉東,李坤,等. 基于吸水剖面資料的油藏層間平均剩余油飽和度計算方法[J].油氣地質與采收率,2014,21(5):98-100.
[12] 吳軍來,劉月田,楊海寧.基于 3D 流線模擬的水驅油藏動態(tài)評價新方法[J].西安石油大學學報(自然科學版),2011,26(2):43-47.
[13] 翟亮.基于流線的示蹤劑技術研究薄層底水油藏開發(fā)規(guī)律[J].特種油氣藏,2011,18(1):79-82.
[14] 劉慧卿.油藏數(shù)值模擬方法專題[M].東營:石油大學出版社,2001:20-25.
[15] 孫致學,楚鵬,張青菇,等.微觀驅油實驗數(shù)值模擬方法[J].大慶石油地質與開發(fā),2015,34(3):124-128.
[16] 呂琦.基于流線的油藏數(shù)值模擬研究[D].青島:中國石油大學,2009.
編輯 姜 嶺
20151117;改回日期:20160106
國家科技重大專項“大型油氣田及煤層氣開發(fā)”之“復雜油藏剩余油分布預測”(2011ZX05009-003-03);中國石油天然氣股份有限公司重大科技項目“沈84-安12塊高凝油油藏提高采收率研究”(QW08037005-3)
喻鵬(1987-),男,講師,2008年畢業(yè)于黃山學院數(shù)學與應用數(shù)學專業(yè),2014年畢業(yè)于中國地質大學(武漢)石油與天然氣工程專業(yè),獲博士學位,現(xiàn)從事油氣田開發(fā)方面的教學和科研工作。
10.3969/j.issn.1006-6535.2016.02.019
TE319
A
1006-6535(2016)02-0081-04