朱日房 ,張林曄 ,李 政 ,王 茹 ,張守春 ,張 蕾
(1.中國(guó)石化勝利油田分公司勘探開(kāi)發(fā)研究院,山東東營(yíng)257015;2.國(guó)家能源局頁(yè)巖油研發(fā)中心勝利油田分室,山東東營(yíng)257015;3.中國(guó)石化頁(yè)巖油氣勘探開(kāi)發(fā)重點(diǎn)實(shí)驗(yàn)室,山東東營(yíng)257015;4.西北大學(xué)大陸動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,陜西西安710069)
北美頁(yè)巖油勘探開(kāi)發(fā)成果極大地激勵(lì)了中國(guó)頁(yè)巖油的勘探研究,中國(guó)東部新生代斷陷湖盆發(fā)育多套陸相富有機(jī)質(zhì)泥頁(yè)巖,并有多口探井發(fā)現(xiàn)了頁(yè)巖油氣流,顯示出良好的勘探前景[1-7]。近10 a來(lái),許多研究者在頁(yè)巖的含油性[8-13]、儲(chǔ)集性及流動(dòng)性方面進(jìn)行了探索[14-18]。頁(yè)巖油資源潛力評(píng)價(jià)從最初以滯留油為主的原地資源量評(píng)價(jià)發(fā)展為以游離油為主要賦存方式的游離油量評(píng)價(jià)[19-21]。由于頁(yè)巖油資源評(píng)價(jià)技術(shù)的不斷進(jìn)步、方法不斷完善,以及對(duì)頁(yè)巖中油氣賦存機(jī)理認(rèn)識(shí)的加深,必然對(duì)頁(yè)巖油的資源潛力評(píng)價(jià)提出更高的要求。頁(yè)巖油資源潛力評(píng)價(jià)研究仍然是目前頁(yè)巖油勘探實(shí)踐中重要的研究?jī)?nèi)容之一。為此,根據(jù)近幾年頁(yè)巖含油性定量分析等方面的進(jìn)展,以東營(yíng)凹陷沙三段下亞段為例,探討陸相斷陷盆地頁(yè)巖油資源潛力評(píng)價(jià)。
東營(yíng)凹陷位于濟(jì)陽(yáng)坳陷東南部,其東接青坨子凸起,南部與魯西隆起、廣饒凸起呈超覆接觸,西與惠民凹陷毗鄰,北與陳家莊凸起和濱縣凸起相鄰,東西長(zhǎng)度約為90 km,南北寬度約為62 km,總面積為5 850 km2(圖1)。從大地構(gòu)造區(qū)劃分上,東營(yíng)凹陷屬于中國(guó)東部渤海灣盆地的一個(gè)次級(jí)構(gòu)造單元,是在古生界基巖古地形背景上,經(jīng)構(gòu)造運(yùn)動(dòng)發(fā)育形成的中、新生代斷陷-坳陷盆地;由于基巖起伏以及中央隆起帶的分隔,凹陷內(nèi)又分為4個(gè)次級(jí)洼陷,它們既是盆地的沉積中心,也是富含有機(jī)質(zhì)泥頁(yè)巖發(fā)育區(qū)。凹陷內(nèi)主要發(fā)育沙四段上亞段和沙三段下亞段2套成熟生烴的泥頁(yè)巖[22],在110口井中的泥頁(yè)巖段發(fā)現(xiàn)油氣顯示,試油37口井,14口井獲得了工業(yè)油氣流,其中對(duì)河54井沙三段下亞段泥頁(yè)巖發(fā)育段試油,日產(chǎn)油量為 91.4 t/d[23],累積產(chǎn)油量為27 896 t,顯示出良好的頁(yè)巖油勘探前景。
圖1 東營(yíng)凹陷區(qū)域構(gòu)造位置Fig.1 Regional tectonic location of Dongying Sag
頁(yè)巖油資源潛力評(píng)價(jià)方法主要分為動(dòng)態(tài)法和靜態(tài)法2大類(lèi)。動(dòng)態(tài)法主要根據(jù)頁(yè)巖油在開(kāi)發(fā)過(guò)程中的動(dòng)態(tài)資料,通過(guò)一定的數(shù)學(xué)模型定量計(jì)算資源量。靜態(tài)法主要應(yīng)用靜態(tài)參數(shù)計(jì)算資源量,按計(jì)算方式的不同細(xì)分為統(tǒng)計(jì)法、類(lèi)比法和成因法3種。統(tǒng)計(jì)法需要大量典型樣本,適用于中-高勘探程度階段;類(lèi)比法需要有可類(lèi)比的刻度區(qū),一般適用于低勘探程度區(qū);成因法重視物質(zhì)平衡原理的應(yīng)用,更多是確定性評(píng)價(jià),適用于盆地勘探的各個(gè)階段。目前由于缺少系統(tǒng)的頁(yè)巖油開(kāi)采動(dòng)態(tài)資料和典型樣本,難以應(yīng)用動(dòng)態(tài)法和統(tǒng)計(jì)法,同時(shí)中國(guó)東部陸相盆地的地質(zhì)情況與北美具有較大差別,難以直接將北美成熟頁(yè)巖油開(kāi)采區(qū)作為中國(guó)頁(yè)巖油資源潛力評(píng)價(jià)類(lèi)比刻度區(qū)。因此,成因法是中國(guó)頁(yè)巖油資源潛力評(píng)價(jià)最常用的方法。
體積法是成因法中最常用的一種方法,其計(jì)算頁(yè)巖油資源量的原理簡(jiǎn)單明確,將泥頁(yè)巖的質(zhì)量(頁(yè)巖體積與密度乘積)與泥頁(yè)巖中頁(yè)巖油含量相乘即可獲得頁(yè)巖油資源量,其計(jì)算式為:
計(jì)算頁(yè)巖油資源量所需的參數(shù)主要包括泥頁(yè)巖的面積、有效厚度、密度和頁(yè)巖油含量。泥頁(yè)巖的面積和有效厚度主要通過(guò)頁(yè)巖厚度等值線(xiàn)圖獲取,密度由實(shí)測(cè)密度剖面獲取。由于頁(yè)巖油主要以游離油的形式賦存于泥頁(yè)巖中,因此,泥頁(yè)巖中的頁(yè)巖油含量用游離油量來(lái)近似代替,而游離油量的定量分析是頁(yè)巖油資源量計(jì)算中的難點(diǎn)之一。
游離油量的定量分析主要有2種方法,一種是孔隙度和含油飽和度法,另一種是地化參數(shù)法。對(duì)于常規(guī)砂巖儲(chǔ)層來(lái)說(shuō),孔隙度和含油飽和度法表征游離油量簡(jiǎn)單明確,然而對(duì)于非均質(zhì)性強(qiáng)、孔隙結(jié)構(gòu)復(fù)雜、以納米孔為主的泥頁(yè)巖來(lái)說(shuō),其孔隙度和含油飽和度的測(cè)量非常困難,數(shù)據(jù)點(diǎn)少,難以用于宏觀(guān)的頁(yè)巖油資源量計(jì)算;因此,地化參數(shù)法成為最常用的游離油量定量分析方法。
應(yīng)用地化參數(shù)法很難直接獲取游離油量和吸附油量,但能夠確定頁(yè)巖中的滯留油量和巖石對(duì)油的吸附潛量。在飽和吸附時(shí),頁(yè)巖對(duì)油的吸附潛量近似于頁(yè)巖的吸附量,頁(yè)巖中的滯留油量與頁(yè)巖的吸附潛量之間的差值即為頁(yè)巖的游離油量。
表征頁(yè)巖滯留油量最常用的參數(shù)為氯仿瀝青“A”含量和熱解參數(shù),相應(yīng)的計(jì)算式分別為:
頁(yè)巖對(duì)油的吸附能力與巖石的結(jié)構(gòu)、礦物組分和有機(jī)質(zhì)豐度、類(lèi)型及演化程度等因素密切相關(guān)[22-28]。不同的礦物組分對(duì)油的吸附能力明顯不同,有機(jī)質(zhì)對(duì)油的吸附能力遠(yuǎn)大于無(wú)機(jī)礦物對(duì)油的吸附能力[27-28],不同類(lèi)型的有機(jī)質(zhì)以及熱演化程度不同的有機(jī)質(zhì)對(duì)油的吸附能力也有明顯區(qū)別[24]??紫督Y(jié)構(gòu)影響無(wú)機(jī)礦物與油的接觸面積,從而影響無(wú)機(jī)礦物對(duì)油的吸附能力[25]。因此,可以分別計(jì)算干酪根及不同孔隙度條件下各礦物組分的吸附潛量來(lái)表征整個(gè)巖石對(duì)油的吸附能力,同時(shí),在計(jì)算無(wú)機(jī)礦物吸附潛量時(shí)加入與孔隙結(jié)構(gòu)有關(guān)的泥頁(yè)巖比表面積相關(guān)系數(shù)。頁(yè)巖對(duì)油吸附潛量的計(jì)算式為:
其中
在獲取滯留油量和頁(yè)巖對(duì)油的吸附潛量后,就可以計(jì)算頁(yè)巖的吸附油量和游離油量。當(dāng)頁(yè)巖中滯留油量小于頁(yè)巖對(duì)油的吸附潛量時(shí),為欠飽和吸附,吸附油量就等于滯留油量,此時(shí)沒(méi)有游離油;當(dāng)頁(yè)巖中滯留油量大于或等于頁(yè)巖對(duì)油的吸附潛量時(shí),為飽和吸附,吸附油量就等于頁(yè)巖對(duì)油吸附潛量,游離油量則為滯留油量與吸附油量的差值。
可動(dòng)油主要指在地層條件下,采用一定的工程措施能夠采出的原油。可動(dòng)油資源量對(duì)頁(yè)巖油的勘探開(kāi)發(fā)具有重要意義,是頁(yè)巖油資源潛力最直接的衡量指標(biāo)。可動(dòng)油資源量表達(dá)式為:
頁(yè)巖油可動(dòng)率與泥頁(yè)巖的骨架、原油、地層水性質(zhì)等有關(guān),還與地層流體壓力有關(guān),巖石骨架、原油和水的壓縮性越大、地層流體壓力越高、能量越大,頁(yè)巖油可動(dòng)率就越高。
東營(yíng)凹陷沙三段下亞段為半咸水沉積,主要為灰色、灰褐色頁(yè)巖與深灰色泥巖組成的不等厚互層,其厚度一般為100~300 m。地化分析資料表明,研究區(qū)沙三段下亞段總有機(jī)碳含量(TOC)分布范圍較寬,為0.5%~18.6%,平均值為4.9%;氯仿瀝青“A”含量為0.11%~2.94%,平均值為0.92%;有機(jī)質(zhì)類(lèi)型以Ⅰ型干酪根為主,少部分為Ⅱ1型干酪根,埋深為1 500~4 000 m,各洼陷內(nèi)基本上處于成熟生烴階段。沙三段下亞段泥頁(yè)巖有機(jī)質(zhì)豐度高、類(lèi)型好、主體處于成熟生烴階段,是一套生烴潛力大、分布穩(wěn)定的主力烴源巖。
3.2.1 頁(yè)巖含油性分級(jí)
后綴法:“詞根+后綴”。 如:-meter“計(jì)量?jī)x器”,可后綴派生出 barometer氣壓計(jì),telemeter測(cè)距儀,spectrometer分光儀,oscilloscope示波器等多個(gè)新詞。
湖相泥頁(yè)巖的含油性在平面上和縱向上均具有很強(qiáng)的非均質(zhì)性,頁(yè)巖含油性分級(jí)描述能有效地消除頁(yè)巖非均質(zhì)性的影響[29-31]。按照頁(yè)巖含油性參數(shù)(TOC、氯仿瀝青“A”含量和S1)及其相關(guān)性,結(jié)合東營(yíng)凹陷沙三段下亞段烴源巖具體特征,將頁(yè)巖含油性分為Ⅰ,Ⅱ和Ⅲ共3個(gè)級(jí)別。具體劃分方案見(jiàn)表1。
表1 東營(yíng)凹陷頁(yè)巖含油性分級(jí)劃分方案Table1 Division scheme for the oil-bearing abundance of shale in Dongying Sag
采用系統(tǒng)取心井與測(cè)井曲線(xiàn)相結(jié)合,應(yīng)用改進(jìn)的Δlog R模型[32]建立TOC、氯仿瀝青“A”含量和S1的解釋模型,進(jìn)而計(jì)算頁(yè)巖含油性參數(shù)。根據(jù)計(jì)算結(jié)果和表1劃分方案,可以對(duì)研究區(qū)各單井的頁(yè)巖含油性進(jìn)行分級(jí)、統(tǒng)計(jì)各級(jí)厚度及含油性參數(shù),編制各級(jí)的厚度等值線(xiàn)和地化參數(shù)等值線(xiàn)(圖2,圖3)。從圖2可以看出,東營(yíng)凹陷Ⅰ級(jí)頁(yè)巖的厚度中心為各洼陷中心,厚度均在150 m以上,其中利津、民豐達(dá)到200 m;各含油性參數(shù)也具有類(lèi)似的分布特征,利津洼陷中心處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í)頁(yè)巖與Ⅰ級(jí)頁(yè)巖具有相似的特征,牛莊、利津洼陷中心頁(yè)巖厚度達(dá)到了200 m,博興、民豐洼陷為100 m,各洼陷中心的TOC、氯仿瀝青“A”含量和S1值分別為2.0%,0.4%和2.0 mg/g。
圖2 東營(yíng)凹陷沙三段下亞段Ⅰ級(jí)頁(yè)巖含油性參數(shù)等值線(xiàn)分布Fig.2 Contour map of oil content parameters of GradeⅠshale of Lower Es3Member in Dongying Sag
3.2.2 頁(yè)巖密度
根據(jù)實(shí)測(cè)的密度資料(圖4)發(fā)現(xiàn),東營(yíng)凹陷泥頁(yè)巖密度表現(xiàn)出隨埋深增加、壓實(shí)程度增強(qiáng)而增大的特征,但在不同埋深段密度隨埋深的變化曲線(xiàn)呈不同的斜率。在同一埋深下,烴源巖密度與TOC呈一定程度的負(fù)相關(guān)關(guān)系,即TOC越高,密度越小。由此,泥頁(yè)巖的密度可以由不同TOC在不同埋深下的密度-埋深關(guān)系曲線(xiàn)(圖4)來(lái)獲取。
3.2.3 輕烴恢復(fù)系數(shù)及重?zé)N比例系數(shù)
輕烴恢復(fù)系數(shù)及重?zé)N比例系數(shù)采用文獻(xiàn)[8]的資料,各系數(shù)隨埋深的增加而增大。在主要的生油窗范圍內(nèi),東營(yíng)凹陷氯仿瀝青“A”恢復(fù)系數(shù)為1.0~1.5,熱解輕烴散失系數(shù)為0.1~0.4,熱解S2中重?zé)N比例系數(shù)則為 0.1~0.8[8]。
根據(jù)東營(yíng)凹陷沙三段下亞段實(shí)際樣品的地化資料、巖礦資料和物性資料,以及不同組分對(duì)油的吸附能力,應(yīng)用(4)式就可計(jì)算不同樣品對(duì)油的吸附潛量。在生油窗范圍內(nèi),東營(yíng)凹陷沙三段下亞段干酪根對(duì)油的單位TOC吸附潛量為50~120 mg/g[28],無(wú)機(jī)礦物中粘土礦物、石英、碳酸鹽和黃鐵礦對(duì)油的吸附潛量分別為18.0,3.0,1.8和10.0 mg/g[25],泥頁(yè)巖的初始孔隙度為 50%,現(xiàn)今孔隙度為5%~15%。東營(yíng)凹陷沙三段下亞段頁(yè)巖油吸附潛量計(jì)算結(jié)果(圖5)表明,東營(yíng)凹陷沙三段下亞段泥頁(yè)巖在主要生烴初始階段頁(yè)巖油單位TOC吸附潛量為200~300 mg/g,略大于國(guó)外的 200 mg/g[29],與應(yīng)用A/TOC確定的230 mg/g結(jié)果較一致[19,28]。在埋深為4 000 m左右,頁(yè)巖油單位TOC吸附潛量約為100~150 mg/g。
圖3 東營(yíng)凹陷沙三段下亞段Ⅱ級(jí)頁(yè)巖含油性參數(shù)等值線(xiàn)分布Fig.3 Contour map of oil content parameters of GradeⅡshale of Lower Es3Member in Dongying Sag
圖4 東營(yíng)凹陷泥頁(yè)巖密度隨深度變化曲線(xiàn)Fig.4 Variation curves of shale density with depth in Dongying Sag
圖5 東營(yíng)凹陷沙三段下亞段頁(yè)巖油吸附潛量演化剖面Fig.5 Evolution of adsorption oil capacity of shale in Lower Es3Member in Dongying Sag
在確定不同深度下的氯仿瀝青“A”輕烴恢復(fù)系數(shù)、熱解輕烴散失系數(shù)、熱解S2中重?zé)N比例系數(shù)、頁(yè)巖油吸附潛量以及泥頁(yè)巖密度等參數(shù)后,根據(jù)研究區(qū)不同級(jí)別的泥頁(yè)巖總有機(jī)碳含量、氯仿瀝青“A”含量、S1、有效厚度、埋深等參數(shù),計(jì)算濟(jì)陽(yáng)坳陷不同級(jí)別泥頁(yè)巖的滯留油量和游離油量。計(jì)算過(guò)程中,為求精確,將不同級(jí)別泥頁(yè)巖分布區(qū)在平面上均分為500 m×500 m的網(wǎng)格區(qū),分別應(yīng)用氯仿瀝青“A”參數(shù)法和熱解參數(shù)法計(jì)算各網(wǎng)格區(qū)的滯留油量和游離油量,所有網(wǎng)格的滯留油量之和為總滯留油量,游離油量之和為總頁(yè)巖油量即頁(yè)巖油資源量。
研究結(jié)果表明,2種方法的計(jì)算結(jié)果(表2)非常接近,采用氯仿瀝青“A”參數(shù)法和熱解參數(shù)法,東營(yíng)凹陷沙三段下亞段總滯留油量分別為65.44×108和64.26×108t;頁(yè)巖油資源量分別為11.39×108和11.57×108t。其中利津洼陷頁(yè)巖油最豐富,資源量超過(guò)6.00×108t,占比超過(guò)50%,其次為牛莊洼陷。應(yīng)用張林曄等建立的“彈性驅(qū)動(dòng)+溶解氣驅(qū)動(dòng)可動(dòng)油率曲線(xiàn)”[28]作為頁(yè)巖油可動(dòng)率來(lái)計(jì)算頁(yè)巖油可動(dòng)資源量,得到東營(yíng)凹陷沙三段下亞段頁(yè)巖油可動(dòng)資源量為1.50×108t。
表2 東營(yíng)凹陷沙三段下亞段頁(yè)巖油資源量計(jì)算結(jié)果Table2 Calculation results of shale oil resources of Lower Es3Member in Dongying Sag
從東營(yíng)凹陷沙三段下亞段頁(yè)巖油資源豐度的平面分布特征看(圖6),其頁(yè)巖油資源量主要分布在各洼陷帶,從洼陷中心向邊部逐漸減少,洼陷中心的資源豐度以利津洼陷最高,可達(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%,反映出深部頁(yè)巖具有更高的可動(dòng)油資源量。從頁(yè)巖油資源豐度和埋深大于3 500 m頁(yè)巖分布范圍疊合圖發(fā)現(xiàn)(圖6),利津洼陷的洼陷帶沙三段下亞段頁(yè)巖主體埋深大于3 500 m,具有高頁(yè)巖油資源豐度和高可動(dòng)率特征,是沙三段下亞段有利的頁(yè)巖油勘探區(qū)。
圖6 東營(yíng)凹陷沙三段下亞段頁(yè)巖油資源豐度等值線(xiàn)分布Fig.6 Contour map of shale oil abundance of Lower Es3Member in Dongying Sag
應(yīng)用體積法能很好地計(jì)算陸相斷陷盆地頁(yè)巖油資源量,其中滯留油定量分析可以選用氯仿瀝青“A”參數(shù)法或熱解參數(shù)法,2種方法計(jì)算的滯留油量有較好的一致性,具有較高的可信度。吸附油量可以通過(guò)吸附潛量進(jìn)行估算,頁(yè)巖中各組分的吸附潛量通過(guò)一定的模型能很好地表征頁(yè)巖的吸附能力,進(jìn)而進(jìn)行頁(yè)巖油資源量的計(jì)算。
東營(yíng)凹陷沙三段下亞段頁(yè)巖油主要分布在各洼陷中心,資源豐度最高可達(dá)250×104t/km2,主要分布在埋深為3 000~3 500 m,但在3 500~4 000 m具有更高的可動(dòng)頁(yè)巖油資源豐度。利津洼陷帶沙三段下亞段頁(yè)巖油資源豐度和可動(dòng)率均相對(duì)較高,為研究區(qū)沙三段下亞段頁(yè)巖油有利的勘探區(qū)。
符號(hào)解釋?zhuān)?/p>
Q資源量——頁(yè)巖油資源量,108t;S——泥頁(yè)巖面積,km2;h——泥頁(yè)巖有效厚度,m;ρ——泥頁(yè)巖密度,t/m3;a——泥頁(yè)巖中頁(yè)巖油含量,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——頁(yè)巖油吸附潛量,mg/g;po——有機(jī)質(zhì)在泥頁(yè)巖中所占比例;xo——有機(jī)質(zhì)吸附潛量,mg/g;γ——泥頁(yè)巖比表面積相關(guān)系數(shù);pi——第i種礦物在泥頁(yè)巖中所占比例;xi——第i種礦物的吸附潛量,mg/g;S0,S——孔隙壓實(shí)前和壓實(shí)后的比表面積,m2/g;?,?0——泥頁(yè)巖現(xiàn)今孔隙度和初始孔隙度;Q可動(dòng)油資源量——可動(dòng)油資源量,108t;Er——頁(yè)巖油可動(dòng)率,%。