[发明专利]一种叠后地震数据绕射波属性提取方法

权利要求书

1.一种叠后地震数据绕射波属性提取方法，其特征在于，包括以下步骤：

1)地震数据预处理；

2)根据实际数据建立正演模型，通过计算，确定主成分分析技术提取绕射波时，最适参数及主成分的维度；

3)将本征矢量计算公式改进为适用于三维地震数据计算式，并应用步骤2)中的最适参数及主成分维度，提取实际地震数据中的绕射波场；

根据三维地震数据特点，通过下式提取空间绕射波场数据：

式中,是计算矩阵中位置为i,j,k的采样点上第r个正交分量；s_i+n,j+m,k+l是计算矩阵中位置为i,j,k的地震数据采样点值；n,m,l是i,j,k方向上的延迟距；N,M,L为最大延迟距；s_avr为计算矩阵中样点振幅的平均值；s_max,s_min分别为计算矩阵中样点振幅的最大值和最小值；分别为矩阵特征值和特征向量。

2.如权利要求1所述方法，其特征在于：所述步骤1)中，预处理包括以下步骤：

1.1)定义如下变量：A＝{A₁,A₂,...,A_n}是一个m×n的矩阵，n为输入的地震道数，每个地震道为m维向量，代表m个采样点；A_n表示第n个地震道的地震数据采集序列；进行地震数据体的标准化预处理：

式中，x_pq为标准化后的地震数据采集序列；为样本均值：A_pq为原始采集的地震数据采集序列；S_p为样本标准差：p＝1,…,n；q＝1,…,m；

1.2)处理后得到标准化地震数据体为X＝{x₁,x₂,...,x_n}；其中，X是一个m×n的矩阵，n为输入的地震道数，每个地震道为m维向量，代表m个采样点，x_n表示第n个地震道标准化后的地震数据采集序列。

3.如权利要求1所述方法，其特征在于：所述步骤2)中，最适参数及主成分的维度确定方法包括以下步骤：

2.1)通过下式将输入数据矢量x_g变换成新的矢量s_g：

s_g＝U^Tx_g

其中，U是一个n×n正交矩阵，它的第g列U_g是样本协方差矩阵；g＝1，…，n；

2.2)计算第g个本征矢量：

λ_gu_g＝Cu_g，g＝1,...,n (1)

其中，C为一个n阶方阵；λ_g是C的一个本征值，u_g是对应的本征矢量；

2.3)通过上述步骤对地震模型中绕射波场进行提取，当达到压制反射波同相轴同时充分保留小尺度异常体绕射波信息的效果时，确定最适参数及主成分维度。