[發(fā)明專利]一種疊后地震信號波形分類方法

有效

• 申請?zhí)枺?/td>: 201710154208.2
• 申請日: 2017-03-15
• 公開（公告）號: CN106707335B
• 公開（公告）日: 2018-06-19
• 發(fā)明（設(shè)計）人: 石好果;張曰靜;徐佑德;許濤;商豐凱;楊東根;陳林;王林;劉子超;韓梓航;于騰飛;張新見;張奉
• 申請（專利權(quán)）人: 中國石油化工股份有限公司勝利油田分公司勘探開發(fā)研究院西部分院
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京科家知識產(chǎn)權(quán)代理事務(wù)所(普通合伙) 11427
• 代理人: 莫文新

1.一種疊后地震信號波形分類方法，包括如下步驟：

(1)提取待分析的地震數(shù)據(jù)，并進(jìn)行預(yù)處理；

(2)使用稀疏編碼方法提取地震數(shù)據(jù)中的特征數(shù)據(jù)x，并對特征數(shù)據(jù)x進(jìn)行降維；

(3)將步驟(2)中得到的特征數(shù)據(jù)x利用聚類算法進(jìn)行分類；

其特征在于，

所述步驟(3)包括下述具體步驟：

(3-1)初始化

(3-1-1)給定用于確定截斷距離d_c的參數(shù)t∈(0,1)；

(3-1-2)計算任意兩樣本數(shù)據(jù)間的距離

(3-1-3)將升序排列，計算截斷距離

(3-1-4)計算樣本數(shù)據(jù)的局部密度使局部密度降序排列的下標(biāo)x_i排列生成即

(3-1-5)計算其中，表示所有局部密度比x_i大的數(shù)據(jù)點x_j中與x_i距離最近的數(shù)據(jù)點的編號，

(3-2)確定聚類中心并初始化數(shù)據(jù)點歸類屬性標(biāo)記

(3-3)按照ρ從大到小遍歷數(shù)據(jù)對非聚類中心數(shù)據(jù)點歸類；

(3-4)若簇的個數(shù)N_c>1，則將每個簇中數(shù)據(jù)點進(jìn)一步分為core和halo，core的局部密度大，對應(yīng)簇核心部分，halo的局部密度小，對應(yīng)簇邊緣部分，設(shè)halo部分為若h_i＝0，則x_i屬于core部分，反之x_i屬于halo部分。

2.根據(jù)權(quán)利要求1所述的疊后地震信號波形分類方法，其特征在于，所述步驟(1)的具體步驟包括：

(1-1)對算法中輸入的寬方位角地震數(shù)據(jù)進(jìn)行降噪；

(1-2)提取待分析層位段地震數(shù)據(jù)，每個層位段地震數(shù)據(jù)的樣本點數(shù)t，組成訓(xùn)練樣本：

X＝{x₁,...x_n}，

其中，x_i為維度t的向量。

3.根據(jù)權(quán)利要求1所述的疊后地震信號波形分類方法，其特征在于，所述步驟(2)中：

輸入：字典D，目標(biāo)信號y，稀疏度K，初始化殘差r₀＝y(tǒng)，索引集Λ＝φ，t＝1；

輸出：y在字典D中的K稀疏逼近所對應(yīng)的稀疏系數(shù)；

具體步驟包括：

(2-1)找出殘差與字典矩陣中列內(nèi)積最大的項，其列號記為λ，即：

(2-2)更新索引集Λ＝Λ∪{λ}，并根據(jù)更新后的索引集構(gòu)建新的重構(gòu)向量集合

(2-3)對向量Φ_t和x_t進(jìn)行處理，通過最小二乘得到：

x_t＝argmin||y-Φ_tx_t||₂ 式(2-2)；

(2-4)更新殘差r_t＝y(tǒng)-Φ_tx_t，t＝t+1；

(2-5)重復(fù)步驟(2-1)-(2-4)，直到滿足停止條件：t＞k。

4.根據(jù)權(quán)利要求1所述的疊后地震信號波形分類方法，其特征在于，所述步驟(3-4)的步驟包括：

(3-4-1)初始化h_i＝0；(3-4-2)為每一個簇生成一個平均局部密度上界

(3-4-3)標(biāo)記簇halo部分。