[發明專利]井控稀疏表征大規模譜聚類地震相劃分方法

有效

• 申請號: 201811010363.8
• 申請日: 2018-08-31
• 公開（公告）號: CN109272029B
• 公開（公告）日: 2022-02-22
• 發明（設計）人: 王峣鈞;李昊陽;胡光岷
• 申請（專利權）人: 電子科技大學
• 主分類號: G06K9/62
• 分類號: G06K9/62;G06F17/16
• 代理公司: 成都虹盛匯泉專利代理有限公司 51268
• 代理人: 王偉

1.井控稀疏表征大規模譜聚類地震相劃分方法，其特征在于，包括以下步驟：

S1、輸入測井數據和測井數據所對應的巖相標簽；

S2、輸入地震數據所得屬性體，對該數據按照地震數據采樣方式進行網格化，將網格化坐標進行隨機采樣，得到k個采樣點及其對應的屬性曲線，隨機采樣點不包含井點所在位置；

S3、以測井數據作為標簽數據，采用支持向量機或神經網絡方法訓練分類器，并將該分類器用于對前述步驟隨機選擇點位置提取出的曲線進行分類操作；

S4、以隨機選擇點和測井數據點作為地標點，進行基于地標點的譜聚類操作，可得到井控譜聚類地震相劃分結果；具體實現方法為：選擇一定數量代表性數據點作為地標，并將原始數據點表示為這些地標的稀疏線性組合；然后利用基于地標的稀疏表示來高效地計算譜聚類中的特征向量從而加速譜聚類；具體包括以下子步驟：

S41、設計具有如下屬性的鄰接矩陣：W＝Z^TZ，其中Z∈R^p×n是稀疏的且p遠小于n；

S42、導出Z，在O(p³+p²n)中計算圖拉普拉斯矩陣的特征向量，O(p³+p²n)表示總的計算時間；Z的每一列都是原始輸入相對于新基向量組U的p維表示；因此，W＝Z^TZ是用新的表示法表示的數據X＝UZ的相似矩陣；

S43、對于任一數據點,近似表示成

其中u_j是U矩陣的第j個列向量，z_ji是矩陣Z的第j行i列的元素；用下式計算z_ji,

上式中K是一個計算相似度的核函數，i表示U的子矩陣，它包含了離數據點x_i最近的r個地標向量；

S44、通過SVD方法分解矩陣：

Z＝A∑B^T

中間矩陣為p乘p對角矩陣，其組成元素為被分解矩陣的p個奇異值，即特征值的算術平方根；A是p*p的左奇異矩陣，其列向量稱為左奇異向量，B是n*p的右奇異矩陣，其列向量為右奇異向量；根據數學定理可知B的列向量是W的特征向量；A的列向量是ZZ^T的特征向量；ZZ^T的維度是p乘p，計算A只需O(p³)的時間；

有了A之后，B通過下式被計算出來：

B^T＝∑^-1A^TZ

通過B求解W的特征向量。