朱日房 ，張林曄 ，李 政 ，王 茹 ，張守春 ，張 蕾

（1.中國石化勝利油田分公司勘探開發(fā)研究院，山東東營257015；2.國家能源局頁巖油研發(fā)中心勝利油田分室，山東東營257015；3.中國石化頁巖油氣勘探開發(fā)重點(diǎn)實(shí)驗(yàn)室，山東東營257015；4.西北大學(xué)大陸動(dòng)力學(xué)國家重點(diǎn)實(shí)驗(yàn)室，陜西西安710069）

北美頁巖油勘探開發(fā)成果極大地激勵(lì)了中國頁巖油的勘探研究，中國東部新生代斷陷湖盆發(fā)育多套陸相富有機(jī)質(zhì)泥頁巖，并有多口探井發(fā)現(xiàn)了頁巖油氣流，顯示出良好的勘探前景［1-7］。近10 a來，許多研究者在頁巖的含油性［8-13］、儲(chǔ)集性及流動(dòng)性方面進(jìn)行了探索［14-18］。頁巖油資源潛力評(píng)價(jià)從最初以滯留油為主的原地資源量評(píng)價(jià)發(fā)展為以游離油為主要賦存方式的游離油量評(píng)價(jià)［19-21］。由于頁巖油資源評(píng)價(jià)技術(shù)的不斷進(jìn)步、方法不斷完善，以及對(duì)頁巖中油氣賦存機(jī)理認(rèn)識(shí)的加深，必然對(duì)頁巖油的資源潛力評(píng)價(jià)提出更高的要求。頁巖油資源潛力評(píng)價(jià)研究仍然是目前頁巖油勘探實(shí)踐中重要的研究?jī)?nèi)容之一。為此，根據(jù)近幾年頁巖含油性定量分析等方面的進(jìn)展，以東營凹陷沙三段下亞段為例，探討陸相斷陷盆地頁巖油資源潛力評(píng)價(jià)。

1 地質(zhì)概況

東營凹陷位于濟(jì)陽坳陷東南部，其東接青坨子凸起，南部與魯西隆起、廣饒凸起呈超覆接觸，西與惠民凹陷毗鄰，北與陳家莊凸起和濱縣凸起相鄰，東西長(zhǎng)度約為90 km，南北寬度約為62 km，總面積為5 850 km2（圖1）。從大地構(gòu)造區(qū)劃分上，東營凹陷屬于中國東部渤海灣盆地的一個(gè)次級(jí)構(gòu)造單元，是在古生界基巖古地形背景上，經(jīng)構(gòu)造運(yùn)動(dòng)發(fā)育形成的中、新生代斷陷-坳陷盆地；由于基巖起伏以及中央隆起帶的分隔，凹陷內(nèi)又分為4個(gè)次級(jí)洼陷，它們既是盆地的沉積中心，也是富含有機(jī)質(zhì)泥頁巖發(fā)育區(qū)。凹陷內(nèi)主要發(fā)育沙四段上亞段和沙三段下亞段2套成熟生烴的泥頁巖［22］，在110口井中的泥頁巖段發(fā)現(xiàn)油氣顯示，試油37口井，14口井獲得了工業(yè)油氣流，其中對(duì)河54井沙三段下亞段泥頁巖發(fā)育段試油，日產(chǎn)油量為 91.4 t/d［23］，累積產(chǎn)油量為27 896 t，顯示出良好的頁巖油勘探前景。

圖1 東營凹陷區(qū)域構(gòu)造位置Fig.1 Regional tectonic location of Dongying Sag

2 頁巖油資源潛力評(píng)價(jià)方法

2.1 頁巖油資源量計(jì)算方法

頁巖油資源潛力評(píng)價(jià)方法主要分為動(dòng)態(tài)法和靜態(tài)法2大類。動(dòng)態(tài)法主要根據(jù)頁巖油在開發(fā)過程中的動(dòng)態(tài)資料，通過一定的數(shù)學(xué)模型定量計(jì)算資源量。靜態(tài)法主要應(yīng)用靜態(tài)參數(shù)計(jì)算資源量，按計(jì)算方式的不同細(xì)分為統(tǒng)計(jì)法、類比法和成因法3種。統(tǒng)計(jì)法需要大量典型樣本，適用于中-高勘探程度階段；類比法需要有可類比的刻度區(qū)，一般適用于低勘探程度區(qū)；成因法重視物質(zhì)平衡原理的應(yīng)用，更多是確定性評(píng)價(jià)，適用于盆地勘探的各個(gè)階段。目前由于缺少系統(tǒng)的頁巖油開采動(dòng)態(tài)資料和典型樣本，難以應(yīng)用動(dòng)態(tài)法和統(tǒng)計(jì)法，同時(shí)中國東部陸相盆地的地質(zhì)情況與北美具有較大差別，難以直接將北美成熟頁巖油開采區(qū)作為中國頁巖油資源潛力評(píng)價(jià)類比刻度區(qū)。因此，成因法是中國頁巖油資源潛力評(píng)價(jià)最常用的方法。

體積法是成因法中最常用的一種方法，其計(jì)算頁巖油資源量的原理簡(jiǎn)單明確，將泥頁巖的質(zhì)量（頁巖體積與密度乘積）與泥頁巖中頁巖油含量相乘即可獲得頁巖油資源量，其計(jì)算式為：

計(jì)算頁巖油資源量所需的參數(shù)主要包括泥頁巖的面積、有效厚度、密度和頁巖油含量。泥頁巖的面積和有效厚度主要通過頁巖厚度等值線圖獲取，密度由實(shí)測(cè)密度剖面獲取。由于頁巖油主要以游離油的形式賦存于泥頁巖中，因此，泥頁巖中的頁巖油含量用游離油量來近似代替，而游離油量的定量分析是頁巖油資源量計(jì)算中的難點(diǎn)之一。

2.2 滯留油量與吸附油量計(jì)算方法

游離油量的定量分析主要有2種方法，一種是孔隙度和含油飽和度法，另一種是地化參數(shù)法。對(duì)于常規(guī)砂巖儲(chǔ)層來說，孔隙度和含油飽和度法表征游離油量簡(jiǎn)單明確，然而對(duì)于非均質(zhì)性強(qiáng)、孔隙結(jié)構(gòu)復(fù)雜、以納米孔為主的泥頁巖來說，其孔隙度和含油飽和度的測(cè)量非常困難，數(shù)據(jù)點(diǎn)少，難以用于宏觀的頁巖油資源量計(jì)算；因此，地化參數(shù)法成為最常用的游離油量定量分析方法。

應(yīng)用地化參數(shù)法很難直接獲取游離油量和吸附油量，但能夠確定頁巖中的滯留油量和巖石對(duì)油的吸附潛量。在飽和吸附時(shí)，頁巖對(duì)油的吸附潛量近似于頁巖的吸附量，頁巖中的滯留油量與頁巖的吸附潛量之間的差值即為頁巖的游離油量。

表征頁巖滯留油量最常用的參數(shù)為氯仿瀝青“A”含量和熱解參數(shù)，相應(yīng)的計(jì)算式分別為：

頁巖對(duì)油的吸附能力與巖石的結(jié)構(gòu)、礦物組分和有機(jī)質(zhì)豐度、類型及演化程度等因素密切相關(guān)［22-28］。不同的礦物組分對(duì)油的吸附能力明顯不同，有機(jī)質(zhì)對(duì)油的吸附能力遠(yuǎn)大于無機(jī)礦物對(duì)油的吸附能力［27-28］，不同類型的有機(jī)質(zhì)以及熱演化程度不同的有機(jī)質(zhì)對(duì)油的吸附能力也有明顯區(qū)別［24］?？紫督Y(jié)構(gòu)影響無機(jī)礦物與油的接觸面積，從而影響無機(jī)礦物對(duì)油的吸附能力［25］。因此，可以分別計(jì)算干酪根及不同孔隙度條件下各礦物組分的吸附潛量來表征整個(gè)巖石對(duì)油的吸附能力，同時(shí)，在計(jì)算無機(jī)礦物吸附潛量時(shí)加入與孔隙結(jié)構(gòu)有關(guān)的泥頁巖比表面積相關(guān)系數(shù)。頁巖對(duì)油吸附潛量的計(jì)算式為：

其中

在獲取滯留油量和頁巖對(duì)油的吸附潛量后，就可以計(jì)算頁巖的吸附油量和游離油量。當(dāng)頁巖中滯留油量小于頁巖對(duì)油的吸附潛量時(shí)，為欠飽和吸附，吸附油量就等于滯留油量，此時(shí)沒有游離油；當(dāng)頁巖中滯留油量大于或等于頁巖對(duì)油的吸附潛量時(shí)，為飽和吸附，吸附油量就等于頁巖對(duì)油吸附潛量，游離油量則為滯留油量與吸附油量的差值。

2.3 可動(dòng)油資源量計(jì)算方法

可動(dòng)油主要指在地層條件下，采用一定的工程措施能夠采出的原油?？蓜?dòng)油資源量對(duì)頁巖油的勘探開發(fā)具有重要意義，是頁巖油資源潛力最直接的衡量指標(biāo)?？蓜?dòng)油資源量表達(dá)式為：

頁巖油可動(dòng)率與泥頁巖的骨架、原油、地層水性質(zhì)等有關(guān)，還與地層流體壓力有關(guān)，巖石骨架、原油和水的壓縮性越大、地層流體壓力越高、能量越大，頁巖油可動(dòng)率就越高。

3 東營凹陷沙三段下亞段頁巖油資源潛力評(píng)價(jià)

3.1 頁巖油地質(zhì)條件

東營凹陷沙三段下亞段為半咸水沉積，主要為灰色、灰褐色頁巖與深灰色泥巖組成的不等厚互層，其厚度一般為100～300 m。地化分析資料表明，研究區(qū)沙三段下亞段總有機(jī)碳含量（TOC）分布范圍較寬，為0.5%～18.6%，平均值為4.9%；氯仿瀝青“A”含量為0.11%～2.94%，平均值為0.92%；有機(jī)質(zhì)類型以Ⅰ型干酪根為主，少部分為Ⅱ1型干酪根，埋深為1 500～4 000 m，各洼陷內(nèi)基本上處于成熟生烴階段。沙三段下亞段泥頁巖有機(jī)質(zhì)豐度高、類型好、主體處于成熟生烴階段，是一套生烴潛力大、分布穩(wěn)定的主力烴源巖。

3.2 頁巖油資源潛力評(píng)價(jià)參數(shù)

3.2.1 頁巖含油性分級(jí)

后綴法：“詞根+后綴”。 如：-meter“計(jì)量?jī)x器”，可后綴派生出 barometer氣壓計(jì),telemeter測(cè)距儀,spectrometer分光儀,oscilloscope示波器等多個(gè)新詞。

湖相泥頁巖的含油性在平面上和縱向上均具有很強(qiáng)的非均質(zhì)性，頁巖含油性分級(jí)描述能有效地消除頁巖非均質(zhì)性的影響［29-31］。按照頁巖含油性參數(shù)（TOC、氯仿瀝青“A”含量和S1）及其相關(guān)性，結(jié)合東營凹陷沙三段下亞段烴源巖具體特征，將頁巖含油性分為Ⅰ，Ⅱ和Ⅲ共3個(gè)級(jí)別。具體劃分方案見表1。

表1 東營凹陷頁巖含油性分級(jí)劃分方案Table1 Division scheme for the oil-bearing abundance of shale in Dongying Sag

采用系統(tǒng)取心井與測(cè)井曲線相結(jié)合，應(yīng)用改進(jìn)的Δlog R模型［32］建立TOC、氯仿瀝青“A”含量和S1的解釋模型，進(jìn)而計(jì)算頁巖含油性參數(shù)。根據(jù)計(jì)算結(jié)果和表1劃分方案，可以對(duì)研究區(qū)各單井的頁巖含油性進(jìn)行分級(jí)、統(tǒng)計(jì)各級(jí)厚度及含油性參數(shù)，編制各級(jí)的厚度等值線和地化參數(shù)等值線（圖2，圖3）。從圖2可以看出，東營凹陷Ⅰ級(jí)頁巖的厚度中心為各洼陷中心，厚度均在150 m以上，其中利津、民豐達(dá)到200 m；各含油性參數(shù)也具有類似的分布特征，利津洼陷中心處TOC、氯仿瀝青“A”含量和S1值分別為6.0%，1.2%和5.0 mg/g，博興洼陷中心處分別為5.0%，1.0%和4.0 mg/g，民豐洼陷中心處分別為4.0%，1.0%和4.0 mg/g。Ⅱ級(jí)頁巖與Ⅰ級(jí)頁巖具有相似的特征，牛莊、利津洼陷中心頁巖厚度達(dá)到了200 m，博興、民豐洼陷為100 m，各洼陷中心的TOC、氯仿瀝青“A”含量和S1值分別為2.0%，0.4%和2.0 mg/g。

圖2 東營凹陷沙三段下亞段Ⅰ級(jí)頁巖含油性參數(shù)等值線分布Fig.2 Contour map of oil content parameters of GradeⅠshale of Lower Es3Member in Dongying Sag

3.2.2 頁巖密度

根據(jù)實(shí)測(cè)的密度資料（圖4）發(fā)現(xiàn)，東營凹陷泥頁巖密度表現(xiàn)出隨埋深增加、壓實(shí)程度增強(qiáng)而增大的特征，但在不同埋深段密度隨埋深的變化曲線呈不同的斜率。在同一埋深下，烴源巖密度與TOC呈一定程度的負(fù)相關(guān)關(guān)系，即TOC越高，密度越小。由此，泥頁巖的密度可以由不同TOC在不同埋深下的密度-埋深關(guān)系曲線（圖4）來獲取。

3.2.3 輕烴恢復(fù)系數(shù)及重?zé)N比例系數(shù)

輕烴恢復(fù)系數(shù)及重?zé)N比例系數(shù)采用文獻(xiàn)［8］的資料，各系數(shù)隨埋深的增加而增大。在主要的生油窗范圍內(nèi)，東營凹陷氯仿瀝青“A”恢復(fù)系數(shù)為1.0～1.5，熱解輕烴散失系數(shù)為0.1～0.4，熱解S2中重?zé)N比例系數(shù)則為 0.1～0.8［8］。

根據(jù)東營凹陷沙三段下亞段實(shí)際樣品的地化資料、巖礦資料和物性資料，以及不同組分對(duì)油的吸附能力，應(yīng)用（4）式就可計(jì)算不同樣品對(duì)油的吸附潛量。在生油窗范圍內(nèi)，東營凹陷沙三段下亞段干酪根對(duì)油的單位TOC吸附潛量為50～120 mg/g［28］，無機(jī)礦物中粘土礦物、石英、碳酸鹽和黃鐵礦對(duì)油的吸附潛量分別為18.0，3.0，1.8和10.0 mg/g［25］，泥頁巖的初始孔隙度為 50%，現(xiàn)今孔隙度為5%～15%。東營凹陷沙三段下亞段頁巖油吸附潛量計(jì)算結(jié)果（圖5）表明，東營凹陷沙三段下亞段泥頁巖在主要生烴初始階段頁巖油單位TOC吸附潛量為200～300 mg/g，略大于國外的 200 mg/g［29］，與應(yīng)用A/TOC確定的230 mg/g結(jié)果較一致［19，28］。在埋深為4 000 m左右，頁巖油單位TOC吸附潛量約為100～150 mg/g。

圖3 東營凹陷沙三段下亞段Ⅱ級(jí)頁巖含油性參數(shù)等值線分布Fig.3 Contour map of oil content parameters of GradeⅡshale of Lower Es3Member in Dongying Sag

圖4 東營凹陷泥頁巖密度隨深度變化曲線Fig.4 Variation curves of shale density with depth in Dongying Sag

圖5 東營凹陷沙三段下亞段頁巖油吸附潛量演化剖面Fig.5 Evolution of adsorption oil capacity of shale in Lower Es3Member in Dongying Sag

3.3 頁巖油資源潛力及分布特征

在確定不同深度下的氯仿瀝青“A”輕烴恢復(fù)系數(shù)、熱解輕烴散失系數(shù)、熱解S2中重?zé)N比例系數(shù)、頁巖油吸附潛量以及泥頁巖密度等參數(shù)后，根據(jù)研究區(qū)不同級(jí)別的泥頁巖總有機(jī)碳含量、氯仿瀝青“A”含量、S1、有效厚度、埋深等參數(shù)，計(jì)算濟(jì)陽坳陷不同級(jí)別泥頁巖的滯留油量和游離油量。計(jì)算過程中，為求精確，將不同級(jí)別泥頁巖分布區(qū)在平面上均分為500 m×500 m的網(wǎng)格區(qū)，分別應(yīng)用氯仿瀝青“A”參數(shù)法和熱解參數(shù)法計(jì)算各網(wǎng)格區(qū)的滯留油量和游離油量，所有網(wǎng)格的滯留油量之和為總滯留油量，游離油量之和為總頁巖油量即頁巖油資源量。

研究結(jié)果表明，2種方法的計(jì)算結(jié)果（表2）非常接近，采用氯仿瀝青“A”參數(shù)法和熱解參數(shù)法，東營凹陷沙三段下亞段總滯留油量分別為65.44×108和64.26×108t；頁巖油資源量分別為11.39×108和11.57×108t。其中利津洼陷頁巖油最豐富，資源量超過6.00×108t，占比超過50%，其次為牛莊洼陷。應(yīng)用張林曄等建立的“彈性驅(qū)動(dòng)+溶解氣驅(qū)動(dòng)可動(dòng)油率曲線”［28］作為頁巖油可動(dòng)率來計(jì)算頁巖油可動(dòng)資源量，得到東營凹陷沙三段下亞段頁巖油可動(dòng)資源量為1.50×108t。

表2 東營凹陷沙三段下亞段頁巖油資源量計(jì)算結(jié)果Table2 Calculation results of shale oil resources of Lower Es3Member in Dongying Sag

從東營凹陷沙三段下亞段頁巖油資源豐度的平面分布特征看（圖6），其頁巖油資源量主要分布在各洼陷帶，從洼陷中心向邊部逐漸減少，洼陷中心的資源豐度以利津洼陷最高，可達(dá)250×104t/km2，其次是牛莊洼陷和博興洼陷，分別為150×104和100×104t/km2，民豐洼陷最低?？v向上，埋深為3 000～3 500 m最高，為7.01×108t，明顯高于埋深為3 500～4 000 m的4.37×108t；可動(dòng)油資源量則2個(gè)埋深段相當(dāng)，3 000～3 500 m 為 0.75×108t，可動(dòng)率為10.75%，3 500～4 000 m 為 0.73×108t，可動(dòng)率為16.70%，反映出深部頁巖具有更高的可動(dòng)油資源量。從頁巖油資源豐度和埋深大于3 500 m頁巖分布范圍疊合圖發(fā)現(xiàn)（圖6），利津洼陷的洼陷帶沙三段下亞段頁巖主體埋深大于3 500 m，具有高頁巖油資源豐度和高可動(dòng)率特征，是沙三段下亞段有利的頁巖油勘探區(qū)。

圖6 東營凹陷沙三段下亞段頁巖油資源豐度等值線分布Fig.6 Contour map of shale oil abundance of Lower Es3Member in Dongying Sag

4 結(jié)論

應(yīng)用體積法能很好地計(jì)算陸相斷陷盆地頁巖油資源量，其中滯留油定量分析可以選用氯仿瀝青“A”參數(shù)法或熱解參數(shù)法，2種方法計(jì)算的滯留油量有較好的一致性，具有較高的可信度。吸附油量可以通過吸附潛量進(jìn)行估算，頁巖中各組分的吸附潛量通過一定的模型能很好地表征頁巖的吸附能力，進(jìn)而進(jìn)行頁巖油資源量的計(jì)算。

東營凹陷沙三段下亞段頁巖油主要分布在各洼陷中心，資源豐度最高可達(dá)250×104t/km2，主要分布在埋深為3 000～3 500 m，但在3 500～4 000 m具有更高的可動(dòng)頁巖油資源豐度。利津洼陷帶沙三段下亞段頁巖油資源豐度和可動(dòng)率均相對(duì)較高，為研究區(qū)沙三段下亞段頁巖油有利的勘探區(qū)。

符號(hào)解釋：

Q資源量——頁巖油資源量，108t；S——泥頁巖面積，km2；h——泥頁巖有效厚度，m；ρ——泥頁巖密度，t/m3；a——泥頁巖中頁巖油含量，t/t；Qa——利用氯仿瀝青“A”含量計(jì)算的滯留油質(zhì)量分?jǐn)?shù)，%；Ka——氯仿瀝青“A”輕烴恢復(fù)系數(shù)；A——氯仿瀝青“A”測(cè)量值，%；Qs——利用熱解參數(shù)計(jì)算的滯留油質(zhì)量分?jǐn)?shù)，%；KS1——熱解輕烴散失系數(shù)；S1，S2——熱解烴量，mg/g；KS2——熱解S2中重?zé)N比例系數(shù)；Sp——頁巖油吸附潛量，mg/g；po——有機(jī)質(zhì)在泥頁巖中所占比例；xo——有機(jī)質(zhì)吸附潛量，mg/g；γ——泥頁巖比表面積相關(guān)系數(shù)；pi——第i種礦物在泥頁巖中所占比例；xi——第i種礦物的吸附潛量，mg/g；S0，S——孔隙壓實(shí)前和壓實(shí)后的比表面積，m2/g；?，?0——泥頁巖現(xiàn)今孔隙度和初始孔隙度；Q可動(dòng)油資源量——可動(dòng)油資源量，108t；Er——頁巖油可動(dòng)率，%。