張學娟，張 雷，馬 磊，袁子龍

(1.東北石油大學地球科學學院，黑龍江大慶 163318;2.東北石油大學非常規(guī)油氣成藏與開發(fā)省部共建重點實驗室，黑龍江大慶 163318)

ΔlogR技術(shù)［1］是由 Q.R.Passey等通過大量數(shù)據(jù)擬合及實驗分析得出的烴源巖評價方法。該技術(shù)利用電阻率和聲波時差兩條測井曲線對烴源巖的差異響應(yīng)刻畫有機質(zhì)豐度。ΔlogR技術(shù)的優(yōu)點是聲波時差和電阻率曲線都對孔隙度敏感，兩條曲線隨孔隙度的變化成一定比例，只要比例選取適當，孔隙度變化使得兩條曲線產(chǎn)生的偏移幅度相同，從而抵消了孔隙度對有機碳測井響應(yīng)的影響。但該方法的缺陷是需人為確定聲波時差和電阻率曲線的基線值，確定基線值的操作較為繁瑣且主觀性較大。另外，預先給定疊合系數(shù)為0.02的合理性需要進一步驗證。

ΔlogR技術(shù)引入國內(nèi)后，得到推廣應(yīng)用。朱光有等［2－3］在對東營凹陷沙四段烴源巖有機質(zhì)豐度非均質(zhì)性研究過程中采用ΔlogR技術(shù)進行評價，取得了較好的評價效果。劉超［4］對ΔlogR技術(shù)進行了一定的改進，總體方法是以貧有機質(zhì)泥巖段為節(jié)點，將目標層段細分，然后取細分層段內(nèi)貧有機碳泥巖ΔlogR為零，動態(tài)選取疊合系數(shù)，對測井曲線動態(tài)疊合，降低了人為確定基線帶來的主觀性。胡慧婷等［5］在其基礎(chǔ)上，通過逐步回歸法建立了改進ΔlogR模型，在松遼盆地及其外圍盆地烴源巖有機質(zhì)豐度的評價中取得了良好效果。湯麗娜［6］對比分析了多元統(tǒng)計回歸定量表征烴源巖和ΔlogR敏感曲線直接定量計算有機碳含量 (TOC)這兩種方法，發(fā)現(xiàn)利用ΔlogR敏感曲線計算的TOC精度高、誤差較小，對湖相烴源巖識別效果可靠。韓雙彪等［7］在評價吐哈盆地侏羅系水西溝群的烴源巖時，依據(jù)不同的沉積環(huán)境分別構(gòu)建ΔlogR的多元線性回歸模型，建立了基于烴源巖成因分類的有機質(zhì)豐度評價方法。高陽［8］在評價東營凹陷沙四下亞段鹽湖相烴源巖時發(fā)現(xiàn)，鹽湖相烴源巖受烴源巖中膏鹽質(zhì)成分、含量及有機質(zhì)豐度的共同影響，測井響應(yīng)特征相當復雜，提出了依據(jù)自然伽馬絕對值的大小 (反映泥質(zhì)含量值)分段ΔlogR擬合計算TOC的方法，實際應(yīng)用效果較好。袁彩萍等［9］在對珠江口盆地文昌組烴源巖TOC評價過程中建立了TOC與電阻率、中子孔隙度、聲波時差及自然伽馬等多元線性回歸方程，取得了較好的效果。

綜上所述，前人的研究對象多以深湖—半深湖相泥質(zhì)烴源巖為主，大民屯凹陷高含有機質(zhì)的油頁巖形成環(huán)境與深湖—半深湖存在一定差異，其主要形成于干旱氣候、湖平面頻繁升降的淺水穩(wěn)定局部臺地構(gòu)造背景，沉積環(huán)境為封閉、半封閉狀態(tài)的弱還原、弱氧化環(huán)境。本文采用改進的ΔlogR技術(shù)對大民屯凹陷油頁巖有機質(zhì)豐度進行評價，以期建立起一套適用于此類成因油頁巖的有機質(zhì)豐度評價方法。

1 大民屯凹陷沙四下亞段烴源巖發(fā)育基本特征

沙四下亞段 (Es42)烴源巖、沙四上亞段(Es41)和沙三段四亞段 (Es34)構(gòu)成了大民屯凹陷的主要烴源巖系統(tǒng)［10－12］。大民屯凹陷中部的安福屯—靜安堡周邊地區(qū)沙四下亞段烴源巖為一套油頁巖和鈣鎂質(zhì)油頁巖 (圖1)，所生成的油氣以高含蠟質(zhì)為主要特征。該套烴源巖含砂率極低(僅為2%～5%)，分布面積較廣 (達250km2)，平均厚度為120～250m，在東西兩側(cè)的安福屯斷階帶和勝東洼陷烴源巖最為發(fā)育，最大厚度達350m，凹陷中央的臺地烴源巖相對欠發(fā)育，烴源巖厚度一般為80～120m。沙四下亞段烴源巖內(nèi)部有機質(zhì)豐度的非均質(zhì)性直接影響生排烴能力的差異，進而影響資源的分布。

圖1 大民屯凹陷區(qū)域構(gòu)造簡圖Fig.1 Sketch of regional structure in Damintun sag

大民屯凹陷安福屯—靜安堡周邊地區(qū)沙四段油頁巖分為3段，其中上段和下段油頁巖中礦物成分大都是氧化硅和黏土礦物，也可能含少量碳酸鹽，常與碳酸鹽巖、硅質(zhì)巖類和含磷酸鹽的巖類伴生，表現(xiàn)為微細的水平紋層油頁巖和具透鏡狀內(nèi)部結(jié)構(gòu)的油頁巖。沙四下亞段中部主要發(fā)育與蒸發(fā)鹽巖共生的灰質(zhì)油頁巖或白云質(zhì)油頁巖，在相對濕潤時期，湖平面上升，沉積范圍擴大，底棲的藍綠藻和真菌形成了絮凝狀的富有機質(zhì)的軟泥;干旱時期，湖平面下降，湖區(qū)范圍縮小，富含有機質(zhì)的軟泥濃縮，與砂和粉砂級的白云質(zhì)微晶內(nèi)碎屑伴生;湖平面再次上升，就形成透鏡狀 “內(nèi)碎屑”油頁巖［13－14］。

2 ΔlogR技術(shù)及數(shù)據(jù)歸位方法

2.1 模型的建立

據(jù) Q.R.Passey［1］提出的 ΔlogR技術(shù)公式為:

式中 AC——實測的聲波時差，μs/ft;

Rt——巖石的實測電阻率，Ω·m。

Rt基線——非烴源巖層段對應(yīng)電阻率測井值;

AC基線——非烴源巖層段對應(yīng)聲波時差;

0.02——疊合系數(shù)，即每 50μs/ft(164μs/m)的聲波時差。

在一定的成熟度條件下，ΔlogR與有機碳呈線性相關(guān)，由ΔlogR計算有機碳的經(jīng)驗公式［1］為:

式中 TOC——計算的有機碳含量，%;

LOM——有機質(zhì)成熟度;

ΔTOC——有機碳含量背景值。

LOM可根據(jù)區(qū)域大量的成熟度指標分析 (如熱變指數(shù)、Tmax分析、鏡質(zhì)組反射率等)得到。

該方法需要成熟度參數(shù)，且確定巖性基線位置的誤差較大，操作繁瑣。對于同一區(qū)塊相同層位烴源巖來說，烴源巖的成熟度基本一致，而AC基線和Rt基線為常量，即得到文獻 ［4］中的公式:

式中 a、b、c——ΔlogR系數(shù)。

改進模型在無需LOM、ΔTOC及人為讀取基線值的條件下便可計算出有機碳含量。

大民屯凹陷鉆遇沙四段的眾多井中，選取實測TOC點較多且具有連續(xù)巖性剖面和完整測井資料的井為模型建立標準井 (沈352井)，建模前首先對聲波時差與電阻率曲線進行標準化，使其量綱一致，將標準化后的曲線與建模標準井的實測TOC進行擬合優(yōu)化，得到系數(shù)a、b和c。

沈352井沙四下亞段在建模過程中，需剔除非烴源巖層段;同時消除非地層響應(yīng)的異常測井段的影響，包括井徑變化、儀器穩(wěn)定性等影響。標準化后，對沙四下亞段電阻率測井值求對數(shù)，再對聲波時差和電阻率對數(shù)進行歸一化處理，使得電阻率的對數(shù)和聲波時差在相同的數(shù)量級范圍內(nèi)。

2.2 數(shù)據(jù)歸位

實測有機質(zhì)樣品的深度受鉆井取心作業(yè)、鉆具測量長度等誤差的影響，巖心深度不準;因此，在整理巖心資料時，一般會參照測井曲線進行歸位，這種歸位受人為主觀影響較大，效果不理想。此外，分析取樣的位置也會帶來一定的誤差［15］。巖屑樣品深度受鉆井液回返遲到時間計算影響，誤差更為顯著。如何對巖心和巖屑樣品的深度進行高精度歸位是進行測井評價烴源巖的必要前提。

劉超［4］認為巖屑的最大歸位幅度為1m，即巖屑樣品原始位置上下各浮動0.5m較為合理，通過實測TOC與計算TOC的相關(guān)系數(shù)最大為目標函數(shù)進行擬合。主觀確定最大歸位幅度存在的問題是過大幅度浮動樣品深度會使得樣品歸位失真，而過小幅度浮動樣品深度又會使得歸位不徹底，所以需要進一步合理地確定歸位深度。在沈352井模型建立過程中，由于巖屑深度的誤差較大，選擇巖屑的最大歸位深度為巖心的2倍，將巖心的最大歸位幅度作為一個變量X(即巖心樣品原始位置上下各浮動X/2)，對應(yīng)巖屑樣品的最大歸位幅度為2X(即巖屑樣品原始位置上下各浮動X)，尋求樣品的合理歸位幅度。

2.2.1 確定歸位步長及最大歸位深度范圍

選取樣品歸位步長d為0.125m(測井數(shù)據(jù)采樣點間距為0.125m)，歸位過程中以淺部樣品的埋深不超過深部樣品埋深為原則，從上至下進行如下歸位。

(1)對于巖心樣品。

①當i=1時:

S1=X，H2－H1≥X/2

S1=X/2+H2－H1，H2－H1≤X/2

②當i=2，3，…，n－1時:

Si=X，Hi+1－ Hi≥ X/2，Hi－ Hi－1'≥ X/2

Si=Hi+1－Hi－1'，Hi+1－Hi＜ X/2，Hi－Hi－1'＜ X/2

Si=X/2+Hi－ Hi－1'，Hi+1－ Hi≥ X/2，Hi－Hi－1' ＜ X/2

Si=X/2+Hi+1－ Hi，Hi+l－ Hi＜ X/2，Hi－Hi－1' ≥ X/2

③當i=n時:

Sn=X，Hn－ Hn－1'≥ X/2

Sn=X/2+Hn－ Hn－1'，Hn－ Hn－1'≤ X/2

式中 Si——第i個烴源巖樣品的最大歸位幅度;

Hi——第i個烴源巖樣品的初始深度;

Hi+1——第i+1個烴源巖樣品的初始深度;

Hi－1'——第 i－1 個烴源巖樣品歸位后的深度;

n——樣品總數(shù)。

(2)對于巖屑樣品。

①當i=1時:

S1=2X，H2－H1≥X

S1=X+H2－H1，H2－H1≤X

②當i=2，3，…，n－1時:

Si=2X，Hi+1－ Hi≥ X，Hi－ Hi－1'≥ X

Si=Hi+1－ Hi－1'，Hi+l－ Hi＜ X，Hi－ Hi－1' ＜ X

Si=X+Hi－ Hi－1'，Hi+l－ Hi≥X，Hi－ Hi－1' ＜ X

Si=X+Hi+1－ Hi，Hi+l－ Hi＜ X，Hi－ Hi－1'≥X

Sn=2X，Hn－ Hn－1'≥ X

Sn=X+Hn－ Hn－1'，Hn－ Hn－1'≤ X

2.2.2 對各樣品點深度逐步歸位校正

由淺至深，從第1個烴源巖樣品到第n個烴源巖樣品，對烴源巖樣品逐個逐步長地進行歸位，即第i個巖心或巖屑樣品歸位后的深度為Hi':

式中 Hi——第i個烴源巖巖心或者巖屑樣品的記錄深度;

j——歸位次數(shù)，第i個樣品的歸位次數(shù)為2Si/d。

總的歸位校正組合數(shù)P可以表示為:

2.2.3 構(gòu)建目標函數(shù)，評價各樣品歸位校正深度

以實測TOC與歸位后計算TOC之差的平方和為目標函數(shù)，對各樣品歸位校正深度組合效果進行評價，當目標函數(shù)值達到最小值時，擬合效果最優(yōu)，即:

式中 Ri'——第i個樣品校正深度后所對應(yīng)的電阻率值;

ACi'——第i個樣品校正深度后所對應(yīng)的聲波時差值;

TOCi——第i個樣品的有機碳含量。

通過最小二乘法 (當目標函數(shù)F達到極小值時，目標函數(shù)對a、b和c偏導值為0)來求取a、b和c，再將求取的a、b和c代入式 (6)中，求得目標函數(shù)的最小值Fmin。從所有歸位校正組合中選擇Fmin的最小值及其對應(yīng)的a、b和c，即該歸位組合為目標歸位組合，記錄Fmin最小值及其所對應(yīng)的a、b和c，以及校正后各樣品的深度Hi'。

2.2.4 確定最大歸位深度范圍

為了進一步探索最大歸位幅度與TOC擬合效果之間的關(guān)系，將巖心最大歸位幅度確定為0～1m，計算步長為0.125m;巖屑的最大歸位幅度確定為0～2m，計算步長為0.25m。對比計算發(fā)現(xiàn)，當巖心樣品的最大歸位幅度小于0.375m(巖屑的最大歸位幅度確定為0.75m)時，TOC的計算值和實測值之間的相關(guān)系數(shù)隨巖心樣品最大歸位幅度的增大而顯著增大，但是當巖心樣品最大歸位幅度大于0.375m(巖屑的最大歸位幅度為0.75m)后，其擬合相關(guān)系數(shù)提升效果有限，且浮動距離過大將導致模型失真，脫離了地質(zhì)實際。所以確定巖心樣品最大歸位幅度為0.375m(巖屑的最大歸位幅度為0.75m)為該模型樣品的歸位校正幅度 (圖2)。

通過計算對比歸位前后TOC與實測TOC的相關(guān)關(guān)系，發(fā)現(xiàn)歸位后計算TOC和實測TOC的相關(guān)性得到顯著提升，其擬合相關(guān)系數(shù)R2由歸位前的0.59提升到0.79(圖2)，預測效果良好。

圖2 沈352井沙四下亞段油頁巖巖心樣品實測TOC最大歸位幅度與擬合相關(guān)系數(shù)的關(guān)系圖Fig.2 Relationship between the maximum homing amplitude and the fitting correlation coefficient of measure TOC of core oil shale samples in the submermber of the fourth member of Shahejie Formation in Well Shen352

沙四下亞段分為3個層段，這3個層段的巖性存在一定差異，特別是中段以白云質(zhì)泥巖和灰質(zhì)泥巖為主，而上段和下段的巖性以泥質(zhì)油頁巖為主。在進行TOC預測時可能需要分層段擬合。將沙四下亞段分成3段分別擬合各層段TOC的相對含量，從效果來看，分段擬合的效果有所提升，擬合相關(guān)系數(shù)R2由0.75提高到0.767，但效果提升有限，且分段擬合和整體擬合的擬合相關(guān)系數(shù)相當 (表1、圖3)，再加上分段擬合工作量更大，所以分段擬合就變得沒有必要。

表1 不同方法擬合的沈352井沙四下亞段TOC效果及ΔlogR系數(shù)對比表Table 1 Comparison of TOC effect and ΔlogR coefficient of Well Shen352 with different fitting methods

圖3 不同方法擬合沈352井TOC與其實測TOC對比圖Fig.3 Comparison between TOC and measure TOC of Well Shen352 by different fitting methods

3 改進ΔlogR模型的應(yīng)用

3.1 模型驗證

通過建立的歸位整段擬合參數(shù)對沈166井和安1井進行了計算TOC與模擬TOC的相關(guān)性分析。將擬合參數(shù)a、b及c帶入標準化和歸一化的聲波時差和電阻率測井數(shù)據(jù)中，計算了這兩口井沙四下亞段的TOC值，擬合效果發(fā)現(xiàn)，計算TOC值在油頁巖段呈現(xiàn)高值，碳酸鹽巖段的TOC值相對較低，基本符合客觀實際;并且計算TOC與實測TOC值高度相關(guān)，其中沈166井沙四段計算TOC與實測TOC的相關(guān)系數(shù)R2為0.759，而安1井的計算TOC與實測TOC的相關(guān)系數(shù)R2高達0.8191(圖4)。

圖4 實測TOC與擬合TOC對比圖Fig.4 Comparison between measured TOC and fitted TOC

3.2 結(jié)果分析

預測TOC與實測TOC高度相關(guān)，說明模型高度可靠。從圖4來看，沙四下亞段下部 (PS1)有機質(zhì)豐度震蕩變化，為氣候變化比較頻繁的表現(xiàn)，以半干旱—半濕潤氣候為主，沉積水體比較動蕩，在弱氧化—弱還原環(huán)境間頻繁交互 (還原環(huán)境占主導地位)，水體鹽度較高，以微咸水—半咸水環(huán)境為主，沉積暗色油頁巖、泥巖 (微咸水還原環(huán)境產(chǎn)物)，沈352井、沈166井和安1井PS1段有機質(zhì)豐度呈現(xiàn)較強的非均質(zhì)性。沙四下亞段中部 (PS2)TOC含量與上下表現(xiàn)出較大的差異，主要原因可能為水體變淺，呈現(xiàn)出弱氧化環(huán)境特征，蒸發(fā)環(huán)境使水體進一步咸化，以微咸水—半咸水環(huán)境為主，氣候干旱、較強水動力條件和微咸水—半咸水的弱氧化蒸發(fā)環(huán)境條件不利于油頁巖的形成，處于安福屯斷階帶的沈166井此時位于構(gòu)造低洼處，水體相對較深，形成了整體低有機質(zhì)豐度背景下的相對高有機質(zhì)豐度特征。沙四下亞段上部 (PS3)沉積時期整體為相對暖濕氣候，沉積水體加深，有機質(zhì)豐度最高，其TOC平均值高達4.3%。

4 結(jié)束語

(1)對ΔlogR技術(shù)提出了改進，通過TOC=a·logR+b·AC+c的預測模型，以實測TOC與計算TOC之差的平方和作為目標函數(shù)，由巖心樣品在一定范圍內(nèi)自由浮動進行擬合，確定該模型的巖心樣品最大浮動距離為0.375m(巖屑樣品最大浮動距離為0.75m)時最為合理，既保證了歸位徹底，又不會造成數(shù)據(jù)過度歸位而失真。

(2)改進ΔlogR技術(shù)在大民屯凹陷沙四下亞段油頁巖的有機質(zhì)豐度非均質(zhì)性評價中應(yīng)用效果良好，沈166井和安1井沙四下亞段油頁巖計算TOC與實測 TOC值高度相關(guān)，相關(guān)系數(shù) R2分別為0.759和0.8191。

(3)沙四下亞段下部 (PS1)有機質(zhì)豐度表現(xiàn)為較強的非均質(zhì)性;沙四下亞段中部 (PS2)沉積時期，臺地中央氣候干燥炎熱，水體大量蒸發(fā)，水體變淺，氧化性增強，陸源有機質(zhì)供給有限，有機質(zhì)豐度較低;沙四下亞段上部 (PS3)沉積時期轉(zhuǎn)為相對暖濕氣候，水體加深，有機質(zhì)豐度最高，其TOC平均值高達4.3%。

［1］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］朱光有，金強.烴源巖的非均質(zhì)性及其研究——以東營凹陷牛38井為例 ［J］.石油學報，2002，23(5):34－39.

［3］朱光有，金強，張林曄.用測井信息獲取烴源巖的地球化學參數(shù)研究［J］.測井技術(shù)，2003，27(2):104.

［4］劉超.測井資料評價烴源巖方法改進及作用 ［D］.東北石油大學，2011:15－21.

［5］胡慧婷，盧雙舫，劉超，等.測井資料計算源巖有機碳含量模型對比及分析 ［J］.沉積學報，2011，29(6):1199－1204.

［6］湯麗娜.烴源巖測井識別在渤海灣石油開發(fā)中的應(yīng)用［J］.科技導報，2011，29(21):51－54.

［7］韓雙彪，張金川，黃衛(wèi)東，等.吐哈盆地烴源巖有機碳質(zhì)量分數(shù)測井定量預測模型 ［J］.斷塊油氣田，2012，19(3):406－408.

［8］高陽.利用測井信息評價鹽湖相烴源巖 ［J］.沉積學報，2013，31(4):730－737.

［9］袁彩萍，徐思煌，薛羅.珠江口盆地惠州凹陷主力烴源巖測井預測及評價 ［J］.石油實驗地質(zhì)，2014，36(1):110－116.

［10］李建華，宋兵，龐雄奇，等.大民屯凹陷源控油氣作用及有利勘探區(qū)帶預測 ［J］.科技導報，2011，29(30):18－23.

［11］朱華鈺，喻梓靚.大民屯凹陷中部沙四段沉積體系研究 ［J］.石油天然氣學報，2013，35(4):11－16.

［12］謝文彥，姜建群，張占文，等.大民屯凹陷高蠟油成因及分布規(guī)律 ［J］.石油學報，2007，28(2):57－61.

［13］黃海平，李虹，馬刊創(chuàng)，等.大民屯凹陷高蠟油的形成條件［J］.石油與天然氣地質(zhì)，2001，22(1):64－71.

［14］謝文彥，姜建群，張占文，等.大民屯凹陷油氣系統(tǒng)研究［J］.石油勘探與開發(fā)，2004，3(2):38－42.

［15］劉新穎，鄧宏文，邸永香，等.海拉爾盆地烏爾遜凹陷南屯組優(yōu)質(zhì)烴源巖發(fā)育特征 ［J］.石油實驗地質(zhì)，2009(1):68－73.