頁巖基質(zhì)孔隙油微觀賦存及可動性定量表征
——以東營凹陷沙河街組為例

李俊乾，宋兆京，王民，張鵬飛，蔡建超

1 中國石油大學(xué)(華東)深層油氣全國重點(diǎn)實(shí)驗(yàn)室，青島 266580

2 中國石油大學(xué)(華東)地球科學(xué)與技術(shù)學(xué)院，青島 266580

3 山東科技大學(xué)地球科學(xué)與工程學(xué)院，青島 266590

4 中國石油大學(xué)(北京)油氣資源與探測國家重點(diǎn)實(shí)驗(yàn)室，北京 102249

0 引言

隨著常規(guī)油氣資源開采難度日益增大，包括頁巖油氣在內(nèi)的非常規(guī)源巖層系油氣逐漸成為我國油氣增儲上產(chǎn)的主陣地和戰(zhàn)略性接替能源[1-5]。近年來，中國陸相頁巖油在松遼盆地古龍凹陷白堊系青山口組、渤海灣盆地滄東凹陷古近系孔店組、渤海灣盆地濟(jì)陽坳陷古近系沙河街組、鄂爾多斯盆地三疊系延長組、準(zhǔn)噶爾盆地吉木薩爾凹陷二疊系蘆草溝組、江漢盆地潛江凹陷古近系潛江組以及準(zhǔn)噶爾盆地瑪湖凹陷二疊系風(fēng)城組等多套層系獲得頁巖油工業(yè)油流[6-13]，展現(xiàn)了廣闊的頁巖油開發(fā)前景；國家能源局先后設(shè)立了新疆吉木薩爾國家級陸相頁巖油示范區(qū)(2020 年)、大慶古龍陸相頁巖油國家級示范區(qū)(2021 年)和勝利濟(jì)陽陸相斷陷湖盆頁巖油國家級示范區(qū)(2022 年)，推動了陸相頁巖油的規(guī)?；_發(fā)進(jìn)程。但中國陸相頁巖油的勘探開發(fā)也面臨著諸多挑戰(zhàn)。相比于北美穩(wěn)定克拉通、前陸盆地的海相頁巖油，中國陸相頁巖油發(fā)育于斷陷、內(nèi)陸坳陷盆地，頁巖油儲層具有巖性/巖相變化快、分布不穩(wěn)定、有機(jī)/無機(jī)非均質(zhì)性強(qiáng)、礦物成分多樣、孔隙結(jié)構(gòu)和類型復(fù)雜等特征[1,6,14-15]，確定頁巖油“甜點(diǎn)”段/區(qū)難度較大。因此，找準(zhǔn)頁巖油富集高產(chǎn)的“甜點(diǎn)”位置是中國陸相頁巖油高效開發(fā)的關(guān)鍵環(huán)節(jié)。

頁巖油儲集于頁巖層系地層的多級孔-縫網(wǎng)絡(luò)系統(tǒng)內(nèi)，其中裂縫系統(tǒng)(頁理縫、層理縫、構(gòu)造縫等)中的頁巖油可動能力強(qiáng)，可采性較好，是頁巖油井初期產(chǎn)量的主要來源[13]。但研究結(jié)果顯示中國陸相頁巖油水平井分段壓裂后產(chǎn)量遞減快，一次采收率普遍小于10%[13,16]，這說明頁巖油井排采過程中仍有大量的頁巖油(大于90%)被滯留于頁巖基質(zhì)內(nèi)無法排出，如何高效動用頁巖基質(zhì)孔隙系統(tǒng)中的石油資源(即頁巖基質(zhì)孔隙油)是當(dāng)前頁巖油開采面臨的重要挑戰(zhàn)。因此，頁巖基質(zhì)含油性和頁巖油可動性長期以來作為頁巖油“甜點(diǎn)”評價的重頭戲而受到業(yè)界的高度重視。目前，在該領(lǐng)域已形成了眾多實(shí)驗(yàn)評價方法，主要包括：①含油飽和度指數(shù)法(OSI)[17]、②多溫階熱解法[18-21]、③溶脹法[22]、④T1-T2二維核磁法[23]、⑤溶劑分步萃取法[24-25]、⑥鏡下觀察法[26-27]、⑦飽和-離心-核磁聯(lián)合法[28]、⑧抽提前后對比法[29]等，為確定不同地區(qū)/凹陷的頁巖油吸附量、游離量及可動量提供了技術(shù)支撐。此外，分子模擬技術(shù)在頁巖油多相態(tài)賦存機(jī)理分析方面發(fā)揮了重要作用，揭示了吸附油呈多層狀吸附于頁巖基質(zhì)孔隙表面，游離油呈自由態(tài)分布于孔隙內(nèi)部[30]；吸附油、游離油的賦存特征受到孔隙大小及形態(tài)、溫度及壓力條件、油組成、巖石表面性質(zhì)等多種因素的綜合影響[31-33]。然而，針對頁巖基質(zhì)微納尺度孔隙中不同相態(tài)油的含量、比例、分布及可動性的定量化數(shù)學(xué)描述還比較薄弱，是亟待加強(qiáng)的一個研究方向，也是當(dāng)前頁巖油地質(zhì)評價領(lǐng)域的一大挑戰(zhàn)。

很多學(xué)者從統(tǒng)計學(xué)角度認(rèn)識了巖石物性、孔喉微觀結(jié)構(gòu)、巖石組成、流體性質(zhì)對頁巖油賦存及可動性的影響[24,29,34-40]，但是它們之間的定量耦合關(guān)系還不清楚，無據(jù)可依的情況下，增加了統(tǒng)計分析結(jié)果的多解性。前期筆者通過考慮多層吸附、孔壁表面性質(zhì)和孔喉微觀結(jié)構(gòu)特征，初步構(gòu)建了頁巖微納尺度孔隙中不飽和吸附狀態(tài)(0＜相對壓力P/P0＜1)下單組分烴的吸附量與游離量評價模型[41]；飽和吸附狀態(tài)(P/P0≥1)下單組分烴和混合烴(殘余油)的吸附量與游離量評價模型[34,42]，為頁巖基質(zhì)孔隙油的微觀賦存與可動性定量化研究奠定了理論基礎(chǔ)。在此基礎(chǔ)上，本文針對頁巖基質(zhì)孔隙內(nèi)不同相態(tài)油的含量、比例、分布及可動性開展系統(tǒng)研究。以油-巖相互作用為紐帶，建立了頁巖含油性、儲集性(宏觀物性、孔喉微觀結(jié)構(gòu))及頁巖油可動性之間的耦合關(guān)系，并從數(shù)學(xué)角度揭示了它們之間的內(nèi)在聯(lián)系；結(jié)合物理模擬實(shí)驗(yàn)(離心-核磁等)定量評價了東營凹陷陸相頁巖油微觀賦存及可動性特征。研究成果深化了對陸相頁巖油微觀賦存機(jī)制的認(rèn)識，有助于頁巖油藏精細(xì)化地質(zhì)“甜點(diǎn)”評價。

1 頁巖基質(zhì)孔隙油含量實(shí)驗(yàn)評價

1.1 理論與方法

按照賦存狀態(tài)，頁巖油可分為吸附油、游離油及溶解油[30,42-44]。吸附油是以吸附態(tài)(或類固態(tài))吸著于干酪根、無機(jī)礦物孔隙內(nèi)表面，具有弱可動性，流動較為困難；游離油是以自由狀態(tài)分布于有機(jī)/無機(jī)孔隙、裂縫內(nèi)部，具有強(qiáng)可動性，是現(xiàn)實(shí)可采的資源；溶解油是以互溶態(tài)吸收于干酪根大分子內(nèi)部，賦存空間顯著受限，幾乎不可流動[18,38,42,45]。由上可見，頁巖基質(zhì)孔隙內(nèi)吸附油與游離油應(yīng)是關(guān)注的重點(diǎn)，特別是游離油的含量與分布評價對于頁巖油開采至關(guān)重要。

頁巖基質(zhì)孔隙油的可動性與其賦存狀態(tài)密切相關(guān)[30]。筆者基于頁巖油賦存狀態(tài)及可動性內(nèi)在聯(lián)系，厘定了頁巖基質(zhì)孔隙油分類方案：在一定溫度、壓力條件下，石油在頁巖基質(zhì)孔隙內(nèi)由吸附油和游離油兩部分組成；其中游離油進(jìn)一步分為毛細(xì)管束縛油和可動油。在外界作用力(如離心力、驅(qū)替力)下，吸附油幾乎不可流動，游離油可流動，其可流動能力受頁巖孔喉微觀結(jié)構(gòu)的影響；可流動的游離油稱為可動油，不可流動的那部分游離油稱為毛細(xì)管束縛油。隨外界條件的改變，毛細(xì)管束縛油可轉(zhuǎn)化為可動油，且理論最大可動油量等于游離油量，吸附油量是最小的束縛油量。在頁巖基質(zhì)孔隙油有效分類的基礎(chǔ)上，基于飽和-離心實(shí)驗(yàn)，建立了一套完整的吸附油/游離油含量實(shí)驗(yàn)評價方法及分析流程。

頁巖基質(zhì)孔隙油在離心(排出)過程中受到毛細(xì)管力(Pc)和離心壓力差(△P)的雙重作用。對于某一油-巖體系，毛細(xì)管力是孔隙油排出的阻力，其與孔徑成反比，孔徑越小，阻力越大，油越不容易排出；離心壓力差是促使孔隙油排出的動力，其值越大，油越容易排出。根據(jù)受力分析，孔隙系統(tǒng)中的油存在以下2種情況：①當(dāng)△P＞Pc時，油被排出頁巖孔隙，排出油即為可動油；②當(dāng)△P≤Pc時，油滯留于頁巖孔隙內(nèi)，滯留油即為毛細(xì)管束縛油。然而，在某一離心壓力差條件下，可動用孔隙系統(tǒng)(三維連通網(wǎng)絡(luò))中的油并非瞬時完全排出，而是需要一定的時間由大孔至小孔、由外部至內(nèi)部陸續(xù)排出，即可動油量是隨時間變化的，也就是說孔隙油的排出是個動態(tài)過程。定義某一時刻的可動油比例為巖心累積排出的可動油量占巖心中賦存游離油量的比例，表示如下：

式中，f為可動油比例，無量綱；Qm為累積可動油量，mg/g；Qf為游離油量，即最大可動油量，mg/g。

在離心過程中孔隙油被排出的概率是離心壓力差的函數(shù)，即隨離心壓力差增加，孔隙油被排出的概率增大。結(jié)合動力學(xué)原理，可得在非平衡動態(tài)過程中可動油比例隨時間的變化速率為：

式中，t為時間，s；km為孔隙油排出速率常數(shù)，1/(s·MPa)；kim為孔隙油滯留速率常數(shù)，1/s；△P為離心壓力差，MPa。

式(2)中km(1-f) △P代表孔隙油排出速率；kimf代表孔隙油滯留速率。孔隙油排出速率與離心壓力差呈正相關(guān)、與可動油比例呈負(fù)相關(guān)，即離心壓力差越大、可動油比例越低，孔隙油越易于排出、排出速率越高；反之亦然。在某一不變離心壓力差條件下，隨時間增加，可動油比例增大，孔隙油排出速率逐漸降低；相反，滯留速率與可動油比例呈正相關(guān)，隨可動油比例增加，孔隙系統(tǒng)中滯留的油受阻力越大、越難以排出，滯留速率越高。因此，在孔隙油動態(tài)排出-滯留過程中，排出速率逐漸降低、滯留速率逐漸增大，當(dāng)孔隙油的排出速率等于滯留速率，達(dá)到動態(tài)平衡，可動油比例不再發(fā)生變化，即進(jìn)一步可得：

由式(3)可得：

令△PL=kim/km，將式(1)代入式(4)可得平衡條件下可動油量與離心壓力差的關(guān)系式[46]，如下：

式中，△PL為中值壓力差，即可動油量達(dá)到最大可動油量的一半時對應(yīng)的離心壓力差，MPa?！鱌L反映了孔隙油排出的難易程度，其值越小，代表孔隙油排出速率常數(shù)越大、滯留速率常數(shù)越小，孔隙油更易于排出孔隙。進(jìn)一步，可將式(5)轉(zhuǎn)化為：

采用離心實(shí)驗(yàn)測試的數(shù)據(jù)，建立1/Qm與1/△P之間的線性關(guān)系，可以估算出Qf與△PL值。同時，根據(jù)測試的總含油量，可計算出吸附油量Qa：

式中，Qt為總含油量，mg/g；Qa為吸附油量，mg/g。

1.2 實(shí)例驗(yàn)證與結(jié)果

本次采用的14 個泥頁巖樣品來源于渤海灣盆地濟(jì)陽坳陷東營凹陷內(nèi)9 口油井的主力頁巖層段(沙三上亞段、沙四下亞段)，樣品信息如表1 所示，取心井位分布見文獻(xiàn)[28]。根據(jù)有機(jī)質(zhì)豐度、無機(jī)礦物組分特征(圖1)，將這些樣品分成兩類：I類以富有機(jī)質(zhì)的鈣質(zhì)頁巖為主(含1 個富泥質(zhì)頁巖樣品)，發(fā)育紋層結(jié)構(gòu)；TOC含量1.64 %～2.43 %(平均2.03 %)，黏土礦物、石英、方解石和白云石平均含量分別為38.0 %、21.0 %、22.8 %和9.0 %，以黏土礦物、方解石和石英為主。II類以含有機(jī)質(zhì)的富硅質(zhì)、硅質(zhì)泥巖為主，發(fā)育塊狀結(jié)構(gòu)；另有1 個樣品為含有機(jī)質(zhì)富泥質(zhì)頁巖，發(fā)育紋層結(jié)構(gòu)；TOC含量0.13 %～1.56 %(平均0.71 %)，黏土礦物、石英、方解石和白云石平均含量分別為30.0 %、34.2 %、6.9 %和9.9 %，以石英和黏土礦物為主。為便于敘述，本文將泥頁巖統(tǒng)稱為頁巖。

圖1 頁巖有機(jī)/無機(jī)組分特征Fig. 1 Organic and inorganic composition characteristics of the shales

表1 頁巖樣品基本信息Table 1 Basic information of shale samples

將頁巖樣品制備成巖心柱，經(jīng)過抽真空干燥后，進(jìn)行抽真空加壓飽和輕質(zhì)油(以正十二烷為例)，然后針對飽和巖心柱開展不同轉(zhuǎn)速的離心測試，在每個離心平衡時取出巖心并測試質(zhì)量，對部分樣品同時測試了核磁共振信號。離心實(shí)驗(yàn)的溫度設(shè)定為20 ℃(低溫可抑制烴揮發(fā))、壓力為常壓條件，離心平衡時間為4 h。詳細(xì)實(shí)驗(yàn)流程見文獻(xiàn)[28]。圖2a、2b中展示了14 個巖心的測試結(jié)果，顯示隨離心壓力差增大，可動油量逐漸增加，其變化趨勢可通過式(5)所示的方程描述。通過式(6)擬合實(shí)驗(yàn)數(shù)據(jù)(圖2c、2d)，依據(jù)線性關(guān)系的截距和斜率可得到巖心中游離油量和中值壓力差。進(jìn)而，依據(jù)飽和油量計算出吸附油含量。分析結(jié)果見表2。本次實(shí)驗(yàn)頁巖中吸附油量介于5～30 mg/g之間、游離油量介于5～25 mg/g之間，吸附油占比普遍大于50 wt%，吸/游比(吸附油與游離油質(zhì)量比)主體介于1～2 之間，顯示了頁巖基質(zhì)孔隙中吸附油占優(yōu)勢的基本特點(diǎn)。總體上，富有機(jī)質(zhì)頁巖中吸/游比略高于含有機(jī)質(zhì)頁巖。富有機(jī)質(zhì)頁巖中吸/游比介于1.13～1.97 之間，平均1.50；含有機(jī)質(zhì)頁巖中吸/游比介于0.61～2.36 之間，平均1.43。此外，富有機(jī)質(zhì)頁巖的吸附油、游離油含量普遍大于含有機(jī)質(zhì)頁巖，屬于頁巖油富集的優(yōu)勢巖相類型。根據(jù)吸附油、游離量含量及吸/游比，可將頁巖含油性進(jìn)行分類(圖3)，后文中將結(jié)合可動性進(jìn)一步討論。

圖2 可動油量與離心壓力差的關(guān)系(a，富有機(jī)質(zhì)頁巖；b，含有機(jī)質(zhì)頁巖)和可動油量倒數(shù)與離心壓力差倒數(shù)之間的關(guān)系(c，富有機(jī)質(zhì)頁巖；d，含有機(jī)質(zhì)頁巖)(離心數(shù)據(jù)來源于文獻(xiàn)[28])Fig. 2 Relationship between mobile oil amount and centrifugal pressure difference (a, organic-rich shales; b, organic-bearing shales) and Relationship between reciprocal mobile oil amount and reciprocal centrifugal pressure difference (c, organic-rich shales; d, organic-bearing shales) (centrifugal data from literature [28])

圖3 頁巖基質(zhì)孔隙油賦存量評價圖版Fig. 3 Evaluation chart on the oil storage amount in shale matrix pores

表2 離心實(shí)驗(yàn)測試分析結(jié)果Table 2 Results obtained by the centrifugation test

東營凹陷頁巖油儲層發(fā)育多種儲集空間類型，主要包括碎屑粒間孔、有機(jī)質(zhì)和黏土礦物收縮孔、黏土礦物晶間孔、碳酸鹽晶間孔和溶蝕孔[47]。在巖石組分-孔喉微觀結(jié)構(gòu)匹配關(guān)系分析的基礎(chǔ)上，通過吸附油、游離油含量與巖石組分含量之間的相關(guān)性統(tǒng)計(圖4)，可以揭示頁巖油賦存空間。結(jié)果顯示，隨TOC含量增加，吸附油含量線性增加(圖4a)，有機(jī)質(zhì)含量對所有類型頁巖的吸附油賦存具有相似的控制作用。含有機(jī)質(zhì)頁巖中吸附油量與黏土含量呈一定正相關(guān)性(圖4b)、與方解石含量無明顯相關(guān)性(圖4c)，這是因?yàn)檠芯繀^(qū)具有“有機(jī)質(zhì)-黏土礦物共生、有機(jī)質(zhì)-碳酸鹽礦物分異”的沉積特點(diǎn)[48]，因此在黏土礦物相關(guān)孔隙中富集了一定量的吸附油。圖5a中展示了黏土含量和TOC含量呈正相關(guān)性，也可佐證該結(jié)論。在富有機(jī)質(zhì)頁巖中吸附油主要賦存于有機(jī)質(zhì)孔內(nèi)部，可能有少量吸附油賦存于方解石相關(guān)孔隙內(nèi)。對于游離油而言，富有機(jī)質(zhì)頁巖中游離油主要賦存于方解石、有機(jī)質(zhì)相關(guān)的孔隙中(圖4d、4e)；含有機(jī)質(zhì)頁巖中游離油主要賦存于有機(jī)質(zhì)、石英礦物相關(guān)孔隙中(圖4d、4f)。圖4d中兩個數(shù)據(jù)點(diǎn)偏離趨勢線，是因?yàn)樵搩蓚€樣品具有較高的石英礦物含量，石英相關(guān)孔隙中賦存了較多的游離油，彌補(bǔ)了有機(jī)質(zhì)貧化的負(fù)面影響。同理，圖4e中所示的具有低方解石含量的頁巖中游離油量偏高，是因?yàn)樵擁搸r較高的有機(jī)質(zhì)含量(TOC=2.43%)的貢獻(xiàn)。基于石英含量與TOC含量的負(fù)相關(guān)性(圖5b)，富有機(jī)質(zhì)頁巖中石英含量與游離油量呈現(xiàn)負(fù)相關(guān)性(圖4f)是因?yàn)槭⒑吭黾訉?yīng)有機(jī)質(zhì)含量降低，進(jìn)而降低了游離油量；在含有機(jī)質(zhì)頁巖中隨石英含量增加，盡管TOC含量降低，但碎屑粒間孔提供了一部分游離油儲集空間，使得游離油量增加。

圖4 吸附油量、游離油量與巖石組分含量的關(guān)系 (a) 吸附油量與TOC含量；(b) 吸附油量與黏土礦物含量；(c) 吸附油量與方解石含量；(d) 游離油量與TOC含量；(e) 游離油量與方解石含量；(f) 游離油量與石英含量Fig. 4 Relationships of adsorbed and free amounts of shale oil with rock composition contents (a) adsorbed oil amount and TOC content; (b) adsorbed oil amount and clay mineral content; (c) adsorbed oil amount and calcite content; (d) free oil amount and TOC content; (e) free oil amount and calcite content; (f) free oil amount and quartz content

根據(jù)上述分析，建立了研究區(qū)頁巖基質(zhì)孔隙油微觀賦存模式(圖6)。富有機(jī)質(zhì)頁巖呈現(xiàn)富有機(jī)質(zhì)黏土紋層與富碳酸鹽紋層疊置的“雙層結(jié)構(gòu)”特征，在富有機(jī)質(zhì)黏土紋層的有機(jī)質(zhì)孔內(nèi)賦存吸附油和部分游離油，而在富碳酸鹽紋層的方解石晶間孔和溶蝕孔內(nèi)賦存大量游離油及少量吸附油(圖6a～c)。正是由于在不同紋層結(jié)構(gòu)中頁巖油賦存狀態(tài)及方式的這種差異性，構(gòu)成了有利的“生烴-儲集”紋層組合[10,49]。含有機(jī)質(zhì)頁巖不發(fā)育紋層，吸附油、游離油遍布基質(zhì)孔隙內(nèi)，在有機(jī)質(zhì)孔和黏土礦物晶間孔內(nèi)賦存吸附油，游離油則賦存于有機(jī)質(zhì)孔和顆粒碎屑粒間孔內(nèi)(圖6d～f)。

圖6 頁巖基質(zhì)孔隙油微觀賦存模式 (a) LX884-1 樣品巖心柱；(b) LX884-1 樣品薄片，正交光，×50，發(fā)育黏土、石英、方解石礦物；(c)富有機(jī)質(zhì)頁巖孔隙油微觀賦存模式；(d) N5-3 樣品巖心柱；(e) N5-3 樣品薄片，正交光，×50，發(fā)育石英、斜長石、方解石礦物；(f)含有機(jī)質(zhì)頁巖孔隙油微觀賦存模式Fig. 6 Shale matrix pore oil microscopic occurrence pattern (a) LX884-1 core column; (b) LX884-1 thin slice, orthogonal light, ×50, developed clay, quartz and calcite minerals; (c) microscopic occurrence pattern of pore oil in organic-rich shale; (d)N5-3 core column; (e) N5-3 thin slice, orthogonal light, ×50, developed quartz, plagioclase and calcite minerals; (f) microscopic occurrence pattern of pore oil in organic-bearing shale

2 頁巖基質(zhì)孔隙油含量數(shù)學(xué)描述

頁巖基質(zhì)孔隙中吸附油受范德華力、庫侖力等作用呈多層狀吸附于有機(jī)質(zhì)和無機(jī)礦物孔隙表面，游離油呈自由狀態(tài)分布于孔隙內(nèi)部[32-33]。理論上，頁巖吸附油量受控于含油孔隙比表面積、吸附相密度和吸附相厚度，游離油量則受控于含油孔隙體積與游離相密度[42]，可數(shù)學(xué)表示為：

式中，St為含油孔隙比表面積，m2/g；Vt為含油孔隙體積，10-3cm3/g；ρ1和ρ2分別為吸附油和游離油的平均密度，g/cm3；H為平均吸附厚度，nm。

考慮到現(xiàn)有實(shí)驗(yàn)方法的精度和探測范圍，含油孔隙的體積(Vt)和比表面積(St)很難測試出準(zhǔn)確數(shù)值。比如常見的低溫氮?dú)馕椒?、高壓壓汞法、覆壓孔滲法等測試的巖心孔體積和比表面積，由于使用的探測流體與頁巖油(本次研究為正十二烷)的性質(zhì)不同，實(shí)驗(yàn)探測范圍和真實(shí)情況存在一個系統(tǒng)誤差。此外，濕重法和核磁共振法未考慮頁巖油賦存狀態(tài)，在測試孔體積時也存在一定誤差。假定系統(tǒng)誤差為定值，可將實(shí)際值和測試值進(jìn)行線性轉(zhuǎn)換，如下：

式中，Smt為含油孔隙比表面積實(shí)驗(yàn)測試值，m2/g；Vmt為含油孔隙體積實(shí)驗(yàn)測試值，10-3cm3/g；C1和C2分別為比表面積和孔體積測試的系統(tǒng)誤差，單位分別為m2/g 和10-3cm3/g。

進(jìn)而，可得：

變形可得：

通過建立Qa～Smt、(Vmt-Qf/ρ2) ～Qa之間的線性擬合關(guān)系，可以確定吸附油平均密度和平均厚度。如圖7 所示，本研究中富有機(jī)質(zhì)頁巖的吸附油平均密度為0.8331 g/cm3、吸附油平均厚度為1.7475 nm，含有機(jī)質(zhì)頁巖的吸附油平均密度為0.8067 g/cm3、吸附油平均厚度為0.5734 nm。兩種類型頁巖吸附油的能力差異較大，富有機(jī)質(zhì)頁巖具有更強(qiáng)的油-巖相互作用，這也是相同比表面積情況下富有機(jī)質(zhì)頁巖吸附油含量更高的原因。

圖7 關(guān)鍵吸附參數(shù)評估交匯圖Fig. 7 Cross plots of evaluating the key adsorption parameters

進(jìn)一步根據(jù)孔體積、比表面積、孔隙大小和孔隙形態(tài)之間的關(guān)系[42]，可得到頁巖基質(zhì)孔隙中吸附油/游離油含量與油-巖相互作用、孔喉微觀結(jié)構(gòu)和巖石物性之間的耦合關(guān)系，如下所示：

式中，φ為巖石孔隙度，分?jǐn)?shù)；So為巖石含油飽和度，分?jǐn)?shù)；ρa(bǔ)為巖石視密度，g/cm3；dm為巖石平均孔隙直徑，nm；F為孔隙形狀因子，無量綱，平行板狀孔F=2，柱狀孔F=4，球形孔F=6。

從公式上可以看出，頁巖吸附油與游離油含量是頁巖孔喉微觀結(jié)構(gòu)、巖石物性與油-巖相互作用耦合的結(jié)果。頁巖油吸附量和游離量均受吸附油厚度、孔隙大小及形態(tài)、巖石孔隙度、含油飽和度和巖石視密度的綜合影響；吸附量與吸附油平均密度相關(guān)，游離量與游離油平均密度相關(guān)。以東營凹陷頁巖油儲層為例，可以得到以下幾點(diǎn)認(rèn)識：①不同巖相類型頁巖的油-巖相互作用具有明顯差異性，集中體現(xiàn)于吸附油平均厚度、密度的不同，進(jìn)而直接影響了不同巖相類型頁巖的吸附油含量。富有機(jī)質(zhì)頁巖比含有機(jī)質(zhì)頁巖的吸附油平均密度略高，但富有機(jī)質(zhì)頁巖中吸附油厚度是含有機(jī)質(zhì)頁巖的3 倍，因此富有機(jī)質(zhì)頁巖中油-巖相互作用更強(qiáng)、吸附油量更多。②采用表3 所示的數(shù)據(jù)，分析了在吸附參數(shù)一定時，孔喉微觀結(jié)構(gòu)對含油量的控制作用。隨平均孔隙直徑增加吸附油量逐漸降低、游離油量逐漸增加；孔隙形狀因子越大，吸附油量越高、游離油量越低(圖8)。不同巖相類型頁巖的變化規(guī)律差異較大，富有機(jī)質(zhì)頁巖的孔隙直徑大于約100 nm時，孔隙內(nèi)以游離油為主，吸附油含量較低，且不受孔隙形態(tài)的影響(圖8a、8c)；含有機(jī)質(zhì)頁巖開始大量富集游離油的孔隙直徑拐點(diǎn)約為30 nm(圖8b、8d)。③孔隙度和含油飽和度對含油量也有直接的影響，隨孔隙度、含油飽和度增加，吸附油、游離油含量均增大。可將孔隙度與含油飽和度的乘積定義為含油孔隙度，并用來反映含油性。本研究中含油飽和度為1，孔隙度與含油飽和度的乘積等于孔隙度，從圖9 上可以看出吸附油量和游離油量均隨孔隙度增加而呈線性增加，驗(yàn)證了公式16 和17 的正確性。因此，含油孔隙度是一個識別頁巖含油性“甜點(diǎn)”有效的參數(shù)。

圖8 平均孔隙直徑對吸附油量、游離油量的影響(a, c為富有機(jī)質(zhì)頁巖；b, d為含有機(jī)質(zhì)頁巖)Fig. 8 Influence of average pore diameter on adsorbed and free oil amounts (a and c, organic-rich shale; b and d, organicbearing shale)

圖9 孔隙度與吸附油量、游離油量的關(guān)系Fig. 9 Relationships of porosity with adsorbed and free oil

表3 東營凹陷沙河街組頁巖參數(shù)Table 3 Shale parameters of the Shahejie Formation in Dongying sag

明確游離油賦存的孔徑下限(即可動油理論孔徑下限)對于頁巖油“甜點(diǎn)”預(yù)測具有重要的指導(dǎo)作用。前人通過實(shí)驗(yàn)確定了不同地區(qū)頁巖油的可動孔徑下限，例如劉惠民等[47]認(rèn)為東營凹陷泥頁巖中游離油賦存孔徑(直徑)下限以及游離油大量富集的孔徑門檻值分別為10 nm和30 nm；王民等[29]認(rèn)為濟(jì)陽坳陷頁巖中游離油賦存孔徑下限為5 nm；黨偉等[39]等認(rèn)為鄂爾多斯盆地延安地區(qū)延長組7段3亞段陸相頁巖孔徑大于3 nm時，孔隙內(nèi)主要為游離油。本次研究提出了游離油賦存孔隙直徑下限(dmin)的理論數(shù)值，應(yīng)為吸附油平均厚度(H)與孔隙形狀因子(F)的乘積(FH)。直徑≤dmin的孔隙內(nèi)全部為吸附油，不含游離油。根據(jù)吸附油平均厚度和F值范圍(2～6)，可知富有機(jī)質(zhì)頁巖的dmin介于3.5～10.5 nm之間，數(shù)值大小與孔隙形態(tài)有關(guān)；含有機(jī)質(zhì)頁巖的dmin介于1.1～3.4 nm之間，低于富有機(jī)質(zhì)頁巖的孔隙直徑下限。

3 頁巖基質(zhì)孔隙油微觀分布規(guī)律

3.1 理論與方法

基于前述吸附油量、游離油量數(shù)學(xué)模型(式16 和17)，可進(jìn)一步推導(dǎo)出頁巖基質(zhì)孔隙中吸附油量占比的數(shù)學(xué)表達(dá)式，即吸附比例方程[51-52]。該方程是用來描述一定溫度、壓力條件下，液體在多孔介質(zhì)孔隙內(nèi)以吸附與游離兩種狀態(tài)共存時，孔隙內(nèi)吸附流體所占總量的質(zhì)量比。該方程不僅適用于頁巖油[28,34]，也適用于頁巖孔隙水[46,51]、煤巖孔隙水[53]，因此具有普遍性意義。吸附比例方程表示如下：

根據(jù)公式可知，頁巖基質(zhì)孔隙中吸附油占比(ra)與吸附油/游離油密度比(ρ1/ρ2)、孔徑/吸附厚度比(dm/H)及孔隙形態(tài)(F)相關(guān)。對于一定的孔隙形態(tài)，隨dm/H增加，ra逐漸減??；隨ρ1/ρ2增加，ra逐漸增加；F值從2 至6 逐漸增加時，ra值逐漸增加；當(dāng)dm=FH時，吸附比例等于1(圖10)。通過吸附比例方程，可認(rèn)識頁巖孔隙系統(tǒng)中吸附油占比及分布特征。圖11 中展示了東營凹陷2 類頁巖中吸附油占比的分布特征，隨孔徑增加，吸附油占比降低、游離油占比增加。富有機(jī)質(zhì)頁巖由于具有較高的吸附油平均厚度及密度，同等尺度孔隙中頁巖油吸附比例要高于含有機(jī)質(zhì)頁巖。根據(jù)吸附油與游離油含量占比隨孔隙直徑的變化趨勢(圖11)，可以看到富有機(jī)質(zhì)頁巖＞100 nm的孔隙中吸附油占比＜30 %、游離油占比＞70 %；含有機(jī)質(zhì)頁巖＞30 nm的孔隙中吸附油占比＜30 %、游離油占比＞70%。通過吸附/游離比例分布特征，可初步判斷在頁巖基質(zhì)連通孔隙系統(tǒng)中吸附油、游離油的賦存孔隙空間。

圖10 頁巖油吸附比例的影響因素(a, d, F=6; b, e, F=4; c, f, F=2)Fig. 10 Factors influencing the adsorption ratio of shale oil (a, d, F=6; b, e, F=4; c, f, F=2)

圖11 不同尺度孔隙中吸附油與游離油比例分布 (a)富有機(jī)質(zhì)頁巖；(b)含有機(jī)質(zhì)頁巖Fig. 11 Distributions of adsorbed and free oil ratios in pores with different size (a) organic-rich shale; (b) organic-bearing shale

為了揭示頁巖基質(zhì)孔隙油的微觀分布，僅僅確定了不同尺度孔隙中吸附油、游離油的占比是不夠的，難以確定出不同尺度孔隙中吸附油、游離油的含量。根據(jù)吸附量、游離量數(shù)學(xué)模型，可以獲得一個巖心中吸附油、游離油的含量，但這些油在復(fù)雜的孔隙系統(tǒng)中是怎么分布的？頁巖基質(zhì)孔隙中吸附油、游離油微觀分布的確定還是一個巨大挑戰(zhàn)。通常，核磁共振T2譜可以反映孔隙中流體的分布，但是難以分解出吸附態(tài)和游離態(tài)流體的分布。吸附比例方程的提出為解決該問題提供了一個新思路。本研究將吸附比例方程與經(jīng)典的核磁共振理論相結(jié)合，建立了不同含油率頁巖中吸附油/游離油微觀分布的定量評價方法。

首先，根據(jù)孔體積、比表面積、孔隙大小和孔隙形態(tài)之間的關(guān)系，將吸附比例方程(式18)改為如下形式[34]：

公式(19)中的孔體積與比表面積參數(shù)在核磁共振橫向弛豫時間(T2)的計算公式中也有出現(xiàn)，這為采用核磁共振技術(shù)快速、無損的檢測頁巖孔隙油分布提供了可能。核磁共振橫向馳豫時間T2的計算公式通常簡寫為：

式中，T2為橫向馳豫時間，ms；ρ為表面弛豫率，nm/ms。

進(jìn)一步，結(jié)合公式(19)和(20)可得[51]：

基于核磁共振T2譜，結(jié)合公式(21)，理論上可計算出核磁共振T2譜上每個T2值對應(yīng)的吸附比例[51]，即為：

式中，T2i為核磁共振T2譜上第i個T2值；rai為T2i值對應(yīng)的吸附比例(0＜rai≤1，當(dāng)計算值rai＞1 時，取rai=1)，分?jǐn)?shù)。

基于飽和油頁巖的核磁共振T2譜，分析不同表面馳豫率時核磁共振計算吸附比例與離心測試吸附比例之差的絕對值(即δ值)，當(dāng)δ值=0 時可得頁巖表面馳豫率。在確定表面馳豫率參數(shù)之后，根據(jù)式(22)可計算頁巖孔隙內(nèi)吸附油微觀分布特征?？傂盘柡臀接托盘柕牟钪悼傻糜坞x油信號分布。此外，根據(jù)頁巖含油量和核磁共振信號幅度之間的線性關(guān)系，亦可將吸附油與游離油的信號分布轉(zhuǎn)換成含油量分布。根據(jù)式(22)可以看出，T2值反映了孔隙尺寸，T2值越大，吸附油占比越小。吸附參數(shù)和表面馳豫率反映了油-巖相互作用，孔徑相同時，密度比(ρ1/ρ2)和吸附厚度H越大，反映油-巖相互作用力越強(qiáng)，吸附油占比越大；表面馳豫率ρ越小，反映了油-巖相互作用力越強(qiáng)，吸附油占比越大。

3.2 實(shí)例分析

對于同一巖相類型的頁巖，吸附參數(shù)相同，但每個巖石的表面馳豫率是不同的。因此，每個巖石具有特定的吸附油、游離油分布特征。以含有機(jī)質(zhì)的L752-2 樣品和富有機(jī)質(zhì)的Y556-1 樣品為例，可得表面馳豫率分別為2.574 nm/ms和1.499 nm/ms。在確定了表面馳豫率之后，采用式(22)，可計算出飽和頁巖中吸附油信號分布，總信號和吸附油信號的差值可得游離油信號分布。圖12a和圖12d展示了兩個頁巖樣品飽和油、吸附油和游離油的分布特征，吸附油主要分布于較小T2值范圍內(nèi)，游離油主要分布于較大T2值范圍，但兩者有一定的重疊區(qū)間。此外，研究發(fā)現(xiàn)吸附油和游離油均有可能出現(xiàn)雙峰分布，并非完全是呈高斯分布的特征。通過離心過程可獲得不同含油率(本文定義為含油量和飽和油量的比值)頁巖，根據(jù)不同含油率頁巖的核磁共振信號幅度可計算出頁巖含油率。在離心過程中，L752-2 樣品中吸附油分布變化較小(圖12e)，Y556-1 樣品中吸附油分布發(fā)生了微調(diào)，由雙峰變成了單峰，但信號幅度變化不大(圖12b)；而游離油在離心過程中含量逐漸減少，譜峰逐漸降低、左移(圖12c和12f)。這說明在離心過程中，吸附油很少被排出，也進(jìn)一步驗(yàn)證了第1.1 小節(jié)中的結(jié)論。

圖12 不同賦存狀態(tài)頁巖油T2 譜分布 (a) Y556-1 樣品吸附油與游離油分布；(b) Y556-1 樣品在不同含油率下的吸附油分布；(c) Y556-1 樣品在不同含油率下的游離油分布；(d) L752-2 樣品吸附油與游離油分布；(e) L752-2 樣品在不同含油率下的吸附油分布；(f) L752-2 樣品在不同含油率下的游離油分布Fig. 12 T2 spectral distributions of shale oil in different occurrence states (a) adsorbed and free oil distributions of sample Y556-1; (b) adsorbed oil distribution of sample Y556-1 under different oil-bearing rate; (c) free oil distribution of sample Y556-1 under different oil-bearing rate; (d) adsorbed and free oil distributions of sample L752-2; (e) adsorbed oil distribution of sample L752-2 under different oil-bearing rate; (f) free oil distribution of sample L752-2 under different oil-bearing rate

根據(jù)孔徑(d)和T2值之間的線性關(guān)系(d=FρT2)，F(xiàn)取平均值4，可將T2譜分布轉(zhuǎn)換為孔徑分布。如圖13所示，Y556-1 樣品吸附油主要分布于＜100 nm的孔隙中、游離油主要分布于＞7 nm的孔隙中；L752-2 樣品吸附油主要分布于＜30 nm的孔隙中、游離油主要分布于＞3 nm的孔隙中。游離油分布主峰孔徑明顯大于吸附油分布主峰孔徑。相比于L752-2 樣品，Y556-1 樣品具有更高的信號幅度，因此含油量較高；Y556-1樣品含油分布的主峰孔徑為20 nm，高于L752-2 樣品的主峰孔徑2 nm，因此L752-2 樣品具有更高的吸附比例。

圖13 不同賦存狀態(tài)頁巖油賦存空間分布(a) Y556-1 樣品吸附油與游離油分布；(b) Y556-1 樣品在不同含油率下的吸附油分布；(c) Y556-1 樣品在不同含油率下的游離油分布；(d) L752-2 樣品吸附油與游離油分布；(e) L752-2 樣品在不同含油率下的吸附油分布；(f) L752-2 樣品在不同含油率下的游離油分布Fig. 13 Spatial distributions of shale oil in different occurrence states (a) adsorbed and free oil distributions of sample Y556-1; (b) adsorbed oil distribution of sample Y556-1 under different oil-bearing rate; (c) free oil distribution of sample Y556-1 under different oil-bearing rate; (d) adsorbed and free oil distributions of sample L752-2; (e) adsorbed oil distribution of sample L752-2 under different oil-bearing rate; (f) free oil distribution of sample L752-2 under different oil-bearing rate

4 頁巖基質(zhì)孔隙油可動性表征

頁巖油可動性與頁巖油賦存特征(賦存狀態(tài)、不同狀態(tài)含量及比例)及賦存的儲集空間(孔喉大小及分布)密不可分，同時也受儲層流體條件(油組分、溫度、壓力)的影響，因此頁巖油可動性是一個綜合的儲層屬性。趙文智院士等認(rèn)為可從熱成熟度、氣油比、原油密度、烴組分構(gòu)成、保存條件等方面評價滯留烴流動性[15,54]。也有學(xué)者采用可動系數(shù)(輕質(zhì)游離油量與滯留油量之比)[20]、可動流體飽和度[36,40,55]、OSI[17]、游離油量及可動油量[28]等參數(shù)來反映頁巖油可動性，尚未形成統(tǒng)一的認(rèn)識。筆者認(rèn)為頁巖油可動性表征存在以下難點(diǎn)：①可動性內(nèi)涵不清晰，沒有形成行業(yè)統(tǒng)一的定量表征參數(shù)，如何量化可動性存在較大爭議；②可動性不是一個獨(dú)立的頁巖油屬性，它依賴于頁巖油的賦存和頁巖孔喉結(jié)構(gòu)，如何體現(xiàn)它們?nèi)咧g的內(nèi)在聯(lián)系，難度較大。

本研究提出頁巖油可動性至少包括2 個方面的內(nèi)涵：可動潛力和可動能力?？蓜訚摿κ侵冈谝欢ㄩ_采條件下頁巖基質(zhì)孔隙中有多少頁巖油可以流動，體現(xiàn)了頁巖油現(xiàn)實(shí)可采資源，可通過最大可動量來表示。最大可動量即為游離量，反映了頁巖油賦存特征，其含量越大，頁巖油可采資源潛力越大，越有利于頁巖油開發(fā)?？蓜幽芰κ侵疙搸r油排出基質(zhì)孔隙的難易程度，體現(xiàn)了頁巖油開采的難度，可通過中值壓力差來表示。中值壓力差反映了頁巖孔喉微觀結(jié)構(gòu)，其值越低，代表頁巖基質(zhì)孔隙中的油越容易排出，反映頁巖油開采過程中需要較小的壓力降即可采出大量的頁巖油。若頁巖中值壓力差越高，頁巖油排出難度越大，需要較高的壓力降才能采出頁巖油。較大的壓力降可能會導(dǎo)致儲層較強(qiáng)的應(yīng)力敏感負(fù)效應(yīng)，對儲層造成較大的傷害。綜上，頁巖油可動性好應(yīng)該是可動潛力大、可動能力強(qiáng)，即具有較高的游離油含量、較低的中值壓力差。如圖14 所示，展示了頁巖油在外力(如離心力)作用下排出基質(zhì)孔隙的過程。曲線A和曲線B具有相同的可動潛力，即最大可動量相同(QfA=QfB)，但曲線B的可動能力要差于曲線A(△PLB＞△PLA)。曲線B和曲線C相比，兩者有相同的可動能力，即中值壓力差相同(△PLB=△PLC)，但曲線C的可動潛力差于曲線B(QfB＜QfA)。綜合來看，曲線A反映了最好的可動性，既具有較高的可動潛力，又具有較大的可動能力。

圖14 可動油含量與離心力關(guān)系模式圖Fig. 14 Model diagram of the relationship between mobile oil amount and centrifugal force

在第1.1 小節(jié)中已介紹了基于離心實(shí)驗(yàn)獲得游離油量(Qf)和中值壓力差(△PL)的相關(guān)理論和分析方法。根據(jù)游離油量和中值壓力差交匯圖，游離油量越高、中值壓力差越小的區(qū)域代表頁巖油具有較好的可動性。然而，實(shí)驗(yàn)結(jié)果顯示游離油量和中值壓力差呈一定正相關(guān)性(圖15)，也就是說游離油量高(可動潛力大)的頁巖中值壓力差也大(可動能力小)、游離油量低(可動潛力小)的頁巖中值壓力差反而小(可動能力大)。這為判斷頁巖油可動性增加了難度。為了定量化頁巖油可動性，本研究提出了可動性指數(shù)(Im)，其數(shù)值上等于游離油量除以中值壓力差，依據(jù)該參數(shù)來反映頁巖基質(zhì)孔隙油的可動性。

圖15 頁巖基質(zhì)孔隙油可動性評價圖版Fig. 15 Evaluation chart on the oil mobility in shale matrix pores

評價頁巖油可動性影響因素也即是考察游離油量和中值壓力差的影響因素。游離油量的影響因素前文已經(jīng)論述；影響中值壓力差的因素主要為頁巖孔喉大小分布(孔徑分布)的復(fù)雜性，體現(xiàn)于孔隙-喉道(孔-喉)的配置關(guān)系。在不同尺度連通孔隙組成的系統(tǒng)中，就兩個連通孔隙比較而言，較小孔可視為喉道，較大孔則為孔隙。如圖16a所示的3 個頁巖樣品的孔徑分布特征，主峰位于10～30 nm，但從孔徑范圍來看，B172-3 樣品孔徑分布較為集中，幾乎不發(fā)育30 nm以上孔隙；LX884-2S樣品發(fā)育30～100 nm的孔隙；Y556-5S樣品發(fā)育30～300 nm的孔隙，且在100～200 nm之間出現(xiàn)了次峰和少量3～5 μm的較大孔隙。因此，頁巖孔徑分布復(fù)雜性由大到小依次為Y556-5S樣品、LX884-2S樣品和B172-3 樣品。相應(yīng)的，頁巖油可動能力由高到低依次出現(xiàn)于B172-3 樣品、LX884-2S樣品和Y556-5S樣品。此外，孔徑呈多峰分布的頁巖中孔隙油的可動能力差于單峰分布的頁巖，如L76-1 樣品和LX884-2S樣品。而且，頁巖孔徑分布越不連續(xù)，孔隙油可動能力越差，如L76-1樣品和Y556-1 樣品(圖16b)。這些現(xiàn)象說明了頁巖孔喉微觀結(jié)構(gòu)的復(fù)雜性對頁巖油可動能力有較大影響，其本質(zhì)是孔-喉配置關(guān)系發(fā)揮了作用。對于孔徑較為均一的頁巖，其孔-喉配置關(guān)系較好，具有較低的孔喉比(孔隙和喉道的尺寸比)，有利于頁巖油流動；而孔徑分布較為復(fù)雜的頁巖，其孔-喉配置關(guān)系相對較差，孔喉比較高，存在孔隙屏蔽效應(yīng)，部分可動空間受阻，不利于頁巖油流動(圖16c)。從頁巖吸附油、游離油微觀分布角度也可論證這一觀點(diǎn)。上文中展示了Y556-1 樣品和L752-2 樣品中吸附油、游離油微觀分布。Y556-1 樣品具有較高的游離油含量(19.7629 mg/g)，L752-2 樣品具有較低的游離油含量(5.9952 mg/g)，說明Y556-1 樣品具有較大的可動潛力；但Y556-1 樣品的中值壓力差(2.08 MPa)大于L752-2 樣品(1.35 MPa)，可動能力較差。從分布譜圖上可以看出，Y556-1 樣品中游離油賦存空間大于L752-2 樣品，但是其分布呈現(xiàn)孤立雙峰型，相比于L752-2 樣品的連續(xù)單峰型，Y556-1 樣品的孔徑分布更復(fù)雜，這也是該樣品可動能力較差的主要原因。

圖16 頁巖高壓壓汞的孔徑分布(a，b)和孔-喉配置模式(c)Fig. 16 Pore size distribution based on high pressure mercury injection test (a, b) and pattern diagram of pore-throat configuration (c)

統(tǒng)計結(jié)果顯示吸附比例與中值壓力差呈現(xiàn)負(fù)相關(guān)性(圖17a)。如前文所述，吸附比例和平均孔隙直徑呈現(xiàn)非線性負(fù)相關(guān)性，即孔徑越大、吸附比例越低。也就是說，平均孔徑較大的頁巖中由于出現(xiàn)了較多的大孔隙，導(dǎo)致較差的孔-喉配置關(guān)系，在大孔隙中賦存的游離油排出受阻，導(dǎo)致頁巖油可動能力降低。因此，對于微納尺度孔喉系統(tǒng)發(fā)育的頁巖而言，吸附比例較低的頁巖，大孔隙占比較高，平均孔徑較大，但是孔-喉配置關(guān)系較差，不利于頁巖油流動；相反，吸附比例較高的頁巖，小孔隙占比較高，平均孔徑較小，孔-喉配置關(guān)系較好，有利于頁巖油流動。此外，統(tǒng)計結(jié)果也顯示隨頁巖油吸附比例增加，頁巖中游離油量逐漸降低(圖17b)。這就導(dǎo)致隨吸附比例增加，一方面頁巖油可動能力增加，另一方面頁巖油可動潛力降低。這是兩個完全相反的趨勢。根據(jù)吸附比例與可動性指數(shù)的關(guān)系(圖17c)，可以看出在吸附比例約為60 %時，可動性指數(shù)出現(xiàn)極大值。因此可以推斷當(dāng)頁巖油吸附比例＜60 %時，頁巖油可動能力增強(qiáng)起主導(dǎo)作用，使得頁巖油可動性指數(shù)增加；當(dāng)頁巖油吸附比例＞60 %時，頁巖油可動潛力降低起主導(dǎo)作用，使得頁巖油可動性指數(shù)降低。根據(jù)這個結(jié)論反觀吸附油量-游離油量圖版(圖3)，可知孔隙油吸/游比約等于1.5 時可動性較好。相對而言，富有機(jī)質(zhì)頁巖中孔隙油可動性指數(shù)(平均6.24 mg·g-1·MPa-1)高于含有機(jī)質(zhì)頁巖(平均5.20 mg·g-1·MPa-1)，富有機(jī)質(zhì)頁巖中孔隙油可動性相對較強(qiáng)。

圖17 頁巖油吸附比例與中值壓力差(a)、游離油量(b)、可動性指數(shù)(c)的關(guān)系Fig. 17 Relationships between adsorption ratio of shale oil with median pressure difference (a), free oil amount (b) and mobility index (c)

5 頁巖油儲層屬性耦合關(guān)系

頁巖油各儲層屬性(包括頁巖儲集性、含油性及頁巖油可動性)不是獨(dú)立存在的，而是具有明確的內(nèi)在聯(lián)系[56]，通過為儲層屬性賦予特征參數(shù)后，可以從數(shù)學(xué)角度描述它們之間的內(nèi)在聯(lián)系(圖18)，以更深層次的理解頁巖油儲層屬性耦合關(guān)系。在儲層屬性耦合關(guān)系中，油-巖相互作用扮演著至關(guān)重要的作用，是銜接頁巖儲集性(孔喉微觀結(jié)構(gòu)、宏觀物性)、含油性和頁巖油可動性的紐帶。對于一個巖石，其由巖石骨架和儲集空間耦合構(gòu)成，頁巖油賦存于儲集空間，若無油-巖相互作用發(fā)生，則所有孔隙油都處于單一游離態(tài)。頁巖含油量通過孔隙度、含油飽和度、游離油密度等參數(shù)即可計算獲得。此時，頁巖油可動性僅受巖石孔喉微觀結(jié)構(gòu)的影響。若存在油-巖相互作用，頁巖儲集空間中吸附油與游離油共存，在孔喉微觀結(jié)構(gòu)(孔喉大小、分布和形態(tài))的疊加影響下，不同相態(tài)孔隙油賦存量及其分布更加復(fù)雜，頁巖油可動性影響機(jī)制也更加復(fù)雜，將受到頁巖孔喉微觀結(jié)構(gòu)和油賦存狀態(tài)的綜合影響。因此，總結(jié)成一句話，即為頁巖儲集空間是一個由不同尺度孔隙連接構(gòu)成的三維網(wǎng)絡(luò)，在孔隙網(wǎng)絡(luò)中發(fā)生油-巖相互作用形成不同相態(tài)孔隙油(吸附油、游離油)，吸附油與游離油在三維孔隙網(wǎng)絡(luò)中的含量、比例和分布，制約了孔隙油的可動性及可采性。由于巖性/巖相與巖石骨架組成有關(guān)，表現(xiàn)出不同的儲集空間發(fā)育特征和油-巖相互作用，因此具有不同的含油性和頁巖油可動性。盡管不同盆地/區(qū)塊的頁巖油儲集層的地質(zhì)條件(如巖性、巖相、成熟度、溫壓系統(tǒng)、流體性質(zhì)、巖石組成、儲集空間等)千差萬別，但控制頁巖油賦存和可動的本質(zhì)相同。通過從數(shù)學(xué)上建立儲層屬性之間的耦合關(guān)系，不僅具有預(yù)測功能，也便于不同盆地之間差異性對比，對于尋找頁巖油富集、可動的“甜點(diǎn)”位置具有積極作用，同時為后期的頁巖油提高采收率提供精細(xì)的儲層描述。具體表現(xiàn)在3 個方面：其一，可明確頁巖基質(zhì)孔隙中游離油含量及其分布空間，落實(shí)現(xiàn)實(shí)可采的頁巖油資源；其二，可明確頁巖基質(zhì)孔隙中吸附油含量及其分布空間，鎖定頁巖油開采后期挖潛目標(biāo)；其三，可揭示頁巖油可動性主控因素及機(jī)制，助力頁巖油高效開采。

圖18 頁巖油儲層評價參數(shù)耦合關(guān)系Fig. 18 Coupling relationships among shale oil reservoir evaluation parameters

6 結(jié)論

本文構(gòu)建了頁巖油儲層屬性(頁巖含油性、儲集性及頁巖油可動性)之間的耦合關(guān)系，建立了頁巖基質(zhì)孔隙中不同相態(tài)油(吸附油、游離油)的含量、比例、分布及可動性數(shù)學(xué)評價模型和實(shí)驗(yàn)分析方法?；陲柡?離心-核磁共振聯(lián)合實(shí)驗(yàn)，以渤海灣盆地濟(jì)陽坳陷東營凹陷沙河街組沙三下亞段和沙四上亞段頁巖油儲層為例，分析了不同巖相類型頁巖基質(zhì)孔隙中輕質(zhì)油(正十二烷)的吸附量、游離量及其比例、分布、可動性特征，主要結(jié)論如下：

(1)頁巖吸附油量介于5～30 mg/g、游離油量介于5～25 mg/g，吸/游比(吸附油與游離油質(zhì)量比)主體介于1～2。富有機(jī)質(zhì)頁巖的吸附油、游離油含量普遍大于含有機(jī)質(zhì)頁巖，且富有機(jī)質(zhì)頁巖中吸/游比略高于含有機(jī)質(zhì)頁巖。富有機(jī)質(zhì)頁巖中吸附油主要賦存于有機(jī)質(zhì)孔內(nèi)，含有機(jī)質(zhì)頁巖中吸附油賦存于有機(jī)質(zhì)孔和黏土孔內(nèi)；富有機(jī)質(zhì)頁巖中游離油主要賦存于方解石孔、有機(jī)質(zhì)孔內(nèi)，含有機(jī)質(zhì)頁巖中游離油主要賦存于有機(jī)質(zhì)孔、石英礦物孔內(nèi)。

(2)富有機(jī)質(zhì)頁巖吸附油平均密度為0.8331 g/cm3、吸附油平均厚度為1.7475 nm，含有機(jī)質(zhì)頁巖吸附油平均密度為0.8067 g/cm3、吸附油平均厚度為0.5734 nm，顯示富有機(jī)質(zhì)頁巖具有更強(qiáng)的油-巖相互作用。富有機(jī)質(zhì)頁巖的孔隙直徑大于約100 nm時，孔隙內(nèi)以游離油為主(＞70 %)，吸附油含量較低(＜30 %)；含有機(jī)質(zhì)頁巖主要富集游離油的孔隙直徑拐點(diǎn)約為30 nm。富有機(jī)質(zhì)頁巖的游離油賦存孔隙直徑下限(dmin)介于3.5～10.5 nm，含有機(jī)質(zhì)頁巖的dmin介于1.1～3.4 nm，數(shù)值大小與孔隙形態(tài)有關(guān)。

(3)頁巖基質(zhì)孔隙油可動性受控于油微觀賦存和頁巖孔徑分布復(fù)雜性，通過可動性指數(shù)(即游離油量與中值壓力差之比)可有效表征頁巖基質(zhì)孔隙油的可動性。富有機(jī)質(zhì)頁巖的孔隙油可動性指數(shù)(平均6.24 mg·g-1·MPa-1)高于含有機(jī)質(zhì)頁巖(平均5.20 mg·g-1·MPa-1)，反映富有機(jī)質(zhì)頁巖中孔隙油可動性相對較強(qiáng)。孔隙油吸/游比約為1.5 時可動性普遍較好。

下一步建議在以下幾個方面加強(qiáng)理論研究：①確定原位儲層條件下頁巖基質(zhì)孔隙內(nèi)吸附油、游離油含量及分布，重點(diǎn)攻關(guān)吸附油-游離油動態(tài)轉(zhuǎn)化規(guī)律；②明確頁巖基質(zhì)孔隙中不同相態(tài)頁巖油的微觀流動特征及機(jī)制；③通過頁巖巖石表面改性、改變儲層溫壓環(huán)境等，促使吸附油向游離油轉(zhuǎn)化，提高頁巖油可動潛力；④通過改善頁巖孔喉微觀結(jié)構(gòu)，降低孔喉比，增強(qiáng)頁巖油可動能力。