馮笑含,趙萬(wàn)春,王婷婷
(東北石油大學(xué),黑龍江 大慶 163318)
目前,非常規(guī)油氣資源的開采需通過(guò)大規(guī)模體積壓裂方式制造壓裂縫網(wǎng)對(duì)儲(chǔ)層進(jìn)行改造。在非均質(zhì)儲(chǔ)層壓裂效果評(píng)價(jià)中,巖石脆性是天然裂縫生成及縫網(wǎng)形態(tài)的重要影響因素之一[1]。因此,在壓裂方案中,需對(duì)儲(chǔ)層進(jìn)行脆性評(píng)價(jià)分析,判斷縫網(wǎng)形態(tài)及連通范圍[2]。目前,巖石脆性評(píng)價(jià)的方法多達(dá)數(shù)十種,各自具有其優(yōu)勢(shì)和局限性[3-13],而現(xiàn)場(chǎng)普遍接受并采用主要是彈性參數(shù)法[14]及礦物成分分析法[15]。彈性參數(shù)法認(rèn)為巖石彈性模量的升高及泊松比的降低均能增強(qiáng)巖石的脆性,巖石礦物成分分析法認(rèn)為巖石中石英等礦物的含量越高,巖石脆性指數(shù)越大。Burnaman等[16-17]認(rèn)為脆性礦物是導(dǎo)致壓裂過(guò)程中產(chǎn)生不同縫網(wǎng)形態(tài)的重要原因,并優(yōu)化了礦物成分分析方法。丁文龍等[18]認(rèn)為石英、長(zhǎng)石和碳酸鹽巖含量增高均會(huì)使巖石脆性增大,儲(chǔ)層易產(chǎn)生天然裂縫?,F(xiàn)有的研究手段都是從巖心的整體失穩(wěn)進(jìn)行描述,在微觀破裂領(lǐng)域無(wú)法進(jìn)行有效分析,而礦物組分法等類似方法將巖石礦物定性為脆性、非脆性礦物,與實(shí)際情況有較大誤差。因此,需要建立新的脆性評(píng)價(jià)方法,用以在微觀領(lǐng)域進(jìn)行研究。
采用逐步回歸分析方法,對(duì)礦物成分進(jìn)行篩選分析,確定重點(diǎn)影響巖石彈性模量及泊松比的參數(shù)。當(dāng)被解釋變量同時(shí)受多個(gè)因素影響時(shí),利用求解求逆緊湊變換法和雙向檢驗(yàn)法,逐個(gè)分析解釋變量對(duì)被解釋變量的貢獻(xiàn)程度,建立最優(yōu)回歸方程。該方法通過(guò)反復(fù)檢驗(yàn),可剔除具有多重共線性的解釋變量,并逐步引入貢獻(xiàn)最大的解釋變量,使回歸方程更全面和精準(zhǔn)。模型包含的被解釋變量為巖石的彈性模量和泊松比,解釋變量為巖石礦物成分含量測(cè)定值,建立的回歸模型為:
(1)
通過(guò)實(shí)驗(yàn)得到n組實(shí)驗(yàn)數(shù)據(jù),根據(jù)實(shí)驗(yàn)數(shù)據(jù)計(jì)算解釋變量與被解釋變量間的相關(guān)系數(shù)。并根據(jù)不同變量間的相關(guān)系數(shù)建立回歸分析初始相關(guān)系數(shù)矩陣:
(2)
式中:rij為變量i、j(同時(shí)包含解釋變量與被解釋變量)間的相關(guān)系數(shù);R為初始相關(guān)系數(shù)矩陣。
根據(jù)初始相關(guān)系數(shù)矩陣,計(jì)算不同解釋變量與被解釋變量間的偏回歸平方和,表征不同解釋變量在回歸方程中的貢獻(xiàn)度:
(3)
式中:ui為回歸變量貢獻(xiàn)度。
選擇最大偏回歸平方和uimax對(duì)應(yīng)的解釋變量作為引入變量,經(jīng)F檢驗(yàn)法進(jìn)行分析,當(dāng)檢驗(yàn)結(jié)果大于經(jīng)驗(yàn)性確定的臨界閾值,該變量可引入回歸方程作為回歸變量,反之,該變量無(wú)法引入回歸方程,需要剔除。
引入成功后,需將所引用的解釋變量作為主元素,通過(guò)求解求逆緊湊變換法剔除該變量產(chǎn)生的貢獻(xiàn)度,計(jì)算新的相關(guān)系數(shù)矩陣,并通過(guò)新相關(guān)系數(shù)矩陣重新計(jì)算各解釋變量偏回歸平方和,求出剩余變量貢獻(xiàn)度,從而通過(guò)F檢驗(yàn)確定新引入的回歸方程變量,直至回歸方程無(wú)法引用與剔除任何變量,此時(shí),所選中解釋變量即為回歸方程的引入自變量。
選定回歸方程解釋變量后,可通過(guò)式(4)計(jì)算各變量對(duì)應(yīng)的回歸系數(shù),未被選中的變量的回歸系數(shù)均為0,計(jì)算完成后將實(shí)驗(yàn)數(shù)據(jù)平均值代入回歸方程計(jì)算回歸系數(shù)。
(4)
式中:σi為解釋變量i的實(shí)驗(yàn)數(shù)據(jù)方差;σY為被解釋變量Y的實(shí)驗(yàn)數(shù)據(jù)方差;l為回歸方程最終選擇變量的數(shù)量。
將計(jì)算得到的回歸系數(shù)帶入逐步回歸方程中,即可得到巖石彈性模量及泊松比的回歸模型。通過(guò)該模型,可由巖石礦物成分計(jì)算得到不同區(qū)域內(nèi)巖石的彈性模量及泊松比。
在微觀研究中,巖石整體可離散出無(wú)數(shù)微元體(圖1),不同微元體中巖石礦物含量具有差異性,根據(jù)離散后不同微元體礦物含量差異即可計(jì)算巖石局部脆性。
圖1 巖心離散化示意圖
此時(shí),微元體方位可由(a,b,c)表征,根據(jù)礦物微區(qū)衍射實(shí)驗(yàn)測(cè)定該單元礦物成分,并通過(guò)回歸方程計(jì)算該單元的彈性模量及泊松比,記為E(a,b,c)與v(a,b,c)。結(jié)合彈性參數(shù)法表征巖石脆性分布,可得離散化巖石脆性分布規(guī)律:
(5)
式中:BRIT(a,b,c)為巖石微元體(a,b,c)處的脆性;Emax為巖石彈性模量最大值,GPa;Emin為巖石彈性模量最小值,GPa;vmax為巖石泊松比最大值;vmin為巖石泊松比最小值。
該模型從微觀層面建立了巖石脆性與巖石力學(xué)性質(zhì)之間的聯(lián)系,計(jì)算離散化巖石局部脆性分布規(guī)律,定量、精準(zhǔn)地表征巖石礦物成分對(duì)微元體脆性的影響,精確判斷巖石壓裂過(guò)程中的受力薄弱點(diǎn),確定巖石易產(chǎn)生微觀破裂的方位,從而根據(jù)微裂縫萌生及擴(kuò)展規(guī)律判斷巖石失穩(wěn)狀態(tài)下的縫網(wǎng)連通擴(kuò)展形態(tài),為巖石從整體失穩(wěn)研究過(guò)度到微觀劣化研究提供了科學(xué)探索手段及研究基礎(chǔ)。
實(shí)例巖心取自遼河油田A區(qū)塊同一目的層,共計(jì)20塊標(biāo)準(zhǔn)巖心。通過(guò)真三軸抗壓實(shí)驗(yàn)測(cè)量巖心彈性模量及泊松比,并選取實(shí)驗(yàn)后的巖心碎塊制成薄片,通過(guò)X射線衍射儀進(jìn)行微區(qū)礦物成分衍射實(shí)驗(yàn),測(cè)量不同巖心礦物組成成分(表1)。
表1 巖心實(shí)驗(yàn)測(cè)定結(jié)果
由表1可知:巖心的彈性模量測(cè)量結(jié)果為20.61~30.98 GPa,泊松比測(cè)量結(jié)果為0.279~0.323;巖心主要由9種礦物成分構(gòu)成,分別為石英、鉀長(zhǎng)石、斜長(zhǎng)石、方解石、白云石、方沸石、菱鐵礦、黃鐵礦以及黏土礦物。將實(shí)驗(yàn)數(shù)據(jù)帶入計(jì)算模型中,建立初始相關(guān)矩陣,巖心彈性模量及泊松比的解釋變量經(jīng)多次引入及剔除后,通過(guò)式(4)計(jì)算得到彈性模量及泊松比的回歸系數(shù)(表2)。
表2 回歸系數(shù)計(jì)算結(jié)果
將逐步分析回歸系數(shù)結(jié)果代入基礎(chǔ)模型中,可以建立巖石礦物成分對(duì)力學(xué)性質(zhì)影響的回歸方程:
E=37.6374+0.2020x石英-0.3856x方解石+
0.0958x方沸石-0.6009x菱鐵礦-0.6816x黏土
(6)
v=0.2472-0.0007x石英+0.0009x方解石+
0.0008x白云石+0.0015x黏土
(7)
分析式(6)、(7)可知:在模擬遼河油田地層實(shí)際儲(chǔ)層參數(shù)情況下,剔除非顯著性因素后,巖石彈性模量主控因素包括石英、方沸石、菱鐵礦、方解石以及黏土含量,其數(shù)值隨石英、方沸石的含量升高而增加,隨菱鐵礦、方解石以及黏土的含量升高而降低;巖石泊松比主控因素包括石英、白云石、方解石以及黏土含量,其數(shù)值隨白云石、方解石以及黏土的含量升高而增加,隨石英含量的升高而降低。
對(duì)巖體微觀力學(xué)性質(zhì)與礦物成分間逐步回歸模型進(jìn)行顯著性檢驗(yàn)分析,結(jié)合彈性參數(shù)法,將式(6)、(7)帶入式(5),建立了新的巖體微觀脆性評(píng)價(jià)方法:
(8)
式(8)為遼河油田A區(qū)塊巖石微觀脆性評(píng)價(jià)模型,該模型從微觀礦物成分分析入手,建立了微觀礦物成分與巖體脆性間的定量表征關(guān)系。通過(guò)該方法計(jì)算巖體脆性分布,忽略了巖體宏觀力學(xué)特征分析,從微觀層面表征儲(chǔ)層壓裂過(guò)程中易破裂區(qū)域,可為壓裂選層提供精細(xì)量化技術(shù)指導(dǎo)。
將巖心實(shí)驗(yàn)測(cè)定結(jié)果帶入新脆性評(píng)價(jià)模型,計(jì)算該層段巖石脆性分布狀態(tài),并與礦物成分分析法、彈性參數(shù)法計(jì)算結(jié)果及現(xiàn)場(chǎng)測(cè)定脆性對(duì)比分析(圖2)。由圖2可知:①巖心脆性的整體變化趨勢(shì)是相同的;②巖石礦物成分分析法計(jì)算的脆性指數(shù)與現(xiàn)場(chǎng)實(shí)際情況相比嚴(yán)重偏低,主要原因是該方法僅對(duì)礦物成分進(jìn)行了定性分析,因此,通過(guò)該結(jié)果對(duì)壓裂后縫網(wǎng)形態(tài)判斷會(huì)產(chǎn)生較大偏差;③彈性參數(shù)法計(jì)算的脆性指數(shù)與現(xiàn)場(chǎng)測(cè)定結(jié)果對(duì)比差值較大,存在一定的計(jì)算誤差,主要原因是彈性參數(shù)法所需的輸入?yún)?shù)——彈性模量是通過(guò)室內(nèi)實(shí)驗(yàn)分析測(cè)定的,與現(xiàn)場(chǎng)聲波測(cè)定結(jié)果存在明顯差異;④新方法計(jì)算的脆性指數(shù)變化規(guī)律與現(xiàn)場(chǎng)測(cè)定結(jié)果誤差最小,可有效代替彈性參數(shù)法進(jìn)行脆性分析,且該方法結(jié)合礦物成分對(duì)儲(chǔ)層脆性的影響,從微觀層面研究?jī)?chǔ)層脆性變化規(guī)律,能夠更容易判斷裂縫破裂區(qū)域,預(yù)測(cè)壓裂縫網(wǎng)形態(tài),為現(xiàn)場(chǎng)施工提供更科學(xué)有效的判別依據(jù)。
圖2 不同計(jì)算方法脆性指數(shù)結(jié)果
(1) 針對(duì)遼河油田實(shí)際儲(chǔ)層情況,采用逐步回歸分析方法,針對(duì)不同礦物成分篩選影響巖石彈性模量、泊松比的主控因素,剔除弱影響參數(shù),建立了巖石力學(xué)參數(shù)逐步回歸模型。通過(guò)彈性參數(shù)法,并引入巖石力學(xué)參數(shù)逐步回歸模型,建立了巖石脆性計(jì)算新方法。該方法可以從微觀領(lǐng)域分析巖心脆性變化規(guī)律,并表征巖心從局部劣化至整體失穩(wěn)的變化特征,實(shí)現(xiàn)巖心離散化后局部脆性表征。
(2) 巖石彈性模量主要由石英、方解石、方沸石、菱鐵礦以及黏土的含量影響,巖石泊松比主要由石英、方解石、白云石以及黏土的含量影響。
(3) 多種巖石脆性計(jì)算方法結(jié)果與實(shí)際測(cè)定值進(jìn)行了對(duì)比,新方法從微觀表征巖石脆性變化規(guī)律,計(jì)算結(jié)果與現(xiàn)場(chǎng)測(cè)定結(jié)果誤差最小,同時(shí)精細(xì)量化了不同礦物成分對(duì)儲(chǔ)層脆性的影響,修正了礦物分析法計(jì)算結(jié)果與實(shí)際情況偏差較大的現(xiàn)象,可精確判斷巖石壓裂過(guò)程中的失穩(wěn)破裂方位,對(duì)現(xiàn)場(chǎng)壓裂施工效果分析更具有參考價(jià)值。