[发明专利]地震相干属性计算方法及系统

权利要求书

1.一种地震相干属性计算方法，包括：

步骤1：根据地震资料，获得工区内的采样点分布图；

步骤2：针对每个采样点，重复执行子步骤201-203：

子步骤201：针对所述采样点设定采样点坐标系和菱形边界图，其中所述采样点作为所述采样点坐标系的原点和所述菱形边界图的长短轴交点，所述菱形边界图的长轴与所述采样点坐标系的X轴重合，所述菱形边界图的短轴与所述采样点坐标系的Y轴重合；

子步骤202：以所述采样点为旋转中心，对所述菱形边界图进行旋转，获得所述菱形边界图内的有效点；

子步骤203：设定时窗，根据所述菱形边界图内的有效点的时窗内的数据，计算每个旋转角度对应的相干值，进而确定所述采样点的最大相干值；

步骤3：根据每个采样点的最大相干值，获得所述工区的相干属性图。

2.根据权利要求1所述的地震相干属性计算方法，其中，获得所述菱形边界图内的有效点包括：

将当所述菱形边界图旋转角度θ时处于所述菱形边界图内的所有其他采样点标记为有效点。

3.根据权利要求2所述的地震相干属性计算方法，还包括：

设置所述菱形边界图的外接圆的外接正方形为扫描范围。

4.根据权利要求3所述的地震相干属性计算方法，其中，获得所述菱形边界图内的有效点包括：

在所述扫描范围内，将当所述菱形边界图旋转角度θ时处于所述菱形边界图内的所有其他采样点标记为有效点，并记录所述有效点的编号。

5.根据权利要求4所述的地震相干属性计算方法，其中，子步骤203包括：

记录每个旋转角度的有效点的编号，根据每个旋转角度的有效点的编号，获取所述有效点的时窗内的数据，进而计算每个旋转角度对应的相干值。

6.根据权利要求1所述的地震相干属性计算方法，其中，根据以下公式(1)计算相干值：

其中，C₃为相干值，C为有J道每道有N个样点组成的地震数据体，D_N*J为协方差矩阵，其中D_N*J为参与计算的由J道N个样点构成的地震子体，是D_N*J的转置，Tr(C)为矩阵的迹，代表协方差矩阵的能量，λ₁是最大特征值，代表占优的能量，c_jj是矩阵对角线上的元素，λ_j是矩阵非负特征值，J为参与计算的地震道数，j为参与计算的地震道编号。

7.根据权利要求1所述的地震相干属性计算方法，其中，所述菱形边界图的长短轴之比为5:2。

8.一种地震相干属性计算系统，其特征在于，该系统包括：

存储器，存储有计算机可执行指令；

处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：

步骤1：根据地震资料，获得工区内的采样点分布图；

步骤2：针对每个采样点，重复执行子步骤201-203：

子步骤201：针对所述采样点设定采样点坐标系和菱形边界图，其中所述采样点作为所述采样点坐标系的原点和所述菱形边界图的长短轴交点，所述菱形边界图的长轴与所述采样点坐标系的X轴重合，所述菱形边界图的短轴与所述采样点坐标系的Y轴重合；

子步骤202：以所述采样点为旋转中心，对所述菱形边界图进行旋转，获得所述菱形边界图内的有效点；

子步骤203：设定时窗，根据所述菱形边界图内的有效点的时窗内的数据，计算每个旋转角度对应的相干值，进而确定所述采样点的最大相干值；

步骤3：根据每个采样点的最大相干值，获得所述工区的相干属性图。

9.根据权利要求8所述的地震相干属性计算系统，其中，获得所述菱形边界图内的有效点包括：

将当所述菱形边界图旋转角度θ时处于所述菱形边界图内的所有其他采样点标记为有效点。

10.根据权利要求8所述的地震相干属性计算系统，其中，根据以下公式(1)计算相干值：

其中，C₃为相干值，C为有J道每道有N个样点组成的地震数据体，D_N*J为协方差矩阵，其中D_N*J为参与计算的由J道N个样点构成的地震子体，是D_N*J的转置，Tr(C)为矩阵的迹，代表协方差矩阵的能量，λ₁是最大特征值，代表占优的能量，c_jj是矩阵对角线上的元素，λ_j是矩阵非负特征值，J为参与计算的地震道数，j为参与计算的地震道编号。