宋光超

摘 ?要：儲層參數(shù)的變化引起地震屬性的變化，但是儲層參數(shù)與地震屬性的變化并沒有一一對應(yīng)的關(guān)系，為解決此問題采用地震多屬性聚類分析技術(shù)。通過屬性優(yōu)選，將弧長、上半周持續(xù)時(shí)間、帶寬、最小振幅、瞬時(shí)頻率五種對沙四段1砂組儲層敏感且相互獨(dú)立的地震屬性參與聚類分析，取得了較好的應(yīng)用效果。

關(guān)鍵詞：八面河;地震屬性;屬性優(yōu)選;聚類分析

中圖分類號：P631.4 ? ? ? ?文獻(xiàn)標(biāo)志碼：A ? ? ? ? 文章編號：2095-2945（2019）18-0186-03

Abstract： The change of reservoir parameters leads to the change of seismic attributes， but there is no one-to-one corresponding relationship between reservoir parameters and seismic attributes. Seismic multi-attribute clustering analysis is used to solve this problem. Through attribute optimization， five kinds of seismic attributes， such as arc length， duration of the first half cycle， bandwidth， minimum amplitude and instantaneous frequency， which are sensitive and independent to the reservoir of ES4 member 1 sand group， are involved in cluster analysis， and good application results are obtained.

Keywords： Bamianhe; seismic attribute; attribute optimization; cluster analysis

八面河西區(qū)位于八面河鼻狀斷裂帶西段，總體構(gòu)造格局是一個(gè)南高北低的受中生界基底控制的繼承性發(fā)育的西緩東陡、北東-南西走向的寬緩鼻狀構(gòu)造。受繼承性古隆起控制，該區(qū)沙四段主要發(fā)育灘壩砂體沉積，砂體橫向分布不穩(wěn)定，橫向預(yù)測難。筆者利用地震多屬性聚類分析技術(shù)，將地震數(shù)據(jù)與地質(zhì)信息之間的非線性關(guān)系正確的反映出來，建立八面河西區(qū)沙四段地震數(shù)據(jù)與儲層分布之間的關(guān)系，有效的預(yù)測了灘壩砂體的發(fā)育和分布規(guī)律。

1 地震屬性提取

地震屬性是指地震數(shù)據(jù)中所包含的幾何學(xué)、動力學(xué)、運(yùn)動學(xué)或統(tǒng)計(jì)學(xué)特征。在目標(biāo)地區(qū)地震地質(zhì)情況確定的情況下，儲層厚薄發(fā)生變化，地震數(shù)據(jù)會有所反映，表現(xiàn)為波形、能量、頻率和相位等一系列變化。

本次研究利用Geoframe軟件中的SATK屬性提取模塊，針對八面河西區(qū)沙四段1砂組提取了計(jì)算地震屬性類（包括均方根振幅、瞬時(shí)相位、弧長、最大強(qiáng)度等）、軌跡屬性類（包括上半周持續(xù)時(shí)間、下半周持續(xù)時(shí)間、上半周面積、下半周面積等）、持續(xù)時(shí)間屬性類（包括正半周平均持續(xù)時(shí)間、最大周期持續(xù)時(shí)間、持續(xù)時(shí)間的標(biāo)準(zhǔn)方差等）、層段屬性類（包括振幅標(biāo)準(zhǔn)差、平均正振幅、波峰平均值等）、頻譜屬性類（包括頻帶寬度、總能量、主頻等）共計(jì)四大類40多種地震屬性。

2 地震屬性優(yōu)選

為了尋找對儲層厚度較為敏感的地震屬性，本次統(tǒng)計(jì)了工區(qū)沙四段1砂組共計(jì)28口井的砂巖垂直厚度（表1）導(dǎo)入Geoframe軟件中，并在Seisclass模塊中將井點(diǎn)處不同地震屬性分別與砂巖垂直厚度進(jìn)行二維交會分析，計(jì)算其相關(guān)系數(shù)，對相關(guān)系數(shù)在0.6以下的地震屬性剔除掉，保證所選地震屬性反映儲層厚度的可靠性。經(jīng)過分析，共優(yōu)選出最大強(qiáng)度、瞬時(shí)頻率、最大振幅、最小振幅、平均波峰值、弧長、均方根振幅、正負(fù)比率、帶寬、上半周持續(xù)時(shí)間、振幅標(biāo)準(zhǔn)差等對儲層厚度敏感的地震屬性。

對初步優(yōu)選出來的地震屬性進(jìn)行屬性相關(guān)性分析（表2），對相似程度較高的兩種或幾種屬性進(jìn)行篩選，以減少屬性空間的數(shù)據(jù)冗余而提高預(yù)測的精度。通過篩選，共篩選出弧長、上半周持續(xù)時(shí)間、帶寬、最小振幅、瞬時(shí)頻率五種相互獨(dú)立但對本區(qū)儲層敏感的地震屬性。

3 聚類分析

3.1 聚類分析方法介紹

3.1.1 K-Means聚類方法

K-Means聚類是一種統(tǒng)計(jì)學(xué)算法，首先在屬性空間內(nèi)根據(jù)用戶定義的分類數(shù)目任意的為每一類選擇一個(gè)聚類中心。由于開始聚類中心選擇的任意性，因此該計(jì)算不能總是收斂于一個(gè)解決方法，但對于輸入數(shù)據(jù)相同的多個(gè)操作，其最終結(jié)果是非常相似的。隨著計(jì)算的運(yùn)行，按數(shù)據(jù)點(diǎn)與聚類中心的最近平均值，所有數(shù)據(jù)點(diǎn)都被指定到某一類別。每一次迭代會減小數(shù)據(jù)點(diǎn)與聚類中心的平均距離。迭代過程時(shí)，當(dāng)每個(gè)點(diǎn)被指定到某一類別后，其與聚類中心的距離將被更新。迭代不斷進(jìn)行，直到每個(gè)點(diǎn)都分配到某一類別，且進(jìn)一步迭代不會改變數(shù)據(jù)點(diǎn)的類別。聚類中心的最初位置是隨機(jī)選擇的。Error Limit是可供選擇的唯一的一個(gè)控制參數(shù)，當(dāng)聚類中心間的平均值小于該參數(shù)時(shí)，聚類完成。當(dāng)將Error Limit參數(shù)設(shè)置較低時(shí)，會提高聚類的精度，但會增加計(jì)算的時(shí)間。

3.1.2 Competitive Learning聚類方法

Competitive Learning聚類是一種神經(jīng)網(wǎng)絡(luò)算法。與K-Means算法類似，該算法不需要對數(shù)據(jù)有先驗(yàn)認(rèn)識，用戶只需定義需要?jiǎng)?chuàng)建的類別數(shù)目。使用Competitive這個(gè)詞是因?yàn)樵谠撍惴ㄖ芯垲愔行男枰獙λ鶎贁?shù)據(jù)點(diǎn)進(jìn)行競爭。聚類中心的最初位置并不像K-Means算法一樣是任意給定的，而是通過算法來選擇的。隨著計(jì)算的進(jìn)行，數(shù)據(jù)點(diǎn)不斷的提供給聚類中心。數(shù)據(jù)點(diǎn)與聚類中心的最近距離不斷更新，本次距離最近的數(shù)據(jù)點(diǎn)在下一次的迭代中比其它數(shù)據(jù)點(diǎn)（在屬性空間中）具有更大的競爭優(yōu)勢。實(shí)際上，在神經(jīng)網(wǎng)絡(luò)術(shù)語中，對于提供給網(wǎng)絡(luò)的每一個(gè)數(shù)據(jù)點(diǎn)，神經(jīng)網(wǎng)絡(luò)彼此間通過競爭來決定誰是勝者。與提供數(shù)據(jù)點(diǎn)的權(quán)重系數(shù)最接近的節(jié)點(diǎn)是競爭的勝者。當(dāng)?shù)螖?shù)達(dá)到預(yù)定值，或者聚類中心間距離變化的平均值小于定義的值時(shí)，聚類完成。

3.2 聚類分析應(yīng)用

利用篩選后的弧長、上半周持續(xù)時(shí)間、帶寬、最小振幅、瞬時(shí)頻率五種地震屬性在Geoframe軟件的Seisclass模塊中進(jìn)行聚類分析運(yùn)算后得到沙四段1砂組砂巖厚度預(yù)測圖（圖1），實(shí)鉆井?dāng)?shù)據(jù)與預(yù)測數(shù)據(jù)吻合度較高，說明聚類分析可以對砂體進(jìn)行很好的預(yù)測，據(jù)此實(shí)施的4口探井均獲得了成功，說明我們所用的方法是正確的。

4 結(jié)束語

（1）弧長、上半周持續(xù)時(shí)間、帶寬、最小振幅、瞬時(shí)頻率五種地震屬性對沙四段1砂組儲層敏感且相互獨(dú)立。

（2）多屬性聚類分析技術(shù)是針對灘壩砂儲藏預(yù)測的一種有效手段，不僅可以避免單屬性的多解性問題，而且可以準(zhǔn)確預(yù)測出研究區(qū)儲層的展布情況（包括平面厚薄關(guān)系和真實(shí)儲層厚度），該方法進(jìn)一步完善了儲層預(yù)測技術(shù)系列，提高了隱蔽油藏儲層預(yù)測技術(shù)水平，降低了勘探風(fēng)險(xiǎn)。

參考文獻(xiàn)：

[1]趙政璋，趙賢正，王英民，等.儲層地震預(yù)測理論與實(shí)踐[M].北京：科學(xué)出版社，2005.

[2]王鑫，王洪國，等.聚類分析方法及工具應(yīng)用研究[J].計(jì)算機(jī)科學(xué)，2006，33（2）：197-200.

[3]徐懷大，王石鳳，陳開遠(yuǎn).地震地層學(xué)解釋基礎(chǔ)[M].北京：中國地質(zhì)大學(xué)出版社，1990.

[4]戈國華，肖海波，張敏.基于FCM的數(shù)據(jù)聚類分析及Matlab實(shí)現(xiàn)[J].福建電腦，2007，04.