李彥婧
(中石化華東分公司物探研究院,江蘇南京210001)
利用測井資料定量計算頁巖層有機碳含量
李彥婧*
(中石化華東分公司物探研究院,江蘇南京210001)
在地球化學分析的基礎上,利用常規測井資料,采用多元統計分析法、?lgR法,建立了測井參數與有機碳含量之間的定量關系,獲取了川東南地區X井龍馬溪組頁巖層連續的有機碳含量值。結果表明,多元統計分析法計算結果與實測值有更好的相關性,為該區的頁巖評價提供了依據。
測井資料;有機碳含量;多元統計分析法;?lgR法
有機碳含量(TOC)是頁巖氣選區評價中的一項關鍵地質參數[1-2],但一般很難獲取其在縱向上的連續分布值[2-5]。常規測井曲線對有機質有特殊的響應特征[5-7],通過建立測井參數與有機碳含量之間的定量關系,可以計算頁巖層連續的有機碳含量值。本文主要應用多元統計分析法、?lgR法來擬合計算頁巖層TOC,為尋找優質頁巖段提供參考。
川東南地區X井龍馬溪組頁巖層有65個巖芯樣品的實測TOC值。通過繪制該井實測TOC值與對應深度的測井參數相關性圖發現(圖1),5種測井參數與TOC回歸方程的斜率不同,回歸方程均選用線性回歸,TOC與自然伽馬、電阻率、聲波時差3種測井參數呈正相關,與中子孔隙度、體積密度2種測井參數呈負相關。綜合以上5種測井參數,多元統計分析的結果將會更準確。
以實測TOC值為因變量,對應深度的自然伽馬、電阻率、聲波時差、中子孔隙度、體積密度等5種測井參數為自變量,進行線性、非線性回歸,建立回歸方程。在多元統計分析過程中,采用逐步回歸的方法,依次排除在可信度區間外的自變量組合類型,建立最優回歸方程。回歸結果見表1。
從表1中結果來看,隨著測井參數種類的增加,相關系數逐漸增大,殘差(實測值與回歸值之間的差異)之和逐漸減小。四參數回歸方程的相關系數為0.578,殘差之和為22.16,屬于該井有機碳含量定量預測的最優回歸方程。
應用多元統計分析法計算X井連續的TOC曲線(圖2),從圖2可以看出,計算TOC值與實測TOC值一致性良好,優質頁巖(TOC>2%)厚度為18.1m。在實測TOC值底部峰值處,多元統計分析方法存在誤差。原因在于,多元統計分析方法計算得到的是平均值,會使高值計算過低,因此反映不明顯。
將算數坐標下的聲波曲線和對數坐標下的電阻率曲線在非烴源巖處疊合在一起(每2個數量級的對數刻度電阻率對應聲波時差間隔為328μs/m),確定基線位置?;€位置確定之后,兩條曲線之間的間距就是?lgR。?lgR計算公式為:

式中:?lgR——聲波曲線與電阻率曲線疊合的幅度差;
R——實測的深側向電阻率,Ω·m;
?t——聲波時差,μs/m;
R基線——基線處的深側向電阻率,Ω·m;
t基線——基線處的的聲波時差,μs/m。
Passey模型[8]的有機碳含量計算公式為:

式中:LOM——熱變指數,反映有機質成熟度;
?TOC——該區有機碳含量背景值。將公式(1)代入公式(2),即可計算有機碳含量值。川東南地區X井龍馬溪組有機碳含量計算公式為TOC=2.756×?lgR+0.22,應用?lgR法計算該井連續的TOC曲線(圖2)。結果顯示:在龍馬溪組上段,計算TOC值與實測值一致性高,受LOM、?TOC取值影響,求得的TOC值總體上均小于2%,比實測值低。?lgR法計算的TOC值反映了有機碳含量在縱向上的相對變化趨勢,在縱向上有2個峰值,有機質豐度較高,優質頁巖段位于龍馬溪組底部。

圖1 川東南地區X井測井參數與TOC相關性圖

表1 川東南地區X井多元統計分析回歸方程
分別將2種方法計算TOC值和實測TOC值擬合,繪制相關性擬合圖(圖3)。從圖3可以看出,計算TOC值與實測TOC值基本呈線性趨勢,多元統計分析法相關系數為0.59、?lgR法相關系數為0.15,2種方法的計算結果存在差異。多元統計分析法優于?lgR法,擬合的效果具有可操作性,誤差在可接受的范圍內。
(1)自然伽馬、體積密度、中子孔隙度、電阻率組合的四參數回歸方程計算TOC值與實測值相關系數達到0.578,是X井龍馬溪組頁巖層有機碳含量預測的最優回歸方程;
(2)?lgR法求得的TOC值總體上小于實測值,縱向上的峰值位于龍馬溪組底部;
(3)在現有利用測井資料計算有機碳含量的方法中,X井的應用表明:多元統計分析法優于?lgR法。
[1]陳新軍,包書景,侯讀杰,等.頁巖氣資源評價方法與關鍵參數探討[J].石油勘探與開發,2012,39(5):566-571.
[2]徐思煌,朱義清.烴源巖有機碳含量的測井響應特征與定量預測模型——以珠江口盆地文昌組烴源巖為例[J].石油實驗地質,2010,32(3):290-300.
[3]劉新穎,鄧宏文,邸永香,等.海拉爾盆地烏爾遜凹陷南屯組優質烴源巖發育特征[J].石油實驗地質,2009,31(1):68-73.
[4]郭龍,陳踐發,苗忠英.一種新的TOC含量擬合方法研究與應用[J].天然氣地球科學,2009,20(6):951-956.
[5]王鸝,唐成磊,王飛.?lgR技術在烴源巖評價中的應用[J].斷塊油氣田,2009,16(2):52-54.
[6]王貴文,朱振宇,朱廣宇.烴源巖測井識別與評價方法研究[J].石油勘探與開發,2002,29(4):50-53.
[7]郭澤清,孫平,劉衛紅.利用△lgR技術計算柴達木盆地三湖地區第四系有機碳[J].地球物理學進展,2012,27(2):626-633.
[8]Passey Q R,Creaney S,Kulla J B.A Practical Model for Organic Richness from Porosity and Resistivity Logs[J].AAPG Bulletin,1990,74(5):1777-1794.

圖2 川東南地區X井龍馬溪組TOC計算值與實測值對比圖

圖3 相關性分析
TE122.1
A
1004-5716(2015)01-0140-03
2014-06-12
2014-06-16
李彥婧(1987-),女(漢族),湖北鄂州人,助理工程師,現從事頁巖氣地質方面工作。