[發明專利]一種基于疊前地震數據的地震相分析方法

有效

• 申請號: 201410438523.4
• 申請日: 2014-08-31
• 公開（公告）號: CN104181597B
• 公開（公告）日: 2017-05-17
• 發明（設計）人: 錢峰;王鑫;胡光岷
• 申請（專利權）人: 電子科技大學
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 成都宏順專利代理事務所(普通合伙)51227
• 代理人: 周永宏

1.一種基于疊前地震數據的地震相分析方法，其特征在于：包括以下步驟：

S1：輸入疊前地震數據；

S2：對輸入的疊前地震數據進行拆分，得到N個“類疊后”數據體，其中，N為地震信號道數；對輸入數據進行拆分的方法為：將M組疊前地震數據記為{D₁,D₂,…,D_M}，每組地震數據含有N道分別編號為1,2,…,N的信號，記為D_ij，其中，1≤i≤M,1≤j≤N，將{D₁,D₂,…,D_M}以{D₁₁,D₂₁,…,D_M1}，…，{D_1N,D_2N,…,D_MN}方式進行分組，共計分為N組；

S3：對每個“類疊后”數據體分別進行SOM聚類分析；

S4：利用聚類分析結果對疊前數據進行相似度分析，統計生成相似度矩陣，構建無向圖；

S5：利用廣度優先遍歷對無向圖進行聚類處理。

2.根據權利要求1所述的地震相分析方法，其特征在于：所述的步驟S3中進行SOM聚類分析的方法為：對N組數據分別進行聚類，定義在第k小組內的聚類為：

3.根據權利要求2所述的地震相分析方法，其特征在于：所述的步驟S4的相似度分析包括以下子步驟：

S41：定義相似度，定義數據D_i與D_j的相似度為：

$P(D_i,D_j)=\frac{1}{N}\underset{k=1}{\overset{k=N}{\Σ}}N(D_{ik},D_{jk})$

S42：定義連通性，0≤P(D_i,D_j)≤1，定義連通性為：

對于情況D_i,D_j連通和D_i,D_j斷開兩種情況，引入離群點的概念進行計算；

S43：定義離群點：若1≤j≤N且j≠i，則D_i為離群點，

S44：將D_i視為離群點，σ為人工設置的離群點闕值，σ越大，產生的離群點越多，聚類簇越緊密，反之離群點越少，聚類簇內部相對疏松，對于離群點的處理，將其單獨聚為一類或將其加入最近的一個聚類簇中；

S45：定義離群點與已有聚類簇的相似度為：

P(D_i,A_k)＝Max{P(D_i,D_j)|D_j∈A_k}，其中A_k為第k個聚類簇。

4.根據權利要求3所述的地震相分析方法，其特征在于：所述的步驟S5中利用廣度優先遍歷對無向圖進行聚類處理包括以下子步驟：

S51：從候選集合中，選取任一數據進行廣度優先遍歷，遍歷過程中遇到離群點，標記后暫不處理，直到遍歷完畢，將遍歷過的非離群點踢出候選集合，保存為新的聚類簇；

S52：判斷候選集合中是否還存在非離群點，如果集合中還有非離群點，則返回步驟S51，否則進行下一步；

S53：處理離群點：遍歷候選集合中的所有點，此時候選集合中沒有非離群點，將離群點加入到相似度最高的已有聚類簇中或將離群點單獨成一類；

S54：判斷候選集合是否為空，若為空則結束，否則返回步驟S53。