多屬性聚類分析技術(shù)在八面河西區(qū)中的應(yīng)用

宋光超

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

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

中圖分類號：P631.4 ? ? ? ?文獻標志碼：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)造格局是一個南高北低的受中生界基底控制的繼承性發(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é)、運動學(xué)或統(tǒng)計學(xué)特征。在目標地區(qū)地震地質(zhì)情況確定的情況下，儲層厚薄發(fā)生變化，地震數(shù)據(jù)會有所反映，表現(xiàn)為波形、能量、頻率和相位等一系列變化。

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

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

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

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

3 聚類分析

3.1 聚類分析方法介紹

3.1.1 K-Means聚類方法

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

3.1.2 Competitive Learning聚類方法

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

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

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

4 結(jié)束語

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

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

參考文獻：

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

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

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

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