曹丹平, 韓金鑫*, 肖竣夫, 劉強, 付飛琪, 劉迪
1 中國石油大學(華東)地球科學與技術學院, 青島 266580 2 中國石油集團東方地球物理勘探有限責任公司物探技術研究中心, 河北涿州 072750
隨著社會經(jīng)濟的發(fā)展以及環(huán)境壓力的不斷提高,全國能源行業(yè)肩負著“碳達峰、碳中和”的任務,頁巖氣逐漸成為我國重要的低碳能源(劉合等,2021;田新等,2023).而頁巖本身在自然狀態(tài)下具有低孔低滲的特性,難以直接獲得頁巖油氣資源,必須通過壓裂的方式來形成復雜的裂縫網(wǎng)絡,從而增加頁巖地層的滲透能力來實現(xiàn)頁巖儲層的有效開發(fā)(賈長貴等,2012,劉喜武等,2022).其中,脆性是開展頁巖可壓裂性評價和壓裂方案設計的重要指標,提高頁巖儲層的脆性評價精度對頁巖油氣開發(fā)具有重要意義.
常用的脆性評價方法主要包括彈性參數(shù)脆性評價方法和礦物組分脆性評價方法等.其中,彈性參數(shù)脆性評價方法主要利用楊氏模量和泊松比參數(shù)來計算脆性指數(shù),通常楊氏模量高、泊松比低的巖石具有較高的脆性.基于礦物組分的脆性評價方法主要利用脆性礦物成分含量與總礦物含量的比值來計算脆性指數(shù),石英等脆性礦物含量越高,巖石脆性越強,在頁巖儲層壓裂時越容易形成復雜的裂縫網(wǎng)絡.兩種方法各有特點且相對獨立,并均在實際應用中取得了很好的成效(Jin et al., 2014;Harbert et al.,2020).
采用頁巖中礦物成分含量來評價頁巖儲層脆性時,其關鍵在于確定頁巖儲層中的脆性礦物成分及其含量,但巖石礦物組分多種多樣,不同區(qū)域、不同儲層類型往往發(fā)育不同的脆性礦物,在巖性比較簡單的區(qū)域可以只考慮石英礦物的影響,而在復雜區(qū)域計算脆性礦物指數(shù)時則需要充分考慮石英、長石、方解石、白云石等更多礦物成分的影響(Zhang et al., 2017;張宜,2019).基于礦物含量的脆性評價方法在計算巖石脆性指數(shù)時,每種礦物往往都被賦予相同的權重系數(shù),也就是將各礦物對脆性指數(shù)的影響等同對待.實際上,不同礦物的脆性是不一樣的,對巖石整體脆性特征的影響程度也具有差異性(胡明毅等,2015;白佳佳,2022),在實際應用中應該充分考慮不同礦物成分對巖石整體脆性的差異性影響,從而有效提高巖石脆性計算結果的合理性及準確性.
如何標定不同脆性礦物對巖石脆性整體特征的定量影響程度是有效提高脆性評價精度的關鍵.不同脆性礦物成分的變化一定會改變巖石的彈性特征,導致楊氏模量、泊松比及其所反映的巖石力學性質(zhì)發(fā)生改變.通過三軸壓縮巖石物理實驗有助于定量分析不同礦物成分對巖石整體脆性在貢獻程度方面的差異性,并可以進一步采用礦物楊氏模量與泊松比的比值來表征不同礦物對巖石脆性特征的差異性影響.但不同礦物成分對巖石整體脆性特征的影響是一個復雜的非線性問題,傳統(tǒng)的實驗難以精確標定不同礦物及其含量對巖石脆性的定量影響關系.近年來,數(shù)字巖心技術在計算巖石彈性參數(shù)等方面發(fā)展迅速(Hou and Cao, 2022),通過對不同礦物組成的巖石圖像進行數(shù)值模擬(Zhu et al.,2017;Das et al., 2019;Cao et al., 2022),能夠直觀反映不同礦物成分及其不同組合方式下的彈性特征(Haddad and Sepehrnoori,2016),通過彈性參數(shù)計算的脆性指數(shù)能夠更加直接地反映出不同礦物組分對巖石脆性的影響,為定量評價不同脆性礦物成分對巖石脆性的影響提供了優(yōu)越的實驗條件,也為采用基于彈性參數(shù)的巖石脆性計算方法來定量標定不同的脆性礦物對巖石整體脆性的差異性影響程度提供了契機.
本文聯(lián)合彈性參數(shù)和礦物組分兩種脆性評價方法,通過開展數(shù)字巖石物理實驗來定量刻畫脆性礦物成分相對于黏土礦物對巖石整體脆性的貢獻程度,在彈性特征的約束下求取不同礦物成分與脆性指數(shù)的定量關系及權重系數(shù),提出基于礦物成分貢獻程度自適應變化的脆性評價方法,避免了不同礦物平均化對巖石脆性指數(shù)計算的影響,更加準確合理地體現(xiàn)了不同礦物組分對頁巖脆性的貢獻.
本文通過固體波動方程,模擬了巖石物理應力應變實驗.固體運動由波動方程給出:
(1)
其中ρ為固體密度,u為位移場,S為壓力場,F為體力項.
將每種礦物都假設為線性彈性的材料,利用胡克定律來模擬固體礦物的本構關系:
σij=λδijεαα+2μεij,
(2)
式中λ和μ是拉梅常數(shù);δij和εij分別是應力張量和應變張量的元素.
在頁巖頂部施加垂直于巖石頂部和平行于巖石頂部的正弦位移:
u=A*sin(2πtf),
(3)
其中u是巖石頂部的位移;f是正弦位移的頻率;t是正弦位移的時間,A為振幅,大小為10-9m.
為了討論礦物成分對巖石脆性的定量影響關系,本文構建了一個頁巖數(shù)字巖石物理模型,模型大小為5 mm×5 mm,來模擬微觀下頁巖巖心的實際狀態(tài)(Wu et al., 2020),采用有限元方法對模型進行網(wǎng)格剖分,以滿足模型求解時的穩(wěn)定性.為了更好地模擬實際頁巖巖心中礦物的分布情況,在數(shù)字模型中礦物隨機分布,不同的顏色代表不同的礦物,包含石英、黏土、有機質(zhì)、長石、方解石、白云石、黃鐵礦七種礦物.模型邊界條件設置是數(shù)字巖石物理實驗在模擬縱橫波傳播時效果好壞的關鍵,為滿足模型的穩(wěn)定性,更加符合數(shù)字巖石物理實驗的條件,本文在模擬縱波傳播時,巖石兩側(cè)的邊界固體位移在X方向上為0(Quintal et al., 2016),使得應變只能發(fā)生在縱向方向;在模擬橫波傳播時,巖石兩側(cè)的邊界為周期性邊界條件(朱偉等,2020),在模擬縱橫波傳播時,底部邊界固體的位移皆為0(圖1).
圖1 頁巖數(shù)字巖石物理模型及數(shù)值模擬邊界條件設置 (a) 縱波傳播時邊界條件; (b) 橫波傳播時邊界條件.Fig.1 Digital petrophysical model of shale and setting of boundary conditions for numerical simulation(a) Boundary conditions of longitudinal wave propagation; (b) Boundary conditions of shear wave propagation.
在模擬縱波傳播求取縱波速度時,在巖石頂部施加一個垂直于巖石頂部的正弦振動,X方向為水平方向,Y方向為豎直方向.計算Y方向的平均應力和平均應變,Y方向的平均應變是通過應變除以Y方向的巖石長度來計算的,有效P波模量Re(M)是由巖石頂部平均應力應變計算得到的,其計算模型為
(4)
(5)
其中,Re(M)表示有效P波模量;σP表示峰值處應力振幅;εP表示峰值處的應變振幅;t0表示峰值處的時間;VP為縱波速度;ρ為礦物的密度,f為100 Hz.
在求取橫波速度時,在巖石頂部施加一個平行于巖石頂部的正弦振動,計算X方向的平均應力和平均應變.剪切模量Re(G)計算方法與有效P波模量的計算方法類似,皆是由巖石頂部平均應力應變比計算得到的,其計算模型具體為
(6)
(7)
其中,Re(G)表示剪切模量;σS表示峰值處的應力振幅,εS表示峰值處的應變振幅;t1表示峰值處的時間;VS為橫波速度,f為100 Hz.
因此,通過數(shù)字巖石物理實驗可以直觀地獲得不同礦物組分條件下的縱橫波速度等彈性參數(shù),為基于彈性參數(shù)和礦物成分開展脆性評價及聯(lián)合分析提供了橋梁.
根據(jù)數(shù)字巖石物理實驗獲得的縱橫波速度可以計算出動態(tài)楊氏模量、體積模量及泊松比(Fj?r,2019),得到不同礦物組成下巖石的彈性參數(shù),其計算模型為
(8)
(9)
(10)
其中Kdyn為動態(tài)體積模量;Edyn為動態(tài)楊氏模量;ν為泊松比.
利用Herbert等人利用的彈性參數(shù)脆性評價方法,對楊氏模量和泊松比做處理后計算得到巖石脆性指數(shù)B1,其計算模型為
(11)
(12)
(13)
其中υ為泊松比,E是楊氏模量,單位為MPSI.BRITE代表楊氏模量對脆性指數(shù)的影響;BRITυ代表泊松比對脆性指數(shù)的影響.楊氏模量的范圍界定在1~8 MPSI,泊松比范圍界定在0.15~0.45之間,當E=8 MPSI且υ=0.15時,認為頁巖將“100%”易碎,可壓裂性好;當E=1 MPSI,υ=0.45時,巖石被視為“0%”脆性的巖石,其可壓裂性差.
本文也采用了Guo等(2013)利用的彈性參數(shù)脆性評價方法B2,利用楊氏模量與泊松比的比值來計算脆性大小:
(14)
(15)
基于礦物成分的脆性指數(shù)評價是對頁巖儲層開展脆性評價的重要方法(黃銳等,2014;Ye et al., 2020),它們的共同特點是將脆性礦物與總礦物含量的比值作為巖石脆性指數(shù),但在具體呈現(xiàn)方式上具有差異性.可以將石英、長石、方解石、白云石等多種脆性礦物的共同影響全部列出,如公式(16)所示(陳吉和肖賢明,2013);也可以將石英與碳酸鹽礦物的影響合并在一起考慮來作為對巖石脆性的影響,如公式(17)所示(李鉅源,2013);也可以僅考慮石英礦物的影響,如公式(18)所示(才博等,2012);還可以進一步考慮楊氏模量與泊松比在礦物脆性評價中對巖石脆性的影響,如公式(19)所示(刁海燕,2013).公式(16)—(19)為
(16)
(17)
(18)
B6=
(19)
其中B3、B4、B5、B6代表各種脆性指數(shù)計算方法;W代表礦物在總礦物含量中所占比重;YM代表礦物的楊氏模量;PR代表礦物的泊松比.礦物脆性評價方法B3、B4將每種礦物對巖石脆性的貢獻視為等同的,例如B3將白云石、方解石、長石、石英礦物對巖石脆性的影響視為等同,而實際上,由于每種礦物的化學成分不同,其彈性特征也會存在差異,對巖石脆性的貢獻也會存在很大差異.B6將礦物楊氏模量與泊松比的比值作為約束進行巖石的脆性評價,但B6與B3都只將石英礦物作為脆性礦物來開展巖石的脆性評價.
彈性參數(shù)脆性評價方法與礦物脆性評價方法都具有其特點且相對獨立,彈性參數(shù)脆性評價方法能夠從巖石的彈性特征上計算巖石整體的脆性特征,礦物脆性評價方法可以通過巖石的礦物組成對其進行脆性評價,因此本文綜合考慮兩種方法的優(yōu)勢,設計了基于數(shù)字巖石物理的彈性參數(shù)和礦物組分脆性評價方法聯(lián)合實驗流程(圖2),構建基于礦物成分貢獻程度自適應變化的脆性評價方法.
圖2 數(shù)字巖石物理實驗流程圖Fig.2 Flow chart of digital rock physics experiment
首先,提取不同礦物成分的楊氏模量與泊松比,明確礦物彈性特征的差異性;其次,在設計數(shù)字巖石物理實驗時,利用楊氏模量、泊松比及密度來表征不同礦物成分的彈性性質(zhì),通過控制不同礦物成分含量的變化,研究縱橫波速度隨著礦物含量的變化規(guī)律,明確楊氏模量、泊松比與礦物成分間的關系,從礦物的彈性特征上綜合分析不同礦物對頁巖脆性的影響,進而明確不同礦物成分與脆性指數(shù)的定量關系;最后,通過提取不同礦物成分與脆性指數(shù)關系曲線的斜率,作為礦物脆性指數(shù)敏感因子(礦物彈性特征標定的自適應系數(shù)),根據(jù)該聯(lián)合實驗流程提出礦物成分貢獻程度自適應變化的脆性評價方法:
(20)
其中B7為構建的基于礦物組分貢獻程度自適應變化的礦物脆性評價方法;i為礦物種類;Wi為地層中礦物的體積百分含量;W總礦物為地層中礦物總含量;Ai為礦物彈性特征標定的自適應系數(shù).不同礦物組分的彈性特征不同,對巖石脆性的影響程度也不同,因此需要針對不同的礦物成分的差異來構建脆性評價方法,本文開展數(shù)字巖石物理的目的就是精確標定不同礦物成分的系數(shù),從而在脆性評價中體現(xiàn)出不同脆性礦物的貢獻程度的差異.
采用巖石物理手冊等資料(Carmichael,1984)中石英、黃鐵礦、方解石、白云石、長石、黏土的楊氏模量、泊松比及密度信息(表1),以此來表征不同的礦物.巖石中礦物組成及其礦物彈性特征的差異都會引起巖石整體彈性特征的變化,從而造成巖石整體脆性發(fā)生改變.在本次數(shù)字巖石物理實驗中,設定了模型中礦物組成的初始值,其中黏土含量40%、長石含量15%,方解石含量10%,石英含量20%、黃鐵礦含量3%,白云石含量10%和有機質(zhì)含量2%(圖3),其中黏土礦物、石英、碳酸巖含量較多,而黃鐵礦與有機質(zhì)含量較少.
表1 不同礦物彈性及密度參數(shù)Table 1 Elasticity and density parameters of different mineral
圖3 模型中礦物組成的初始值Fig.3 Initial values of mineral composition in the model
同時也設定了礦物在頁巖模型中的合理變化范圍(圖4),由于不同工區(qū)對脆性礦物的定義不同,而黏土礦物作為典型的塑性礦物是普遍認可的,可以作為標尺來衡量其他礦物相對于黏土對脆性的貢獻情況.因此本文將黏土礦物對巖石脆性的貢獻假設為0,在控制礦物含量變化時,黏土和其他礦物之間存在礦物替代.即在模型的初始條件下增加或減少其中的一種礦物成分含量,相對應的,黏土含量減少或增加對應的含量,其他礦物含量保持不變,以此來控制礦物含量的變化,研究不同礦物成分相對于黏土礦物對巖石脆性的相對影響.
圖4 不同礦物含量變化范圍Fig.4 Variation range of different mineral contents
頁巖儲層中礦物組成不同,其彈性特征也會存在明顯的差異.圖5是不同礦物成分組成的數(shù)字巖石模型,根據(jù)公式(9)—(10),可以計算出不同礦物成分與動態(tài)楊氏模量和泊松比的關系(圖5),得到了不同礦物成分與模型的縱橫波速度、動態(tài)楊氏模量及泊松比的定量影響關系.隨著石英、黃鐵礦、方解石、白云石、長石礦物含量的增加,模型的縱橫波速度也在不斷增加;而有機質(zhì)含量越高,縱橫波速度越小.其中泊松比反映了巖石的抗破壞能力,楊氏模量反映了巖石在壓裂后保持裂縫穩(wěn)定的能力(Rickman et al., 2008),兩者綜合反映了巖石的力學性質(zhì),楊氏模量越大,泊松比越小,巖石的脆性特征越明顯.隨著石英、黃鐵礦、白云石、長石礦物含量的增加,動態(tài)楊氏模量變大,泊松比減小,礦物與動態(tài)楊氏模量及泊松比關系曲線在圖像上呈交叉形態(tài),可以判定石英、黃鐵礦、白云石、長石,相對于黏土礦物來說,對巖石脆性具有的積極的貢獻;隨著方解石的含量增加,楊氏模量變大,泊松比也隨之增高,由于楊氏模量增加、泊松比減小是判斷礦物是否為脆性礦物的標準,所以并不能從礦物與模型動態(tài)楊氏模量與泊松比的關系上去判斷方解石、有機質(zhì)對巖石脆性的貢獻是積極的還是消極的,還需要進一步的討論.
圖5 不同礦物對彈性參數(shù)的定量影響關系Fig.5 Quantitative relationships between the effects of different minerals on elastic parameters
圖5所得結論表明,僅僅通過楊氏模量和泊松比是無法準確地判斷礦物對巖石脆性的貢獻程度的,因此通過脆性評價方法B1、B2、B3、B4、B5、B6對不同礦物組成的頁巖模型開展脆性評價,研究礦物在不同脆性指數(shù)評價方法下對巖石脆性影響的差異性,結果如圖6所示.
圖6 不同礦物成分對頁巖脆性指數(shù)的影響關系Fig.6 Relationship between different mineral compositions and shale brittleness index
可以發(fā)現(xiàn)隨著礦物含量的增加,不同的脆性評價方法在相同的礦物組成下存在很大的差異,其中脆性指數(shù)數(shù)值大小與變化趨勢發(fā)生比較明顯的變化.對比礦物脆性評價方法B3、B4、B5、B6的結果,B6的脆性指數(shù)的數(shù)值明顯大于B3、B4、B5所得的脆性指數(shù),與B3、B4、B5方法所得結果作比較,B6所得結果中方解石對礦物脆性指數(shù)的影響得到了相反的趨勢.B5、B6僅把單一礦物作為脆性礦物進行脆性指數(shù)分析,B3、B4則把不同脆性礦物對頁巖脆性指數(shù)的影響平均化,導致脆性評價結果難以體現(xiàn)不同礦物成分的差異.結果表明,礦物成分的差異對脆性指數(shù)的數(shù)值影響很大,需要針對脆性礦物的差異來體現(xiàn)不同礦物成分對脆性的作用程度.
彈性參數(shù)脆性評價方法B1、B2所得結果與礦物脆性評價方法B3、B4、B5、B6所得結果進行比較,可以發(fā)現(xiàn),隨著石英、黃鐵礦、白云石、方解石、長石礦物含量的增加,彈性參數(shù)脆性評價結果數(shù)值上大于礦物脆性評價的結果,且礦物含量的增加對兩種彈性參數(shù)脆性評價方法B1、B2所得結果的影響趨勢是一致的,其中石英、黃鐵礦、白云石、方解石、長石都對巖石脆性有積極的響應;有機質(zhì)含量的增加,則會使巖石的脆性減小,而B3、B4、B5、B6方法所得結果沒有體現(xiàn)出有機質(zhì)含量對頁巖脆性指數(shù)的影響.
(21)
(22)
(23)
圖7 礦物與彈性特征表征的脆性指數(shù)關系曲線斜率歸一化結果Fig.7 The slope normalization result of the brittle index relationship curve characterized by minerals and elastic characteristics
(24)
(25)
表2 礦物脆性評價結果與彈性參數(shù)脆性評價結果相關性Table 2 Correlation between mineral brittleness evaluation results and elastic parameter brittleness evaluation results
本文以下寺灣地區(qū)的頁巖數(shù)據(jù)為例開展脆性評價測試,該工區(qū)各礦物的體積分數(shù)如圖8所示(秦曉艷等,2016),對其進行頁巖脆性評價.其中,在15塊頁巖試樣中,黏土礦物相對含量最高,范圍在34.2%~68.6%之間,平均為51.4%;黃鐵礦的相對含量最少,范圍大致分布在0.4%~3.4%之間,平均為1.9%.長石的相對含量在14.4%~42.8%,平均為28.6%,石英的相對含量在7.8%~22.4%,平均為15.1%;碳酸鹽礦物在總礦物含量的比重范圍在0.9%~14.4%之間,平均為7.65%,將碳酸鹽礦物劃分為方解石和白云石,方解石和白云石在總礦物含量的所占比重皆為3.8%,我們?nèi)∶糠N礦物含量的平均值來計算脆性指數(shù).對其應用多種礦物脆性評價方法開展脆性評價,并與本文提出的礦物彈性特征約束下的脆性評價方法計算結果進行比較(圖9).
圖8 下寺灣地區(qū)試樣的礦物組成對比Fig.8 Comparison of mineral composition of samples in Xiasiwan area
圖9 下寺灣地區(qū)試樣不同脆性評價結果對比Fig.9 Comparison of different brittleness evaluation results of specimens in the Xiasiwan area
結果表明,不同的脆性評價方法之間所得的脆性指數(shù)存在著很大的差異.利用B3進行頁巖脆性指數(shù)評價時,雖然將長石、石英、碳酸鹽作為脆性礦物,但是也忽略了黃鐵礦的影響,同時該方法將每種礦物對巖石脆性的貢獻視為等同的;利用B4進行頁巖脆性評價時,該方法將碳酸鹽和石英作為脆性礦物,同樣將石英和碳酸鹽對脆性指數(shù)的貢獻視為等同,且忽略了長石與黃鐵礦等礦物對脆性指數(shù)的影響;應用B5脆性評價方法開展脆性評價時,只將石英作為脆性礦物,而忽略了碳酸鹽等礦物對頁巖脆性指數(shù)的影響;而通過礦物脆性評價方法B6開展脆性評價時,雖然考慮了礦物彈性特征對脆性指數(shù)的影響,但僅把石英作為脆性礦物.而本文構建的B7是各礦物彈性特征約束下的礦物脆性評價方法,能夠體現(xiàn)各礦物對巖石脆性貢獻的差異性,能夠更合理地反映不同脆性礦物對巖石整體脆性特征的影響程度,避免了平均化各礦物對脆性指數(shù)的影響,能夠提高頁巖脆性評價及可壓裂性評價時的合理性.
頁巖儲層的脆性評價是勘探開發(fā)過程中重要的評價內(nèi)容.本文在收集、調(diào)研前人研究成果的基礎上,針對我國頁巖儲層礦物含量復雜以及礦物脆性評價方法存在局限性等問題,通過數(shù)字巖石物理實驗,開展了礦物成分與頁巖脆性關系的研究,聯(lián)合彈性參數(shù)脆性評價方法,提出了根據(jù)礦物成分貢獻程度自適應變化的頁巖脆性礦物指數(shù)計算方法.本文主要取得以下認識:
(1) 刻畫了礦物對頁巖彈性特征影響的差異性.通過數(shù)字巖石物理實驗的方式,得到了礦物含量變化與縱橫波速度的關系,隨著石英、黃鐵礦、長石、方解石、白云石礦物含量的增加,頁巖模型整體的縱橫波速度也在不斷增加.隨著方解石含量增加,楊氏模量變大,泊松比也隨之增高;而隨著石英、黃鐵礦、白云石、長石礦物含量的增加,楊氏模量變大,泊松比減小.
(2) 構建了根據(jù)礦物成分貢獻程度自適應變化的頁巖脆性礦物指數(shù)計算方法.通過數(shù)字巖石物理實驗的方式,明確了不同礦物與脆性指數(shù)的定量關系,基于礦物成分貢獻程度自適應的變化,構建了頁巖脆性礦物指數(shù)計算方法,避免了將各礦物對脆性指數(shù)的影響平均化,準確地表征了不同礦物對巖石脆性貢獻的差異性.將不同礦物脆性評價結果與彈性參數(shù)脆性評價結果進行相關性分析,發(fā)現(xiàn)該方法所得結果與彈性參數(shù)脆性評價方法所得結果間存在較好的相關性,從而驗證了該方法的合理性.
(3) 以下寺灣地區(qū)試樣的礦物組成為例,研究了不同脆性評價方法的差異性.結果表明,利用B3、B4對其進行脆性評價時,平均化了各礦物對脆性指數(shù)影響;B5只考慮了石英作為脆性礦物;B6雖然考慮了礦物彈性特征對脆性指數(shù)的影響,但僅把石英作為脆性礦物,忽略了長石、白云石、黃鐵礦對脆性指數(shù)的影響.本文構建的B7是在基于不同礦物的彈性特征約束下的礦物組分脆性評價方法,綜合考慮了礦物對脆性指數(shù)的貢獻的差異性,能夠更合理地反映不同脆性礦物對巖石整體脆性特征的影響程度,能夠更精確地指導頁巖可壓裂性評價.