鄭永建,李家駿,郭海敏,曾桃

(1.中海石油(中國)有限公司湛江分公司,廣東湛江524057;2.長江大學(xué)油氣資源與勘探技術(shù)教育部重點實驗室,湖北武漢430100;3.長江大學(xué)地球物理與石油資源學(xué)院,湖北武漢430100)

0 引 言

碳氧比測井技術(shù)可以確定套管井剩余油飽和度,是油田油藏動態(tài)監(jiān)測的重要手段,并對水淹層的判別及裸眼井中低電阻率油層的識別等問題具有一定指導(dǎo)意義[1]。由于海上油田生產(chǎn)管柱結(jié)構(gòu)復(fù)雜,以斜井為主,碳氧比測井解釋結(jié)果常常存在較大的誤差。同時海上油田缺乏完備的RPM碳氧比測井資料處理解釋軟件,解釋和處理一般由服務(wù)公司完成后提交處理成果,制約了該技術(shù)與其他測井資料的同步結(jié)合,給地質(zhì)、油藏研究帶來了極大不便。因此,本文開展了RPM碳氧比測井資料解釋軟件的研發(fā)工作,并在AECOLog測井解釋軟件平臺上使用Qt開發(fā)框架實現(xiàn)了工業(yè)級的數(shù)據(jù)處理解釋軟件。

Qt是基于C++的跨平臺圖形界面程序開發(fā)軟件,具有跨系統(tǒng)平臺移植性,可提供“一次編寫,隨處編譯”的開發(fā)框架[2-3]?？商岣邔ο筇幚淼撵`活性與多樣性,形成的軟件系統(tǒng)也具有良好的可移植性。因此該解釋系統(tǒng)具有操作方便,處理速度快等特點,在提供準(zhǔn)確的處理成果的同時,也為實際生產(chǎn)使用提供了極大的便利。

1 RPM碳氧比測井資料解釋流程

1.1 曲線預(yù)處理

深度校正。以原始地層自然伽馬曲線為準(zhǔn),RPM-C/O測井前應(yīng)進(jìn)行RPM-PNC測量,用SGFC與原始自然伽馬校深,進(jìn)而使RPM-C/O測量的CCL深度與RPM-PNC測量的CCL深度一致以保證深度的準(zhǔn)確性。

合并曲線。對于碳氧比測井來說,因為測量過程中存在一定的計數(shù)統(tǒng)計漲落誤差,故使得最終測量得到的曲線存在較大的誤差和噪聲。RPM碳氧比測井是多趟多次測量,為了合理的利用數(shù)據(jù)通常是需要對多次測量的曲線進(jìn)行合并濾波處理。

針對不同的曲線質(zhì)量情況,所研發(fā)的RPM測井解釋軟件為用戶提供3種合并處理方法。當(dāng)多趟測量曲線質(zhì)量均較好時,可以使用常見的均值法進(jìn)行曲線的合并處理;當(dāng)多趟測量曲線中有1趟測量或2趟曲線異常時,使用中值法進(jìn)行曲線的合并處理;若所測量的多趟曲線中僅有1趟或2趟曲線質(zhì)量較好,則可以選用測量質(zhì)量好的一條曲線為基準(zhǔn)曲線,使用皮爾遜相關(guān)法根據(jù)閾值選取其他曲線數(shù)據(jù)進(jìn)行均值合并處理?？梢绘I處理長短源距的所有曲線合并處理,提高了數(shù)據(jù)處理的效率。

1.2 影響因素校正

碳氧比測井儀器所發(fā)射的14 MeV的快中子需要穿透儀器的外壁、油管、套管和水泥環(huán)等介質(zhì)才能到達(dá)地層[4]。會受到套管尺寸、水泥環(huán)厚度、泥質(zhì)含量、礦化度、環(huán)空流體和鉆井液侵入等對測量結(jié)果的影響。在實際進(jìn)行解釋處理時,需要對其環(huán)境影響因素進(jìn)行校正。該系統(tǒng)提供了6種影響因素的校正方法。

1.2.1 套管尺寸校正

套管往往影響非彈性散射碳氧比值、俘獲硅鈣比值和非彈性散射鈣硅比值等,從而影響所計算的含水飽和度。為了校正管柱對含水飽和度的影響,在歸一化解釋模型的基礎(chǔ)上得到了含水飽和度計算公式

(1)

式中,C/Oc為經(jīng)校正后非彈性散射碳氧比值,無量綱;C/Ow為100%含水地層非彈性散射碳氧比值,無量綱;ΔC/O為100%含油地層與100%含水地層非彈性散射碳氧比差值,無量綱;M為ΔC/O乘法因子,表征碳氧比曲線對地層中碳與氧原子數(shù)比的靈敏度,隨著套管尺寸變大而變小。

1.2.2 水泥環(huán)厚度校正

水泥環(huán)的厚度往往也會影響含水飽和度的計算,根據(jù)Atlas實驗室刻度6.625 in[注]非法定計量單位,1 ft=12 in=0.304 8 m,下同,套管充淡水得到視含油飽和度Soa與真含油飽和度Sot的圖版,通過線性插值得到不同水泥環(huán)厚度下Soa與Sot的關(guān)系進(jìn)行校正[5]。

圖1 不同水泥環(huán)厚度下視含油飽和度Soa與 真實含油飽和度Sot關(guān)系圖版(Atlas)

如圖1所示,隨著水泥環(huán)厚度的增加,視含油飽和度乘法因子隨之增加。根據(jù)圖版中的函數(shù)關(guān)系式,可以得到視含油飽和度和真實含油飽和度存在關(guān)系式

Sot=k×Soa

(2)

k=0.01091h2+0.076448h+1

(3)

式中,Soa為視含油飽和度,%;Sot為真含油飽和度(經(jīng)水泥環(huán)校正后的含油飽和度),%;k為校正系數(shù)即斜率,無綱量;h為水泥環(huán)的厚度,in。

1.2.3 泥質(zhì)含量校正

在砂泥巖地層剖面中,泥質(zhì)含量中的有機(jī)碳和束縛水中所含的氫元素往往都會對測量的碳氧比值均有影響[6],所研發(fā)的RPM測井解釋軟件提供了2種泥質(zhì)校正的方法。

(1)在分析砂泥巖骨架的碳氧比值和泥質(zhì)含量的關(guān)系時,根據(jù)節(jié)選出的純泥巖層數(shù)據(jù),進(jìn)行多元線性回歸分析得到砂泥巖骨架的碳氧比值和泥質(zhì)含量的關(guān)系[7]

(4)

式中,C/Oma為砂泥巖骨架的非彈性散射碳氧比值,無量綱;Vsh為地層的泥質(zhì)含量值,小數(shù);a為擬合所得到的二次項系數(shù),無量綱;b為擬合所得到的一次項系數(shù),無量綱;c為擬合所得到的常數(shù)項系數(shù),無量綱。

當(dāng)Vsh=0時,即為純砂巖層時,C/Oma=c。

一般當(dāng)Vsh<-b/2a時,碳氧比值均低于純砂巖層的碳氧比值,所以

(5)

將測量的碳氧比值恢復(fù)到純巖性骨架點上,即對碳氧比值進(jìn)行泥質(zhì)含量校正的公式為

C/Oc=C/O+ΔC/O

(6)

式中,C/O為實測的非彈性散射碳氧比曲線值,無量綱。

(2)圖2為改進(jìn)體積模型(雙水模型:可動水和束縛水),考慮到束縛水的影響,得到改進(jìn)體積模型見式(7)

(7)

式中,φwb為束縛水孔隙度;m、n均為校正系數(shù),由實際資料擬合得到。

圖2 儲層組成雙水模型示意圖

1.2.4 礦化度校正

當(dāng)?shù)V化度較高時,俘獲硅鈣比曲線受礦化度影響嚴(yán)重。通過Atlas實驗室確定的校正圖版(見圖3),通過插值計算可以得到不同礦化度下俘獲硅鈣比與孔隙度的關(guān)系,對俘獲硅鈣比進(jìn)行校正[5]。

圖3 不同礦化度下硅鈣比值與孔隙度的關(guān)系圖版(Atlas)

根據(jù)圖版可得硅鈣曲線和孔隙度到在礦化度為0和230 000 mg/L時的關(guān)系

0 mg/L: Si/Ca=0.0014φ+1.309

(8)

230 000 mg/L: Si/Ca=-0.0065φ+1.309

(9)

由此可以刻度得到不同礦化度情況下的礦化度校正公式

(10)

式中,Si/Ca為俘獲鈣硅比值,無量綱;Si/Cac為校正后的俘獲硅鈣比值,無量綱;Ct為地層水的礦化度,×104mg/L;φ為孔隙度,小數(shù)。

1.2.5 環(huán)空流體校正

當(dāng)套管和油管之間的環(huán)空空間中存在環(huán)空流體油時,在測井曲線上會導(dǎo)致碳氧比值呈現(xiàn)“臺階式”增大的現(xiàn)象,碳?xì)浔惹€所受影響較小?；谔佳醣鹊牧黧w靈敏度應(yīng)大于碳?xì)浔?提出利用二者的靈敏度差異反映井眼環(huán)空持率的變化。

(11)

由式(11)推導(dǎo)可得出

ΔS(C/Omax-C/Omin)+C/Omin-

(12)

式中,ΔS為二者的靈敏度差異值,無量綱;C/H為實測非彈性散射的碳計數(shù)率與俘獲的氫計數(shù)率的比值,無量綱;C/Omax為非彈性散射碳氧比曲線的最大值,無量綱;C/Omin為非彈性散射碳氧比曲線的最小值,無量綱;C/Hmax為碳?xì)浔惹€的最大值,無量綱;C/Hmin為碳?xì)浔惹€的最小值,無量綱。

由式(12)可知,碳?xì)浔群吞佳醣惹€存在比較好的線性關(guān)系,假定其如式(13)所示

C/O=a×C/H+b

(13)

則針對于未受環(huán)空流體影響層段不同的深度點應(yīng)存在關(guān)系

(14)

式中,C/On為不同深度點的實測非彈性散射碳氧比值,無量綱;C/Hn為不同深度點的實測非彈性散射的碳計數(shù)率與俘獲的氫計數(shù)率的比值,無量綱;an為不同深度點擬合線性關(guān)系的斜率,無量綱;bn為不同深度點擬合線性關(guān)系的截距,無量綱。

通過解超定方程組得到最優(yōu)解的a和b的值,輸入到解釋參數(shù)中,可通過實測碳?xì)浔惹€計算出一條經(jīng)過環(huán)空校正后的碳氧比曲線,進(jìn)行含水飽和度的計算。

1.2.6 鉆井液侵入校正

由于鉆井過程中存在壓差導(dǎo)致鉆井液濾液侵入地層,改變了周圍井眼的物理特性,從而導(dǎo)致鉆井液侵入前后測得的響應(yīng)值存在差異。

一般SVD法在這里給出了矩陣條件數(shù)小于等于SNR的限制,避免了解的不穩(wěn)定性。根據(jù)體積模型得到侵入前后含水飽和度關(guān)系式并對含水飽和度直接進(jìn)行校正[8],如式(15)所示

(15)

式中,Swc為鉆井液侵入前的地層含水飽和度,%;Sw為鉆井液侵入后的所計算的地層含水飽和度,%;r為侵入深度(可由電阻率曲線校正圖版或者滲流資料計算得到)[9-10],m。

1.3 飽和度計算

經(jīng)過環(huán)境校正后所得到的非彈性散射碳氧比值、俘獲硅鈣比值和非彈性散射鈣硅比值可以計算儲層的含水飽和度(含油飽和度)信息。考慮到環(huán)境影響、地層孔隙度和巖性的變化,該系統(tǒng)提供了3種含水飽和度的計算方法,分別是Si/Ca—C/O交會圖法、Ca/Si—C/O交會圖法和扇形圖版。

1.3.1 Si/Ca—C/O交會圖法

Si/Ca-C/O交會圖法是根據(jù)實驗室的非彈性散射碳氧比值和俘獲硅鈣比值作交會圖版所得到的含油飽和度的解釋公式[5],如式(16)所示

(16)

式中,So為計算的含油飽和度,%;φ為地層的孔隙度,小數(shù);Si/Ca為俘獲硅鈣比值,無量綱;C/O為非彈性散射碳氧比值,無量綱;kw、Lw為分別非彈性散射碳氧比值和俘獲硅鈣比值交會水線的斜率和截距,無量綱。

1.3.2 Ca/Si—C/O交會圖法

Ca/Si—C/O交會圖法同樣也是根據(jù)實驗室的交會圖版所得到的含油飽和度的解釋公式[5],如式(17)所示

(17)

式中,Ca/Si為非彈性散射鈣硅比值,無量綱;kwcs、Lwcs分別為非彈性散射碳氧比值和鈣硅比值交會水線的斜率和截距,無量綱。

1.3.3 扇形圖法

扇形圖版法是通過孔隙度與碳氧比值關(guān)系模擬圖版建立碳氧比解釋的模型,通過100%含油和100%含水2種臨界條件,得到純油地層和純水地層的碳氧比差值ΔC/O。根據(jù)簡單線性解釋方法含水飽和度的歸一化對Sw的計算公式進(jìn)行了改進(jìn)[見式(1)] 。

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

RPM解釋軟件整體分為碳氧比測井預(yù)處理、影響因素校正和碳氧比測井含水飽和度解釋3個部分。其中,預(yù)處理部分完成對測井?dāng)?shù)據(jù)的曲線深度校正、合并濾波等功能;碳氧比測井含水(油)飽和度解釋主要提供環(huán)境因素校正和3種含水(油)飽和度解釋方法。

圖4 碳氧比測井解釋處理成果

圖4是某區(qū)塊A井的碳氧比解釋處理與取心資料對比結(jié)果圖。其中紅色的飽和度曲線(SWCO_S)為利用該RPM軟件解釋處理所得,與Express軟件處理的含水飽和度曲線(SOCO)基本一致,絕對誤差均小于8%,相對誤差小于0.5%。

從取心資料分析,所計算的含水飽和度值與取心巖樣分析基本吻合,符合率達(dá)85%以上。

其中,第1道是井徑、自然伽馬曲線;第2道是電阻率曲線;第3道是經(jīng)本文解釋軟件預(yù)處理后得到的碳氧比和硅鈣比曲線;第4道是深度;第5道是飽和度曲線,其中紅色曲線為該RPM解釋軟件所計算的含水飽和度曲線,綠色曲線為Express處理的含水飽和度曲線,藍(lán)色為裸眼井計算的含水飽和度曲線;第6道是該RPM解釋軟件的含水飽和度曲線和密閉取心資料對比;第7道為碳氧比測井解釋計算的純油純水線;第8道為巖性剖面曲線。

3 結(jié) 論

(1)基于Qt框架編制的RPM碳氧比測井資料解釋系統(tǒng)包括曲線預(yù)處理、影響因素校正和碳氧比測井飽和度計算等模塊,功能齊全,可以滿足油田現(xiàn)場的應(yīng)用需求。

(2)該系統(tǒng)集成于AECOLog測井解釋平臺上,流程簡單、操作方便,便于用戶掌握。處理效果與國外同類處理軟件水平一致,達(dá)到了油田實際生產(chǎn)解釋處理的需要。