鹽池地區(qū)延長(zhǎng)組長(zhǎng)8段致密油儲(chǔ)層脆性指數(shù)測(cè)井評(píng)價(jià)

李坪東,蘇幽雅,邵曉巖,徐 寧,張 雁,楊國(guó)斌

(中國(guó)石油長(zhǎng)慶油田分公司第三采油廠(chǎng),寧夏銀川 750006)

隨著全球能源行業(yè)的日益發(fā)展,能源需求逐漸增加,非常規(guī)油氣資源逐漸成為新的接替能源。鄂爾多斯盆地三疊系延長(zhǎng)組發(fā)育陸相致密油儲(chǔ)層,該類(lèi)儲(chǔ)層非均質(zhì)性強(qiáng),滲透率通常小于0.1×10-3μm2,孔隙度通常小于10%[1-2],在自然條件下無(wú)法獲得產(chǎn)能或工業(yè)產(chǎn)能,只能通過(guò)水力壓裂實(shí)現(xiàn)工業(yè)開(kāi)發(fā)[3]。

脆性指數(shù)是地層可壓裂性評(píng)價(jià)的重要參數(shù),對(duì)儲(chǔ)層改造起到關(guān)鍵作用[4-5]。本文以鄂爾多斯盆地鹽池地區(qū)三疊系延長(zhǎng)組長(zhǎng)8段為實(shí)例,利用巖石礦物成分測(cè)試、陣列聲波測(cè)井、常規(guī)測(cè)井等資料,建立基于測(cè)井資料的致密砂巖儲(chǔ)層脆性指數(shù)計(jì)算模型,該方法可為該區(qū)及地質(zhì)特征相似區(qū)域的致密砂巖油藏勘探開(kāi)發(fā)提供參考和借鑒。

1 地質(zhì)背景

1.1 勘探開(kāi)發(fā)現(xiàn)狀

鹽池地區(qū)位于寧夏回族自治區(qū)東部,北起古峰莊、南至馮地坑、西抵大水坑、東達(dá)馬坊。構(gòu)造位置處于鄂爾多斯盆地西緣的天環(huán)坳陷,構(gòu)造相對(duì)簡(jiǎn)單,每千米起伏約為10 m。該地區(qū)為一多套層系疊合的油氣富集區(qū)域,其中,中生界三疊系延長(zhǎng)組致密砂巖油藏是近年來(lái)非常規(guī)資源產(chǎn)能建設(shè)的重點(diǎn)目標(biāo)[6]。但目前針對(duì)該區(qū)延長(zhǎng)組致密砂巖油藏研究主要集中在儲(chǔ)層沉積演化規(guī)律、砂體構(gòu)型、孔隙結(jié)構(gòu)等方面。關(guān)于儲(chǔ)層脆性、可壓裂性評(píng)價(jià)等工作的開(kāi)展相對(duì)滯后,一定程度上導(dǎo)致了產(chǎn)能建設(shè)“甜點(diǎn)區(qū)”優(yōu)選存在盲目性,在產(chǎn)能建設(shè)、開(kāi)發(fā)管理的過(guò)程中,存在試采液量低、單井產(chǎn)量低等問(wèn)題。本文圍繞該地區(qū)延長(zhǎng)組長(zhǎng)8段致密砂巖儲(chǔ)層,開(kāi)展脆性指數(shù)的測(cè)井預(yù)測(cè)方法研究,旨在加快該地區(qū)勘探開(kāi)發(fā)一體化整體推進(jìn)。

1.2 基礎(chǔ)地質(zhì)特征

鄂爾多斯盆地鹽池地區(qū)三疊系延長(zhǎng)組自下而上可劃分為長(zhǎng)10至長(zhǎng)1共十段,長(zhǎng)8段的巖性主要為灰色長(zhǎng)石巖屑砂巖、巖屑長(zhǎng)石砂巖及黑色泥巖。陸相砂泥巖頻繁互層,多期砂體縱向疊置,沉積微相包括分流河道、分流間灣。目的層砂巖粒度細(xì),分選程度中等-差,膠結(jié)類(lèi)型多為孔隙式膠結(jié)。長(zhǎng)8段屬于三角洲平原沉積環(huán)境,分流河道沿北西至南東方向呈“條帶狀”展布,儲(chǔ)層平均厚度18.3 m,孔隙類(lèi)型主要包括原生粒間孔、粒間溶孔,儲(chǔ)層平均孔隙度為7.1%、平均滲透率為0.212×10-3μm2。長(zhǎng)7段暗色泥巖厚度大、分布范圍廣,含大量有機(jī)質(zhì),生烴潛量大,為優(yōu)質(zhì)烴源巖。烴類(lèi)在剩余壓力及擴(kuò)散作用下,以斷裂系統(tǒng)為運(yùn)移通道,油氣自上而下運(yùn)移至下伏長(zhǎng)8段致密砂巖儲(chǔ)層中聚集[7]。

2 脆性指數(shù)定量計(jì)算

巖石脆性是指巖石發(fā)生破裂前的瞬間變化難易程度[8]。巖石脆性的評(píng)價(jià)指標(biāo)為脆性指數(shù),目前存在多種脆性指數(shù)計(jì)算方法,如基于礦物組分的評(píng)價(jià)方法、基于巖石力學(xué)性質(zhì)的評(píng)價(jià)方法等。

關(guān)于基于礦物組分的評(píng)價(jià)方法,主要是利用脆性礦物含量計(jì)算脆性指數(shù)。但是,針對(duì)不同類(lèi)型的儲(chǔ)層,對(duì)脆性礦物的界定往往存在差異。對(duì)于頁(yè)巖,將石英、長(zhǎng)石、碳酸鹽視為脆性礦物[5,9];而對(duì)于致密砂巖,多數(shù)情況下將石英、碳酸鹽視為脆性礦物[8-9]。此外,對(duì)于不同地區(qū)、不同層位的致密砂巖,由于成巖作用等條件的差異,脆性礦物的界定也并不完全一致[8-9]。相反,基于巖石力學(xué)性質(zhì)的脆性指數(shù)評(píng)價(jià)方法,對(duì)于不同類(lèi)型的巖石,均采用同一套實(shí)驗(yàn)標(biāo)準(zhǔn)和工作流程。

利用巖石力學(xué)實(shí)驗(yàn)獲取相關(guān)參數(shù)計(jì)算脆性指數(shù)是目前最常用的巖石脆性定量評(píng)價(jià)方法[9-13],可直接測(cè)定致密儲(chǔ)層的力學(xué)參數(shù),如楊氏模量、泊松比等。楊氏模量是巖體剛性的度量,可表征巖石形變難易程度,與脆性指數(shù)呈正相關(guān)關(guān)系;而泊松比為巖體橫向與縱向應(yīng)變的比值,常被用來(lái)表征巖體的脆性或塑性特征,與脆性指數(shù)呈負(fù)相關(guān)關(guān)系[12-13]?；诹W(xué)參數(shù)的脆性評(píng)價(jià)方法通常最能代表復(fù)雜礦物組合致密油儲(chǔ)層的的脆性特征[10-13]。

2.1 計(jì)算方法

當(dāng)聲波透過(guò)巖石介質(zhì)時(shí),巖石的剪切及壓縮變形性質(zhì)可以通過(guò)物理方程進(jìn)行定量表征,巖石可以傳播剪切波(橫波),也可以傳播壓縮波(縱波)[14-16]。對(duì)于不同的巖石而言,聲波在傳播時(shí)的速度、頻率、幅度、衰減等聲學(xué)特性存在差異。

聲波監(jiān)測(cè)技術(shù)常被用來(lái)探測(cè)巖石動(dòng)態(tài)巖石力學(xué)屬性,通過(guò)陣列聲波測(cè)井、偶極聲波測(cè)井、巖石聲波特性室內(nèi)測(cè)試,可以獲得橫波時(shí)差;進(jìn)而,根據(jù)縱波時(shí)差、橫波時(shí)差、巖石密度,計(jì)算楊氏模量為:

(1)

式中:E為楊氏模量,GPa;ρb為巖石密度,g/cm3;Δts、Δtc分別為橫波時(shí)差及縱波時(shí)差,μs/m。

目的層巖石泊松比的表達(dá)方程為:

(2)

式中:μ為泊松比,小數(shù)。

在以往針對(duì)致密砂巖、頁(yè)巖儲(chǔ)層的研究中,利用歸一化之后的楊氏模量和泊松比來(lái)計(jì)算致密儲(chǔ)層的脆性指數(shù):

(3)

式中:Enorm為歸一化的楊氏模量,%;Emax、Emin為楊氏模量最大值、最小值,GPa。

(4)

式中:μnorm為歸一化的泊松比,%;μmax、μmin為泊松比最大值、最小值,小數(shù)。

(5)

式中:BI為巖石脆性指數(shù),%。

2.2 計(jì)算結(jié)果

鹽池地區(qū)內(nèi)Y162井長(zhǎng)8段巖石聲波特性測(cè)試數(shù)據(jù)如表1所示。Y162井位于鹽池地區(qū)中部,該井長(zhǎng)8砂巖層段進(jìn)行了鉆井取心,針對(duì)長(zhǎng)8段的17塊巖心樣品,進(jìn)行了巖石聲波特性室內(nèi)實(shí)驗(yàn)測(cè)試(表1)。根據(jù)測(cè)試結(jié)果,橫波時(shí)差分布范圍為339.21～395.71 μs/m,平均值為354.65 μs/m;縱波時(shí)差分布范圍為181.25～214.56 μs/m,平均值為202.15 μs/m;巖石密度主要分布在2.55～2.66 g/cm3,平均值為2.60 g/cm3。最終,計(jì)算得出目的層的楊氏模量分布范圍為43.9～56.2 GPa,平均值為52.1 GPa;泊松比的分布范圍為0.205～0.314,平均值為0.257;脆性指數(shù)分布范圍為36.5%～57.2%,且平均值為 50.1%。

表1 Y162井長(zhǎng)8段巖石聲波特性測(cè)試數(shù)據(jù)

F70井與Y162井相距約1.5 km,對(duì)該井進(jìn)行了陣列聲波測(cè)井現(xiàn)場(chǎng)施工,針對(duì)該井的長(zhǎng)8段,分別利用聲波時(shí)差及巖石密度測(cè)井資料及公式(1)、(2)對(duì)目的層段的楊氏模量和泊松比進(jìn)行了精細(xì)測(cè)井解釋,解釋結(jié)果如圖1所示。該井長(zhǎng)8段楊氏模量最大值為52.9 GPa,最小值為24.0 GPa,平均值為 38.8 GPa;對(duì)于巖石泊松比,最大值為0.340,最小值為0.187,平均值為0.269。最終,將目的層段的楊氏模量和泊松比的測(cè)井計(jì)算結(jié)果進(jìn)行了歸一化處理,并利用公式(5)計(jì)算了脆性指數(shù)并建立了相應(yīng)的測(cè)井曲線(xiàn)(圖1)。計(jì)算結(jié)果顯示,該井長(zhǎng)8段脆性指數(shù)最大值為64.1%,最小值為30.0%,平均值為46.5%。

圖1 F70井長(zhǎng)8段脆性指數(shù)測(cè)井解釋成果

3 脆性指數(shù)測(cè)井評(píng)價(jià)

巖石的脆性指數(shù)可以通過(guò)上述方法進(jìn)行計(jì)算,但問(wèn)題在于這些方法需要在鉆井取心、巖樣采集或開(kāi)展過(guò)陣列聲波測(cè)井、偶極聲波測(cè)井的基礎(chǔ)上應(yīng)用。因此,利用上述方法,只能針對(duì)計(jì)算巖心樣品或進(jìn)行過(guò)陣列聲波、偶極聲波測(cè)井的井位進(jìn)行計(jì)算,無(wú)法應(yīng)用于巖石脆性縱向連續(xù)評(píng)價(jià)和平面精細(xì)刻畫(huà)等地質(zhì)研究。

因此,充分利用常規(guī)測(cè)井資料,建立儲(chǔ)層脆性指數(shù)測(cè)井評(píng)價(jià)模型,對(duì)于致密砂巖油氣資源的精準(zhǔn)勘探和高效開(kāi)發(fā)顯得尤為重要,而通過(guò)上述試驗(yàn)方法計(jì)算的脆性指數(shù),則可以作為約束條件,用來(lái)標(biāo)定常規(guī)測(cè)井預(yù)測(cè)的脆性指數(shù),從而使得測(cè)井預(yù)測(cè)精度更為可靠。

3.1 構(gòu)建巖性指數(shù)預(yù)測(cè)脆性指數(shù)

陸相致密油儲(chǔ)層實(shí)質(zhì)上由多種復(fù)雜的微觀(guān)骨架礦物及基質(zhì)礦物組成,其中,對(duì)儲(chǔ)層脆性影響最直接的組分為脆性礦物,當(dāng)儲(chǔ)層中脆性礦物組分含量大于某個(gè)臨界值,則可以達(dá)到較好的壓裂改造效果,有利于提高初期單井產(chǎn)能。借鑒以往關(guān)于鄂爾多斯盆地三疊系致密砂巖儲(chǔ)層脆性評(píng)價(jià)的研究,將石英和碳酸鹽類(lèi)視為脆性礦物。

自然伽馬測(cè)井屬于放射性測(cè)井的一種方法,自然伽馬曲線(xiàn)常被用來(lái)對(duì)儲(chǔ)層巖性進(jìn)行判別。黏土巖類(lèi)含有較高的放射性元素,且黏土礦物對(duì)陽(yáng)離子的吸附能力較強(qiáng),可吸附周邊放射性物質(zhì)。因此,當(dāng)儲(chǔ)層中表現(xiàn)為高泥質(zhì)含量時(shí),自然伽馬值通常偏高;相反,當(dāng)儲(chǔ)層中含有大量石英、碳酸鹽類(lèi)等非放射性物質(zhì)時(shí),自然伽馬表現(xiàn)為低值特征。

當(dāng)?shù)湍芰抠ゑR射線(xiàn)與地層中的元素發(fā)生光電效應(yīng)時(shí),伽馬射線(xiàn)會(huì)被原子吸收并釋放出光電子,因此,光電吸收截面可以表征地層中元素對(duì)伽馬射線(xiàn)吸收能力的強(qiáng)弱。地層中不同的巖石礦物成分對(duì)應(yīng)了不同的化學(xué)元素,這些元素組成的差異表現(xiàn)為光電吸收截面指數(shù)的差異。

本文利用常規(guī)測(cè)井資料來(lái)定義巖性指數(shù),即自然伽馬與光電吸收截面的乘積。巖性指數(shù)表示單位光電吸收截面指數(shù)條件下地層自然伽馬測(cè)井定量響應(yīng)特征。

LI=GR·PE

(6)

式中:LI為巖性指數(shù),API·b/e;GR為自然伽馬測(cè)井值,API;PE為光電吸收截面指數(shù)測(cè)井值,b/e。

該式為數(shù)據(jù)驅(qū)動(dòng)模型,物理量的單位不參與運(yùn)算。

針對(duì)F70井,繪制巖性指數(shù)與脆性指數(shù)之間對(duì)應(yīng)關(guān)系散點(diǎn)圖(圖2),兩者之間存在明顯的負(fù)相關(guān)關(guān)系,隨著巖性指數(shù)的增加,脆性指數(shù)也呈現(xiàn)出降低的趨勢(shì)。建立兩者之間回歸公式,得出了利用巖性指數(shù)計(jì)算巖石脆性指數(shù)的數(shù)學(xué)表達(dá)式:

圖2 LI與BI對(duì)應(yīng)關(guān)系

BIa1=2 374.5LI-0.537 6

(7)

式中:BIa1為利用巖性指數(shù)計(jì)算的脆性指數(shù),%。

3.2 構(gòu)建物性指數(shù)預(yù)測(cè)脆性指數(shù)

一般情況下,縱波時(shí)差、密度均可通過(guò)常規(guī)測(cè)井資料獲得;受作業(yè)成本等因素的限制,僅少數(shù)井實(shí)施陣列聲波或偶極聲波測(cè)井,因此,橫波時(shí)差測(cè)井資料較為有限。

根據(jù)以往的研究結(jié)果,認(rèn)為脆性指數(shù)與儲(chǔ)層物性存在負(fù)相關(guān)關(guān)系。脆性較高的巖石受到外界載荷作用,在應(yīng)變很小的情況下發(fā)生破裂;而脆性較低的巖石則會(huì)產(chǎn)生較大的應(yīng)變。孔隙度是巖石受力后產(chǎn)生應(yīng)變的一個(gè)重要控制因素,孔隙度大,巖石受力時(shí)產(chǎn)生的應(yīng)變則更為明顯;孔隙度小,應(yīng)變則不明顯。因此,孔隙度與巖石脆性指數(shù)負(fù)相關(guān)。

在常規(guī)測(cè)井資料解釋過(guò)程中,孔隙度由縱波時(shí)差、密度、補(bǔ)償中子等資料計(jì)算得出。孔隙度與上述測(cè)井?dāng)?shù)據(jù)之間存在明顯的對(duì)應(yīng)關(guān)系,即孔隙度與縱波時(shí)差、補(bǔ)償中子正相關(guān),與密度負(fù)相關(guān)。因此,縱波時(shí)差、補(bǔ)償中子與巖石脆性指數(shù)負(fù)相關(guān),密度與巖石脆性指數(shù)正相關(guān)。

利用常規(guī)測(cè)井資料,將縱波時(shí)差、補(bǔ)償中子的乘積與密度的比值定義為物性指數(shù),表示單位密度下地層縱波時(shí)差與補(bǔ)償中子乘積的定量測(cè)井響應(yīng)特征。

PI=AC·CNL/DEN

(8)

式中:PI為物性指數(shù),無(wú)量綱;AC為縱波時(shí)差測(cè)井值,μs/m;CNL為補(bǔ)償中子測(cè)井值,%;DEN為密度測(cè)井值,g/cm3。

針對(duì)F70井,繪制物性指數(shù)與脆性指數(shù)之間對(duì)應(yīng)關(guān)系散點(diǎn)圖(圖3),兩者之間存在明顯的負(fù)相關(guān)關(guān)系,隨著物性指數(shù)的增加,脆性指數(shù)呈現(xiàn)出減小的趨勢(shì)。建立兩者之間回歸公式,得出了利用物性指數(shù)計(jì)算巖石脆性指數(shù)的數(shù)學(xué)表達(dá)式:

圖3 PI與BI對(duì)應(yīng)關(guān)系

BIa2=1 158.7PI-0.573 2

(9)

該式為數(shù)據(jù)驅(qū)動(dòng)模型,物理量的單位不參與運(yùn)算。

式中:BIa2為利用物性指數(shù)計(jì)算的脆性指數(shù),%。

3.3 脆性指數(shù)多元回歸預(yù)測(cè)模型

公式(7)、(9)分別通過(guò)巖石的巖性和物性進(jìn)行脆性指數(shù)預(yù)測(cè)。根據(jù)相關(guān)性分析,巖性指數(shù)與脆性指數(shù)的擬合優(yōu)度為0.733,物性指數(shù)與脆性指數(shù)的擬合優(yōu)度為0.719。為了進(jìn)一步提高脆性指數(shù)常規(guī)測(cè)井計(jì)算模型的精度,利用BI對(duì)LI、PI進(jìn)行多元回歸,并得到相應(yīng)的函數(shù)關(guān)系,擬合優(yōu)度為0.789。通過(guò)多元回歸,實(shí)現(xiàn)了常規(guī)測(cè)井信息的全面提取,建立了脆性指數(shù)綜合測(cè)井預(yù)測(cè)模型,記為BIa。

BIa=-0.057 87LI-0.007 81PI+76.288 73

(10)

4 應(yīng)用實(shí)例

利用公式(10),針對(duì)Y162井長(zhǎng)8段2 610.8～2 624.6 m計(jì)算巖石脆性指數(shù),并繪制脆性指數(shù)測(cè)井解釋成果圖(圖4)。連續(xù)的曲線(xiàn)為利用常規(guī)測(cè)井資料計(jì)算的脆性指數(shù)BIa,離散桿狀圖為通過(guò)巖石聲波特性室內(nèi)測(cè)試得到的脆性指數(shù)BI,兩者在趨勢(shì)上存在明顯的一致性。通過(guò)繪制兩者對(duì)應(yīng)關(guān)系散點(diǎn)圖(圖5),可以看出,散點(diǎn)數(shù)據(jù)存在明顯的相關(guān)性,相關(guān)系數(shù)為0.936。根據(jù)該實(shí)例,測(cè)井預(yù)測(cè)脆性指數(shù)與巖心測(cè)試結(jié)果匹配程度較高。

圖4 Y162井長(zhǎng)8段測(cè)井解釋綜合成果

圖5 BIa與BI對(duì)應(yīng)關(guān)系

對(duì)于石油、天然氣勘探開(kāi)發(fā)而言,當(dāng)脆性指數(shù)大于60%時(shí),儲(chǔ)層壓裂改造過(guò)程中容易形成規(guī)模較大的縫網(wǎng)系統(tǒng),從而為石油、天然氣在儲(chǔ)層內(nèi)的滲流提供有利通道[5];相反,當(dāng)脆性指數(shù)小于40%時(shí),地層難以在壓裂過(guò)程中形成復(fù)雜縫網(wǎng),不利于油氣田的高效勘探開(kāi)發(fā)。針對(duì)Y162井,選取脆性指數(shù)較高的2 613～2 617 m進(jìn)行射孔,采用體積壓裂的儲(chǔ)層改造方式,加陶粒40 m3、砂比12.6%、排量5 m3/min,試油日產(chǎn)純油12.60 t、試采初期日產(chǎn)油1.81 t。該井脆性指數(shù)、單井產(chǎn)能適中,生產(chǎn)效果與脆性指數(shù)評(píng)價(jià)結(jié)果相符合。

5 結(jié)論

1)在聲波振動(dòng)作用下,機(jī)械振動(dòng)使得巖石產(chǎn)生剪切變形和壓縮變形,進(jìn)而形成了在巖石中傳播的剪切波和壓縮波即橫波和縱波。利用陣列聲波測(cè)井、偶極聲波測(cè)井、巖石聲波特性室內(nèi)測(cè)試,可獲取橫波時(shí)差、縱波時(shí)差等參數(shù),用以求取模量、泊松比等巖石彈性力學(xué)參數(shù),并進(jìn)而計(jì)算巖石的性指數(shù)。

2)Y162井長(zhǎng)8段的楊氏模量平均值為52.14 GPa;泊松比平均值為0.257;脆性指數(shù)平均值為50.13%, F70井長(zhǎng)8段楊氏模量楊氏模量平均值為38.77 GPa,泊松比平均值為0.269,脆性指數(shù)平均值為46.46%。

3)針對(duì)F70井,構(gòu)建了巖性系數(shù),通過(guò)擬合回歸,建立了基于巖性指數(shù)的脆性指數(shù)測(cè)井預(yù)測(cè)模型,擬合優(yōu)度為0.733;根據(jù)孔隙度與脆性指數(shù)的負(fù)相關(guān)性,利用縱波時(shí)差、密度、補(bǔ)償中子,構(gòu)建了物性指數(shù),并擬合出基于物性指數(shù)的脆性指數(shù)測(cè)井預(yù)測(cè)模型,擬合優(yōu)度為0.719;利用脆性指數(shù)對(duì)新構(gòu)建的巖性、物性指數(shù)進(jìn)行了多元回歸,并得到相應(yīng)的函數(shù)關(guān)系,擬合優(yōu)度為0.789。

4)計(jì)算了Y162井長(zhǎng)8段脆性指數(shù),并利用該井巖石聲波特性室內(nèi)實(shí)驗(yàn)測(cè)試數(shù)據(jù)對(duì)計(jì)算結(jié)果進(jìn)行了驗(yàn)證,預(yù)測(cè)值與測(cè)試值變化趨勢(shì)一致,存在明顯的對(duì)應(yīng)關(guān)系,相關(guān)系數(shù)為0.936。