易覺(jué)非

(長(zhǎng)江大學(xué)信息與數(shù)學(xué)學(xué)院，湖北 荊州 434023)

基于常規(guī)測(cè)井曲線頻譜分析的多尺度地質(zhì)分層

易覺(jué)非

(長(zhǎng)江大學(xué)信息與數(shù)學(xué)學(xué)院，湖北 荊州 434023)

測(cè)井曲線的頻譜具有其特殊的地質(zhì)意義，即測(cè)井曲線的不同頻率與不同尺度地層相對(duì)應(yīng)。首先討論了測(cè)井曲線頻率與地層的關(guān)系，進(jìn)一步給出了用帶通濾波進(jìn)行多尺度地質(zhì)分層的方法，最后，通過(guò)一個(gè)實(shí)際分層實(shí)例，說(shuō)明了方法的正確性和有效性。

測(cè)井曲線;頻譜;帶通濾波

在信號(hào)分析中，有時(shí)域分析方法和頻域分析方法2類(lèi)，有些信號(hào)特征在時(shí)域中難以獲得，在頻域中則容易獲得，反之亦然。因此，時(shí)域分析方法和頻域分析方法兩者缺一不可，互為補(bǔ)充。測(cè)井曲線的分析也不例外，可以在深度域(即時(shí)域)中，分析其依深度變化的規(guī)律，當(dāng)將它轉(zhuǎn)換到頻域時(shí)，則可以分析其結(jié)構(gòu)特征，發(fā)現(xiàn)其頻率對(duì)應(yīng)的地質(zhì)意義。比如用頻譜分析方法來(lái)確定沉積旋回和層序地層[1,2]。在利用測(cè)井資料進(jìn)行地質(zhì)分層方面，筆者在文獻(xiàn)[3]中詳細(xì)闡述了活度分層法的原理和算法，該方法能夠快速有效的確定地層界面，但不能反映地層的尺度變化規(guī)律。下面，筆者從頻譜分析的基本原理出發(fā)，結(jié)合測(cè)井曲線的頻譜特征及其地質(zhì)意義，給出了利用測(cè)井曲線頻譜進(jìn)行多尺度地質(zhì)分層的方法。

1 基本原理

常規(guī)測(cè)井曲線是一系列不同厚度地層的數(shù)字化，因此，可將它理解為一系列寬度不同的矩形脈沖的疊加，圖1為塔里木盆地哈得地區(qū)某井的自然伽瑪曲線，體現(xiàn)了這一解釋的合理性。

圖1 近似矩形脈沖的自然伽瑪曲線

依據(jù)矩形脈沖及其頻譜響應(yīng)，地層與常規(guī)測(cè)井曲線的頻譜應(yīng)具有如下關(guān)系：

1)厚地層譜幅高，頻率衰減快，薄地層則相反；

2)頻率與層厚成反比例關(guān)系，厚地層頻率低，薄地層則相反。

因此，可以認(rèn)為不同厚度地層對(duì)應(yīng)不同頻率范圍，一般來(lái)說(shuō)，較厚地層占據(jù)低頻段，較薄地層占據(jù)高頻段。依據(jù)這一原理，筆者提出了采用測(cè)井曲線帶通濾波進(jìn)行多尺度地質(zhì)分層的方法。

2 主要算法

實(shí)現(xiàn)測(cè)井曲線多尺度地質(zhì)分層的主要算法流程如圖2所示。

圖2 算法主要流程

1)傅里葉變換 輸入為常規(guī)測(cè)井曲線，輸出為常規(guī)測(cè)井曲線的頻譜，可采用基2快速傅里葉變換(FFT:Fast Fourier Transform)算法[4]。

2)帶通濾波 輸入為測(cè)井曲線頻譜，輸出為經(jīng)過(guò)濾波器濾波后的測(cè)井曲線頻譜。這一步是實(shí)現(xiàn)分層的關(guān)鍵，濾波器設(shè)計(jì)[5]不合理，則不能體現(xiàn)地層特征，影響分層結(jié)果。主要從以下2個(gè)方面考慮濾波器的設(shè)計(jì)：濾波器的頻帶范圍[flow,fhigh]。在此引入空間頻率的概念，物理意義為單位深度內(nèi)常規(guī)測(cè)井曲線波形變化次數(shù)，單位為次/m，記為1/m。其中，flow為低截頻；fhigh為高截頻，該頻帶既不能太寬(以防止不同層級(jí)的地層頻率相互混淆),又不能太窄(以防止同一尺度地層頻率遺漏)。②[flow,fhigh]之間的優(yōu)勢(shì)頻率fext。它是同尺度地層的主要頻率。

圖3 三角形帶通濾波器

3)傅里葉反變換 輸入為經(jīng)過(guò)濾波器濾波后的測(cè)井曲線頻譜，輸出為濾波后的測(cè)井曲線。

4)確定地層界面 輸入為濾波后的測(cè)井曲線，輸出為地層界面點(diǎn)。地層界面點(diǎn)可以采用活度分層法[3]計(jì)算。

5)分層結(jié)束 分層結(jié)束條件設(shè)定為測(cè)井曲線的截止頻率fmax，若深度采樣間隔τ=0.125m，截至頻率可設(shè)為fmax=4(1/m)。如果濾波器頻帶范圍[flow,fhigh]中fhighlt;fmax，則重復(fù)2)、3)、4)步，最終得到不同尺度的地層界面。

3 分層實(shí)例

圖4 某井自然伽瑪曲線頻譜

常規(guī)測(cè)井曲線中，自然伽瑪曲線最能反映泥質(zhì)含量的變化，筆者選用塔里木盆地哈得地區(qū)某井4970～5098m一段自然伽瑪曲線進(jìn)行分層。 采樣點(diǎn)數(shù)N=1024，其深度采樣間隔τ=0.125m，截至頻率為4(1/m)。

首先，對(duì)該段測(cè)井曲線進(jìn)行傅里葉變換,得到其頻譜,如圖4所示。其頻譜曲線具有以下特點(diǎn):①低頻部分譜幅大,高頻部分譜幅小,表明測(cè)井曲線中低頻部分對(duì)應(yīng)的厚層信息較豐富;②低頻段譜幅衰減快,而隨頻率的增高譜幅衰減速度呈下降趨勢(shì)。筆者通過(guò)對(duì)同一地區(qū)多口井自然伽瑪曲線進(jìn)行頻譜分析,發(fā)現(xiàn)其頻譜曲線具有較大的相似性，這說(shuō)明在該地區(qū)地層分布比較穩(wěn)定。

對(duì)上述自然伽瑪曲線用2個(gè)圖3所示三角形帶通濾波器濾波，濾波器參數(shù)分別取flow=0.002，fmax=0.02和flow=0.02,fmax=0.2，得到圖5所示2條帶通濾波曲線，其中頻率0.002～0.02(1/m)的濾波曲線對(duì)應(yīng)厚地層，頻率0.02～0. 2(1/m)的濾波曲線對(duì)應(yīng)薄地層，進(jìn)一步采用活度分層法分別對(duì)這2條曲線進(jìn)行地層界面提取，提取結(jié)果如圖5所示。分層結(jié)果與實(shí)際結(jié)果基本一致，體現(xiàn)了計(jì)算的有效性。

圖5 對(duì)帶通濾波曲線進(jìn)行分層

4 結(jié) 語(yǔ)

利用測(cè)井曲線頻譜分析對(duì)地層進(jìn)行研究，符合地層周期性沉積規(guī)律。通過(guò)帶通濾波器濾波，可以得到不同尺度的地層旋回特征，從而實(shí)現(xiàn)多尺度地質(zhì)分層。實(shí)例研究表明，采用該方法產(chǎn)生的分層結(jié)果與實(shí)際情況基本吻合。

[1]王貴文,鄧清平,唐為清.測(cè)井曲線譜分析方法及其在沉積旋回研究中的應(yīng)用[J].石油勘探與開(kāi)發(fā), 2002,29(1):93～95.

[2]王紅罡,呂炳全,孫小剛.長(zhǎng)慶氣田儲(chǔ)層的高頻層序地層學(xué)分析[J].同濟(jì)大學(xué)學(xué)報(bào),2000, 28(5): 543～546.

[3]易覺(jué)非.利用活度分層法實(shí)現(xiàn)測(cè)井自動(dòng)地質(zhì)分層[J].石油天然氣學(xué)報(bào)(江漢石油學(xué)院學(xué)報(bào)),2007,29(1):78～80.

[4]胡廣書(shū).數(shù)字信號(hào)處理[M].北京：清華大學(xué)出版社,1997.55～174.

[5]Sophocles J.Orfanidis.Introduction to Signal Processing[M]. 北京：清華大學(xué)出版社,1999.541～572.

[編輯] 易國(guó)華

P631.84

A

1673-1409(2009)03-N060-03

2009-06-04

易覺(jué)非(1966-)，男，1987年大學(xué)畢業(yè)，講師，現(xiàn)主要從事應(yīng)用數(shù)學(xué)和軟件開(kāi)發(fā)方面的教學(xué)和研究工作。