——以東營凹陷為例"/>
侯慶杰,劉顯太,韓宏偉,劉浩杰,魏國華,陳雨茂,于文政,王奇韻
1.中國石油化工股份有限公司勝利油田分公司物探研究院,山東東營 257022
2.勝利石油管理局博士后科研工作站,山東東營 257000
3.中國石油化工股份有限公司勝利油田分公司,山東東營 257000
致密油是繼頁巖氣之后全球非常規(guī)油氣勘探開發(fā)的新熱點[1-3],被石油工業(yè)界譽為“黑金”[4]。美國Bakken 致密油已實現(xiàn)規(guī)?;_發(fā),扭轉(zhuǎn)了北美石油產(chǎn)量遞減趨勢,改變了全球能源格局,成為全球重要油氣接替能源[5]。中國致密油勘探近兩年在渤海灣、鄂爾多斯、準噶爾和塔里木等盆地相繼取得了一定突破,以勝利油田為例,致密油控制儲量高達2.7×108t,是常規(guī)油氣藏接替能源[6],其中,致密灘壩砂儲量為1.7×108t,占致密油藏儲量的63%,為最主要的致密油儲層類型。灘壩砂體是發(fā)育于陸相湖盆濱淺湖地區(qū)非常重要的一種沉積砂體,多分布在湖泊邊緣、湖灣、湖中局部隆起的緩坡一側(cè),是波浪和湖流將各處帶來的砂質(zhì)組分再改造后在有利場所堆積而成[7]。沉積物粒度細,一般為細砂巖、粉砂巖、泥質(zhì)粉砂巖等[8],總體具有薄互層特征。由于非均質(zhì)性強,儲量動用難度較大,動用率較低,如何實現(xiàn)在致密灘壩砂儲層中“差中選優(yōu)”是提高致密油儲量動用率的關(guān)鍵問題之一。從其他地區(qū)致密油藏開發(fā)經(jīng)驗來看[9-10],孔隙結(jié)構(gòu)不僅決定了儲集層儲集和滲流能力,也影響著最終可動油氣總量[11]。
中高滲透儲層通常物性較好,開發(fā)難度低,儲層評價參數(shù)通常以孔隙度、滲透率及豐度為主,評價參數(shù)簡單。與常規(guī)砂巖儲層相比,致密砂巖儲層具有孔喉尺度小、非均質(zhì)性強和孔隙結(jié)構(gòu)復(fù)雜的特點[12-14],評價時需要加入孔喉半徑、分選系數(shù)和退汞效率等孔隙結(jié)構(gòu)參數(shù)。然而,根據(jù)致密儲層的特點,僅僅使用有限的樣品實驗參數(shù)表示整個評價層段和區(qū)域,忽略致密儲層非均質(zhì)性,使評價結(jié)果與實際情況有很大的偏差;如果使用測井資料預(yù)測各個評價參數(shù),參數(shù)種類繁多導(dǎo)致工作量巨大,均無法綜合全面地評價儲層優(yōu)劣。孔隙分形概念的出現(xiàn)在一定程度上解決了這些困難。前人研究表明,砂巖的微觀孔隙結(jié)構(gòu)具有很好的統(tǒng)計自相似性[15-16]。與經(jīng)典的歐式幾何相比,利用分形幾何能夠?qū)紫督Y(jié)構(gòu)的非均質(zhì)性和復(fù)雜性進行更真實地表征,而分形維數(shù)即為定量表征孔隙分形程度的參數(shù)??紫斗中闻c儲層物性、孔隙結(jié)構(gòu)的關(guān)系已有相關(guān)討論,但是孔隙分形如何影響開發(fā)效果,如何利用孔隙分形評價儲層優(yōu)劣還少見討論。此外,目前對于孔隙分形維數(shù)的計算多停留在單點實驗測試分析的層面上,無法形成單井或者平面上連續(xù)性的預(yù)測。因此,針對以上問題,本文以東營凹陷沙四上純下次亞段的致密灘壩砂儲集層為例,利用薄片、物性及壓汞等相關(guān)測試數(shù)據(jù),計算致密灘壩砂孔喉分形維數(shù),并探討分形維數(shù)與儲層物性、孔隙結(jié)構(gòu)參數(shù)之間的關(guān)系,然后優(yōu)選測井數(shù)據(jù),建立了一種致密灘壩砂儲層孔隙分形維數(shù)的測井預(yù)測模型,對東營凹陷西部地區(qū)致密灘壩砂儲集層分形維數(shù)的平面分布進行了預(yù)測?;诜中尉S數(shù)儲層評價標準,對典型區(qū)的致密灘壩砂儲層進行了分類與評價,優(yōu)選出有利儲層分布區(qū)域,為確定開發(fā)區(qū)域動用次序提供依據(jù)。
研究區(qū)東營凹陷為渤海灣盆地中受陳南斷裂控制的一個二級構(gòu)造單元,是總體呈“北斷南超”的中新生代箕狀斷陷盆地,南接廣饒凸起,北部與陳家莊、林樊家和濱縣凸起相接,東部為青坨子凸起,西部跨過青城凸起與惠民洼陷相接,面積約為5 850 km2(圖1)[17-20]。古近系沙四段,主要分為沙四上亞段和沙四下亞段,其中,沙四上亞段又可進一步分為純上和純下兩個次亞段。根據(jù)前人研究成果,致密灘壩砂儲層主要分布于純下次亞段,其中,壩砂單層厚度相對較大,一般大于2 m,灘砂單層厚度通常小于2 m,壩砂整體物性優(yōu)于灘砂,分布區(qū)域如圖1所示。
圖1 東營凹陷沙四上純下次亞段灘壩沉積體系分布圖(據(jù)唐東[17])Fig.1 Map of EsU4 beach-bar deposition system (after Tang[17])
本次研究選取了12塊東營凹陷沙四上純下次亞段致密灘壩砂儲集層樣品進行巖石鑄體薄片觀察、物性和高壓壓汞測試,并基于壓汞實驗測試數(shù)據(jù)進行孔喉分形計算。根據(jù)薄片統(tǒng)計顯示(表1),灘砂與壩砂的礦物碎屑組成近似,無明顯差別。12 個樣品巖性主要為巖屑質(zhì)長石砂巖。石英含量為37%~51%,平均值為46.25%;長石含量為30%~38%,平均值為34%;巖屑類型主要為變質(zhì)巖、沉積巖和巖漿巖為主,含量為16%~25%,平均值為19.67%。填隙物以泥質(zhì)雜基為主,含量為2%~8%,平均值為5.5%。膠結(jié)物含量為4%~21%(平均值為9.42%),主要為方解石和白云石(平均含量分別為5.08%和3.42%),其次為硅質(zhì)膠結(jié)(平均含量為0.92%)。砂巖中碎屑顆粒粒徑為0.06~0.32 mm,平均值為0.15 mm,主要為粉砂和細砂,成分成熟度中等,分選中等,碎屑磨圓度主要為次棱角狀。
表1 致密灘壩砂樣品巖石學(xué)特征Table 1 Petrological properties of tight beach-bar sand samples
灘壩砂儲集層孔隙度和滲透率的統(tǒng)計結(jié)果表明(圖2):壩砂儲集層孔隙度為1.0%~30.1%(主體分布在5%~18%,平均值為11.8%);滲透率主體分布在0.1×10-3μm2~10×10-3μm2(平均值為0.92×10-3μm2);灘砂儲集層孔隙度為0.8%~26%(主體分布在3%~16%,平均值為8.9%),滲透率主體分布在0.01×10-3μm2~1×10-3μm2(平均值為0.67×10-3μm2),其孔隙度和滲透率均小于壩砂儲集層。從整體上看,73%的樣品孔隙度小于15%,滲透率小于10×10-3μm2,孔隙度平均值小于12%,滲透率平均值小于1×10-3μm2,為典型的致密砂巖儲集層。
圖2 灘壩砂儲集層孔隙度和滲透率分布直方圖Fig.2 Porosity and permeability histograms for beach-bar sand reservoir
12個樣品的壓汞曲線特征如圖3所示,其中,樣品編號1、2、4、9、10和11為壩砂樣品,樣品編號3、5、6、7、8和12為灘砂樣品。壓汞參數(shù)如表2所示,壩砂樣品整體體現(xiàn)為粗歪度、分選較好的特點,排驅(qū)壓力0.131~0.881 MPa(平均值為0.375 MPa),最大孔喉半徑0.834~5.611μm(平均值為2.954μm),飽和度中值壓力0.83~4.59 MPa(平均值為1.87 MPa),飽和度中值喉道半徑0.16~0.886μm(平均值為0.541μm);灘砂樣品整體體現(xiàn)為細歪度、分選較差的特點,排驅(qū)壓力0.334~3.562 MPa(平均值為1.753 MPa),最大孔喉半徑0.206~2.201μm(平均值為0.752μm),飽和度中值壓力1.95~21.09 MPa(平均值為11.82 MPa),飽和度中值喉道半徑0.035~0.377 μm(平均值為0.117 μm)。灘砂樣品孔喉半徑整體小于壩砂樣品。
表2 樣品壓汞參數(shù)表Table 2 Mercury intrusion parameters for 12 samples
圖3 研究樣品壓汞曲線特征Fig.3 Mercury intrusion curves for 12 samples
根據(jù)壓汞曲線可以獲得孔喉半徑頻率分布圖(圖4),如圖所示,孔喉半徑多呈單峰分布,主要分布區(qū)間為0.1~1μm,但壩砂與灘砂的孔喉峰值半徑(分布頻率最高對應(yīng)的孔喉半徑)有一定的區(qū)別,其中,壩砂孔喉峰值半徑為0.183~1.529 μm,平均值為0.8μm;灘砂孔喉峰值半徑明顯小于壩砂,為0.059~0.189 μm,平均值為0.119 μm。隨孔喉峰值半徑增加,孔隙度相應(yīng)增大,當孔喉峰值半徑小于0.4 μm時,孔隙度均小于10%。
圖4 不同孔隙度灘壩砂樣品的孔喉半徑分布圖Fig.4 Pore throat radius graphs for beach-bar sand samples with different porosity
分形通常被定義為“一個粗糙或零碎的幾何形狀,可以分成數(shù)個部分,且每一部分都(或近似)是整體縮小后的形狀”,即具有“自相似或標度不變性”[21]。前人的研究表明,致密砂巖的孔隙結(jié)構(gòu)存在分形特征,分形理論可用于表征致密砂巖儲集層的結(jié)構(gòu)特征及其非均質(zhì)性[22-24];多孔巖石的分形維數(shù)通常介于2~3[25-27],分形維數(shù)越小,儲層孔喉分布越均勻,均質(zhì)性越強[28-29]。
根據(jù)前人的研究成果[30],儲層孔隙分布表示式為:
根據(jù)油層物理學(xué)中對于飽和度的定義,孔徑小于r的累積體積分數(shù)s即為潤濕相飽和度[31],且對于致密儲層,孔喉半徑的非均質(zhì)性強,因此,rmin<<rmax,則有下式:
式中:Sw為潤濕相飽和度,%;rmax為最大孔喉半徑,μm。
對式6兩邊取對數(shù)可得:
式中:SHg為累積進汞飽和度,%。
根據(jù)12 塊樣品的lg(1-SHg)-lgr關(guān)系圖(圖5),可以計算出樣品的分形維數(shù)。從圖中可以看出,12 塊致密灘壩砂樣品的lg(1-SHg)-lgr關(guān)系圖具有明顯兩段性,將lg(1-SHg)-lgr關(guān)系圖中分段點對應(yīng)的孔喉半徑命名為孔喉轉(zhuǎn)折半徑,大于孔喉轉(zhuǎn)折半徑對應(yīng)的孔隙為大孔喉,用?1表示,小于孔喉轉(zhuǎn)折半徑對應(yīng)的孔隙為小孔喉,用?2表示。根據(jù)各孔隙空間的孔隙度加權(quán)平均得到整個孔喉的總分形維數(shù)D[27,32],為
圖5 12 塊樣品的lg(1-SHg)-lgr 關(guān)系圖Fig.5 Diagrams of log(1-SHg) vs. logr for 12 samples
根據(jù)計算結(jié)果,兩段分形維數(shù)D1、D2以及總分形維數(shù)D均在2~3范圍內(nèi)(表3),對壩砂和灘砂的分形維數(shù)分別統(tǒng)計,壩砂分形維數(shù)為2.222 4~2.531 9,平均值為2.441 2,而灘砂分形維數(shù)為2.585 0~2.742 4,平均值為2.650 9,明顯高于壩砂。
表3 物性參數(shù)與分形維數(shù)Table 3 Physical properties and fractal dimensions
2.2.1 分形維數(shù)與儲層物性關(guān)系
分別建立分形維數(shù)與儲層孔隙度和滲透率的關(guān)系。根據(jù)樣品物性統(tǒng)計數(shù)據(jù)(表3),總分形維數(shù)D與儲層孔隙度?和滲透率k之間呈負向關(guān)關(guān)系(圖6),R2分別為0.819和0.902,相關(guān)程度高。說明分形維數(shù)可以有效表征儲層孔滲性能,即隨著分形維數(shù)增大,儲層的孔隙度與滲透率均變小,儲層的孔滲性變差。
圖6 分形維數(shù)D-孔隙度? 和分形維數(shù)D-滲透率k 交會圖Fig.6 Plots of (a) fractal dimension D vs. porosity ? and (b) D vs. permeability k
2.2.2 分形維數(shù)與孔隙結(jié)構(gòu)參數(shù)關(guān)系
孔隙結(jié)構(gòu)包含孔喉大小、分選性和連通性等信息,排驅(qū)壓力PT、最大孔喉半徑rmax、飽和度中值壓力Pc50及孔喉中值半徑r50表示孔喉大小,分選系數(shù)Sp表示孔喉分選性,退汞效率WE可表示孔喉連通性。根據(jù)分形維數(shù)與孔隙結(jié)構(gòu)參數(shù)數(shù)據(jù)統(tǒng)計發(fā)現(xiàn)(圖7),分形維數(shù)D與PT、rmax、Pc50、r50、Pc50、Sp和WE均具有較強相關(guān)性,其中,D與PT和Pc50呈正相關(guān)關(guān)系,與rmax和r50呈負相關(guān)關(guān)系,R2分別為0.566、0.519、0.823 和0.785,表明隨著分形維數(shù)升高,孔喉半徑有減小的趨勢;D與Sp呈正相關(guān)關(guān)系,R2為0.732,表明隨著分形維數(shù)升高,孔喉的分選性變差;D與WE呈負相關(guān)關(guān)系,R2為0.656,表明隨著分形維數(shù)升高,孔喉的連通性變差。因此,分形維數(shù)可以作為綜合反映儲層孔隙結(jié)構(gòu)特征的指標,分形維數(shù)越大,儲層的儲集條件越差。
圖7 分形維數(shù)與孔隙結(jié)構(gòu)參數(shù)交會圖Fig.7 Intersection graphs of fractal dimension D and pore structure parameters
致密灘壩砂儲層非均質(zhì)性強,連續(xù)系統(tǒng)的鉆井取心和實驗分析費用高昂且難以操作,以有限的樣品取其平均值代表整套砂體的做法不能反映細小差別,得到的結(jié)論很可能誤導(dǎo)開發(fā)決策。測井信息可以間接地反映地層巖性、物性及其流體等特征,部分測井曲線與孔隙結(jié)構(gòu)具有很強的相關(guān)性,因此,可嘗試建立一種準確定量的、易于操作的孔隙分形維數(shù)測井預(yù)測模型,預(yù)測孔隙分形維數(shù)的單井縱向分布,解決實驗資料縱向連續(xù)性差的問題,提高預(yù)測精度。
2.3.1 測井參數(shù)優(yōu)選
通過對致密灘壩砂樣品的分形維數(shù)與常用測井曲線進行相關(guān)性分析(表4),發(fā)現(xiàn)微電位與微梯度差值的絕對值(|RLML-RNML|)、聲波時差(AC)、自然伽馬(GR)和電阻率(RT)相關(guān)程度較高,相關(guān)系數(shù)分別為-0.907、-0.896、0.877和-0.867,隨著分形維數(shù)D增大,|RLML-RNML|值減小、AC 值減小、GR 值增大、RT值減小。
表4 致密灘壩砂樣品分形維數(shù)與測井類型相關(guān)系數(shù)Table 4 Correlation coefficient between fractal dimension and logging type for samples of dense beach-bar sand
(1)微電位與微梯度差值絕對值(|RLMLRNML|)和聲波時差(AC)主要反映物性變化,D值增大,孔隙度和滲透率均呈減小的趨勢(圖6),導(dǎo)致|RLML-RNML|和AC值減小。
(2)自然伽馬(GR)主要反映巖石的放射性。D為大值的樣品多為灘砂樣品,粒度較細,且泥質(zhì)含量較D為小值的壩砂樣品大,導(dǎo)致放射性增強,GR 值增大。
(3)電阻率(RT)主要反映含油性,隨著D升高,孔喉半徑減小,分選性變差(圖7),油氣運移至儲層困難加大,導(dǎo)致含油量減小,RT值降低。
2.3.2 測井預(yù)測模型的建立與檢驗
建立東營凹陷致密灘壩砂測井預(yù)測模型的具體方法為:首先根據(jù)灘壩砂沉積相特點,剔除非灘壩砂層段;為了抵消各類測井參數(shù)數(shù)量級不同帶來的影響,將優(yōu)選出的測井曲線(|RLML-RNML|、AC、GR 和RT),進行標準歸一化處理,標準歸一化公式為:
式中:測井樣本數(shù)據(jù)為x,xmax=max{x},xmin=min{x},測井樣本歸一化數(shù)據(jù)為X;然后,將歸一化后的|RLMLRNML|、AC、GR和RT值對D進行回歸擬合,得到D的測井預(yù)測公式:
式中:D為分形維數(shù);|RLML-RNML|為微電極曲線差值絕對值,Ω·m;AC 為聲波時差值,μs/m;GR為自然伽馬值,API;RT為電阻率值,Ω·m。
將通過測井預(yù)測模型預(yù)測的分形維數(shù)與通過壓汞數(shù)據(jù)計算得到的分形維數(shù)進行比對(圖8),R2值達0.936,基本滿足實際生產(chǎn)和科研的需要。
圖8 分形維數(shù)D 測井預(yù)測值與壓汞數(shù)據(jù)計算值相關(guān)圖Fig.8 Correlation of fractal dimension D log predicted value and value calculated from mercury intrusion data
分形維數(shù)D可以有效表征儲層孔滲性能,綜合反映儲層孔隙結(jié)構(gòu)特征,因此,推測分形維數(shù)D值大小與油井產(chǎn)能具有一定的相關(guān)性。為了揭示分形維數(shù)D對儲層開發(fā)效果的影響,本次研究選取已動用區(qū)塊無壓裂措施下的日產(chǎn)油,累積產(chǎn)油與分形維數(shù)D之間的關(guān)系(圖9),建立了基于分形維數(shù)D的東營凹陷致密灘壩砂儲層評價標準(表5)。
圖9 分形維數(shù)D 值與日產(chǎn)油、累積產(chǎn)油關(guān)系Fig.9 Relationship between fractal dimension D and (a) daily oil production and (b) cumulative oil production
根據(jù)日產(chǎn)油、累積產(chǎn)油與分形維數(shù)D的數(shù)據(jù)統(tǒng)計結(jié)果顯示,以日產(chǎn)油量為主要分類指標,當2<D<2.35 時,可以認定為相對優(yōu)質(zhì)的儲層,此時,對應(yīng)層段的日產(chǎn)油量大于10 t/d,累積產(chǎn)油量大于90 t,孔隙度大于14%,滲透率大于0.9×10-3μm2,PT小于0.18 MPa,rmax大于4.1 μm,Pc50小于1.05 MPa,r50大于0.7μm,Sp小于0.19,WE大于33%;當2.35<D<2.55 時,認定為相對一般的儲層,日產(chǎn)油量為5~10 t/d,累積產(chǎn)油量為50~90 t,孔隙度為10%~14%,滲透率為0.4×10-3μm2~0.9×10-3μm2,PT介于0.18~0.37 MPa,rmax為2.0~4.1μm,Pc50為1.05~2.3 MPa,r50為0.32~0.7μm,Sp為0.19~0.39,WE為26%~33%;當2.55<D<3 時,認定為相對較差的儲層,日產(chǎn)油量小于5 t,累積產(chǎn)油小于50 t,孔隙度小于10%,滲透率小于0.4×10-3μm2,PT大于0.37 MPa,rmax小于2.0μm,Pc50大于2.3 MPa,r50小于0.32μm,Sp大于0.39,WE小于26%。
本次研究以東營凹陷西部區(qū)塊沙四上純下次亞段的致密灘壩砂分布區(qū)域為例,利用分形維數(shù)D測井預(yù)測結(jié)果對致密灘壩砂體進行了儲層評價。根據(jù)分形維數(shù)D平面分布等值線圖(圖10),結(jié)合東營凹陷致密灘壩砂儲層綜合評價標準(表5),優(yōu)質(zhì)儲層的分形維數(shù)D為2~2.35,孔隙度大于14%,主要集中在兩個區(qū)域,分別為B424井區(qū)和C107-C66-B3井區(qū),全部分布于壩砂,該區(qū)域單井日產(chǎn)能高,日產(chǎn)油量10 t以上,為最有利的開發(fā)區(qū)域;一般儲層的分形維數(shù)D為2.35~2.55,孔隙度為10%~14%,主要集中在四個區(qū)域,分別為B424 井區(qū)、F119 井區(qū)、C108-L108-C44井區(qū)和F145-F144-F134井區(qū),主要為壩砂,少量為灘砂,該區(qū)域的單井日產(chǎn)油量為5~10 t,為較好的開發(fā)區(qū)域。利用分形維數(shù)D對致密灘壩砂儲層進行綜合評價,為下一步有利開發(fā)區(qū)域的選取及動用次序的確定提供了理論依據(jù)。
圖10 東營凹陷西部區(qū)塊分形維數(shù)D 值平面分布及儲層分類Fig.10 Plane distribution of fractal dimension D and reservoir classification in western block, Dongying Sag
表5 東營凹陷致密灘壩砂儲層綜合評價標準Table 5 Comprehensive evaluation criterion for tight beach-bar sand reservoirs in Dongying Sag
表5 所列的儲層評價標準適用于東營凹陷灘壩砂型致密儲層,現(xiàn)場應(yīng)用效果顯著。本文所述基于分形維數(shù)的儲層評價方法也在東營凹陷的砂礫巖和濁積巖型致密儲層中進行了推廣和應(yīng)用,但是儲層的優(yōu)劣對應(yīng)的分形維數(shù)范圍與灘壩砂有一定的差異。因此,基于分形維數(shù)的儲層評價方法可以應(yīng)用于致密儲層評價,但是標準的制定還需要根據(jù)致密儲層的實際地質(zhì)特點作相應(yīng)的調(diào)整。
(1)東營凹陷沙四上純下次亞段分布典型的致密灘壩砂體沉積,利用壓汞曲線數(shù)據(jù)計算孔隙分形維數(shù),分析其與表征儲層孔滲性能與孔隙結(jié)構(gòu)參數(shù)的相關(guān)性,發(fā)現(xiàn)分形維數(shù)與總孔隙度?、滲透率k、最大孔喉半徑rmax、中值喉道半徑r50、退汞效率WE呈負相關(guān)關(guān)系,與排驅(qū)壓力PT、中值壓力Pc50、分選系數(shù)Sp呈正相關(guān)關(guān)系,表明隨著分形維數(shù)增大,儲層孔滲性能、孔喉半徑、分選性及連通性均變差。
(2)優(yōu)選微電位與微梯度差值的絕對值(|RLML-RNML|)、聲波時差(AC)、自然伽馬(GR)和電阻率(RT)測井曲線參數(shù),標準歸一化擬合建立分形維數(shù)測井預(yù)測模型,將通過測井預(yù)測模型預(yù)測的分形維數(shù)與通過壓汞數(shù)據(jù)計算得到的分形維數(shù)進行比對,R2值達0.936,基本滿足實際生產(chǎn)和科研的需要。
(3)選取已動用區(qū)塊無壓裂措施下的日產(chǎn)油,累積產(chǎn)油與分形維數(shù)D之間的關(guān)系,建立基于分形維數(shù)D的東營凹陷致密灘壩砂儲層評價標準,當2<D<2.35 時為優(yōu)質(zhì)儲層,當2.35<D<2.55 時為一般儲層,當2.55<D<3時為差儲層。根據(jù)測井預(yù)測模型預(yù)測東營凹陷西部區(qū)塊沙四上純下次亞段的致密灘壩砂有利儲層,認為優(yōu)質(zhì)儲層主要集中在B424 井區(qū)和C107-C66-B3井區(qū),為下一步儲層動用次序的確定提供了理論依據(jù)。