[发明专利]一种基于三维地震资料表征的地下火山通道识别方法

权利要求书

1.一种基于三维地震资料表征的地下火山通道识别方法，其特征在于，包括：

步骤1：对叠后三维地震数据体进行各项异性扩散滤波，获取扩散滤波后的地震数据体；

步骤2：基于所述扩散滤波后的地震数据体，获取相似性数据体；

步骤3：对所述相似性数据体进行二值化处理，获得二值属性数据体；

步骤4：基于所述二值属性数据体，计算目的层以下的二值属性累计值；

步骤5：对所述二值属性累计值进行增益计算，获得增益后属性；

步骤6：搜索所述增益后属性，确定局部最值点集，通过所述局部最值点集验证地下火山通道预测结果，获得地下火山通道坐标集；

步骤7：基于所述地下火山通道坐标集和所述二值属性数据体，获得地下火山通道点集空间分布；

所述步骤7包括：根据以下步骤获得所述地下火山通道坐标集中的每个坐标点的空间分布：

以所述二值属性数据体的任意平面时间切片为基础，采用圆形表征任意平面时间切片上地下火山通道的等效聚合分布形态，以所述地下火山通道坐标点为中心构建搜索窗，在所述搜索窗内计算所述地下火山通道坐标点的数据汇聚中心及所述圆形的半径，获得所述地下火山通道点的空间分布。

2.根据权利要求1所述的基于三维地震资料表征的地下火山通道识别方法，其特征在于，所述步骤2包括：采用纵横比大于第一阈值的时窗计算所述扩散滤波后的地震数据体的相似性，获得相似性数据体。

3.根据权利要求1所述的基于三维地震资料表征的地下火山通道识别方法，其特征在于，所述步骤3包括：采用二值截断法对所述相似性数据体进行二值化处理，获得所述二值属性数据体。

4.根据权利要求1所述的基于三维地震资料表征的地下火山通道识别方法，其特征在于，所述步骤4包括：基于所述二值属性数据体，采用大于或等于第一预设时间的时窗计算目的层以下的二值属性累计值，获得二值体累计属性平面分布。

5.根据权利要求1所述的基于三维地震资料表征的地下火山通道识别方法，其特征在于，所述步骤5包括：采用滑动增强算法对所述二值属性累计值进行增益计算，如以下公式(1)所示：

其中，F(x,y)为增益后属性，L为滑动窗口边长的一半，V(x,y)为t₀至t₁时窗内二值属性累计值，V(x,y,t)为(x，y)道在t时间切片时的二值属性，t为t₀至t₁时间内任一时刻。

6.根据权利要求1所述的基于三维地震资料表征的地下火山通道识别方法，其特征在于，所述搜索所述增益后属性，确定局部最值点集包括：

步骤61：将所述增益后属性按搜索窗大小划分为m×n个网格，其中，m为横向格数，n为纵向格数；

步骤62：分别搜索每个网格的局部最大值点，获得局部最值点列表，所述局部最值点列表包括每个网格的局部最大值点的线道号和增益后属性；

步骤63：分别比较所述局部最值点列表中的每个局部最大值点的增益后属性与预设截断阈值，当所述局部最大值点的增益后属性小于所述预设截断阈值时，则删除所述局部最大值点，否则保留所述局部最大值点，形成局部最值点集。

7.根据权利要求6所述的基于三维地震资料表征的地下火山通道识别方法，其特征在于，在所述步骤62中，根据以下步骤搜索一个所述网格的局部最大值点：

步骤621：针对所述网格，以当前搜索窗在所述网格内进行最大值搜索，获得第一最大值点；

步骤622：以所述第一最大值点为中心构建新的搜索窗，以所述新的搜索窗在所述网格内进行最大值搜索，获得第二最大值点，若所述第一最大值点与所述第二最大值点重合，将所述第二最大值点作为所述网格的局部最大值点，否则以所述第二最大值点作为所述第一最大值点，重复所述步骤622。