[发明专利]一种叠前方位P波各向异性裂缝检测方法及装置

权利要求书

1.一种叠前方位P波各向异性裂缝检测方法，其特征在于，所述裂缝检测方法包括以下步骤：

1)沿目的层提取单元区间内基于方位的叠加道集地震振幅数据；

2)将步骤1)获取的各方位振幅信息在极坐标中进行投影；

3)对步骤2)中投影得到的一组样本点进行椭圆拟合；

4)求解出拟合的椭圆方程的最优椭圆参数解；

5)基于所述最优椭圆参数解确定椭圆几何参数。

2.根据权利要求1所述的叠前方位P波各向异性裂缝检测方法，其特征在于，基于椭圆一般方程Ax²+Bxy+Cy²+Dx+Ey+F＝0定义投影平面上一个样本点(x,y)到椭圆的几何距离为：

d(a,x)＝Ax²+Bxy+Cy²+Dx+Ey+F (1)

其中，a＝[A B C D E F]^T，x＝[x² xy y² x y 1]^T。

3.根据权利要求2所述的叠前方位P波各向异性裂缝检测方法，其特征在于，通过以下公式求得最小二乘意义下的最优椭圆参数解：

4.根据权利要求3所述的叠前方位P波各向异性裂缝检测方法，其特征在于，通过将方程式(1)表示为广义特征值问题的形式：

G^TGa＝λHa (3)

其中，λ为特征值，H为约束矩阵，a为特征向量a＝[A B C D E F]^T，

式中x₁,x₂,...,x_n为第1、第2…第n点的横坐标，

y₁,y₂,...,y_n为第1、第2…第n点的纵坐标；

基于椭圆参数的数值关系约束矩阵求解方程式(3)以得到特征向量a。

5.根据权利要求4所述的叠前方位P波各向异性裂缝检测方法，其特征在于，通过以下公式计算椭圆几何参数：

其中，(x_e,y_e)为椭圆中心相对坐标，a、b为椭圆长轴长度、短轴长度，θ为椭圆长轴与正北方向的夹角。

6.一种叠前方位P波各向异性裂缝检测装置，其特征在于，所述装置包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现以下步骤：

1)沿目的层提取单元区间内基于方位的叠加道集地震振幅数据；

2)将步骤1)获取的各方位振幅信息在极坐标中进行投影；

3)对步骤2)中投影得到的一组样本点进行椭圆拟合；

4)求解出拟合的椭圆方程的最优椭圆参数解；

5)基于所述最优椭圆参数解确定椭圆几何参数。

7.根据权利要求6所述的叠前方位P波各向异性裂缝检测装置，其特征在于，基于椭圆一般方程Ax²+Bxy+Cy²+Dx+Ey+F＝0定义投影平面上一个样本点(x,y)到椭圆的几何距离为：

d(a,x)＝Ax²+Bxy+Cy²+Dx+Ey+F (1)

其中，a＝[A B C D E F]^T，x＝[x² xy y² x y 1]^T。

8.根据权利要求7所述的叠前方位P波各向异性裂缝检测装置，其特征在于，通过以下公式求得最小二乘意义下的最优椭圆参数解：