劉春雷，崔式濤 （中國(guó)石油集團(tuán)測(cè)井有限公司油氣評(píng)價(jià)中心，陜西 西安710077）

景歡 （中國(guó)石油集團(tuán)測(cè)井有限公司長(zhǎng)慶事業(yè)部，陜西 西安710200）

薛猛 （中石油大慶油田有限責(zé)任公司第一采油廠，黑龍江 大慶163001）

曹鵬 （中石油吉林油田分公司地球物理勘探研究院，吉林 松原138000）

胡高賢 （中國(guó)石油集團(tuán)測(cè)井有限公司油氣評(píng)價(jià)中心，陜西 西安710077）

雅布賴盆地位于河西走廊中段，是在前古生界變質(zhì)巖結(jié)晶基底和上古生界褶皺基底上發(fā)育起來的中－新生代斷陷沉積盆地［1］。盆地總體西窄東寬，呈北東東向展布，南以北大山為界與潮水盆地相連，東以巴彥烏拉山為界與巴音浩特盆地相隔，西北以雅布賴山為界與銀額盆地相鄰，面積約為1．5×104km2。侏羅系是主要的烴源巖和勘探目的層［2，3］。

總有機(jī)碳質(zhì)量分?jǐn)?shù)（w （TOC））是雅布賴盆地的烴源巖評(píng)價(jià)研究中的一個(gè)基礎(chǔ)指標(biāo)，可以用來識(shí)別生油巖，指示有機(jī)質(zhì)豐度，判斷生油效率、轉(zhuǎn)化效率和演變程度，估計(jì)生油量，推測(cè)石油初次運(yùn)移方向等［4～7］。通過巖心的有機(jī)碳熱解分析、鏡質(zhì)體反射率等試驗(yàn)分析，可以得到較準(zhǔn)確的數(shù)據(jù)。但巖心分析試驗(yàn)昂貴且費(fèi)時(shí)，因此通過測(cè)井分析獲取適合研究區(qū)w （TOC）經(jīng)驗(yàn)公式是非常必要的。目前，經(jīng)典的w （TOC）算法是ΔlgR 法和多元回歸法，有效彌補(bǔ)了沒有巖心試驗(yàn)分析情況下的w （TOC）確定。

1 經(jīng)典ΔlgR 法確定w （TOC）

ΔlogR 法的基本原理是利用油層、蒸發(fā)巖、火成巖、欠壓實(shí)沉積物和井壁垮塌嚴(yán)重層段在自然伽馬或者自然電位曲線上的不同特點(diǎn)進(jìn)行辨別并剔除，然后選擇適合的刻度范圍，將三孔隙度曲線與電阻率曲線重疊，在含少量有機(jī)質(zhì)的層段，2條曲線相互重合或平行；在含豐富有機(jī)質(zhì)層段，兩者分離［8］。低密度和低速度 （高聲波時(shí)差）的干酪根響應(yīng)是其原因，在未成熟的富含有機(jī)質(zhì)的巖石中還沒有生成油氣，觀測(cè)到的2條曲線之間的差異僅僅是由孔隙度曲線響應(yīng)造成的；而在成熟的烴源巖中，除了孔隙度曲線響應(yīng)外，由于烴類的存在，電阻率增加，使2條曲線產(chǎn)生更大的差異 （或稱間距）［9～11］。ΔlgR 法是一種測(cè)井計(jì)算w （TOC）的經(jīng)典算法，由Passey等［12～16］經(jīng)過試驗(yàn)分析提出：

式中：ΔlgR 為孔隙度曲線與電阻率曲線顯示的包絡(luò)線面積，1；ρp 和ρb 分別為目的層和非烴源巖（基線）電阻率，Ω·m；φp 和φb 分別為目的層和非烴源巖（基線）孔隙度，1；Ilom為反映有機(jī)質(zhì)成熟度的指數(shù)，與鏡質(zhì)體反射率相關(guān)，1；w （TOC）0為非烴源巖的w （TOC），1；K 為互溶刻度比例系數(shù)，1。

該算法的主要缺點(diǎn)包括：①ρb 和φb 難以準(zhǔn)確確定，由于致密油氣儲(chǔ)層往往含有烴源巖，如簡(jiǎn)單地認(rèn)為儲(chǔ)層即為基線容易導(dǎo)致計(jì)算的w （TOC）偏低；②Ilom不易準(zhǔn)確確定，該參數(shù)對(duì)w （TOC）的計(jì)算精度起重要的控制作用，要結(jié)合烴源巖成熟度來確定；③式 （1）中2．297、－0．1688等數(shù)值是特定地區(qū)的參數(shù)，在其他地區(qū)計(jì)算效果不好。

2 多元線性回歸分析確定w （TOC）

w （TOC）常常與多個(gè)因素相聯(lián)系，由多個(gè)自變量的最優(yōu)組合來預(yù)測(cè)或估計(jì)因變量，比只用一個(gè)自變量進(jìn)行預(yù)測(cè)或估計(jì)更有效，更符合實(shí)際。因此，多元線性回歸比一元線性回歸的實(shí)用意義更大。通過對(duì)研究區(qū)的自然伽馬 （qAPI）、自然電位 （Usp）、電阻率 （ρ）、聲波時(shí)差 （Δt）、補(bǔ)償密度 （ρc）、補(bǔ)償中子孔隙度 （φnc）以及ΔlgR 等變量與w （TOC）建立關(guān)系，優(yōu)選出與w （TOC）關(guān)系更密切的qAPI、ρ、Δt、ρc 以及ΔlgR 等變量進(jìn)行詳細(xì)研究。

2．1 ΔlgR 法結(jié)合常規(guī)測(cè)井的w （TOC）計(jì)算模型

經(jīng)過對(duì)測(cè)井資料的相關(guān)性分析，qAPI、Δt及ΔlgR 與w （TOC）的回歸方程為：

如圖1所示，計(jì)算求取的w （TOC）與巖心試驗(yàn)分析得到的w （TOC）相關(guān)性一般：

圖1 巖心w （TOC）與式 （3）計(jì)算w （TOC）交會(huì)圖

式中：x 為巖心w （TOC），1；y 為計(jì)算w （TOC），1；R2為相關(guān)系數(shù)，下同。

2．2 基于常規(guī)測(cè)井的w （TOC）計(jì)算模型

經(jīng)過對(duì)測(cè)井資料的相關(guān)性分析，qAPI、ρ、Δt及ρc 與w （TOC）的回歸方程為：

式中：ρlld2為深雙側(cè)向電阻率，Ω·m。

圖2 巖心w （TOC）與式 （4）計(jì)算w （TOC）交會(huì)圖

如圖2所示，計(jì)算求取的w （TOC）與巖心試驗(yàn)分析得到的w （TOC）相關(guān)性一般：

3 改進(jìn)的ΔlgR 法確定有機(jī)碳含量

針對(duì)經(jīng)典ΔlogR 法的缺點(diǎn)及多元線性回歸法精度不夠的問題，筆者提出了改進(jìn)的ΔlogR 法方法。

1）重新確定基線，不將非烴源巖段測(cè)井值定為基線，而是將電阻率與孔隙度測(cè)井道左刻度值定為基線計(jì)算ΔlgR，有效避免了基線確定不準(zhǔn)確帶來的誤差。

式中：ρr 為電阻率右刻度值，Ω·m；ρl 為電阻率左刻度值，Ω·m；Δtl為聲波時(shí)差左刻度值，μs／m；Δtr為聲波時(shí)差右刻度值，μs／m；Δtb為聲波時(shí)差基線值，μs／m。

2）利用巖心試驗(yàn)分析的w （TOC）計(jì)算得到 （w （TOC）－w （TOC）0）／ΔlgR，再與巖心分析鏡質(zhì)體反射率 （Ro）建立擬合關(guān)系：

3）Ro與深度（d）建立擬合關(guān)系：

4）最后得到研究區(qū)的w （TOC）回歸方程：

如圖3 所示，改進(jìn)的ΔlogR 法計(jì)算求取w （TOC）與巖心試驗(yàn)分析得到的w （TOC）相關(guān)性較高：

圖3 巖心w （TOC）與式 （10）計(jì)算w （TOC）交會(huì)圖

4 實(shí)例分析

由雅探－X 井改進(jìn)的ΔlgR 法與其他方法的效果對(duì)比情況 （圖4）可見，改進(jìn)的ΔlgR 法計(jì)算的w （TOC）精度最高，ΔlgR 法＋常規(guī)測(cè)井共同回歸方法的效果要比常規(guī)測(cè)井的回歸方法要好，經(jīng)典的ΔlgR 法，由于其參數(shù)的不合適，計(jì)算誤差較大。

5 結(jié)論

致密油研究中，w （TOC）定量評(píng)價(jià)是很重要的一部分，但經(jīng)典的w （TOC）計(jì)算方法ΔlgR 法在雅布賴盆地侏羅系計(jì)算效果與巖心差距很大。筆者利用巖心刻度測(cè)井重新構(gòu)建w （TOC）評(píng)價(jià)方法，并與其他多元回歸方法進(jìn)行對(duì)比，改進(jìn)的ΔlgR 法效果最好。在沒有巖心資料情況下，改進(jìn)的ΔlgR 法可以利用深度、電阻率以及聲波時(shí)差等基本測(cè)井資料，快速、準(zhǔn)確地評(píng)價(jià)w （TOC），為雅布賴盆地致密油評(píng)價(jià)提供了良好的基礎(chǔ)。

圖4 巖心w （TOC）與4種方法計(jì)算w （TOC）的檢驗(yàn)測(cè)井圖

［1］吳茂炳，劉春燕，鄭孟．內(nèi)蒙古西部雅布賴盆地侏羅紀(jì)沉積－構(gòu)造演化及油氣勘探方向 ［J］．地質(zhì)通報(bào)，2007，26 （7）：857～863．

［2］甘肅省地質(zhì)礦產(chǎn)局．甘肅省區(qū)域地質(zhì)志 ［M］．北京：地質(zhì)出版社，1989：244～260．

［3］王廷印，王金榮，王士政．阿拉善古生代陸殼的形成和演化 ［M］．蘭州：蘭州大學(xué)出版社，1994：65～71．

［4］周總瑛．烴源巖演化中有機(jī)碳質(zhì)量與含量變化定量分析 ［J］．石油勘探與開發(fā)，2009，3 （4）：463～468．

［5］朱光有，金強(qiáng)，張林曄．用測(cè)井信息獲取烴源巖的地球化學(xué)參數(shù)研究 ［J］．測(cè)井技術(shù)，2003，27 （2）：104～109．

［6］于翔濤．測(cè)井技術(shù)在烴源巖評(píng)價(jià)中的應(yīng)用 ［J］．長(zhǎng)江大學(xué)學(xué)報(bào) （自然科學(xué)版），2009，6 （2）：198～200．

［7］Schmoker J W，Hester T C．Organic carbon in Bakken Formation United States Portion of Williston Basin ［J］．AAPG Bulletin，1983，67 （12）：2165～2174．

［8］賀小林，王振奇，周鳳娟，等．尼日爾三角洲烴源巖評(píng)價(jià)研究 ［J］．長(zhǎng)江大學(xué)學(xué)報(bào) （自然科學(xué)版），2011，8 （3）：38～40．

［9］陳曉娜，文志剛，宋換新．非洲乍得Bongor盆地西部凹陷烴源巖評(píng)價(jià) ［J］．長(zhǎng)江大學(xué)學(xué)報(bào) （自然科學(xué)版），2011，8 （8）：20～23．

［10］張志偉，張龍海．測(cè)井評(píng)價(jià)烴源巖的方法及其應(yīng)用效果 ［J］．石油勘探與開發(fā)，2000，27 （3）：84～87．

［11］王貴文，朱振宇，朱廣宇．烴源巖測(cè)井識(shí)別與評(píng)價(jià)方法研究 ［J］．石油勘探與開發(fā)，2002，29 （4）：50～52．

［12］趙政璋，杜金虎．致密油氣 ［M］．北京：石油工業(yè)出版社，2012：118～120．

［13］王方雄，候英姿，夏季．烴源巖測(cè)井評(píng)價(jià)新進(jìn)展 ［J］．測(cè)井技術(shù)，2002，26 （2）：89～93．

［14］Lang W H．The determination of thermal maturity in potential source rocks using interval transit time interval velocity ［J］．Log Analyst，1994，35 （7）：47～59．

［15］趙彥超，馬正，姚光慶．Waxman－Smtih方程在生油巖評(píng)價(jià)中的應(yīng)用：重疊法和雙孔隙度法 ［J］．地球科學(xué) （中國(guó)地質(zhì)大學(xué)學(xué)報(bào)），1995，20 （3）：306～313．