[发明专利]一种基于快速独立分量分析的地震属性融合方法

摘要

本发明涉及独立分量分析（ICA）技术领域和地震多属性融合领域，提供了一种基于快速独立分量分析（FICA）的地震多属性融合方法。本发明方案是通过把参与融合的属性都分成相同大小和块数的属性块，从这些块里面选出一定数量的块，应用FICA的原理对选出的块做处理，求得分离矩阵和与之互逆的混合矩阵，用分离矩阵把所有的块映射到ICA域。根据融合规则分别对各属性对应的块在ICA域做融合，最后用把ICA域的融合结果映射到空间域得到融合结果。融合结果有助于分析复杂的地层信息，提高储层预测的精度。本方法可广泛应用于地震属性分析、综合解释、地震储层预测、岩性与流体识别。

主权项

一种基于快速独立分量分析的地震属性融合方法，其特征在于包括：（1）对待融合属性数据分块：选取大小为N×N的窗W对各大小均为M1×M2待融合属性切片平滑分块，把分出来的每一个小块存储成N2×1的向量，放到各属性块的集合{Ii}里，得到一个块的集合C，C={{I1},{I2},…,{In}}；（2）选块：根据集合C的维数和各小块的均值，从集合C里面随机的选取P块属性块作为ICA算法的输入来求取ICA变换的变换核T{·}；（3）求取ICA域的变换核：由步骤2得到一个由P块属性块组成的矩阵mixedX,其维数为N2×P，把mixedX作为快速独立分量分析算法的输入，求出mixedX的协方差矩阵CovX，CovX的维数为N2×N2，经过特征分解得到白化矩阵WhiteMattix和去白化矩阵deWhiteMattix，去白化矩阵与白化矩阵是互逆矩阵，通过求出的白化矩阵，求出白化信号Whitesig，把求得的白化信号作为快速独立分量核心算法中的x，求出ICA域的变换核T{·}即分离矩阵W和ICA域的变换核的逆T‑1{·}即混合矩阵A；（4）空间域映射到ICA域：求出各属性对应块的独立分量，即把空域的块映射到ICA域；（5）ICA域融合：由步骤4把各属性的所有块均映射到ICA域里面，用如下的融合规则得到第i块在ICA域的融合结果： Fi ( l ) = Σ k = 1 n | SSki ( l ) | SSki ( l ) Σ k = 1 n | SSki ( l ) | , l = 1,2 , . . . , N 2 - - - ( 2 ) 依照公式(2)求出所有块在ICA域的融合结果；（6）ICA域映射到空域：第i块在空域融合的结果, MFi = A × Fi + 1 n Σ s = 1 n msi - - - ( 3 ) 依照公式(3)得到所有块在空域的融合结果，然后把(M1‑N+1)*(M2‑N+1)*n块属性块的空间域融合结果组合成M1*M2大小的融合属性，对得到的融合结果采用求平均值法分块去重叠处理，得到最终的融合结果。