[發明專利]一種基于快速獨立分量分析的地震屬性融合方法

有效

• 申請號: 201210365999.0
• 申請日: 2012-09-28
• 公開（公告）號: CN102879823A
• 公開（公告）日: 2013-01-16
• 發明（設計）人: 彭真明;普艷香;趙敏;謝成芳;黃振星;李全忠;張萍
• 申請（專利權）人: 電子科技大學
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 成都華典專利事務所(普通合伙) 51223
• 代理人: 徐豐;楊保剛

1.一種基于快速獨立分量分析的地震屬性融合方法，其特征在于包括：

（1）對待融合屬性數據分塊：選取大小為N×N的窗W對各大小均為M1×M2待融合屬性切片平滑分塊，把分出來的每一個小塊存儲成N²×1的向量，放到各屬性塊的集合{Ii}里，得到一個塊的集合C，C={{I1},{I2},…,{In}}；

（2）選塊：根據集合C的維數和各小塊的均值，從集合C里面隨機的選取P塊屬性塊作為ICA算法的輸入來求取ICA變換的變換核T{·}；

（3）求取ICA域的變換核：由步驟2得到一個由P塊屬性塊組成的矩陣mixedX,其維數為N²×P，把mixedX作為快速獨立分量分析算法的輸入，求出mixedX的協方差矩陣CovX，CovX的維數為N²×N²，經過特征分解得到白化矩陣WhiteMattix和去白化矩陣deWhiteMattix，去白化矩陣與白化矩陣是互逆矩陣，通過求出的白化矩陣，求出白化信號Whitesig，把求得的白化信號作為快速獨立分量核心算法中的x，求出ICA域的變換核T{·}即分離矩陣W和ICA域的變換核的逆T^-1{·}即混合矩陣A；

（4）空間域映射到ICA域：求出各屬性對應塊的獨立分量，即把空域的塊映射到ICA域；

（5）ICA域融合：由步驟4把各屬性的所有塊均映射到ICA域里面，用如下的融合規則得到第i塊在ICA域的融合結果：

Fi(l)=Σk=1n|SSki(l)|SSki(l)Σk=1n|SSki(l)|,l=1,2,...,N2---(2)]]>

依照公式(2)求出所有塊在ICA域的融合結果；

（6）ICA域映射到空域：第i塊在空域融合的結果,

MFi=A×Fi+1nΣs=1nmsi---(3)]]>

依照公式(3)得到所有塊在空域的融合結果，然后把(M1-N+1)*(M2-N+1)*n塊屬性塊的空間域融合結果組合成M1*M2大小的融合屬性，對得到的融合結果采用求平均值法分塊去重疊處理，得到最終的融合結果。