本發(fā)明屬于石油地質(zhì)勘探領(lǐng)域，具體涉及一種碎屑巖儲(chǔ)層中地震屬性砂體空間分布描述及評(píng)價(jià)的方法。

背景技術(shù)：

在碎屑巖儲(chǔ)層的油氣勘探開發(fā)中，砂體空間分布描述與定量評(píng)價(jià)至關(guān)重要，特別是當(dāng)前油氣勘探從傳統(tǒng)的構(gòu)造性尤其藏向隱蔽的致密巖性油氣藏轉(zhuǎn)變，對(duì)含油氣儲(chǔ)層砂體的識(shí)別和描述精度要求愈來愈高，必須借助高精度三維地震數(shù)據(jù)才能實(shí)現(xiàn)。一致以來，地震屬性的應(yīng)用以定性描述和分析為主，對(duì)砂體的分布描述和厚度預(yù)測(cè)也以局域性的統(tǒng)計(jì)關(guān)系為主，盡管在實(shí)際生產(chǎn)中廣泛應(yīng)用，但實(shí)際應(yīng)用的可靠性并不好。

隨著地震勘探逐步向高精度、高密度和高分辨率的方法發(fā)展，今年采集的地震數(shù)據(jù)通過新方法的處理與分析已能滿足對(duì)地下薄互儲(chǔ)層砂體識(shí)別和描述的要求。相對(duì)于地震剖面層位解釋，地震屬性切片解釋更多的強(qiáng)調(diào)地震屬性異常的識(shí)別和統(tǒng)計(jì)。在利用地震屬性切片研究薄層結(jié)構(gòu)及其沉積煙花特征時(shí)，國(guó)內(nèi)中石油東方地球物理勘探公司的凌云研究組開展了大量理論和實(shí)際應(yīng)用的研究，提出基于參考標(biāo)準(zhǔn)層的地震屬性連續(xù)提取和動(dòng)態(tài)解釋的思路。在地質(zhì)綜合分析和沉積解釋方面，美國(guó)曾洪流開展了地震屬性切片的地震沉積學(xué)研究，拓展了薄互層地震識(shí)別和儲(chǔ)層預(yù)測(cè)的應(yīng)用領(lǐng)域和深度。

近年來，以地震屬性切片解釋為技術(shù)內(nèi)涵的地震沉積學(xué)引起了地震資料解釋人員的廣泛關(guān)注。它綜合地質(zhì)認(rèn)識(shí)和測(cè)井?dāng)?shù)據(jù)最大限度地挖掘地震數(shù)據(jù)的潛在信息，對(duì)沉積過程和沉積內(nèi)幕進(jìn)行預(yù)測(cè)和分析，綜合了地震屬性切片技術(shù)和沉積學(xué)知識(shí)，使地震屬性切片賦予了沉積演化史的含義。誠(chéng)然，與傳統(tǒng)的地震剖面解釋方法相比，地震沉積學(xué)的技術(shù)內(nèi)涵還比較模糊，地球物理理論基礎(chǔ)還較為薄弱，技術(shù)的適用條件和應(yīng)用對(duì)象尚缺乏比較嚴(yán)密的理論支撐和相對(duì)系統(tǒng)的試驗(yàn)論證，加之應(yīng)用人員知識(shí)結(jié)構(gòu)和實(shí)踐經(jīng)驗(yàn)等因素的限制，在實(shí)際應(yīng)用中，出現(xiàn)了將地震屬性切片解釋技術(shù)簡(jiǎn)單化和程式化的傾向，這種傾向加劇了對(duì)地震切片技術(shù)解決地質(zhì)問題的分歧和爭(zhēng)議。傳統(tǒng)的地震屬性解釋方法并不能準(zhǔn)確識(shí)別和評(píng)價(jià)致密砂巖儲(chǔ)層的高孔滲甜點(diǎn)區(qū)分布和物性參數(shù)，實(shí)際的預(yù)測(cè)效果和精度已不能滿足當(dāng)下油氣勘探的要求。

技術(shù)實(shí)現(xiàn)要素：

本發(fā)明的目的是提供一種碎屑巖儲(chǔ)層中地震屬性砂體空間分布描述及評(píng)價(jià)的方法，從而解決現(xiàn)有的地震屬性解釋方法預(yù)測(cè)精度低，尤其對(duì)薄互層砂體識(shí)別性差的問題。

為了實(shí)現(xiàn)以上目的，本發(fā)明所采用的技術(shù)方案是：

一種碎屑巖儲(chǔ)層中地震屬性砂體空間分布描述及評(píng)價(jià)的方法，包括以下步驟：

1)對(duì)疊后地震數(shù)據(jù)進(jìn)行樣條插值和加密采樣，得到加密采樣的地震數(shù)據(jù)；

2)依據(jù)加密采樣的地震數(shù)據(jù)計(jì)算不同種類的地震屬性；根據(jù)每一類地震屬性與巖性、孔滲性和含油氣性的敏感關(guān)系，優(yōu)化地震屬性組合，得到綜合地震屬性；

3)根據(jù)綜合地震屬性識(shí)別甜點(diǎn)砂體的空間分布；

4)結(jié)合甜點(diǎn)砂體的空間分布，通過測(cè)井孔隙度與含油氣性的標(biāo)定，實(shí)現(xiàn)碎屑巖儲(chǔ)層的砂體識(shí)別和綜合評(píng)價(jià)。

地震屬性在儲(chǔ)層預(yù)測(cè)中起著至關(guān)重要的作用。在具體地震地質(zhì)解釋過程中，通常用地震屬性切片定走向，用地震屬性剖面定傾向，兩者綜合定產(chǎn)狀。在地震屬性切片和剖面屬性分析中，在薄互層預(yù)測(cè)中地震道的時(shí)間采樣率偏大，并不能得到準(zhǔn)確的解釋層位。

步驟1)通過對(duì)地震數(shù)據(jù)進(jìn)行樣條插值和加密采樣，從而消除采樣數(shù)據(jù)誤差，使差值數(shù)據(jù)更好的反應(yīng)原始數(shù)據(jù)。優(yōu)選的，所述樣條差值和加密采樣包括以下步驟：

a)假定有n+1個(gè)數(shù)據(jù)點(diǎn)：(x₀,y₀)，(x₁,y₁)，(x₂,y₂)，……，(x_n,y_n)；計(jì)算步長(zhǎng)h_i＝x_i+1-x_i；(i＝0,1,2,…,n-1)；

b)將數(shù)據(jù)節(jié)點(diǎn)和指定的首端點(diǎn)條件代入矩陣方程：

c)解三對(duì)角矩陣方程，求得二次微分值m_i；

d)計(jì)算樣條曲線的系數(shù)：a_i＝y(tǒng)_i；b_i＝(y_i+1-y_i)/h_i-h_im_i/2-h_i(m_i+1-m_i)/6；c_i＝m_i/2；d_i＝(m_i+1-m_i)/(6h_i)；其中，i＝0,1,2,…,n-1；

e)在每個(gè)子區(qū)間x_i≤x≤x_i+1中，創(chuàng)建方程：

g_i(x)＝a_i+b_i(x-x_i)+c_i(x-x_i)²+d_i(x-x_i)³，即得。

步驟2)中，所述不同種類的地震屬性為振幅類、復(fù)數(shù)道類、頻譜統(tǒng)計(jì)類、序列統(tǒng)計(jì)類和相關(guān)統(tǒng)計(jì)類地震屬性。該步驟為地震屬性的提取和優(yōu)化步驟，通過篩選出對(duì)儲(chǔ)層流體變化具有敏感性的地震屬性進(jìn)行組合，所得綜合地震屬性與砂體的分布、巖性、物性參數(shù)、含油氣性的相關(guān)性好，從而有利于提高薄互層砂體的識(shí)別與預(yù)測(cè)精度。

步驟3)中，甜點(diǎn)砂體的空間分布包括砂體的形態(tài)、平面分布范圍、空間連通特征和縱向分布厚度。所述識(shí)別包括以下步驟：

ⅰ)將綜合地震屬性進(jìn)行褶積模型運(yùn)算，得到巖層或水平地震屬性切片并定位切片的零值點(diǎn)；

ⅱ)根據(jù)零值點(diǎn)之間的間隔，分離和識(shí)別不同深度的砂體及平面分布；

ⅲ)根據(jù)零值點(diǎn)之間的間隔和砂體出現(xiàn)的雙程旅行時(shí)，定量評(píng)價(jià)砂體厚度。

本發(fā)明的碎屑巖儲(chǔ)層中地震屬性砂體空間分布描述及評(píng)價(jià)的方法，綜合考慮測(cè)井資料標(biāo)定的地震屬性數(shù)據(jù)體特征，通過綜合應(yīng)用測(cè)井?dāng)?shù)據(jù)標(biāo)定、地震數(shù)據(jù)樣條插值加密采樣、地震屬性提取、屬性切片零點(diǎn)分析技術(shù)，實(shí)現(xiàn)地震屬性砂體厚度與物性參數(shù)評(píng)價(jià)。

本發(fā)明的碎屑巖儲(chǔ)層中地震屬性砂體空間分布描述及評(píng)價(jià)的方法，具有以下優(yōu)點(diǎn)：

(1)可以實(shí)現(xiàn)不同間隔、深度疊置砂體的識(shí)別與評(píng)價(jià)，所得到的相應(yīng)特征與砂體厚度及疊置關(guān)系的統(tǒng)計(jì)規(guī)律，適用于零相、最小相位子波，提高了復(fù)雜疊置關(guān)系的薄互層砂體的識(shí)別和評(píng)價(jià)；

(2)綜合了多種地層屬性切片，開展了砂體位置關(guān)系與地震波干涉影響的半定量試驗(yàn)分析，所建立的地震屬性砂體空間分布識(shí)別和預(yù)測(cè)方法，與傳統(tǒng)地震解釋方法相比，對(duì)薄互層砂體的識(shí)別與預(yù)測(cè)相比，獲得了更好精度的結(jié)果；

(3)用希爾伯特變換計(jì)算的不同地震屬性，在噪聲背景下可以更好地識(shí)別薄砂體，并對(duì)薄互砂體進(jìn)行定量評(píng)價(jià)，根據(jù)所提出的優(yōu)化地震屬性分類，提高了地震屬性定量分析砂體厚度與物性參數(shù)的精度，對(duì)低孔滲泥質(zhì)砂巖儲(chǔ)層和致密微裂縫儲(chǔ)層有更好的適應(yīng)性，應(yīng)用前景良好；

(4)該方法給出了不同地震屬性切片的零點(diǎn)間隔與砂體厚度的對(duì)應(yīng)關(guān)系，縮小了實(shí)際測(cè)井資料標(biāo)定的地震屬性砂體厚度預(yù)測(cè)結(jié)果的誤差范圍，拓展了地震屬性進(jìn)行薄互砂體預(yù)測(cè)和定量評(píng)價(jià)的應(yīng)用范圍，將會(huì)發(fā)展成為復(fù)雜巖性砂泥巖儲(chǔ)層分布評(píng)價(jià)的有力工具；

(5)該方法適用于薄互層砂體空間分布的預(yù)測(cè)，也可適用于除碎屑儲(chǔ)層外的地層巖性，提高了復(fù)雜巖性儲(chǔ)層預(yù)測(cè)和評(píng)價(jià)的精度。

附圖說明

圖1為厚度統(tǒng)一為5m的砂-泥-砂互層地震響應(yīng)；

圖2為楔形模型圖，其中楔形區(qū)域代表砂巖，其他區(qū)域?yàn)槟鄮r；

圖3為圖2的楔形模型褶積地震記錄；

圖4為2ms采樣數(shù)據(jù)與0.5ms差值后的數(shù)據(jù)對(duì)比；

圖5為主要標(biāo)志層的解釋結(jié)果；

圖6為圖5的地震剖面層位差值結(jié)果；

圖7為地震屬性分析流程示意圖；

圖8為地震記錄(左)與頻譜(右)屬性參數(shù)示意圖；

圖9為模型1的示意圖及x、y方向的剖面圖；其中，左圖為模型1的三維顯示圖，中圖為y＝5m平行于x軸剖面圖，右圖為x＝5m平行于y軸剖面圖；

圖10為模型1中不同的y值的褶積模型結(jié)果；其中(a)y＝5m，第二層厚度5m，(b)y＝25m，第二層厚度10m，(c)y＝40m，第二層厚度15m，(d)y＝55m，第二層厚度20m，(e)y＝70m，第二層厚度25m，(f)y＝85m，第二層厚度30m；

圖11為模型2的示意圖及x、y方向的剖面圖；其中，左圖為模型1的三維顯示圖，中圖為y＝5m平行于x軸剖面圖，右圖為x＝5m平行于y軸剖面圖；

圖12為模型2中不同的y值的褶積模型結(jié)果；其中(a)y＝5m，最左側(cè)地層間隔為5m，(b)y＝25m，最左側(cè)地層間隔為10m，(c)y＝40m，最左側(cè)地層間隔為15m，(d)y＝55m，最左側(cè)地層間隔為20m，(e)y＝70m，最左側(cè)地層間隔為25m，(f)y＝85m，最左側(cè)地層間隔為30m；

圖13為不同子波頻率對(duì)分辨率的影響，其中(a)、(b)、(c)分別為y＝5m處，35HZ、39HZ與30HZ褶積模型結(jié)果；

圖14為本發(fā)明的地震屬性分析風(fēng)閥和砂體識(shí)別和評(píng)價(jià)技術(shù)的流程圖；

圖15為模型4的兩層砂泥巖薄互層三維地質(zhì)模型示意圖；

圖16為模型4的上層砂體整體形態(tài)示意圖(左)和下層砂體整體形態(tài)示意圖(右)；

圖17為模型4的零相位子波230ms-247ms的均方根振幅切片；

圖18為模型4的最小相位子波192ms-203ms的均方根振幅切片；

圖19為模型5的零相位子波的均方根振幅切片；

圖20為模型5的最小相位子波的均方根振幅切片；

圖21為模型5的零相位子波的均方根振幅切片；

圖22為模型5的最小相位子波的均方根振幅切片；

圖23為TK1115井區(qū)T46-s3t0-3零時(shí)窗均方根振幅屬性顯示的砂體分布；

圖24為TK1115井區(qū)T46-s3t0-4零時(shí)窗均方根振幅屬性顯示的砂體分布；

圖25為TK1115井區(qū)T46-s3t0-5零時(shí)窗均方根振幅屬性顯示的砂體分布；

圖26為TK1115井區(qū)T46-s3t0-6零時(shí)窗均方根振幅屬性顯示的砂體分布；

圖27為TK1115井區(qū)T46-s3t0-7零時(shí)窗均方根振幅屬性顯示的砂體分布。

具體實(shí)施方式

下面結(jié)合具體實(shí)施例對(duì)本發(fā)明作具體說明。

1、地震屬性儲(chǔ)層預(yù)測(cè)的基本概念

1.1垂向分辨率

我國(guó)東部的大部分地區(qū)都是以中、新生代陸相含油盆地為主，其沉積特點(diǎn)都是以薄層砂、泥巖沉積為主，地層巖性在縱向上變化較大，多數(shù)儲(chǔ)層厚度遠(yuǎn)小于地震勘探的垂向分辨率。地震勘探中通常定義地層厚度小于λ/4的地層為薄層，λ為地震子波的波長(zhǎng)，即地震波在該層內(nèi)傳播的雙程旅行時(shí)小于波的半個(gè)周期。實(shí)際應(yīng)用中，垂向分辨率定義為：

其中，v是地震波的傳播速度。假設(shè)地層速度為2500m/s，子波的頻率為f＝25Hz，則得△h<λ/4為25m。然而，在東部陸相含油盆地中，常見的含油氣砂層為10m以下的層。由此可見，基于地震記錄處理的剖面，識(shí)別薄互層是十分困難的，以下結(jié)合圖1的模擬地震記錄說明這一困難。

圖1給出了一個(gè)由褶積模型模擬的砂泥巖薄互層地震響應(yīng)。由圖1可見在復(fù)合波一個(gè)視周期范圍內(nèi)，復(fù)合波峰值能量不能對(duì)應(yīng)砂體頂、底反射深度。圖中用豎直實(shí)線指示的是上下層記錄的峰值，從與單層的獨(dú)立響應(yīng)的比較看到，由于薄互層間距太小，子波之間形成強(qiáng)烈干涉，導(dǎo)致復(fù)合波的波形上難以分辨出上下層的差異。

為了說明地震剖面分辨薄層厚度的能力，設(shè)計(jì)了圖2所述的楔形砂體模型，用褶積模型的模擬響應(yīng)來考察地震響應(yīng)與厚度之間的關(guān)系。

圖2中的楔形砂體從左到右由1米增加到40米，用紅色區(qū)域表示，其中，砂巖速度為2865.3m/s，密度為2.29g/cm³，泥巖速度為2816.9m/s，密度為2.35g/cm³，通過計(jì)算波阻抗和反射系數(shù)，轉(zhuǎn)換成雙程旅行時(shí)，選用35Hz的雷克子波進(jìn)行褶積，得到圖3的地震記錄。由圖3的模擬地震記錄可以看到，y為25m處對(duì)應(yīng)的厚度為5m，由于薄互層上下界面的相互干涉導(dǎo)致在圖3的記錄上難以分辨楔形砂體的頂、底界面。

根據(jù)Knapp的定義，垂向分辨率用地震子波脈沖的時(shí)間延續(xù)長(zhǎng)度來定義，并且把這種分辨率稱為厚度分辨率。它與目前廣泛應(yīng)用的一個(gè)不太嚴(yán)格的垂向分辨率(時(shí)間分辨率)公式存在較為復(fù)雜的關(guān)系，實(shí)際上：

T_R＝△t＝1/(2.3f^*) (2)

式中，f^*是地震子波的視頻率(或主頻)。由(2)式出發(fā)，只要剖面上主頻高，就認(rèn)為其時(shí)間分辨率高。若把地震波傳播速度乘上該時(shí)間分辨率便可以得到厚度分辨率：

Z_R＝△z＝(vT_R)/2＝λ^*/4.6 (3)

式中λ^*是視波長(zhǎng)，Z_R稱為可分辨厚度。

(2)和(3)中均假設(shè)地震子波為理想雷克子波。按照瑞利(Rayleigh)準(zhǔn)則，即：“一個(gè)反射波的分辨率的極限是1/4波長(zhǎng)”，若帶通子波通頻帶的下限稱為f₁，上限稱為f₂，則通頻帶的中心頻率f_c決定了視頻率f_p(或稱主頻)，即

f_p＝f_c＝(f₁+f₂)/2

實(shí)際上，主頻就是我們眼睛看到的剖面上同相軸的胖瘦程度。對(duì)于低截頻是f1＝5Hz，高截頻f2＝70Hz的實(shí)際資料，這里主頻f＝35是地震子波為零相位雷克子波時(shí)的主頻，假設(shè)地層的速度為3000米/秒，由Rayleigh準(zhǔn)則：

△H＝vT/4＝v/(4f)＝λ/4

其中，_△H表示可分辨的地層垂直厚度，_T表示地震波的周期，v表示地層的層速度，λ表示地震波的波長(zhǎng)。

可以知道垂向上最高可分辨出25m左右的地層。Knapp劃分了韻律層的分辨率，即對(duì)薄層狀的周期性重復(fù)的沉積韻律層來說，可以用高頻來檢測(cè)薄互層的厚度，這種韻律層產(chǎn)生一種對(duì)頻率高次諧波的調(diào)諧作用。當(dāng)然這只是在理論模型情況下計(jì)算的結(jié)果，實(shí)際沉積剖面是一個(gè)復(fù)雜的地層組合體，沒有足夠的頻寬無法準(zhǔn)確描述它，同時(shí)，實(shí)際地震資料的垂向分辨率還受到地下地層結(jié)構(gòu)、巖性和地震資料頻帶寬度和信噪比的影響；垂向剖面解釋并不能準(zhǔn)確預(yù)測(cè)和評(píng)價(jià)薄互層砂體結(jié)構(gòu)。

1.2地震屬性分析

地震屬性分析是上世紀(jì)八十年代后迅速發(fā)展起來的地震解釋技術(shù)，主要用于儲(chǔ)層的正向預(yù)測(cè)。地震屬性是指由疊前或疊后地震資料，利用現(xiàn)代信號(hào)處理理論和其他數(shù)學(xué)工具而導(dǎo)出的有關(guān)地震波的幾何形態(tài)、動(dòng)力學(xué)特征、運(yùn)動(dòng)學(xué)特征和統(tǒng)計(jì)學(xué)特征。近年來，因其在油氣勘探和開發(fā)中發(fā)揮著越來越重要的作用，因此有關(guān)地震屬性的提取、分析等方法得到了飛速發(fā)展，地震屬性研究成為了油藏地球物理和勘探地球物理研究的重要內(nèi)容。

目前常用的地震屬性主要有疊前處理、疊后處理和瞬時(shí)地震屬性。實(shí)際生產(chǎn)中所用的地震屬性分類主要有五大類，即振幅類、復(fù)數(shù)道類、頻譜統(tǒng)計(jì)類、序列統(tǒng)計(jì)類、相關(guān)統(tǒng)計(jì)類。其中，振幅類屬性有均方根振幅、平均絕對(duì)振幅、平均波峰、波谷振幅、振幅總量、總能量、振幅變化等，它們能識(shí)別振幅異常或?qū)有蛱卣?，有效識(shí)別巖性或含油氣砂巖的變化，預(yù)測(cè)含油氣性。尤其是均方根振幅，即對(duì)于分析時(shí)窗內(nèi)所有點(diǎn)振幅平方的平均值得平方根。因?yàn)樵谄骄郧白髁似椒剑跃礁穹鶎?duì)于檢測(cè)由地層巖性變化或儲(chǔ)層流體引起的振幅橫向變化更為敏感。

復(fù)數(shù)道類屬性是基于希爾伯特變換導(dǎo)出的屬性，如平均瞬時(shí)頻率、平均反射強(qiáng)度等屬性，用這類屬性能追蹤由于含油氣飽和度、斷裂、巖性或地層變化引起的相關(guān)的頻率吸收特征的變化，低值(25-30HZ)常常對(duì)應(yīng)于亮點(diǎn)(高RMS振幅)，指示含油氣性。

頻譜統(tǒng)計(jì)類屬性：對(duì)地震信號(hào)的頻率譜、能量譜進(jìn)行描述，如有效帶寬、波形弧長(zhǎng)、主頻序列(F1、F2、F3)、峰值頻率到最大頻率的斜率等屬性。例如波形弧長(zhǎng)反映了單位時(shí)間內(nèi)地震波的弧線長(zhǎng)度，對(duì)于同是強(qiáng)振幅特征、但有高頻、低頻之分的地層(含油氣砂巖)更有效果。另外一些頻率屬性如主頻序列(F1、F2、F3)、峰值頻率到最大頻率的斜率等屬性，反映油氣層段有明顯的高頻地震波被吸收的現(xiàn)象，揭示吸收衰減引起的頻率趨勢(shì)變化。該類屬性一般需要開辟一定的時(shí)窗才能統(tǒng)計(jì)出來，對(duì)于研究的目標(biāo)層位較薄，該類屬性可能包含上、下層位的信息，從而對(duì)有效信息產(chǎn)生干擾。

序列統(tǒng)計(jì)類屬性：如目前較流行的能量過半、半能量時(shí)間處的斜率能反映垂直地層變化趨勢(shì)，進(jìn)行沉積旋回分析。

2、儲(chǔ)層砂體預(yù)測(cè)中的地震屬性分析方法

地震屬性在儲(chǔ)層預(yù)測(cè)中起著至關(guān)重要的作用。在具體地震地質(zhì)解釋過程中，通常用地震屬性切片定走向，用屬性剖面定傾向，兩者綜合定產(chǎn)狀。在地震屬性切片和剖面屬性分析中，發(fā)現(xiàn)在薄互層預(yù)測(cè)中地震道的時(shí)間采樣率偏大，需要對(duì)疊后地震記錄進(jìn)行插值和加密采樣，具體的采樣率視研究目標(biāo)區(qū)的最小砂體厚度而定。

2.1地震數(shù)據(jù)的樣條差值和加密采樣

地下構(gòu)造在連續(xù)的時(shí)間切片上會(huì)表現(xiàn)為可追蹤的特點(diǎn)，并且可以容易的確定構(gòu)造的走向和構(gòu)造的高點(diǎn)，從而輔助井位的部署。本發(fā)明設(shè)計(jì)了對(duì)2ms的地震道數(shù)據(jù)進(jìn)行時(shí)間域的四點(diǎn)插值，獲得0.5ms采樣的地震數(shù)據(jù)體。然后提取其水平切片，從這些切片中可以清楚的識(shí)別可追蹤的小型構(gòu)造，然后對(duì)插值獲得的小層進(jìn)行修正便可以得到更準(zhǔn)確的解釋層位。

本發(fā)明的樣條差值和加密采樣方法包括以下步驟：

a)假定有n+1個(gè)數(shù)據(jù)點(diǎn)：(x₀,y₀)，(x₁,y₁)，(x₂,y₂)，……，(x_n,y_n)；計(jì)算步長(zhǎng)h_i＝x_i+1-x_i；(i＝0,1,2,…,n-1)；

b)將數(shù)據(jù)節(jié)點(diǎn)和指定的首端點(diǎn)條件代入矩陣方程：

c)解三對(duì)角矩陣方程，求得二次微分值m_i；

d)計(jì)算樣條曲線的系數(shù)：a_i＝y(tǒng)_i；b_i＝(y_i+1-y_i)/h_i-h_im_i/2-h_i(m_i+1-m_i)/6；c_i＝m_i/2；d_i＝(m_i+1-m_i)/(6h_i)；其中，i＝0,1,2,…,n-1；

e)在每個(gè)子區(qū)間x_i≤x≤x_i+1中，創(chuàng)建方程：

g_i(x)＝a_i+b_i(x-x_i)+c_i(x-x_i)²+d_i(x-x_i)³，即得。

圖4顯示了原始的模擬信號(hào)與2ms采樣后及經(jīng)四點(diǎn)0.5ms插值后數(shù)據(jù)的對(duì)比，從圖中可以看出，當(dāng)采樣時(shí)間不在2ms的整數(shù)倍時(shí)，2ms采樣數(shù)據(jù)會(huì)產(chǎn)生較大的誤差，但四點(diǎn)0.5ms插值數(shù)據(jù)較好地恢復(fù)了原始數(shù)據(jù)的特征。

針對(duì)復(fù)雜構(gòu)造和復(fù)雜巖性分布區(qū)，需要選擇標(biāo)準(zhǔn)井，對(duì)全區(qū)的井資料進(jìn)行標(biāo)志層的對(duì)比劃分，經(jīng)過精細(xì)對(duì)比，對(duì)目的層段內(nèi)的同相軸波峰、波谷進(jìn)行追蹤，得到目的層精確的層位，如圖5所示的附帶層位信息的地震剖面。針對(duì)研究區(qū)儲(chǔ)層較薄和砂體規(guī)模小的情況，為了對(duì)細(xì)小砂體進(jìn)行更好的識(shí)別，對(duì)每個(gè)波峰與波谷間采用上述非線性插值的方式獲得加密的層位(如圖6所示)，這樣就得到加密的的地震切片。

2.2地震屬性提取和優(yōu)化

地震數(shù)據(jù)體中包括地下大量的地質(zhì)特征信息，從地震數(shù)據(jù)體中提取屬性，反映的是幾何學(xué)、動(dòng)力學(xué)、運(yùn)動(dòng)學(xué)、統(tǒng)計(jì)學(xué)等方面的特征，因此不同的屬性值可能與某些地質(zhì)參數(shù)具有很大的相關(guān)性。如果能將這些豐富的地震屬性與特定的物理和地質(zhì)現(xiàn)象建立起較好的定性或定量關(guān)系，則可以為地層的橫向變化和物性預(yù)測(cè)提供許多必需的參數(shù)。

圖7是地震屬性分析流程示意圖。地震屬性參數(shù)的提取正越來越多地應(yīng)用在地震勘探、開發(fā)階段上，目前可以從地震數(shù)據(jù)體中提取近百種屬性，大致可分為瞬時(shí)類參數(shù)(如瞬時(shí)相位、瞬時(shí)頻率、瞬時(shí)振幅等)、相關(guān)統(tǒng)計(jì)類參數(shù)、頻(能)譜類參數(shù)、層序統(tǒng)計(jì)類參數(shù)、混沌參數(shù)、突變參數(shù)等，用于幫助識(shí)別巖性、地層層序變化、不整合、斷層、流體的變化、儲(chǔ)層的孔隙率變化、河流、三角洲砂體、某種類型的礁體、地層調(diào)諧效應(yīng)?？傊?，不同的地震屬性，從不同的角度反應(yīng)了地層的物理特性。

對(duì)常用的幾種地震屬性，其簡(jiǎn)單的計(jì)算方法如下。如圖8所示，假設(shè)地震道數(shù)據(jù)為x(t)，有效波的時(shí)窗長(zhǎng)度為T，時(shí)間采樣間隔為△，時(shí)窗內(nèi)樣點(diǎn)數(shù)為N，x(i△)則表示第i個(gè)樣點(diǎn)的振幅值；將時(shí)窗內(nèi)每道地震數(shù)據(jù)分別進(jìn)行傅立葉變換，可得傅立葉譜X(ω)，X(k△ω)表示第k個(gè)樣點(diǎn)頻譜值，△ω表示頻率采樣間隔，ω₁為有效頻帶的低頻，ω₂為有效頻帶的高頻。

a)均方根振幅：對(duì)于分析時(shí)窗內(nèi)所有點(diǎn)振幅平方的平均值得平方根。

其中，a_i為第i個(gè)樣點(diǎn)的振幅值。因?yàn)樵谄骄郧白髁似椒?，所以均方根振幅?duì)于檢測(cè)由地層巖性變化或儲(chǔ)層流體引起的振幅橫向變化更為敏感。

b)振幅變化(振幅標(biāo)準(zhǔn)方差值)：指在分析時(shí)窗內(nèi)，每道內(nèi)每個(gè)樣點(diǎn)的振幅值與所有樣點(diǎn)振幅的平均值之差的平方并求和，再除以時(shí)窗內(nèi)所有的樣點(diǎn)數(shù)。

其中，a_i為第i個(gè)樣點(diǎn)的振幅值，時(shí)窗內(nèi)振幅的平均值。根據(jù)振幅變化的大小，可以用來檢測(cè)地層中流體、裂縫及巖性等變化，一般分析時(shí)窗為20-100ms。

c)弧線長(zhǎng)度：用于計(jì)算時(shí)窗內(nèi)波形的弧線長(zhǎng)度，是一個(gè)綜合了振幅和頻率特性的聯(lián)合屬性。其計(jì)算公式如下：

其中，a_i+1和a_i分別為第i+1和第i個(gè)樣點(diǎn)的振幅值，△T為采樣間隔，N為時(shí)窗內(nèi)采樣點(diǎn)數(shù)。與其它的振幅和頻率屬性聯(lián)合使用，對(duì)砂地比比較適中的地區(qū)反映最敏感，對(duì)物源區(qū)敏感度較低。該屬性可用于區(qū)別同是高振幅特征，但有高頻、低頻之分的地層情況，在砂泥巖互層中可識(shí)別富泥或富砂的地層。

d)平均反射強(qiáng)度：也稱為瞬時(shí)振幅、振幅包絡(luò)，設(shè)原始記錄為_x(i△)，經(jīng)希爾伯特(Hilbert)變換得到_y(i△)，也即復(fù)地震道中的虛部。用來識(shí)別儲(chǔ)層中流體成份、巖性、地層學(xué)變異、油氣聚集而引起的振幅異常。有：

地震屬性提取和優(yōu)化過程中，可有針對(duì)性的開時(shí)窗，提取振幅統(tǒng)計(jì)類、復(fù)地震道統(tǒng)計(jì)類、譜統(tǒng)計(jì)類、層序統(tǒng)計(jì)類、相關(guān)統(tǒng)計(jì)類等地震屬性；在與井對(duì)比分析的基礎(chǔ)上，可優(yōu)選能夠反映地層、巖性、儲(chǔ)層厚度與含油氣性特性的敏感性屬性，作為儲(chǔ)層預(yù)測(cè)及沉積相帶劃分的參考。

針對(duì)主要目的層段的巖性為泥包砂地層，可采用多種沿層地震屬性來分析目的層儲(chǔ)層沉積特征，從中優(yōu)選出適合各區(qū)塊的地震屬性來詳細(xì)分析砂體展布形態(tài)：

對(duì)于上述大量屬性，分別作出每種屬性與巖性、物性及含油性的交匯圖，得到不同屬性對(duì)巖性、物性及含油性的敏感程度，再按屬性對(duì)目標(biāo)體的敏感程度進(jìn)行分類，在所有屬性中選擇對(duì)儲(chǔ)層流體變化敏感性高的屬性制作流體分布圖或進(jìn)行預(yù)測(cè)。首先對(duì)研究區(qū)所提取的屬性進(jìn)行定性對(duì)比，分析屬性敏感程度，選取對(duì)研究對(duì)象(如巖性、物性參數(shù)、含油氣性等)比較敏感的屬性。通過將不同地震屬性與測(cè)井?dāng)?shù)據(jù)進(jìn)行對(duì)比發(fā)現(xiàn)，均方根振幅屬性與砂體的分布相關(guān)性較好。

3、不同厚度薄互層模型的地震屬性特征分析

該部分通過三個(gè)模型分析單層厚度、薄互層之間的間隔對(duì)地震記錄的影響。模型的基本屬性是模型的基本屬性是砂泥巖交互層，砂巖速度為2865.3米每秒，泥巖速度為2816.9米每秒。(速度的選取是通過實(shí)際中的測(cè)井?dāng)?shù)據(jù)選取的)。通過以下3個(gè)模型，分別研究1)在確定砂體間隔時(shí)改變砂體厚度得到的地震響應(yīng)；2)在確定砂體厚度時(shí)改變其間隔得到的地震響應(yīng)；3)在確定砂體厚度和間隔改變其子波頻率得到的地震響應(yīng)。

3.1模型1

圖9所示的模型1由兩層砂體疊置而成。上層砂體由6組平行于y軸厚度隨x坐標(biāo)增加而增加的砂體組成，下層為6組平行于x軸的厚度隨y軸坐標(biāo)增加而增加的砂體組成。這樣就形成了上下層不同厚度的36種組合。上下層不同的厚度分段x和y坐標(biāo)分別為5、25、40、55、70、85。

由圖9可見，上層砂體由6個(gè)平行于y軸的長(zhǎng)方砂體組成，厚度隨x軸增大。下層砂體由平行于x軸的6個(gè)長(zhǎng)方砂體組成。厚度沿y軸增大。上下層砂體之間間隔為5米。用30HZ雷克子波合成褶積模型利用褶積模型生成網(wǎng)格為1m X 1m的三維地震數(shù)據(jù)體。

圖10給出了y＝5、25、40、55、70、85米處的褶積模型結(jié)果。單圖能表示下層厚度不變，上層厚度分別為5、25、40、55、70、85米地震記錄。由圖10的結(jié)果可知，由于上下楔形體的四個(gè)反射界面形成的地震子波會(huì)相互干涉，使得上下砂體難以分辨，如圖10a顯示的最薄的圖層響應(yīng)，只有在上層厚度為20m時(shí)才能將上下層在同相軸上區(qū)別開。而在上下砂體間隔保持5米的時(shí)候，下覆河道厚度固定時(shí)，上覆河道越厚越容易分辨，極限大約是10m；在上覆河道最薄的情況下(5米厚度)，下覆地層達(dá)到20米厚度以上才能明顯分辨出來；由于薄砂層速度較高，當(dāng)下層厚度很小的時(shí)候，會(huì)對(duì)下層同相軸產(chǎn)生下拉作用，然而當(dāng)下層厚度逐漸增加時(shí)，這種現(xiàn)象又逐漸減弱。

3.2模型2

模型2研究砂體之間的不同間隔對(duì)薄層砂體分辨力的影響。模型2如圖11所示，上覆地層由相同厚度六個(gè)長(zhǎng)方砂體組成，與下層砂體最小間隔沿x軸增大，分別為5米、10米、15米、20米、25米、30米，下覆地層由相同厚度六個(gè)長(zhǎng)方砂體組成，與上層砂體最小間隔沿y軸增大，分別為5米、10米、15米、20米、25米、30米。對(duì)比厚度(都為5M)一致但間隔不同的薄層砂體的分辨力提取y＝5、25、40、55、70、85米處的褶積模型結(jié)果，結(jié)果如圖12所示。

由圖12可以看出，上下砂體厚度保持5米的時(shí)候：間隔為5米的時(shí)候，相鄰很近的薄層會(huì)有一些干涉，對(duì)分辨薄層會(huì)有影響。當(dāng)間隔大于20m時(shí)候，薄層之間的干涉會(huì)減弱到基本沒有影響。由此知道：1)當(dāng)兩個(gè)砂體間隔小于10m時(shí)，會(huì)對(duì)復(fù)合波造成影響。大于10m以后基本不會(huì)產(chǎn)生影響；2)由于只用5m厚度的砂體，砂體形成的響應(yīng)復(fù)合波會(huì)很短。

3.3模型3

模型3研究不同子波頻率對(duì)分辨力的影響。將子波頻率更改為35HZ、39HZ，與30HZ的褶積模型結(jié)果進(jìn)行對(duì)比，結(jié)果如圖13所示。

由圖13看出，在間隔相同厚度不同的疊復(fù)砂體的中，頻率越高得到的地震記錄更容易識(shí)別薄層砂體，而且在頻率由30上升至39的過程中，不僅是薄層分辨的更清晰了。厚層的分辨率也一樣得到了等效的提高。以上結(jié)果可以看出薄互層響應(yīng)難以分辨原因還是單層頂?shù)字g的干涉效應(yīng)。

在實(shí)際工作中，剖面地震資料的解釋遇到了瑞雷準(zhǔn)則的分辨率極限，這一極限可以由本發(fā)明的方法加以突破，從而廣泛應(yīng)用于薄層識(shí)別和解釋。

4、地震屬性切片薄紗體識(shí)別技術(shù)與流程

地震屬性種類繁多，它是在地震勘探工作中經(jīng)前人總結(jié)的能對(duì)地震波的幾何、運(yùn)動(dòng)學(xué)、動(dòng)力學(xué)特征進(jìn)行度量的參量，包括瞬時(shí)振幅，瞬時(shí)相位，瞬時(shí)頻率，相干屬性等。地震屬性能等同的反應(yīng)地下介質(zhì)的結(jié)構(gòu)信息。近年來關(guān)于地震振幅屬性與薄層結(jié)構(gòu)識(shí)別的研究越來越多，包括凌云研究組對(duì)地震屬性進(jìn)行了連續(xù)提取，極大地提高了地震解釋技術(shù)識(shí)別薄互層的準(zhǔn)確性。而地震切片技術(shù)在識(shí)別薄層的優(yōu)勢(shì)在于可以通過薄層的構(gòu)造形態(tài)來巧妙地避開地震分辨率的影響。前人在地震屬性分析方面已開展了許多有意義的研究，本專利發(fā)明了一種有別于前人研究成果的地震屬性分析方法和砂體識(shí)別和評(píng)價(jià)技術(shù)。本發(fā)明提出了一種有別于現(xiàn)有技術(shù)的地震屬性分析方法和砂體識(shí)別和評(píng)價(jià)技術(shù)，流程圖如圖14所示。

對(duì)于兩層砂體而言，四個(gè)反射界面所形成的復(fù)合波更為復(fù)雜，理論上更難以從地震剖面上進(jìn)行識(shí)別。下面通過模型4來驗(yàn)證本發(fā)明的通過復(fù)合波零值響應(yīng)的特點(diǎn)在地震屬性切片上進(jìn)行識(shí)別的效果。

模型4的示意圖如圖15和圖16所示。模型4為兩層砂泥巖薄互層三維地質(zhì)模型，其中，上下兩層的厚度都是5m，相距5m，上層深度為270m，下層深度為280m。砂巖的速度為2865.3米每秒，泥巖的速度為2816.9米每秒。

將模型分別用振幅譜相同的最小相位子波和零相位子波進(jìn)行簡(jiǎn)單的褶積模型運(yùn)算，可以得到如圖17所示的地震屬性均方根振幅切片結(jié)果。圖17為230ms-247ms，以1ms為間隔，得出的均方根振幅切片。

圖17中，238ms處，上層砂體地震反射的樣點(diǎn)值為零，地震切片上只出現(xiàn)下層砂體的反射。245ms處則只出現(xiàn)上層砂體的反射。上下層第一次出現(xiàn)零值的時(shí)間分別為238ms和245ms，相差7ms，泥巖速度為2816.9m/s，由于代表雙程旅行時(shí)，計(jì)算后測(cè)得零值之間的距離為9.859m，與砂體中心間距10m一致。

注意到240ms處的切片，可以看到上下兩層的重疊部分的響應(yīng)幾乎為零，但是非重疊部分能表現(xiàn)單層的地方卻依然有響應(yīng)，這一點(diǎn)就充分說明了地震切片技術(shù)的可識(shí)別性。

由于真實(shí)情況下，野外炸藥震源激發(fā)出的地震子波往往是最小相位的，零相位子波一般是在后期處理中給出的，為了貼近實(shí)際情況，圖18研究了采用最小相位子波進(jìn)行褶積得到的結(jié)果。

圖18為將子波換為最小相位子波進(jìn)行褶積得到均方根振幅切片結(jié)果。最小相位第一次出現(xiàn)零值的時(shí)間分別為193ms和200ms，相差7ms，泥巖速度為2816.9m/s，由于代表雙程旅行時(shí)，計(jì)算后測(cè)得零值之間的距離為9.859m，與砂體中心間距10m一致。可以看出當(dāng)子波為零相位和最小相位時(shí)，雖然出現(xiàn)零值時(shí)間不同，但兩層砂體的零值時(shí)間差，均為砂體中心雙程反射時(shí)間之差。

最小相位子波由于能量聚集要比零相位子波快，所以地震振幅切片整體響應(yīng)到達(dá)峰值的時(shí)間要早。

為進(jìn)一步研究砂體厚度對(duì)上述結(jié)論的影響，在該模型的基礎(chǔ)上，將砂體厚度由5m增加到10m，間隔不變(記為模型5)，得到不同子波的均方根振幅切片(圖19、圖20)。

圖19為模型5的零相位子波的均方根振幅切片，由圖可知，上下層第一次出現(xiàn)零值的時(shí)間分別為240ms和250ms，相差10ms，泥巖速度為2816.9m/s，由于代表雙程旅行時(shí)，計(jì)算后測(cè)得零值之間的距離為14.08m，與砂體中心間距15m基本一致。

圖20為模型5的最小相位子波的均方根振幅切片，由圖可知，雖然上下層第一次零值的響應(yīng)為195ms和205ms，之差依然是10ms，但是上層第一次零值出現(xiàn)時(shí)(195ms)，下層還沒有響應(yīng)。如果在現(xiàn)實(shí)的解釋過程當(dāng)中，上層第一次到出現(xiàn)零值，有時(shí)會(huì)將這種現(xiàn)象與單獨(dú)一層出現(xiàn)的響應(yīng)歸到一類中去，所以第一次零值不具備雙層相關(guān)的討論要素。取上下兩層第2次3次零值之間的間距，發(fā)現(xiàn)其與第一次零值間距基本一致。

這里所說的相同次數(shù)的零值之間時(shí)間間隔一致，是因?yàn)榱阒禃r(shí)間建立在上下界面的反射系數(shù)時(shí)間間隔不變的情況下，褶積模型的原理相當(dāng)于將復(fù)合波代替反射系數(shù)，所以復(fù)合波零值時(shí)間間隔應(yīng)與反射系數(shù)間隔一致。

因此，選取第2次零值出現(xiàn)時(shí)間207ms和217ms。計(jì)算可得零值之間距離15.7399m，與砂體中心間距15m一致。

根據(jù)上述分析，分析薄層間隔的影響，無影響的區(qū)域大概在10m左右，所以在應(yīng)用地震屬性切片進(jìn)行識(shí)別時(shí)需要討論10m間隔時(shí)的有效性。在上一個(gè)模型的基礎(chǔ)上增大間隔，考察地震波的干涉響應(yīng)有何種影響。將兩層間隔加大5m，每層厚度仍然為5m(記為模型6)。

圖21為模型6的零相位子波的均方根振幅切片，由圖可知，上下層第一次出現(xiàn)零值的時(shí)間分別為238ms和249ms，相差11ms，泥巖速度為2816.9m/s，由于代表雙程旅行時(shí)，計(jì)算后測(cè)得零值之間的距離為15.7399m，與砂體中心間距15m一致。

圖22為模型6的最小相位子波的均方根振幅切片，由圖可知，最小相位子波的結(jié)果依然是振幅出現(xiàn)的整體時(shí)間提前了，193ms出上層砂體第一次零值出現(xiàn)，但沒有下層響應(yīng)，依然選擇雙層的第二次零值時(shí)間，205ms和216ms，相差11ms，泥巖速度為2816.9m/s，由于代表雙程旅行時(shí)，計(jì)算后測(cè)得零值之間的距離為15.7399m，與砂體中心間距15m一致。

由此可見，在不同的子波作用下，地震振幅屬性切片的零值時(shí)間間隔在識(shí)別薄互層間距的作用上有著極大的穩(wěn)定性。

本發(fā)明的方法從多層界面地震波干涉原理出發(fā)創(chuàng)造性的提出了一套完整的地震屬性分析/提取/優(yōu)化的技術(shù)思路，建議了系列配套的地震屬性處理和應(yīng)用技術(shù)：基于薄互層砂體預(yù)測(cè)與識(shí)別為目的的地震資料插值加密采樣方法，給出了具有普遍意義的算法與處理流程。本發(fā)明的儲(chǔ)層砂體預(yù)測(cè)和識(shí)別評(píng)價(jià)方法，以尋找地層中發(fā)育的砂體為出發(fā)點(diǎn)建立，但并未引入砂體形態(tài)或結(jié)構(gòu)分布的假設(shè)和先驗(yàn)?zāi)Ｐ妥鳛榧s束，屬于數(shù)據(jù)驅(qū)動(dòng)模型的預(yù)測(cè)與識(shí)別方法，對(duì)先驗(yàn)?zāi)Ｐ鸵蕾嚦潭鹊停m用于未知探區(qū)的應(yīng)用。

以下實(shí)施例中，對(duì)中國(guó)西部某油田的一個(gè)區(qū)域進(jìn)行處理和分析，通過目的層上油組頂面和下油組頂面之間的波峰和波谷的詳細(xì)解釋，對(duì)鄰近的波峰和波谷之間采用上述樣條差值和加密采樣方法獲得了三個(gè)加密層位，沿這些層位提取各種屬性，并優(yōu)選利用均方根振幅屬性作為識(shí)別河道砂體的參考；同時(shí)，利用地震屬性的水平時(shí)間切片來對(duì)局部的細(xì)微特征進(jìn)行修正。

實(shí)施例1

本實(shí)施例的碎屑巖儲(chǔ)層中地震屬性砂體空間分布描述及評(píng)價(jià)的方法，針對(duì)TK1115區(qū)塊淺部河道砂體進(jìn)行，包括以下步驟：

1)采用上述2.1節(jié)的方法對(duì)疊后地震數(shù)據(jù)進(jìn)行樣條插值和加密采樣，得到加密采樣的地震數(shù)據(jù)；

2)由上述2.2節(jié)的方法對(duì)加密采樣的地震數(shù)據(jù)計(jì)算不同種類的地震屬性；根據(jù)每一類地震屬性與巖性、孔滲性和含油氣性的敏感關(guān)系，篩選出阻抗、最大值、波形分類、總振幅、平均振幅、均方根振幅、最大波峰振幅屬性為綜合地震屬性，其中，均方根振幅屬性很好的反映了區(qū)塊內(nèi)油氣層的曲流河道砂體形態(tài)；

3)根據(jù)綜合地震屬性識(shí)別甜點(diǎn)砂體的空間分布；

4)結(jié)合甜點(diǎn)砂體的空間分布，通過測(cè)井孔隙度與含油氣性的標(biāo)定，實(shí)現(xiàn)碎屑巖儲(chǔ)層的砂體識(shí)別和綜合評(píng)價(jià)。

本實(shí)施例中，依據(jù)均方根振幅屬性顯示的砂體分布如圖23、圖24所示。由圖可識(shí)別出砂巖厚度大，儲(chǔ)層物性好的區(qū)域，此外，圖中與河道相連接的河口壩砂體也清晰可見；比圖23和圖24可以清楚看到與河道相連的砂體向下分布范圍變小的情況。

實(shí)施例2

本實(shí)施例的碎屑巖儲(chǔ)層中地震屬性砂體空間分布描述及評(píng)價(jià)的方法，針對(duì)TK1115區(qū)塊深部河道砂體進(jìn)行，包括以下步驟：

1)采用上述2.1節(jié)的方法對(duì)疊后地震數(shù)據(jù)進(jìn)行樣條插值和加密采樣，得到加密采樣的地震數(shù)據(jù)；

2)由上述2.2節(jié)的方法對(duì)加密采樣的地震數(shù)據(jù)計(jì)算不同種類的地震屬性；根據(jù)每一類地震屬性與巖性、孔滲性和含油氣性的敏感關(guān)系，篩選出阻抗、最大值、波形分類、總振幅、平均振幅、均方根振幅、最大波峰振幅屬性為綜合地震屬性，其中，均方根振幅屬性很好的反映了區(qū)塊內(nèi)油氣層的曲流河道砂體形態(tài)；

3)根據(jù)綜合地震屬性識(shí)別甜點(diǎn)砂體的空間分布；

4)結(jié)合甜點(diǎn)砂體的空間分布，通過測(cè)井孔隙度與含油氣性的標(biāo)定，實(shí)現(xiàn)碎屑巖儲(chǔ)層的砂體識(shí)別和綜合評(píng)價(jià)。

本實(shí)施例中，依據(jù)均方根振幅屬性顯示的砂體分布如圖25～27所示。其較好的反映了在圖23和圖24上顯示的河道之下并與該河道相交的河道砂體形態(tài)，該河道的屬性與背景值的差異十分明顯，而且與河道相連的朵葉狀砂體也在屬性切片上清晰可見。

由圖可識(shí)別出砂巖厚度大，儲(chǔ)層物性好的區(qū)域，此外，圖中與河道相連接的河口壩砂體也清晰可見；比圖23和圖24可以清楚看到與河道相連的砂體向下分布范圍變小的情況。圖27是該井區(qū)T46-s3t0-7時(shí)刻對(duì)應(yīng)的均方根振幅，其較好的反映了在圖25和圖26上顯示的河道相連的朵葉狀砂體在空間上的變化。圖中屬性顯示，該河道的屬性與背景值的差異十分明顯，而且與河道相連的朵葉狀砂體隨深度加大,分布范圍變大，由于砂體變純，相應(yīng)的地震波振幅也增強(qiáng)。

本發(fā)明的特點(diǎn)是通過綜合測(cè)井信息、多類地震屬性，利用非線性插值和加密采樣的優(yōu)化地震屬性，考慮了不同地震屬性類對(duì)巖性、物性和含油氣性不同的敏感程度，拓展了地震數(shù)據(jù)儲(chǔ)層砂體預(yù)測(cè)和評(píng)價(jià)的應(yīng)用范圍。由于在本發(fā)明實(shí)施例中不僅僅是根據(jù)多類地震屬性，還包括了測(cè)井等信息的標(biāo)定和綜合，不僅提高了儲(chǔ)層砂體預(yù)測(cè)和評(píng)價(jià)精度，也為油氣勘探開發(fā)中精確評(píng)價(jià)砂體的空間變化及預(yù)測(cè)有利油氣勘探開發(fā)遠(yuǎn)景區(qū)提供支撐。

需要說明的是：上述實(shí)施例提供的地震屬性切片分析在識(shí)別砂體的平面分布和預(yù)測(cè)砂體厚度時(shí)，僅以上述各功能模塊的劃分進(jìn)行舉例說明，實(shí)際應(yīng)用中，可以根據(jù)需要而將上述功能分配由不同的功能模塊完成，即將裝置的內(nèi)部結(jié)構(gòu)劃分成不同的功能模塊，以完成以上描述的全部或者部分功能。另外，上述實(shí)施例提供的地震屬性分析進(jìn)行砂體預(yù)測(cè)和評(píng)價(jià)的裝置與河道砂分布預(yù)測(cè)和評(píng)價(jià)的方法實(shí)施例屬于同一構(gòu)思，其具體實(shí)現(xiàn)過程詳見方法實(shí)施例，這里不再贅述。

上述本發(fā)明實(shí)施例序號(hào)僅僅為了描述，不代表實(shí)施例的優(yōu)劣。

本領(lǐng)域普通技術(shù)人員可以理解實(shí)現(xiàn)上述實(shí)施例的全部或部分步驟可以通過硬件來完成，也可以通過程序來指令相關(guān)的硬件完成，所述的程序可以存儲(chǔ)于一種計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)中，上述提到的存儲(chǔ)介質(zhì)可以是只讀存儲(chǔ)器，磁盤或光盤等。

以上所述僅為本發(fā)明的較佳實(shí)施例，并不用以限制本發(fā)明，凡在本發(fā)明的精神和原則之內(nèi)，所作的任何修改、等同替換、改進(jìn)等，均應(yīng)包含在本發(fā)明的保護(hù)范圍之內(nèi)。