徐 慧,林承焰,鄭元超,賈珍臻,丁 圣
(1.中國石油大學,山東 青島 266580;2.克拉瑪依職業技術學院,新疆 克拉瑪依 833600;3.中石化江蘇油田分公司,江蘇 揚州 225009)
水下分流河道為陸上分流河道的水下延伸部分,是淺水三角洲沉積體系中的主力層[1-2],一般由細砂巖、粉砂巖和泥質粉砂巖組成,表現出正韻律或正韻律疊加特征[3]。該類砂體在平面上分布形態較為復雜,對井網及注采系統的優化程度要求較高,且平面挖潛的難度較大[4]。油藏數值模擬技術是油田開發后期認識剩余油分布的科學、實用的技術手段[5]。因此,研究水下分流河道砂體的油藏數值模擬方法,對描述剩余油分布規律、指導下一步挖潛具有重要的意義。
地質模型的精細程度直接影響著油藏數值模擬的精度,目前儲層內部建筑結構劃分技術與流動單元劃分技術[6-8]已趨于成熟,并且流動單元約束的油藏數值模擬模型具有初始擬合程度高、歷史擬合時間短、模型修改程度低、能更精確地反映層內和平面剩余油分布的特點[9]。筆者以大慶油田葡萄花油層葡Ⅰ油組為例,在劃分單河道砂體的基礎上,劃分流動單元,進而建立了精細地質模型。

圖1 S29-J17井巖心分析
綜合巖性、電性特征,以 Miall[10]的層次劃分思想為指導,識別層次界面。其中,Ⅲ級界面為單一分支河道砂體的內部夾層,屬于自旋回細粒沉積。Ⅳ級界面為一期水下分流河道中單河道砂體的邊界面。Ⅴ級界面為一期河道形成的復合砂體邊界面(圖1)。共劃分27個單河道砂體,其中主力3小層劃分為3個單河道砂體。單河道間夾層較薄,平面延伸范圍較小,一般不超過1個井距,少 數可以達2個井距以上。夾層類型主要有3類,分別是泥質夾層、鈣質夾層和物性夾層,以物性夾層最常見,占57.75%,鈣質夾層次之,占25.35%,泥質夾層較少,僅占16.9%。考慮到夾層精細描述的需要,垂向網格長度應小于單河道砂體間的最小夾層厚度,平面網格長度應小于30 m。
大慶油田幾十年注水開發的實踐證實,以沖刷面互相疊置的厚層河道砂巖基本上都是多層多段水淹,垂向上的竄流并不十分嚴重。為避免儲層描述層次的混亂以及動態分析造成的麻煩,平面上流動單元的劃分應當限定在單一河道砂體內部及難以詳細追溯單砂體的復合型河道砂體內部[8]。因此,以單河道砂體為劃分單元,選取流動帶指數(FZI)、孔隙度、滲透率這3種參數為指標,通過聚類分析,劃分了3類流動單元,利用流動單元約束,建立了精細地質模型。Ⅰ類流動單元物性最好,主要分布于寬度大于800 m的河道中心區域;Ⅱ類流動單元物性次之,主要分布于寬度大于800 m的河道側緣及寬度小于600 m的河道中心區域,所占儲量比例最大;III類流動單元物性最差,主要分布于寬度小于600 m的河道邊部,所占儲量比例最小(表1)。

表1 流動單元劃分結果
所測相對滲透率曲線的多個樣品點,按照流動單元劃分的標準進行歸類。對歸屬同一流動單元的多個樣品點,進行歸一化處理,得到反映不同流動單元滲流特征的相對滲透率曲線(圖2)。從圖2中可以看出,各類流動單元的相對滲透率曲線特征區別較大。I類流動單元的兩相共存區最寬,油相相對滲透率下降最慢,水相相對滲透率上升最慢;II類流動單元次之,III類流動單元兩相共存區最窄,油相相對滲透率下降最快,水相相對滲透率上升最快。

圖2 各類流動單元相對滲透率曲線
研究表明[11-12],油藏進入高含水期后,由于黏土礦物的水化、膨脹、分散、遷移及其他地層微粒運移導致儲層物性發生變化。總的來看,孔隙度的變化幅度不大,滲透率向增大方面變化,巖石潤濕性向強親水方向變化。因此,整個開發階段的歷史,應該分低含水期、中含水期、高含水期3個含水階段進行擬合。一般來講,低含水期水驅孔隙體積倍數較小,物性變化不大。如果地質建模精度高,則在油藏數值模擬時,模型中地質儲量應與實際地質儲量吻合,且無需做參數調整,低含水期的各指標即可達到擬合標準。儲層物性與相對滲透率曲線調整主要針對中高含水期。以上一階段的網格壓力和含油飽和度為初始值,通過調整儲層參數和油水相對滲透率曲線,完成該階段的擬合。
關于儲層參數調整,油藏數值模擬工作者一般結合自己的經驗以及對油藏的認識來進行調整,目前尚無確定的標準,因此難度較大。筆者以實際開發數據為基礎,基于系統分析法[13]定量判別井間連通性,進而確定模型中各網格的滲透率變化因子,量化參數調整值,達到歷史擬合標準。


建立如下線性規劃:

如果該注水井有吸水剖面資料,設各單河道砂體的相對吸水量為aj,則可在上述線性規劃中增加約束條件:

利用修正單純形法[14]進行計算,即可得到上述線性規劃的最優解λij,λijkij即為調整后的滲透率。對滲透率發生變化的網格所屬的流動單元重新歸類,并適當調整相應流動單元的油水相對滲透率曲線,即可達到擬合標準。
利用流動單元約束建立精細地質模型,采用上述方法分階段歷史擬合(圖3)。從圖3中可以看出累計產油量、累計產水量、日產油、綜合含水等歷史擬合曲線吻合度高,單井擬合效果好的井數比例達90.5%(表2),碳氧比解釋結果與產液剖面測試結果對比結果誤差小(表3、4)。

圖3 各指標擬合曲線

表2 單井擬合結果

表3 碳氧比解釋結果與數值模擬結果對比

表4 產液剖面測試結果與數值模擬結果對比
(1)地質模型的精細程度應達到單河道砂體級別,垂向網格長度應小于最小夾層厚度,平面網格長度應小于30 m。
(2)在單河道砂體劃分結果的基礎上劃分不同流動單元,不同流動單元采用相應的油水相對滲透率曲線,能夠表征儲層內部油水滲流特征的差異。
(3)通過建立線性規劃,確定儲層參數調整標準,分階段歷史擬合能夠解決注水開發過程中儲層物性變化的影響。
(4)利用該方法進行油藏數值模擬,得到的結果與實際生產結果吻合度高,可用于其他類似油藏的數值模擬研究。
[1]姚光慶,馬正,趙彥超,等.淺水三角洲分流河道砂體儲層特征[J].石油學報,1995,16(1):24-30.
[2]姜在興.沉積學[M].北京:石油工業出版社,2003:384-385.
[3]陸道林,狄幫讓,李向陽,等.松遼盆地南部乾安北地區扶余油層沉積微相分析[J].中國石油大學學報:自然科學版,2011,35(6):8 -13.
[4]房寶財,張玉廣,許洪東,等.窄薄砂體油藏開發調整技術[M].北京:石油工業出版社,2004:1.
[5]魏兆亮,黃尚軍.高含水期油藏數值模擬技術和方法[J].西南石油大學學報:自然科學版,2008,30(1):103-105.
[6]趙翰卿,付志國,呂曉光,等.大型河流-三角洲沉積儲層精細描述方法[J].石油學報,2000,21(4):109-113.
[7]何宇航,于開春.分流平原相復合砂體單一河道識別及效果分析[J].大慶石油地質與開發,2005,24(2):17-19.
[8]趙翰卿.對儲層流動單元研究的認識與建議[J].大慶石油地質與開發,2001,20(3):8-10.
[9]周賢文,湯達幀,張春書.精細油藏數值模擬研究現狀及發展趨勢[J].特種油氣藏,2008,15(4):1-6.
[10]Miall A D.Architecture element analysis:a new method of facies analysis applied to fluvial deposit[J].Earth –Science Review,1985,22(4):261 -308.
[11]鄧玉珍,吳素英,張廣振,等.注水開發過程中儲層物理特征變化規律研究[J].油氣采收率技術,1996,3(4):44-52.
[12]王洪光,蔣明,張繼春,等.高含水期油藏儲集層物性變化特征模擬研究[J].石油學報,2004,2005(6):53-58.
[13]趙輝,李陽,高達,等.基于系統分析方法的油藏井間動態連通性研究[J].石油學報,2010,31(4):633-636.
[14]唐煥文,秦學志.實用最優化方法[M].大連:大連理工大學出版社,2004:41-55.