李家駿,郭海敏,郭常偉,劉珈辰

(1.長江大學(xué)油氣資源與勘探技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室, 湖北 武漢 430100;2.長江大學(xué)地球物理與石油資源學(xué)院, 湖北 武漢 430100;3.中國石油集團(tuán)測井有限公司華北分公司, 河北 任丘 062550)

0 引 言

TNIS熱中子成像測井系統(tǒng)(Thermal Neutron Imaging System),簡稱熱中子成像系統(tǒng),是加拿大GPN公司開發(fā)的一種用于油田新老井油氣飽和度的測井系統(tǒng)[1]。其測量原理是中子發(fā)生器向地層發(fā)射高能的中子脈沖,并由近、遠(yuǎn)熱中子探頭探測快中子經(jīng)過井眼及地層減速以后還沒有被俘獲的熱中子。在熱中子俘獲能力較低的地層仍可以保持相對高的熱中子計(jì)數(shù)率,在一定程度上改進(jìn)了傳統(tǒng)中子壽命測井的局限性[2]。在衰減過程中形成的熱中子衰減譜和熱中子成像譜可直觀快速地定性判定油水性質(zhì)。TNIS不但在中高孔隙度和礦化度地層有很好的適應(yīng)性,在低孔隙度、低礦化度地層(目前大多數(shù)油田生產(chǎn)的難點(diǎn))相對其他測井方式也具有更高的分辨率[3]。

TNIS在中國油田中應(yīng)用廣泛,在吐哈油田多地區(qū)多井次的測試,證實(shí)其在孔隙度大于10%,氯離子濃度大于2.4×104mg/L的低滲透率油藏也可分辨油水層;在吐哈油田雁木西油田,平均孔隙度為22%,總礦化度11.2×104mg/L,Cl-離子濃度(7～9)×104mg/L的地層中,應(yīng)用效果顯著[4]。大慶油田利用該技術(shù)在水驅(qū)層系井網(wǎng)調(diào)整試驗(yàn)區(qū)進(jìn)行了9口井的剩余油監(jiān)測,并成功指導(dǎo)了8項(xiàng)配套措施調(diào)整,取得較好的開發(fā)效果[5]。

然而,中國油田缺乏與之配套的TNIS測井資料處理解釋軟件。室內(nèi)的數(shù)據(jù)處理解釋主要依靠國內(nèi)GPN的服務(wù)公司進(jìn)行解釋處理,制約了油田解釋人員對資料的利用和解釋精度。在分析TNIS熱中子成像測井的原理和處理方法后,在LEAD3.5測井解釋平臺(tái)上開發(fā)了一套TNIS熱中子成像測井的數(shù)據(jù)處理與解釋系統(tǒng)。

該系統(tǒng)的算法采用C++語言編寫,基于LEAD測井系統(tǒng)底層平臺(tái),能夠在主流的Windows平臺(tái)上運(yùn)行。具有操作方便,處理數(shù)據(jù)快速等優(yōu)點(diǎn),與GPN公司所提供的解釋成果和實(shí)際采油資料相對比表明,其處理效果已經(jīng)達(dá)到國外公司軟件的同等水平。

1 TNIS熱中子測井?dāng)?shù)據(jù)處理方法

1.1 TNI格式數(shù)據(jù)解編

GPN公司的TNIS測井儀記錄了一段特定時(shí)間內(nèi)未被俘獲的熱中子數(shù),并由此提取出地層俘獲截面進(jìn)而計(jì)算得到地層含油飽和度。其記錄的數(shù)據(jù)為*.TNI格式,它是一種二進(jìn)制文件,為了更好地應(yīng)用TNIS測井?dāng)?shù)據(jù),需要對其數(shù)據(jù)進(jìn)行數(shù)據(jù)解編和格式轉(zhuǎn)換。

TNIS測井儀器記錄從中子發(fā)生器發(fā)射15 μs后的700 μm時(shí)間內(nèi)的熱中子記數(shù)率,并將熱中子衰減譜記錄分成180道[4],其數(shù)據(jù)的邏輯結(jié)構(gòu)見圖1。

圖1 TNIS測井?dāng)?shù)據(jù)邏輯結(jié)構(gòu)圖

經(jīng)過數(shù)據(jù)解編驗(yàn)證,其數(shù)據(jù)的字節(jié)信息見表1。長、短源距計(jì)數(shù)率、自然伽馬曲線、比率曲線等信息分別以不同字節(jié)寬度共計(jì)256 B存放于一個(gè)數(shù)據(jù)記錄中。

表1 TNIS數(shù)據(jù)格式整體信息

為了準(zhǔn)確求取地層宏觀俘獲截面(Σ)需要準(zhǔn)確提取地層信息,確定地層區(qū)對熱中子計(jì)數(shù)衰減的貢獻(xiàn)。將得到的熱中子成像譜采取了多次滑動(dòng)計(jì)算Σ,最后取滑動(dòng)計(jì)算的Σ平均值作為最終提取的Σ。

圖2 宏觀俘獲截面求取方法

由于熱中子衰減譜是采用各自相應(yīng)的時(shí)間間隔的計(jì)數(shù)方式,根據(jù)計(jì)算點(diǎn)在譜線上滑動(dòng)求平均的方法得到宏觀俘獲截面。一般選取36個(gè)計(jì)數(shù)點(diǎn),時(shí)間間隔為540 μs,把數(shù)據(jù)分成6組,每組6個(gè)點(diǎn)的計(jì)數(shù)分別表示為N0,N1,N2,N3,N4,N5,計(jì)算地層熱中子壽命為

τ=(τ2-t1)/(lnN1-lnN2)

(1)

實(shí)際計(jì)算取作

(2)

每組求得的宏觀俘獲截面為

Σj=4550/τj

(3)

(4)

數(shù)據(jù)解編軟件,可以自動(dòng)提取自然伽馬曲線、熱中子俘獲譜和熱中子衰減譜等信息,且對求取的Σ曲線采用并高斯12點(diǎn)濾波進(jìn)行去噪。根據(jù)操作人員的需要選取合適的曲線導(dǎo)出,并將其轉(zhuǎn)化為LEAD和中國軟件均可以識(shí)別的WIS格式數(shù)據(jù)。

1.2 TNIS數(shù)據(jù)解釋

原始數(shù)據(jù)解編后得到的地層俘獲截面曲線將進(jìn)行解釋處理,處理流程如圖3所示。

圖3 TNIS數(shù)據(jù)解釋處理流程

該解釋系統(tǒng)編寫的解釋方法共5種:簡單圖版法、增強(qiáng)圖版法、Hingle圖版法[6]、體積模型以及雙因子校正模型。解釋方法的選擇需要根據(jù)實(shí)際的井況進(jìn)行選取。簡單圖版法主要適用于非均質(zhì)性弱、孔隙度較高的地層,未考慮泥質(zhì)含量影響;增強(qiáng)圖版法是基于簡單圖版法的改進(jìn);Hingle圖版法需要用地層的電阻率Rt和地層俘獲截面值Σlog作交會(huì)圖,但裸眼井電阻率曲線不能反映當(dāng)前的地層含油飽和度情況,一般只能在檢查井中應(yīng)用。

圖4為儲(chǔ)層的體積模型。組成地層的介質(zhì)泥質(zhì)、骨架和地層的流體分別都有不同的地層俘獲截面[7]。

圖4 儲(chǔ)層體積模型

這些地層介質(zhì)的綜合俘獲截面曲線測井響應(yīng)就是整個(gè)儲(chǔ)層的俘獲截面曲線即實(shí)際所測地層的俘獲截面值,計(jì)算公式為

Σlog=Σma(1-φ-Vsh)+ΣwφSw+

Σhφ(1-Sw)+ΣshVsh

(5)

由公式(1)推導(dǎo)得到含水飽和度計(jì)算的公式

(6)

式中,Σlog、Σma、Σsh、Σw、Σh分別為地層、巖石骨架、泥質(zhì)、地層水及烴的熱中子宏觀俘獲截面,c.u.[注]非法定計(jì)量單位,1 c.u.=10-3 cm-1,下同;φ為地層孔隙度,小數(shù);Vsh為地層泥質(zhì)體積含量,小數(shù);Sw為地層含水飽和度,小數(shù)。

3種圖版法是采用實(shí)際的井?dāng)?shù)據(jù)進(jìn)行數(shù)據(jù)交會(huì),得到Σma、Σw、Σh這3個(gè)基本的解釋參數(shù),為后續(xù)的解釋作準(zhǔn)備。圖5為增強(qiáng)圖版法獲取參數(shù)的界面,其中,x軸為地層俘獲截面值,y軸為地層的孔隙度,數(shù)據(jù)點(diǎn)為不同深度點(diǎn)的數(shù)據(jù),顏色代表含水飽和度的大小,如圖5中的色標(biāo)所示。紅色線表示為油線(即含水飽和度為0),藍(lán)色線為水線(即含水飽和度為100%),黑色線為含水飽和度為50%的線,通過點(diǎn)擊界面下方水、烴類或骨架后,在交會(huì)圖中用鼠標(biāo)拖動(dòng)油線或水線使不同含水飽和度的數(shù)據(jù)點(diǎn)落在相應(yīng)的區(qū)域,以確定水、烴類和骨架的俘獲截面值。

圖5 增強(qiáng)圖版法參數(shù)獲取界面

該模型只適用于地層水礦化度較高的地層,只有地層水和油氣的俘獲截面值有區(qū)分度,才可以保證所計(jì)算的含水飽和度的精度。根據(jù)劉珈辰等[7]TNIS測井資料中參數(shù)敏感性的研究和大量的油井處理表明,在礦化度較低和泥質(zhì)含量較高的區(qū)域,用體積模型法所計(jì)算的含水飽和度偏高。為此,胡冰恒等[2]引入了具有區(qū)域特征的雙因子校正體積模型[見式(7)和式(8)]

Σlog=Σma(1-φ-Vsh)+K1ΣshVsh+

K2ΣwφSw+Σh(1-Sw)φ

(7)

Sw=

(8)

如式(8)所示,為了求取含水飽和度,除了體積模型中提到的6個(gè)參數(shù)外,還需確定K1和K2這2個(gè)區(qū)域特征參數(shù)。其需要針對目標(biāo)井段,結(jié)合裸眼井解釋和實(shí)際地質(zhì)資料,選取至少2段純水層數(shù)據(jù)代入,通過最優(yōu)化解出適合解釋井段的特征參數(shù)。

針對于參數(shù)泥質(zhì)俘獲截面參數(shù)的獲取,可以根據(jù)井所在層位的實(shí)際情況而定,從其所測的純厚層泥巖段來讀其測量曲線值。還可以距離該層位較近或者該研究層位的別的區(qū)塊的分布相對穩(wěn)定的泥巖段進(jìn)行讀值,并進(jìn)行俘獲截面統(tǒng)計(jì)分析,得出該段的泥巖俘獲截面。

在實(shí)際應(yīng)用中還可以利用泥巖段油水線相等的方法確定較為準(zhǔn)確的泥質(zhì)俘獲截面值[8],定義油線為孔隙內(nèi)100%含油,則純油線為

Σo,l=Σma(1-φ-Vsh)+Σhφ+ΣshVsh

(9)

定義水線為孔隙100%含水,則水線為

Σw,l=Σma(1-φ-Vsh)+Σwφ+ΣshVsh

(10)

調(diào)整Σsh,在泥巖段使油線Σo,l、水線Σw,l和測井曲線Σlog基本重合,從而得到泥質(zhì)的俘獲截面值。

2 解釋與應(yīng)用效果

A井所處區(qū)塊的均礦化度為1×104mg/L,平均孔隙度為16%,因此,采用傳統(tǒng)體積模型計(jì)算含水飽和度。TNIS解釋系統(tǒng)處理得到的熱中子俘獲譜、熱中子衰減譜和含水飽和度與GPN公司解釋軟件提供的基本一致。

圖6為2種解釋系統(tǒng)獲得的HB油田X區(qū)塊A井的解釋成果圖。圖6(a)第1道是巖性剖面;第2道是巖性指示曲線;第3道是俘獲截面曲線;第4道是為飽和度曲線;第5道為熱中子俘獲譜;第6道為熱中子衰減譜;第7道為熱中子俘獲成像;第8道為熱中子衰減成像;第9道為井周成像;第10道為俘獲截面曲線,其中紫色為TNIS解釋系統(tǒng)處理所得曲線,紅色為GPN公司處理得到的曲線。

圖6(b)為GPN公司軟件解釋成果圖。第1道是巖性剖面;第2道是校深曲線;第3道是為飽和度曲線;第4道是為孔隙度曲線;第5道為熱中子俘獲譜;第6道為熱中子衰減譜;第7道為熱中子俘獲成像;第8道為井周成像。

2個(gè)系統(tǒng)所計(jì)算的含油飽和度誤差均小于4%(見表2)。

圖6 TNIS解釋系統(tǒng)與GPN公司軟件獲得的解釋成果圖對比

層號頂界/m底界/mGPN含油飽和度/%TNIS含油飽和度/%絕對誤差/%7-12436.02436.83.763.840.087-22436.82441.65.26.10.98-12444.52445.616.9614.832.138-22445.62446.80.10.130.038-32446.82448.419.5717.781.798-42448.42457.10.241.331.0992482.12484.65.24.750.45102489.72491.00.051.31.2511-12495.32499.02.311.870.4412-32531.12532.700012-42532.72534.34.365.941.58142559.82563.610.211.080.8815-22567.62568.600017-12581.62582.600017-22582.62584.415.3213.361.96182586.82591.335.4131.93.51192592.32595.917.1816.460.72202601.72603.435.8539.563.71

3 結(jié) 論

(1)基于LEAD測井解釋平臺(tái)底層開發(fā)的TNIS熱中子成像測井資料系統(tǒng),可以從TNIS測井的原始TNI數(shù)據(jù)解編處理到含油飽和度計(jì)算劃分油水層,功能齊全,流程簡單、操作方便,可以滿足油田現(xiàn)場的應(yīng)用需求。

(2)處理效果與國外處理軟件水平一致,在低礦化度高泥質(zhì)含量地層采用雙因子校正模型處理,與試油結(jié)論吻合,達(dá)到了油田實(shí)際生產(chǎn)解釋處理的需要。