>其中，ΔΤ表示滑動(dòng)時(shí)間窗， tstep表示采樣步長，t。表示層位上樣本點(diǎn)所在的波段中距樣本點(diǎn)最近的波谷或波峰與樣本 點(diǎn)之間的采樣時(shí)間間隔。
[0067] 波段聚類方式可以包括：基于均值的波段聚類方式、基于屬性拼合的波段聚類方 式或者基于差值和的波段聚類方式。
[0068] 下面分別對上述幾種波段聚類方式進(jìn)行說明。
[0069] (1)基于均值的波段聚類方式
[0070] 聚類單元22可以包括：
[0071] 第一確定模塊，用于根據(jù)當(dāng)前層位上的樣本點(diǎn)e所在的波段，按照以下公式確定 樣本點(diǎn)e的屬性向量：
[0072]
[0073] 其中，樣本點(diǎn)e的原始屬性向量為（alyaSg，…a。，…ar〇，aiA_j表示樣本點(diǎn)e的 滑動(dòng)時(shí)間窗AT內(nèi)的點(diǎn)集合A中第j個(gè)點(diǎn)的第i個(gè)屬性值，表示滑動(dòng)時(shí)間窗 為ΔΤ的情況下樣本點(diǎn)e的第i個(gè)屬性的平均屬性值，i= 1，2，…，η;
[0074] 第一計(jì)算模塊，用于按照上述公式對當(dāng)前層位上的其他樣本點(diǎn)進(jìn)行計(jì)算，得到當(dāng) 前層位上各樣本點(diǎn)的屬性向量；
[0075] 第一聚類模塊，用于依次沿著點(diǎn)集合Α中各樣本點(diǎn)所在層位，對各樣本點(diǎn)的屬性 向量進(jìn)行聚類。
[0076] (2)基于屬性拼合的波段聚類方式
[0077] 聚類單元22可以包括：
[0078] 拼合模塊，用于將當(dāng)前層位上的樣本點(diǎn)e的滑動(dòng)時(shí)間窗ΔΤ內(nèi)的點(diǎn)集合A中的各 樣本點(diǎn)的所有屬性進(jìn)行拼合，得到樣本點(diǎn)e的波段屬性向量：(A1全屬性，A2全屬性，…，Aj 全屬性，…），其中，j=l，2，…，ΔΤ，點(diǎn)集合A中樣本點(diǎn)的全屬性為（al，a2，…ai，… an)；
[0079] 第二計(jì)算模塊，用于按照上述步驟對當(dāng)前層位上的其他樣本點(diǎn)進(jìn)行計(jì)算，得到當(dāng) 前層位上各樣本點(diǎn)的波段屬性向量；
[0080] 第二聚類模塊，用于對當(dāng)前層位上的所有樣本點(diǎn)的波段屬性向量進(jìn)行聚類。
[0081] (3)基于差值和的波段聚類方式
[0082] 聚類單元22可以包括：
[0083] 第三計(jì)算模塊，用于按照以下公式分別計(jì)算當(dāng)前層位上的樣本點(diǎn)e的屬性向量與 預(yù)先存儲(chǔ)的多個(gè)聚類中心的距離D
g中，aiA]表示樣本點(diǎn)e的 滑動(dòng)時(shí)間窗AT內(nèi)的點(diǎn)集合A中第j個(gè)點(diǎn)的第i個(gè)屬性值，aiej表示聚類中心的滑動(dòng)時(shí)間 窗AT內(nèi)的點(diǎn)集合C中第j個(gè)點(diǎn)的第i個(gè)屬性值，聚類中心與類別一一對應(yīng)；
[0084] 第二確定模塊，用于確定與樣本點(diǎn)e距離最近的聚類中心對應(yīng)的類別為樣本點(diǎn)e 的類別。
[0085] 在一個(gè)實(shí)施例中，調(diào)整單元23包括：標(biāo)簽?zāi)K，用于將起始層位的各樣本點(diǎn)按照 所屬類別進(jìn)行分類，并按預(yù)設(shè)順序進(jìn)行標(biāo)簽；區(qū)域確定模塊，用于按照相鄰層位的空間相似 度，在與上層層位相鄰的下層層位中，確定與上層層位的對應(yīng)區(qū)域；調(diào)整模塊，用于根據(jù)上 層層位的標(biāo)簽，將下層層位區(qū)域內(nèi)的樣本點(diǎn)的標(biāo)簽調(diào)整為上層層位中對應(yīng)區(qū)域的標(biāo)簽。重 復(fù)利用區(qū)域確定模塊和調(diào)整模塊完成各自的功能，直到完成待研究層段的所有層位的聚類 結(jié)果調(diào)整。
[0086] 在一個(gè)實(shí)施例中，上述裝置還可以包括：預(yù)處理單元，用于在利用滑動(dòng)時(shí)間窗對待 研究層段的地震屬性數(shù)據(jù)進(jìn)行聚類之前，對待研究層段的地震屬性數(shù)據(jù)進(jìn)行以下預(yù)處理： 對地震屬性數(shù)據(jù)中的異常數(shù)據(jù)進(jìn)行過濾處理；以及對過濾后的地震屬性數(shù)據(jù)進(jìn)行歸一化處 理。
[0087] 當(dāng)然，上述模塊劃分只是一種示意劃分，本發(fā)明并不局限于此。只要能實(shí)現(xiàn)本發(fā)明 的目的的模塊劃分，均應(yīng)屬于本發(fā)明的保護(hù)范圍。
[0088] 為了對上述地震屬性聚類裝置及裝置進(jìn)行更為清楚的解釋，下面結(jié)合具體的實(shí)施 例來進(jìn)行說明，然而值得注意的是該實(shí)施例僅是為了更好地說明本發(fā)明，并不構(gòu)成對本發(fā) 明不當(dāng)?shù)南薅ā?br>[0089] 以K-means聚類算法為例，進(jìn)行如下假設(shè)：
[0090] 選擇η種地震屬性，屬性名分別為al，a2，…，ai，…，其中i= 1，2，…，η。
[0091] 如圖3所示，假設(shè)TimeSlice方向的采樣步長為tstep，選取層位線hl，hi向下 增加At得到層位線h2。在hi上取點(diǎn)e，體聚類是沿著層位向下進(jìn)行波段聚類，這里引入 滑動(dòng)窗口ΔΤ，它等于波段中所含采樣點(diǎn)的個(gè)數(shù)，如e點(diǎn)，ΔΤ由e上下最鄰近的波谷決定， 這里基于e點(diǎn)上下波形是對稱的，若e上下波谷距e的時(shí)間不同，取較小值，設(shè)此較小值為 t0。e點(diǎn)的波段內(nèi)共有ΔΤ個(gè)采樣點(diǎn)，將它們按TimeSlice升序排列，組成集合Α。ΔΤ的 計(jì)算公式：
[0092]
[0093] 其中：t0表示e點(diǎn)上下最鄰近波谷距e點(diǎn)的采樣時(shí)間間隔，tstep表示采樣步長。 AT表不滑動(dòng)窗口。
[0094] 波段聚類一般有三種處理方式，分別是基于均值的波段聚類、基于屬性拼合的波 段聚類、基于差值和的波段聚類，如圖4所示。
[0095] (1)基于均值的波段聚類：
[0096] e的原始屬性向量：(ale，a2e，·ι?β，·ιι〇?？紤]e點(diǎn)的波段，取屬性均值作為 e的屬性向量：
[0097]
[0098] 其中：aiAj表示e點(diǎn)的滑動(dòng)時(shí)窗ΔT內(nèi)的點(diǎn)集合A中第j個(gè)點(diǎn)的第i個(gè)屬性值， ^^1/〃'4//471表示當(dāng)滑動(dòng)時(shí)窗為六1'時(shí)4點(diǎn)的第1個(gè)屬性的平均屬性值，1 = 1，2，···，!!。
[0099] 對e點(diǎn)所在層位上其他點(diǎn)進(jìn)行上述處理，依次沿著集合A中采樣點(diǎn)所在層位進(jìn)行 聚類。
[0100] (2)基于屬性拼合的波段聚類：
[0101] 將e點(diǎn)對應(yīng)的集合A中的采樣點(diǎn)的所有屬性拼合，形成波段屬性向量；
[0102] e的原始屬性向量：(ale，a2e，…aie，…ane)。
[0103] e的波段屬性向量：(A1全屬性，A2全屬性，…，Aj全屬性，…），其中j= 1，2,…， ΔΤ〇
[0104] (3)基于差值和的波段聚類
[0105] 類似于K-means聚類，先在層位面產(chǎn)生預(yù)先設(shè)置的k個(gè)聚類中心，而距離計(jì)算不同 于普通的K-means聚類。設(shè)C為層位面hi初始化聚類中心時(shí)產(chǎn)生的一個(gè)聚類中心（參見 圖5):
[0106] e到C的距離計(jì)算如下：
[0107]
[0108] 其中：aiAj表示e點(diǎn)的滑動(dòng)時(shí)窗△T內(nèi)的點(diǎn)集合A中第j個(gè)點(diǎn)的第i個(gè)屬性值， aiCj表示e點(diǎn)所在的聚類中心的滑動(dòng)時(shí)窗△T內(nèi)的點(diǎn)集合C中第j個(gè)點(diǎn)的第i個(gè)屬性值，D 為距離度量，根據(jù)D來判斷e的類別歸屬。
[0109] 聚類標(biāo)簽統(tǒng)一是基于相鄰層位面的聚類結(jié)果類似這個(gè)假設(shè)，利用空間對應(yīng)相似度 來統(tǒng)一相鄰層位面的標(biāo)簽。例如，可以將起始層位的各樣本點(diǎn)按照所屬類別進(jìn)行分類，并按 降序重新標(biāo)簽，記為1，2，···，根據(jù)上層標(biāo)簽順序，分別找出下層層位中空間區(qū)域中最多的標(biāo) 簽，將下層所有該標(biāo)簽轉(zhuǎn)換為上層該區(qū)域的標(biāo)簽，再將下層中對應(yīng)的上層標(biāo)簽轉(zhuǎn)換為已被 上層標(biāo)簽替換的標(biāo)簽，執(zhí)行完所有層位，完成調(diào)整。
[0110] 參見圖6A至6C，圖6A示出了層位面Η的聚類結(jié)果，圖6B中實(shí)線為層位面H1聚類 結(jié)果，虛線為層位面Η的聚類結(jié)果（假定層位面Η1與Η相鄰，Η1位于Η的下方），圖6C示 出了應(yīng)用聚類標(biāo)簽統(tǒng)一算法后，層位面Η1聚類標(biāo)簽的更新結(jié)果。
[0111] 具體地，首先，對地震體屬性數(shù)據(jù)進(jìn)行預(yù)處理，本實(shí)施例中選擇極大極小值法進(jìn)行 預(yù)處理。然后，設(shè)定程序界面輸入?yún)?shù)，確定起始層位（TopHorizon)、層位偏移量（Add TimeSlice)、聚類算法參數(shù)、半時(shí)間窗（HalfTime)。本實(shí)施例中波段聚類方式以基于屬性 拼合的方式為例。執(zhí)行完體聚類算法后，要進(jìn)行不同層位聚類結(jié)果的統(tǒng)一。從起始層位的 下一個(gè)層位開始，每個(gè)層位根據(jù)其上一個(gè)層位的聚類結(jié)果來調(diào)整其聚類結(jié)果。圖7是本發(fā) 明實(shí)施例的體聚類結(jié)果統(tǒng)一示意圖，圖7中，類別號為0~7。圖8是本發(fā)明實(shí)施例的地震 屬性聚類方法的具體實(shí)例的實(shí)施流程圖，即輸入?yún)?shù)，進(jìn)行聚類，聚類結(jié)果（標(biāo)簽）統(tǒng)一，存 儲(chǔ)統(tǒng)一后的標(biāo)簽，然后進(jìn)行可視化，方便用戶查看。