[发明专利]一种基于矩阵降维的特征结构相干体快速实现方法

权利要求书

1.一种基于矩阵降维的特征结构相干体快速实现方法，其特征在于，包括以下步骤：

步骤01：采集得到三维地震数据体，并设定二维空间分析窗大小和一维时间分析窗大小；

步骤02：依照预先设定的二维空间分析窗大小从步骤01采集的三维地震数据体中切取得到一个小三维数据体；

步骤03：在小三维数据体中，计算每一时刻x方向相干度量和y方向相干度量；

步骤04：在小三维数据体内将所有时刻的x方向相干度量和y方向相干度量相乘得到最终的相干度量；

步骤05：判断x及y是否遍历完毕，若完毕输出三维相干度量数据体，否则改变x及y并返回步骤02；

步骤01具体包括：

用f(x,y,t)表示大小为N_x×N_y×N_t的三维地震数据体，x和y表示三维地震数据体在两个空间方向的指标，t表示三维地震数据体在时间方向的指标，其中x和y的变化范围分别为[1,N_x]和[1,N_y]，t的变化范围为[1,N_t]；设定的分析窗大小为M_x×M_y×M_t；其中M_x与M_y为二维空间分析窗的大小，M_t为时间分析窗的大小，且三者均为奇数；由于M_x、M_y及M_t均为奇数，将其表示为M_x＝2W_x+1，M_y＝2W_y+1，M_t＝2W_t+1，其中W_x、W_y及W_t为x、y及t方向的半窗大小；

步骤03中具体包括：

步骤031：依据其前一时刻的x方向协方差矩阵对当前时刻的x方向协方差矩阵进行更新，以获得每一时刻的x方向协方差矩阵；

步骤032：在小三维数据体中对每一时刻，对其对应的x方向协方差矩阵进行特征分解以获得最大特征值，进而获得每一时刻x方向相干度量；

步骤033：在小三维数据体中对每一时刻，依据其前一时刻的y方向协方差矩阵对当前时刻的y方向协方差矩阵进行更新，以获得每一时刻的y方向协方差矩阵；

步骤034：在小三维数据体中在每一时刻，对其对应的y方向协方差矩阵进行特征分解以获得最大特征值，进而获得每一时刻y方向相干度量。

2.根据权利要求1所述的一种基于矩阵降维的特征结构相干体快速实现方法，其特征在于，步骤02中，依照预先设定的空间分析窗尺寸并以x及y为中心从三维地震数据体中切取得到一个小三维数据体，具体包括：

以x及y为中心从三维地震数据体中划取得到一个小三维数据体；其大小为M_x×M_y×N_t，其在原始三维地震数据体中的指标范围为[x-W_x,x+W_x]、[y-W_y,y+W_y]及[1,N_t]；同时假设小三维数据体内沿x方向、y方向的局部指标为mx及my，两者的变化范围为[1,M_x]及[1,M_y]。